Abstract
Background
During the HIV infection several quasispecies of the virus arise, which are able to use different coreceptors, in particular the CCR5 and CXCR4 coreceptors (R5 and X4 phenotypes, respectively). The switch in coreceptor usage has been correlated with a faster progression of the disease to the AIDS phase. As several pharmaceutical companies are starting large phase III trials for R5 and X4 drugs, models are needed to predict the coevolutionary and competitive dynamics of virus strains.
Results
We present a model of HIV early infection which describes the dynamics of R5 quasispecies and a model of HIV late infection which describes the R5 to X4 switch. We report the following findings: after superinfection (multiple infections at different times) or coinfection (simultaneous infection by different strains), quasispecies dynamics has time scales of several months and becomes even slower at low number of CD4+ T cells. Phylogenetic inference of chemokine receptors suggests that viral mutational pathway may generate a large variety of R5 variants able to interact with chemokine receptors different from CXCR4. The decrease of CD4+ T cells, during AIDS late stage, can be described taking into account the X4related Tumor Necrosis Factor dynamics.
Conclusion
The results of this study bridge the gap between the withinpatient and the interpatients (i.e. worldwide) evolutionary processes during HIV infection and may represent a framework relevant for modeling vaccination and therapy.
Background
The relationship between phenotype and survival of the genotype is central to both genetics and evolution. Viruses represent a good sizedcomplexity phenotype to study and they show rapid evolution. Selection pressures mainly depend on the interaction strength with the receptor vital for the entry into the target cell. Human immunodeficiency virus type 1 (HIV1) infection is characterized by the progressive loss of CD4+ T cells. Infection by most strains of HIV1 requires interaction with CD4 and a chemokine receptor, either CXCR4 or CCR5. During early stages of HIV1 infection, viral isolates most often use CCR5 to enter cells and are known as R5 HIV1. Later in the course of HIV1 infection, viruses that use CXCR4 in addition to CCR5 (R5X4) or CXCR4 alone (X4 variants) emerge in about 50% patients (switch virus patients) [1,2]. These strains are syncytiuminducing and are capable of infecting not only memory T lymphocytes but also naive CD4+ T cells and thymocytes through the CXCR4 coreceptor. The switch to using CXCR4 has been linked to an increased virulence and progression to AIDS, probably through the formation of cell syncytia and apoptosis of T cell precursors. X4 HIV strains are rarely, if ever, transmitted, even when the donor predominantly carries X4 virus. Clevestig [3] found that in children, the X4 virus developed from their own R5 population, and was not caused by transmission from the mother. CXCR4 is expressed on a majority of CD4+ T cells and thymocytes, whereas only about 5 to 25% of mature T cells and 1 to 5% of thymocytes express detectable levels of CCR5 on the cell surface [4]. It is noteworthy that X4 HIV strains stimulate the production of cellular factor called Tumor Necrosis Factor (TNF), which is associated with immune hyperstimulation, a state often implicated in Tcell depletion [5].
Tumor necrosis factor (TNF)related apoptosisinducing ligand (TRAIL) is produced mainly by monocytes; it is a type II transmembrane protein belonging to the TNF family, homologous to Fas ligand (FasL), which is a wellcharacterized apoptosisinducing ligand. In particular the Fasassociated domain or a Fasassociated domainlike adaptor molecule leads to activation of the caspase cascade, resulting in apoptotic cell death. TNF seems able to both inhibit the replication of R5 HIV strains while having no effect on X4 HIV and to down regulate the number of CCR5 coreceptors that appear on the surface of Tcells [6]. Plasma TRAIL is elevated in HIV1 infected patients and is decreased by Highly Active Anti Retroviral Therapy (HAART). Thus, when HAART decreased viral load, there is a concomitant decrease in plasma TRAIL, which may be one of the reasons for the efficacy of antiviral therapy [7].
Mathematical modeling
The use of mathematical models is an insightful and essential complement to in vivo and in vitro experimental design and interpretation. Indeed mathematical models of HIV dynamics have proven valuable in understanding the mechanisms of many of the observed features of the progression of the HIV infection, see for example [815].
A powerful concept in understanding HIV variability and its consequences is that of quasispecies, accounting for the result of evolution not being the selection of a single sequence (genotype), rather a distribution of quasiidentical sequences, termed the quasispecies, centered around a master sequence. Quasispecies are the combined result of mutations and recombination, that originate variability, and of coinfection (simultaneous infection), superinfection (delayed secondary infection) and selection, that keep variability low. HIV1infected individuals show heterogeneous viral populations, best described as viral quasispecies [16]. Infact, the infection capacity of mutants may vary and also their speed of replication [17].
Moreover, since the number of targets (the substrate) is limited, fitter clones tend to eliminate less fit mutants, which are subsequently regenerated by the mutation mechanism. Mutations are a key ingredient for exploring the genetic space in the search for the fitness maximum, but are also responsible for the disappearance of the quasispecies when exceeding a critical mutation rate, the error threshold [18].
Taking into account stateoftheart models of HIV infection, we address the issue of studying the coevolutive and competitive dynamics of different strains of HIV1 virus also leading to the R5 to X4 phenotype switching. Ribeiro and colleagues [19] have recently presented a model of R5 to X4 switch based on the hypothesis that X4 and R5 viruses have a preferential tropism for naive and memory T cells, respectively. Here we prefer to follow the mutational hypothesis supported by several experiments [14]. In the next section we introduce two models: a quasispecies model for R5 phase in which several R5 strains appear by mutations, coinfection and superinfection. In the limit of a single quasispecies we are able to find the same values observed experimentally and in other models (most notably Perelson's standard model). We test the model in the scenarios of coinfection and superinfection using parameters derived from biological and medical literature. A second model focuses on the R5 to X4 shift and the hyperstimulation of T cell precursors through TNF. The results of the numerics reproduce well the decreasing dynamics of CD4+ after the appearance of X4 strains and make the model suitable for further investigations on antiretroviral therapies and their effects on disease progression. The R5 to X4 switch is also investigated by using phylogenetic models of the amino acid sequences of the human and mice chemokine receptor families. The analysis of the mutational pathway suggests that the switch from R5 to X4 may allow the HIV to bind to other chemokine receptors, thus likely leading to immune system signaling disfunctions.
Results and discussion
Modeling HIV dynamics and the immune response
The basic model of HIV1 dynamics, first introduced by Perelson and colleagues in 1995 [8], is a class model considering three variables: uninfected cells, infected cells and free viruses. A system of three differential equations describes how these quantities change over time. In the basic model by Perelson uninfected cells encounter free virus and turn into infected cells. The rate of production of infected cells is proportional to the product of the density of uninfected cells times the density of free virions. Free virions are produced by infected cells. Uninfected cells, infected cells and free virions die at fixed rate. Uninfected cells are also assumed to be constantly replaced by the immune system. These assumption define the basic model of virus dynamics.
We extended this basic model of HIV1 dynamics, by including the response by the immune system and the role of B lypmhocytes, see Figure 1. There are several works that discuss the role of B cells in the immune response [20]. Since the progression to AIDS has been found to correlate well with CD4+ T cells decrease, B cells are thought to play a minor role in the immune response to HIV. Note that B cells can act only as predator to the HIV, so their coupling with HIV dynamics is different from that of T cells. Our aim is to present a more general model framework of both T and B immune responses to HIV. We have first considered the following system of differential equations describing the dynamics of a single viral strain:
Figure 1. Representation of the interactions between cells and/or cells and viruses. Representation of the interactions between cells and/or cells and viruses. HIV strains infect CD4+ T cells that become infected and produce new viruses. At the same time HIV peptides are presented via APC cells to THelper cells that become activated. Activated CD4+ T cells trigger B and CD8+ cells reactions: the first release antibodies that bind to the antigen while the latter directly remove infected CD4+ T cells. '+' and '' signs indicate cell/virus production or removal.
This model considers the Thelper (CD4+) cells (T) and HIV virus particles (V); the T cells can become infected (I). With respect to Refs. [8,9], in Eq. (1), we describe how the number of naive T cells (T) which have passed the thymus selection, depends on rate of formation in the bone marrow (λ) and on clonal amplification upon stimulation by infected cells, I (term IT). They decrease with a rate that is the sum of a natural clearance (δ_{T}T) due to cell aging and cell destruction upon virus infection, (VT). The density of T cells is limited by a saturating density/lymphonode capacity factor, K. Following Ref. [14] we have set K = 10^{12}.
The second equation describes the rate of infection, described by β, of naive T cells upon the interaction with the virus (term VT). Infected T cells are cleared out at a fixed rate, δ_{I}, and due to the action of natural killer cells, CD8+ (term TI).
The third equation describes the budding of viruses from infected cells, π. Virus particles are cleared out at rate c (defective viruses) and after immunoglobulin binding and subsequent engulfments by the macrophages (term TV).
The γ parameters have the same meaning of the constant of association in chemistry, or can be thought as a combination of both the probability of interaction and the interaction strengths between cells (γ^{(T)}, γ^{(I)}) or between cells and viruses (γ^{(V)}). Note that in the limit γ → 0 and K → ∞, we recover the pattern of the standard model [13].
The B cell response is modeled using the parameters corresponding to the T cells which activate them by receptor recognition. Here, we have assumed the immunoglobulin concentration, which represent the B cell response, to be linearly correlated to the concentration of activated B cells. These in turn are supposed linearly correlated to the concentration of the CD4+ T cells. Exploratory analysis with different parameters, suggested us that Eq. (3) allows to keep a minimum number of parameters without loss of important details.
A quasispecies framework
A meaningful model to study the evolution of a viral population is that of quasispecies. The quasispecies concept was put forward by Manfred Eigen and Peter Schuster [21] within the framework of physical chemistry. Their model accounts for the description of the process of the Darwinian evolution of selfreplicating entities and shows that the result of evolution is not the selection of a single sequence, rather a distribution of sequences centered around a master sequence. The quasispecies behavior of viruses is nowadays recognised as a key element for the understanding and modeling of viral evolution and disease control [22]. The extreme mutation rate of RNA genetic elements, together with their large population sizes, endows them with the characteristics of quasispecies populations. Most progeny genomes of individual (clonal) RNA viruses generally differ from one another in having one or more mutations scattered randomly along their genomes; thus comprising a quasispecies cloud or swarm of related, but not identical variants.
A quasispecies model of HIV infection can be derived by the model previously introduced. We adopt a quasispecies description of the process by first hypothesizing the evolution of the virus to be characterized by only one slight changing phenotypic trait, e.g. the shape of the envelope and its corresponding binding affinity. This assumption is justified if the phenotypes are determined by few viral protein functional determinants which are both independent and differ only few DNA bases, i.e. few mutations can change one determinant into another. This corresponds to map the different phenotypes of the virus on a linear space (see [23,24] for similar assumptions). In fact this binding does not need to be a perfect match, but has a certain degree of tolerance. Similar phenotypes compete for the same coreceptor and mutations can be the key factor to give origin to a new phenotype able to use a different coreceptor. On the other hand T lymphocytes look forward to kill infected T cells and free virions. Their recognition ability of viral antigens depends on the binding between a T cell receptor (TCR) and a viral epitope.
Based upon the above schematization, the population of T cells can be thought of as composed by different phenotypes corresponding to the different TCRs and the virus populations as characterized by a specific binding ability. Thus the clonal amplification of naïve T cells (of population i) depends on the ability of (i^{th }class) T cell to recognize all the infected T cells (carrying an epitope from the k^{th }class of viruses). The rate of infection of naïve T cells (of class i) depends on all the viruses (containing the epitope k). In the same way, the clonal expansion of infected T cells (of class k, meaning that they have been infected by a virus of class k irrespective of the original T class) depends on the interaction of these viruses with all the T cells (of any class i).
It's noteworthy that this schematization is valid for the quasispecies case and for the undifferentiated one. In the latter case, the indication of strains (in parenthesis) should be neglected.
A model of the early R5 phase evolution
In someone who is newly infected by HIV, several variants of the virus, called R5, are often the only kind of virus that can be found. A meaningful way to model strain mutations, coinfection and superinfection is to extend the model previously introduced by incorporating multi strain evolutive dynamics. Thus the R5 quasispecies dynamics can be described by the following set of equations:
The following cell types are considered: Thelper (CD4+) cells carrying the CCR5 coreceptor responding to virus strain i, (T_{i}); T cells infected by virus strain k, (I_{k}); k strains of R5 virus, (V_{k}). We have thus assumed that viral strain k are identified by just one epitope, which is then displayed on the surface of the T cell of class k, and that a T cell of class i can be activated at least by one CD4+ T cell carrying the epitope k, which is specific of the viral strain k. The indices i (k) range from 1 to N_{i }(N_{k}), and in the following we have used N_{i }= N_{k }= N.
In Equation (4) T cells are generated through two mechanisms: the bonemarrow source (and selection in the thymus) and the duplication of T cell strains activated upon the recognition with an antigen carrying cell that may be even an infected one. We modeled T cells activation as a logistic term mimicking the global carrying capacity (the K parameter) of immune system [25].
The deathrate term is composed by a natural death rate proportional to the population, and by the infection rate of T cells due to any viral strain. The term ∑_{k}β_{k}V_{k}T_{i }and the sum over T_{i }in the I cell birth rate reproduce the infection probability, that is the same irrespective of the T class. As a cell become infected it does no more contribute to the immune response.
Equation (5) describes the infection dynamics. The two death rate parameters account for the decrease of the infected cells due to cellular death and after the action of Tkiller cells (CD8+). Even if there are clear experimental evidences that CD4+ cells decrease during the late HIV infection stages and in the AIDS state, as far as the asymptomatic phase of the infection is concerned, the parameter δ_{I }may be assumed as a constant, medical literature referred, value.
The μ_{kk' }term is responsible of the mutation process affecting the phenotype, essential for the formation of new quasispecies. The choice of a mutation rate of the order of 10^{5 }is based on considering only those nonsynonymous mutations that alter the phenotype (protein structure) [25].
In Equation (6) the virus replicative dynamics is described. The birth rate term is proportional to the virus "budding" numerosity while the viral death rate parameters depend on the rate of natural death and accounts for the recognition of virus by B cells.
It's worth noting that Bcells and Tkiller cells are only implicitly included in the model in order to reduce the dimensionality without loosing too many details. We assume that these responses are fast enough to be at equilibrium and they are just proportional to the abundance of (cognate) T helper cells.
The three γ parameters γ^{(T)}, γ^{(I) }and γ^{(V) }are matrices describing the interactions between cells and/or cells and viruses, i.e. who will interact with whom, in terms of geometry and strength of the interaction. It is thus possible to consider which strains of the virus are recognized and with which accuracy, and the same for the action of Tkiller cells and Bcells. It is also worth noting that is the most important determinant of the viral fitness [25].
Modeling the transition R5 to X4
In about half of the people who develop advanced HIV disease, the virus begins to use another coreceptor called CXCR4 (X4 viral phenotype). The shift to using CXCR4 it is generally accompanied by a dramatic increase in the rate of Tcell depletion. The inability of the thymus to efficiently compensate for even a relatively small loss of naive T cells may be a key factor for CD4+ T cells depletion and AIDS progression. We hypothesize that it may not be exhaustion of homeostatic responses, but rather thymic homeostatic failure along with gradual wasting of T cell supplies through hyperactivation of the immune system that lead to CD4 depletion in HIV1 infection.
We here introduce a modified version of the previous model, that considers only CD4 dynamics. Neither B nor CD8+ T cells are explicitly modeled and the space that is considered is that of different phenotypes of the virus. Under those assumptions we may focus on the appearance of X4 viruses and on their subsequent interaction with R5 strains.
In the equations above, the variables modeled are the pool of immature CD4+ T cells, U, the different strains of uninfected and infected T cells (T and I, respectively), HIV virus, V, and the excess of TNF concentration, F. TNF concentration, in uninfected people, i.e. when X4 = 0, has a low basal value. Here we consider the excess of TNF abundance, by assuming its concentration to be proportional to the X4 concentration. A schematic view of the model is depicted in Figure 2. The value of the parameters introduced with respect to the R5 model are summarized in Table 1.
Table 1. Additional parameters introduced in the R5 to X4 phenotypic switch model.
Figure 2. Schematic description of the model for the switching from R5 to X4 viral phenotype. Schematic description of the model for the switching from R5 to X4 viral phenotype. Naïve Tcells, U, are generated at constant rate N_{U }and removed at rate δ^{U}. They give birth to differentiated, uninfected Tcells, T. These in turn are removed at constant rate δ^{T }and become infected as they interact with the virus. Infected Tcells, I, die at rate δ^{I }and contribute to the budding of viral particles, V, that are cleared out at rate c. As soon as the X4 phenotype arise, the production of the TNF starts, proportional to the X4 concentration and contribute to the clearance of naïve Tcells, via the parameter.
In particular, Equation (7) describes the constant production of immature T cells by the thymus N_{U }and their differentiation into mature T cells at rate δ^{U}. If X4 viruses are present, upon the interaction with TNF, immature Tcells are cleared at fixed rate .
Equation (8) describes how uninfected mature T cells of strain i are produced at fixed rate δ^{U }by the pool of immature T cells. Those cells, upon the interaction with any strain of the virus, V_{k}, become infected at rate β_{k }= β ∀ k. The infectiousness parameter, β, is not constant over time, but depends on the interplay between R5 and X4 viruses. In particular, due to the increase in TNF abundance upon X4 emergence, the infectivity of R5 strains is reduced (β_{R5}(t) = β  k_{R5}F(t)), while the one of X4 viruses increases, with constant of proportionality k_{X4 }(β_{X4}(t) = β + k_{X4 }F(t)), mimicking the cell syncytium effect induced by the TNF molecule.
Equation (9) describes the infection of mature Tcells. Infected Tcells of strain k arise upon the interaction of a virus of strain k with any of the mature Tcell strains. The infected cells, in turn, are cleared out at a rate δ^{I}. When TNF is released, this value increases linearly with constant , δ^{I}(t) = δ^{I }+ F(t).
Equation (10) is close to that in the R5 quasispecies model, a part from different viral phenotypes being here considered.
Finally, in Equation (11), we model the dynamics of accumulation of TNF by assuming the increase in TNF level to be proportional, via the constant k_{F}, to the total concentration of X4 viruses present.
Investigating the mutational pathway from R5 to X4
In our model we represent the different phenotypes by using a linear strain space ordered in terms of phenotype similarity. Although we are aware of the several recent works on HIV mutational dynamics and phylogenetic assessments, we thought that a meaningful way to estimate the mutational pathways between R5 and X4 seen is to use phylogenetic inference on chemokine receptor families. The assessment of phylogenies using likelihood framework depends on the choice of an evolutionary model [26]. We computed the maximum likelihood (ML) analysis of the CRs data set under different models of evolution: [27], JTT [28], WAG [29]. We used these models considering the incorporation of the amino acid frequencies of the chemokine data sets, ('+F '), and the heterogeneity of the rates of evolution, implemented using a gamma distribution ('+Γ') [30,31].
Results
We have first extended Perelson's standard model to incorporate different antigen recognition abilities by the immune system and coexistence dynamics of different R5 strains of HIV virus. Our approach is a mean field one, i.e. we investigate the average quantities of these molecular species [3236]. Then, we have modified the model to describe features of the latent phase of the infection, i.e. the R5 to X4 switch and the hyperstimulation of the T cell precursors through the TNF. Finally we have analysed the results of the phylogenetic inference on chemokine receptors.
Amplification of R5 strains: mutation, coinfection and superinfection
Recent works have shown that HIV quasispecies may compete [22,37] and that persistence of the initial or ancestor quasispecies is a good indicator for disease progression [38]. Burch and Chao [39] have shown that the evolution of an RNA virus is determined by its mutational neighborhood. As the phenotype divergence among viral strains arises from differences in selection pressure, these differences may lead, for instance, to a higher infection rate. Since the competition is through the immune system response and given that the phase space of antigen recognition is not homogeneously covered [40], the HIV high mutation rate allows the quasispecies to find regions with weak immune response. This competition may lead to speciation of viral strains.
If we consider the model of the early phase of the infection, the evolution of T cells abundances in a scenario of quasispecies is shown in Figure 3a. Here we consider five viral phenotypes being differently recognized by the immune system. In particular, the recognition strength, corresponding to the selection pressure on viruses, is reported on the right yaxis in the inset of the figure as a dashed line. The temporal evolution of different T cells strains is characterized by the amplification of the T cells clones responding to viral infection, while asymptotic viral abundance is represented as solid stems in the inset. The asymptotic state of our model is a fixed point, thus the asymptotic distribution is insensitive of the initial conditions and the more abundant strains are that corresponding to the lower pressure exerted by the immune system (see inset of Figure 3b). However, one should consider that this asymptotic state may be reached after such a long time that it may be outside any practical scenario of the progression of a disease. The role of mutations in the transitory regime is quite particular. First of all, starting from the first inoculum at time t = 0 on the zeroth phenotype, mutations are necessary to populate the other strains of the virus, see also Figure 3. Moreover, in the presence of coupling among strains, due to competition or to a global constraints (the K parameter), the specific form of mutations does not play a fundamental role, see also Ref. [41].
Figure 3. Typical time evolution of the T cell abundance in the R5 model. Typical time evolution of the T cell abundance in the R5 model for shortterm (a) and longterm behavior (b). We set μ = 10^{5 }with N = 5. In the inset, the yaxis on the left reports the abundance of virus strains at time t = 50 in plot (a) and t = 2000 in plot (b); the yaxis on the right shows the interaction strength (dashed line) between T cells and virus phenotypes (xaxis).
Figures 4a and 4b show the results of short and long term viral coevolution after superinfection. We considered the first viral inoculum to happen at time t = 0, with the superinfection event occurring at time t = 20 days, when the immune response to the first inoculum has completed and the virus has established a chronic infection. After the second inoculum the model exhibits a short transient, followed by a slow mounting of the second infection. Due to the resulting low dynamics, the time needed by the second quasispecies to reach the same level of the other amounts to several months (Figure 4a) and represents another example of a slow relaxation toward a fixedpoint equilibrium. We may also account for the progression of the disease by considering a weak immune system, for example being characterized by a lower thymic activity, i.e. a lower value of λ, with respect to the previous scenario. In this case the strain corresponding to the second inoculum requires much longer time to reach the same abundances of the first strain.
Figure 4. Viral counts, V, during a superinfection scenario. Viral counts, V, during a superinfection scenario. We set N = 5 and no mutation is considered, thus μ = 0. (a) A slow mounting of the second viral infection (□), having time scale of several months, is observed. In (b) a compromised immune system is considered. The time for the second strain to reach the same abundance of the firstinfecting strain (◇) is greater than in (a).
We have also investigated how the coevolutionary and competitive dynamics of viral strains, mediated by the immune response, may lead to the formation of new viral quasispecies. We have considered a phenotypic space of 25 strains; this is an estimate of the number of different phenotypes (such those targeting different cell receptors or having different binding specifity for the same receptor) while the effective number of strains may be much larger. In Figure 5 the initial inoculum is at phenotype 15 (Figure 5a). In this simulation, since the immune system does not discriminate among similar phenotypes, there is an induced competition among neighboring strains. The result of this induced competition is the separation of the original quasispecies into two clusters (quasispeciation), Figure 5b. However the immune system response continues to change in time (Figures 5c–d), resulting into a complex coevolution with viral populations.
Figure 5. Speciation of virus quasispecies. Speciation of virus quasispecies and uninfected T cells dynamics after competitive superinfection at four different times: t = 0 (a), t = 4.5 (b), t = 5.25 (c) and t = 5.75 (d). Virus strain 15 is present at time t = 0, while strain 5 is inoculated at time t = 1. Mutation rate μ = 10^{4 }and nonuniform interaction strength as in Figure 3. The dashed line represents the abundances of T cells targeting each viral phenotype, represented as vertical stems.
R5 to X4 switch
Research into HIV dynamics has much to gain from investigating the evolution of chemokine coreceptor usage. Although CCR5 and CXCR4 are the major coreceptors used by HIV1 a number of chemokine receptors display coreceptor activities in vitro. Several other chemokine receptors, possibly not present on the T cell membrane, may act as targets. To date, a number of human receptors, specific for these chemokine subfamilies, have been described, though many receptors are still unassigned. Several viruses, for example EpsteinBarr, Cytomegalovirus, and Herpes Samiri, contain functional homologous to human CRs, an indication that such viruses may use these receptors to subvert the effects of host chemokines [42]. Cells different from CD4 and CD8, such as macrophages, express lower levels of CD4, CCR5, and CXCR4 on the cell surface compared with CD4+ T cells [4345], and low levels of these receptors expressed on macaque macrophages can restrict infection of some nonMtropic R5 HIV1 and X4 simian immunodeficiency virus (SIV) strains [46,47].
We studied the coevolutive dynamics leading to X4 strain appearance by successive mutations of the ancestor R5 strain. The stimulated production of TNF regulate the interactions between immune response and the virus and between the different strains of HIV virus. The results of these interactions are a decline in Tcells level, leading to the AIDS phase of the disease, and the decline in levels of viruses using the R5 coreceptor. In Figure 6 we report the result of the numerical simulations of the infection. The early stages of the infection are characterized by R5 strains. As time goes on, mutation accumulates that finally lead to the X4 phenotype appearance. After a short transitory regimes, T cells abundance starts declining. We have finally studied the influence of switching coreceptor usage in superinfection dynamics. In Figure 7 we show Tcells dynamics for different times of the superinfection event. We may observe that if the superinfection occurs after the appearance of the X4, the new R5 strain does not have any effect on Tcells behavior. On the other hand is worth noting that if the new R5 inoculum takes place before the X4 appearance, this may speed up the switching to the X4 phenotype if the new strain is mutationally close to the X4.
Figure 6. R5 to X4 switch. Time evolution of the concentrations of uninfected Tcells (◇) and viruses (+), during R5 to X4 switch, occurring at time t ≈ 900. After the appearance of the X4 phenotype a continuous slow decline in CD4+ Tcells level leads to AIDS phase (CD4 counts below 200 cells/ml). We set μ ≈ 0.001.
Figure 7. CD4+ Tcells concentration during HIV1 superinfection. CD4+ Tcells concentration during HIV1 superinfection by a R5 viral strain. Different signs represent: evolution without superinfection, (◇); superinfection occurring at time t = 100 and 400, (+) and (□), respectively. For a superinfection event occurring after the R5 to X4 switching the dynamics is qualitatively the same as for a single infection, (◇). If the second delayed infection occurs before the R5 to X4 switching, the time of appearance of X4 viruses may be shorter, when the superinfecting strain is closer to the X4 phenotypes, (+, □). Parameters as in Figure 6.
Investigating the mutational pathway between R5 and X4
Fundamental to the evolutionary approach is the representation of the evolution of sequences along lineages of evolutionary trees, as these trees describe the complex patterns of dependence amongst sequences that are caused by their common ancestry [29,48,49].
The ML tree, obtained using the JTT+F+Γ model of evolution, is shown in Figure 8. The topology clearly shows that the CCR family is not homogeneous: CCR6, CCR7, CCR9 and CCR10 are separated from the other CCRs; in particular, CCR10 clusters with CXCRs; CXCR4 and CXCR6 do not cluster with the CXCRs. The tree shows that there are many mutational steps between CCR5 and CXCR4. The phylogeny suggests that the mutations that allow the virus env to cover a wide phenotypic distance from R5 to X4, may also lead to visit other receptors.
Figure 8. Maximum likelihood phylogeny. The maximum likelihood phylogeny under the JTT+F+Γ model of evolution for the set of human and mouse (mouse sequences are labelled with "M") chemokine receptors. We have considered only the external loop regions. The scale bar refers to the branch lengths, measured in expected numbers of amino acid replacements per site.
Conclusion
The worldwide presence of several strains of the HIV virus and their often simultaneous presence within a patient, due to the increased frequency of multiple infections, are the remarkable features of HIV pandemia. For example, HIV1 exists as several groups, subdivided into a growing number of subtypes which are slightly predominant in different geographical regions [50]. HIV1infected CD4+ T cells isolated from the spleens of two individuals were recently shown to harbor anywhere between one and eight proviruses, with an average of three to four proviruses per cell [51]. Mutations, recombination and selection pressure cause the appearance of quasispecies [52].
We first presented a model of the withinpatience persistence of HIV quasispecies, by extending to multiple strain the Perelson's standard model [13]. This approach allows to incorporate coevolutive and competitive dynamics resulting from the different strains of HIV virus and different antigen recognition abilities by the immune system considered. Our model shows that the time evolution of the competition between quasispecies is slow and has time scales of several months.
Phylogenetic inference of chemokine receptors shows that there are several mutational patterns linking CCR5 to several receptors that have the same branch length of that from CCR5 to CXCR4. The evolutive dynamics towards the selection of X4 requires a preferential binding to a CXCR4 receptor. Recent works show that TNF is a prognostic marker for the progression of HIV disease [5,7] and has role in regulating the interactions between the different strains of HIV virus. The second model we have introduced shows that keeping low the concentration of TNF, both the depletion of Tcells precursors repertoire and the X4 dominance over R5 strains slow down.
The observed coevolutionary dynamics of virus and immune response opens the way to the challenging possibility of the introduction or modulation of a quasispecies to be used in therapy against an already present aggressive strain, as experimented by Snell and colleagues [53,54]. Different drug treatments can alter the population of quasispecies. Will R5 blocking drugs cause HIV to start using X4? And will that be worse than letting the R5using virus stay around along at its own, slower, but no less dangerous activity? The presence of a large number of R5 will increase the mutational spectra in R5 strains (late R5) and the probability of getting closer to the binding specificities of other chemokine receptors. Therefore, our model suggests that a therapy such the HAART to decrease the HIV load should be started the sooner the better. The large effect of TNF on T cells dynamics described by our model, suggests the benefit of a TNF buffering therapy. It is known that the dynamics of TNF is related to the dynamics of TNFrelated apoptosisinducing ligand (TRAIL). In this model we consider constant the concentration of TRAIL [7] and we do not consider many other important players such as Rantes.
Our models represent also a general framework to investigate intermittency or switching dominance of strains and the arising of new dominant strains during different phases of therapy; how superinfection will evolve in case of replacement of drugresistant virus with a drugsensitive virus and acquisition of highly divergent viruses of different strains; to investigate whether antiviral treatment may increase susceptibility to superinfection by decreasing antigen load.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
All authors conceived this project and drafted the manuscript. LS carried out most of simulations. All authors have interpreted and discussed the results, read and approved the final manuscript.
Acknowledgements
We thank Stefano Ruffo for insightful suggestions. P.L. would like to thank EU 6th Framework Bioinfogrid project. L.S. research at the Computer Laboratory of University of Cambridge is supported by an EMBO fellowship.
This article has been published as part of BMC Evolutionary Biology Volume 7 Supplement 2, 2007: Second Congress of Italian Evolutionary Biologists (First Congress of the Italian Society for Evolutionary Biology). The full contents of the supplement are available online at http://www.biomedcentral.com/14712148/7?issue=S2
References

Karlsson I, Antonsson L, Shi Y, Oberg M, Karlsson A, Albert J, Olde B, Owman C, Jansson M, Fenyo EM: Coevolution of RANTES sensitivity and mode of CCR5 receptor use by human immunodeficiency virus type 1 of the R5 phenotype.
J Virol 2004, 78:1180711815. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Gorry PR, Churchill M, Crowe SM, Cunningham AL, Gabuzda D: Pathogenesis of macrophage tropic HIV.
Curr HIV Res 2005, 3:5360. PubMed Abstract

Clevestig P, Maljkovic I, Casper C, Carlenor E, Lindgren S, Naver L, Bohlin AB, Fenyo EM, Leitner T, Ehrnst A: The X4 phenotype of HIV type 1 evolves from R5 in two children of mothers, carrying X4, and is not linked to transmission.
AIDS Res Hum Retroviruses 2005, 21:371378. PubMed Abstract  Publisher Full Text

Gray L, Sterjovski J, Churchill M, Ellery P, Nasr N, Lewin SR, Crowe SM, Wesselingh SL, Cunningham AL, Gorry PR: Uncoupling coreceptor usage of human immunodeficiency virus type 1 (HIV1) from macrophage tropism reveals biological properties of CCR5restricted HIV1 isolates from patients with acquired immunodeficiency syndrome.
Virology 2005, 337:38498. PubMed Abstract  Publisher Full Text

Herbeuval JP, Hardy AW, Boasso A, Anderson SA, Dolan MJ, Dy M, Shearer GM: Regulation of TNFrelated apoptosisinducing ligand on primary CD4+ T cells by HIV1: Role of type I IFNproducing plasmacytoid dendritic cells.
Proc Nat Acad Sci USA 2005, 102:1397413979. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Pavlakis GN, Valentin A, Morrow M, Yarchoan R: Differential effects of TNF on HIV1 expression: R5 inhibition and implications for viral evolution.
Int Conf AIDS 2004., 15
(abstract no. MoOrA1048).

Herbeuval JP, Grivel JC, Boasso A, Hardy AW, Chougnet C, Dolan MJ, Yagita H, Lifson JD, Shearer GM: CD4+ Tcell death induced by infectious and noninfectious HIV1: role of type 1 interferondependent, TRAIL/DR5mediated apoptosis.
Blood 2005, 106:35243531. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Ho D, Neumann AU, Perelson AS, Chen W, Leonard JM, Markowitz M: Rapid turnover of plasma virions and CD4 lymphocytes in HIV1 infection.
Nature 1995, 373:123126. PubMed Abstract  Publisher Full Text

Chao L, Davenport MP, Forrest S, Perelson AS: A stochastic model of cytotoxic T cell responses.
J Theor Biol 2004, 228:227240. PubMed Abstract  Publisher Full Text

Celada F, Seiden PE: Affinity maturation and hypermutation in a simulation of the humoral immune response.
Eur J Immunol 1996, 26:13501358. PubMed Abstract

De Boer RJ, Perelson AS: Towards a general function describing Tcell proliferation.
J Theor Biol 1995, 175:567576. PubMed Abstract  Publisher Full Text

Wodarz D, Nowak MA: HIV dynamics and evolution.
BioEssays 2002, 24:11781187. PubMed Abstract  Publisher Full Text

Perelson AS, Neumann AU, Markowitz M, Leonard JMD, Ho D: HIV1 Dynamics in Vivo: Virion Clearance Rate, Infected Cell LifeSpan, and Viral Generation Time.
Science 1996, 271:15821586. PubMed Abstract  Publisher Full Text

Wiegel FW, Perelson AS: Some scaling principles for the immune system.
Immunol Cell Biol 2004, 82:127131. PubMed Abstract  Publisher Full Text

Wei X, Ghosh SK, Taylor ME, Johnson VA, Emini EA, Deutsch P, Lifson JD, Bonhoeffer S, Nowak MA, Hahn BH, Saag MS, Shaw GM: Viral dynamics in human immunodeficiency virus type 1 infection.
Nature 1995, 373:117122. PubMed Abstract  Publisher Full Text

Biebricher CK, Eigen M: The error threshold.
Virus Res 2005, 107:11727. PubMed Abstract  Publisher Full Text

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 Interferonalpha Therapy.
Science 1998, 282:103107. PubMed Abstract  Publisher Full Text

Eigen M, Schuster P: The hypercycle.
Naturwissenschaften 1977, 64:541565. PubMed Abstract

Ribeiro MR, Hazenberg DM, Perelson AS, Davenport PM: Naive and Memory Cell Turnover as Drivers of CCR5toCXCR4 Tropism Switch in Human Immunodeficiency Virus Type 1: Implications for Therapy.
J Virol 2006, 80:802809. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

De Boer RJ, Segel LA, Perelson AS: Pattern formation in one and twodimensional shapespace models of the immune system.
J Theor Biol 1992, 155:295333. PubMed Abstract

Eigen M, Schuster P: The Hypercycle: A Principle of Natural SelfOrganization. Berlin: Springer; 1979.

Vignuzzi M, Stone JK, Arnold JJ, Cameron CE, Andino R: Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population.
Nature 2006, 439:344348. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Gog JR, Grenfell BT: Dynamics and selection of manystrain pathogens.
Proc Nat Acad Sci 2002, 99:1720917214. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Lin J, Andreasen V, Casagrandi R, Levin SA: Traveling waves in a model of influenza A drift.
J Theor Biol 2003, 222:437445. PubMed Abstract  Publisher Full Text

Bagnoli F, Lió P, Sguanci L: Modeling viral coevolution: HIV multiclonal persistence and competition dynamics.

Whelan S, Liò P, Goldman NP: Molecular phylogenetics: stateoftheart methods for looking into the past.
Trends Genet 2001, 17:262272. PubMed Abstract  Publisher Full Text

Dayhoff MO, Schwartz RM, Orcutt BC: A model of evolutionary change in proteins. In Atlas of Protein Sequence and Structure. Volume 5. Edited by Dayhoff MO. National Biomedical Research Foundation; 1978::345352.

Jones DT, Taylor WR, Thornton JM: A mutation data matrix for transmembrane proteinsm.
FEBS Lett 1994, 339:26975. PubMed Abstract  Publisher Full Text

Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum likelihood approach.
Mol Biol Evol 2001, 18:691699. PubMed Abstract  Publisher Full Text

Yang Z: Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods.
J Mol Evol 1994, 39:306314. PubMed Abstract

Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood.
CABIOS 1997, 13:555556. PubMed Abstract

Ramos A, Hu DJ, Nguyen L, Phan KO, Vanichseni S, Promadej N, Choopanya K, Callahan M, Young NL, McNicholl J, Mastro TD, Folks TM, Subbarao S: Intersubtype human immunodeficiency virus type 1 superinfection following seroconversion to primary infection in two injection drug users.
J Virol 2002, 76:74447452. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Levy JA: Is HIV superinfection worrisome?
Lancet 2003, 361:9899. PubMed Abstract  Publisher Full Text

ten Haaft P, Verschoor EJ, Verstrepen B, Niphuis H, Dubbes R, Koornstra W, Bogers W, Rosenwirth B, Heeney JL: Readily acquired secondary infections of human and simian immunodeficiency viruses following single intravenous exposure in nonhuman primates.
J Gen Virol 2004, 85:37353745. PubMed Abstract  Publisher Full Text

Altfeld M, Allen TM, Yu XG, Johnston MN, Agrawal D, Korber BT, Montefiori DC, O'Connor DH, Davis BT, Lee PK, Maier E, Harlow J, Goulder PJR, Brander C, Rosenberg ESB, Walker D: Evolution and transmission of stable CTL escape mutations in HIV infection.

Gerhardt M, Mloka D, Tovanabutra S, SandersBuell E, Hoffmann O, Maboko L, Mmbando D, Birx DL, McCutchan FE, Hoelscher M: Indepth, longitudinal analysis of viral quasispecies from an individual triply infected with latestage human immunodeficiency virus type 1, using a multiple PCR primer approach.
J Virol 2005, 79:82498261. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Costa L, Janini M, Pinto M, Diaz R: Quasispecies analyses of a familial cluster infected with multiple HIV1 subtypes.
Int Conf AIDS 2000., 13
(abstract no. MoPpA1005).

Bello G, Casado C, Sandonis V, AlonsoNieto M, Vicario JL, García S, Hernando V, Rodríguez C, Delomero JR, LópezGalíndez C: A subset of human immunodeficiency virus type 1 longterm nonprogressors is characterized by the unique presence of ancestral sequences in the viral population.
J Gen Virol 2005, 86:355364. PubMed Abstract  Publisher Full Text

Burch CL, Chao L: Evolvability of an RNA virus determined by its mutational neighborhood.
Nature 2000, 406:625628. PubMed Abstract  Publisher Full Text

DeBoer RJ, Segel LA, Perelson AS: Pattern formation in one and twodimensional shapespace models of the immune system.
J Theor Biol 1992, 155:295333. PubMed Abstract

Bagnoli F, Bezzi M: Speciation as Pattern Formation by Competition in a Smooth Fitness Landscape.

Murphy PM: Viral exploitation and subversion of the immune system through chemokine mimicry.
Nature Immunology 2001, 2:116122. PubMed Abstract  Publisher Full Text

Lewin SR, Sonza S, Irving LB, McDonald CF, Mills J, Crowe SM: Surface CD4 is critical to in vitro HIV infection of human alveolar macrophages.
AIDS Res Hum Retroviruses 1996, 12:877883. PubMed Abstract

Ometto L, Zanchetta M, Cabrelle A, Esposito G, Mainardi M, ChiecoBianchi L, De Rossi L: Restriction of HIV type 1 infection in macrophages heterozygous for a deletion in the CCchemokine receptor 5 gene.
AIDS Res Hum Retroviruses 1999, 15:14411452. PubMed Abstract  Publisher Full Text

Wang J, Crawford K, Yuan M, Wang H, Gorry PR, Gabuzda D: Regulation of CC chemokine receptor 5 and CD4 expression and human immunodeficiency virus type 1 replication in human macrophages and microglia by T helper type 2 cytokines.
J Infect Dis 2002, 185:885897. PubMed Abstract  Publisher Full Text

Bannert N, Schenten DS, Craig J, Sodroski J: The level of CD4 expression limits infection of primary rhesus monkey macrophages by a Ttropic simian immunodeficiency virus and macrophagetropic human immunodeficiency virusesm.
J Virol 2000, 74:1098410993. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Mori K, Rosenzweig M, Desrosiers RC: Mechanisms for adaptation of simian immunodeficiency virus to replication in alveolar macrophages.
J Virol 2000, 74:1085210859. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Shields DC: Gene conversions among chemokine genes.
Gene 2000, 246:239245. PubMed Abstract  Publisher Full Text

Goh CS, Bogan AA, Joachimiak M, Walther D, Cohen FE: Coevolution of proteins with their interaction partners.
J Mol Biol 2000, 299:283293. PubMed Abstract  Publisher Full Text

Yang OO, Daar ES, Jamieson BD, Balamurugan A, Smith DM, Pitt JA, Petropoulos CJ, Richman DD, Little SJ, LeighBrown AJ: Human immunodeficiency virus type 1 clade B superinfection: evidence for differential immune containment of distinct clade B strains.
J Virol 2005, 79:860868. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Dixit NM, Perelson AS: HIV dynamics with multiple infections of target cells.
Proc Natl Acad Sci USA 2005, 102:8198203. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Eigen M, McCaskill J, Schuster P: The Molecular QuasiSpecies.

Schnell J, Johnson E, Buonocore L, Rose JK: Construction of a novel virus that targets HIV1infected cells and controls HIV1 infection.
Cell 1997, 90:849857. PubMed Abstract  Publisher Full Text

Revilla T, GarcýaRamos G: Fighting a virus with a virus: a dynamic model for HIV1 therapy.