A theoretical model of the West Nile Virus survival data.

BACKGROUND: In this work, we develop a theoretical model that explains the survival data in West Nile Virus infection.
RESULTS: We build a model based on three cell populations in an infected host; the collateral damage cells, the infected dividing cell, and the infected non-dividing cells. T cell-mediated lysis of each of these populations is dependent on the level of MHC-1 upregulation, which is different in the two infected cell populations, interferon-gamma and free virus levels.
CONCLUSIONS: The model allows us to plot a measure of host health versus time for a range of initial viral doses and from that infer the dependence of minimal health versus viral dose. This inferred functional relationship between the minimal host health and viral dose is very similar to the data that has been collected for WNV survival curves under experimental conditions.

Background

Viruses in the family Flaviviridae are single-stranded, plus-sense RNA viruses, principally transmitted by mosquitoes and ticks. They are either viscerotropic, causing diseases such as dengue and yellow fever, or neurotropic, causing central nervous system disease, like West Nile virus encephalitis, Japanese encephalitis or Saint Louis encephalitis, all of which may be fatal. These viruses are found worldwide and the diseases they cause pose a significant public health burden. A careful discussion of the many facets of flavivirus infections can be found in [1], and we therefore confine our introductory remarks to the features salient to this paper.

Virus-infected cells process virus proteins into peptide fragments in the proteasome and these are bound to class I major histocompatibility complex molecules (MHC-I) in the endoplasmic reticulum. With the transport of MHC molecules to the cell surface, viral peptides displayed in the context of MHC on the infected cell can be recognized by virus-specific cytotoxic T lymphocytes (CTL). Such CTL kill infected cells by a variety of lytic effector molecules prior to the release of mature virus progeny from infected cells, thereby progressively reducing virus numbers and ultimately eradicating virus from surviving hosts (reviewed in [2]).

Paradoxically, flaviviruses such as West Nile (WNV) and others, directly induce increased expression of MHC-I, as well as MHC-II, and several adhesion molecules involved in immune recognition by CTL. This increased expression results in a marked increase in the efficiency of recognition and killing of infected cells by WNV-specific CTL ([3, 4]), because although the affinity of individual T cell receptor (TcR)-MHC-virus peptide interaction is unchanged, the multiple intermolecular interactions increases the avidity of interaction of virus-specific CTL with the infected cell. This increased avidity also enables the functional interaction with MHC , infected cells by CTL clones of low MHC-virus peptide affinity, i.e., clones previously below the recognition threshold. Some of these low-affinity CTL clones are likely to be self-reactive [4] or even able to recognize MHC without peptide specificity [5]. Thus, the increased avidity brought about by high MHC expression enables low-affinity, self-reactive clones, not normally involved in anti-viral immune responses, to lyse both infected and uninfected target cells [6]. Adhesion molecules such as ICAM-1, as non-specific accessory molecules upregulated by WNV infection, also increase the avidity of CTL-target cell interactions, further lowering the affinity threshold for T cell recognition and target cell lysis [7]. In addition, interferon- γ (IFN- γ), released by CTL on recognition of their cognate ligand, strongly increases MHC and ICAM-1 expression on neighboring target cells, further contributing to the progressive increase in avidity of interaction between CTL and target cells [8]. In this context, the stage of the cell cycle in which a cell is infected is also important; cells infected by WNV in G 0 (resting) increase MHC-I expression by 6-10-fold, compared to a 2-3-fold increase observed in cells infected during the cell cycle (G 1, S, G 2+M) [3]. WNV-infected G 0 cells are approximately 10-fold more susceptible to CTL lysis than infected cycling cells exposed to the same CTL [3]. Thus, WNV-infected cycling cells are less easily recognized, while the avidity of interaction between CTL and infected G 0 cells is significantly enhanced by the higher levels of MHC-I and ICAM-1. Notably, WNV replicates significantly better in the poorly recognized cycling cells than in G 0 cells. In vivo, most cells are in G 0, presumably presenting an easy CTL target once infected, but a small population of productively infected cycling cells maintaining a low immunological profile, could substantially increase the probability of virus transmission to the next host [4]. As indicated above, the release of IFN- γ upon target cell recognition by CTL would also increase MHC and ICAM-1 expression on uninfected cells in the vicinity of virus- infected cells, making high MHC-I-expressing (uninfected) cells a potential target for lysis by low affinity virus-specific and/or self-reactive (i.e., cross-reactive) CTL clones. The collateral destruction of uninfected cells by low-affinity clones of this kind would cause substantial additional damage to the brain, with corresponding increases in morbidity and mortality [9]. We have developed a model of the collateral damage caused by a West Nile Virus infection, which is supported by simulation results [10]. A discussion of the underlying simulation code can be found in [11]. We have also used these simulation results to explain the unusual ragged survival data seen in West Nile Virus infections [12]. The previous work, in focusing on very low level details, based on first principles, could be regarded as a micro model. Here we develop a macro model, based on a much higher level approach and ideas. Each model has its pros and cons but we believe each brings complementary insights into a very complicated problem. Our focus here is therefore on developing a macro level theoretical model of how the survival of a host depends on the level of initial viral dose. We believe this approach provides an abstract focus which can also be directed towards more general models of immunopathology and we will briefly mention those connections at the end. These more general models of autoimmune response are the focus of additional work we are doing.

WNV survival data

In a survival experiment, a population of animals such as mice are infected with a common amount of virus and the number of animals that die are counted. This experiment is repeated for a range of virus. With most viruses, the survival experiment gives a classical dose-response curve which progressively and smoothly decays down to a survival of 0 animals at high virus. However, WNV has a peculiar survival curve as shown in [12] but we can easily simulate what such a curve looks like as we have done in Fig. 1. The purpose of this paper is to find a theoretical model that explains this data. The derivations here use standard ideas from advanced calculus and differential equations. The necessary background can be found in any nonlinear modeling text such as [13]. We assume we have a large population of cells which consists of cells which are dividing and are infected, , cells which are not dividing but are infected, , non infected cells, and non infected cells which will be removed due to auto immune action which we call , for collateral damage. We assume all of these variables depend on the three parameters, , the IFN- γ signal; the MHCI upregulation factor and the free virus level . The approach we take here is a common one in such a nonlinear interaction environment. An example of how this is used to develop a model of diabetes detection can be found in [13].

Fig. 1

Simulated survival curve

The model

We assume the dynamics here are

There are then three nonlinear interaction functions F 1, F 2 and F 3 because we know , and depend on each other’s levels in very complicated ways. Usually, we assume the initial dose gives rise to some fraction of infected cells which will differ in both the dividing and nondividing cells. In our previous WNV simulations of collateral damage [10], we assume p 0=1.0e−3 of the initial dose gives rise to infected cells. This gives the number of infected cells to be p 0 which is split into p 1 p 0 infected non dividing cells - i.e. in , and p 2 p 0 dividing infected cells, i.e. in , where p 1+p 2=1. Typically, we let p 1=0.99 and p 2=0.01 but other choices could be made. Thus, the total amount of virus that goes into infected cells is p 0 and the amount of free virus is therefore (1−p 0) . Thus, we could expect =0, =p 2 p 0 and = =p 1 p 0 . However, we will explicitly assume we are starting from a point of equilibrium prior to the administration of the viral dose . We could assume there is always some level of collateral damage, in a host, but we will not do that. We will therefore assume , and have achieved these values =0, =0 and =0 right before the moment of infection. Hence, we don’t expect there to be initial contribution to ′(0), ′(0) and ′(0); i.e. F 1( , , )=0, F 2( , , )=0 and F 3( , , )=0. We are interested in the deviation of , and from their optimal values , and , so let =− , =− and =− . We can then write = +, = + and = + The model can then be rewritten as

or

Next, we do a standard tangent plane approximation on the nonlinear dynamics functions F 1, F 2 and F 3 to derive approximation dynamics.

Linearization details

To develop a useful approximation to the dynamics, we use an appropriate tangent plane approximation on the nonlinear dynamics functions F 1, F 2 and F 3. The mathematics behind this approximation come from multivariate calculus and can easily be reviewed if required. The standard tangent plane expansions is as follows.

But the terms F 1( , , )=0, F 2( , , )=0 and F 3( , , )=0. so we can simplify to

It seems reasonable to assume that since we are so close to ordinary operating conditions, the errors , and will be negligible. Also, to save space, we let () denote that we are evaluating the partial derivatives at the point ( , , ). Thus our model approximation is

and our corresponding nonlinear dynamics approximation is

We can write this approximation In matrix - vector form also:

where we now use a standard subscript scheme to indicate the partials. Now let’s add IFN- γ, upregulation and virus level.

The CDN IFN- γ, upregulation and virus model

We can think each variable , and as depending on the interferon level , the upregulation level and the free virus level . Thus, we have

We assume the dynamics here are then

As before assume , and have achieved the same optimal values =0, =0 and =0 prior to the moment of infection with virus dose . These correspond to the starting values of prior to infection for base values , and . Initially, we don’t expect IFN- γ signals so =0. Eventually, we do expect some level of upregulation due to this initial dose and from experimental data, we expect the upregulation to be proportional to the level of the dose ; we will assume this is a simple scaling factor, i.e. =q 1 for some suitable parameter q 1. Also, once the viral dose is administered, we would expect some fraction of it to remain as free virus which as discussed in [10] is modeled as =(1−p 0) . We will deal with these initial values after infection later in the Section “The IFN- ”. But now, we think of all the initial values as zero; i.e. =0, =0 and =0. We still don’t expect to have any contribution to ′(0), ′(0) and ′(0); i.e. H 1( , , )=0, H 2( , , )=0 and H 3( , , )=0. We are interested in the deviation of , and from their optimal values , and due to the changes , and from the base IFN- γ, upregulation and virus values. So let =− , =− and =− . We can then write = +, = + and = + The model can then be rewritten as

which, as usual, implies

Next, we again perform a tangent plane approximation on the nonlinear dynamics functions H 1, H 2 and H 3 just as we did in Section “Linearization details” for the F 1, F 2 and F 3 we use for the model.

Once again, it seems reasonable to assume that since we are so close to ordinary operating conditions, the tangent plane errors are negligible. Letting () denote that we are evaluating the partial derivatives at the point ( , , ), we find

The corresponding nonlinear dynamics approximation in matrix - vector form is then

where we now use a standard subscript scheme to indicate the partials. We therefore find the nonlinear dynamics approximation is

where we now use a standard subscript scheme to indicate the partials.

The algebraic signs of the linearization matrix

If we hold everything constant except which we increase to +δ , what happens? Increasing the IFN- γ level should increase collateral damage. Hence, . What about the number of dividing cells that are infected? We do not think this should have an effect; hence, too. A similar argument suggests as well. Thus, the coefficient matrix above which we call so far looks like

Now hold everything constant except the upregulation and increase to +δ . What happens? Increasing the upregulation level should increase collateral damage. Hence, . Now infected cells are lysed once their IFN- γ upregulation hits a certain level and the increase in upregulation that makes the infected cells more visible certainly effects this as it will take fewer additional IFN- γ signals to trigger lysis. In general, there is a lag of say and time steps before this happens. With the increase δ , we expect and to both go down. Hence and decrease and we have and for each type of infected cell. Thus, the coefficient matrix looks like

Next hold everything constant except the virus level and increase to +δ . What happens? Let’s think of the virus increase δ as giving rise to an increase in the amount of virus stored inside a dividing cell or a non dividing cell. Now if the amount of virus in the cell goes up, that means when the cell is lysed, there is more virus available to infect cells which means more cells will be infected in later times. An increase in virus means an increase in collateral damage in general, so . At a given time then, A is the virus level. We can write A=A +A +A where A is the free virus, A is the virus inside the dividing cells and A is the virus inside the nondividing cells. So if A goes up, we expect the amount of A and A to remain constant. Hence, if A goes up, A goes down.

The increase in A would mean more virus is released on lysis and hence an increase in . However, that increase in A is a concomitant decrease in . So we should have and . Thus, the coefficient matrix now looks like

Or letting , , , , the coefficient matrix now looks like

Thus, we have the changes in collateral damage and infection levels

Now we need to estimate , and .

The IFN- γ, upregulation and free virus model

The amount of interferon level , the upregulation level and the free virus level depend on the initial amount of virus applied when in the equilibrium state; i.e. this is the amount that causes the initial infection. This is . We assume the dynamics here are

In the model of “The CDN IFN- ” section, we assumed , and depended on the perturbations of , and from a zero state. Now, we want to model the , and deviations from a base state , and which is not zero. As previously discussed, we expect =(1−p 0) , the initial IFN- γ level =0 and the initial upregulation level =q 1 . Let the deviations from these equilibrium values be given by =− , =− and =− . The model can then be rewritten as

or

The approximation to the , and model is handled in a way that is very similar to the previous two expansions which were explained in some detail in “Linearization details” section.

The linearization again involves a tangent plane approximation to nonlinear dynamics functions, which are now the functions G 1, G 2 and G 3. It seems reasonable to assume that since we are so close to ordinary operating conditions, the tangent plane errors as usual will be negligible. Thus, using to denote the partials evaluated at the base point the model approximation, we have

Thus our dynamics approximation in matrix - vector form is

The dynamics approximation is then given by

The analysis of the signs of these partials is next. This is similar to what we did for the previous model. If we hold everything constant except which we increase to +δ , what happens? Increasing the IFN- γ level should increase IFN- γ. Hence, . What about the upregulation level? We do not think this should have an effect; hence, too. A similar argument suggests as well. Thus, the coefficient matrix above which we call so far looks like

Now hold everything constant except the upregulation and increase to +δ . What happens? Increasing the upregulation level should increase the IFN- γ level. Hence, . If upregulation increases, the ′ goes up, hence . What about the level of virus? If upregulation increases, the increase in upregulation makes infected cells more visible as it will take fewer additional IFN- γ signals to trigger lysis. In general, there is a lag of say and time steps before this happens. With the increase δ , we expect and to both go down. Hence, free virus should increase; i.e. . Thus, the coefficient matrix looks like

Now hold everything constant except the virus level and increase to +δ . What happens? Let’s think of the virus increase δ as giving rise to an increase in the amount of virus stored inside a dividing cell or a non dividing cell. Now if the amount of virus in the cell goes up, that means when the cell is lysed, there is more virus available to infect cells which means more cells will be infected in later times. An increase in virus should not effect IFN- γ levels, so . An increase in virus must imply Now reason as we did earlier: at a given time then, A is the free virus level. We can write A=A +A +A where A is the free virus, A is the virus inside the dividing cells and A is the virus inside the nondividing cells. So if A goes up and the total amount of virus at this time is constant, we expect A and A to go down. If the amount of virus inside the cells decreases, we expect the upregulation to decrease. So . Thus, the coefficient matrix now looks like

Or letting , , , , we have the coefficient matrix now looks like

Oscillations in upregulation and free virus

We now use standard results from the theory of linear differential equation systems. This analysis relies on the eigenvalues of our model systems. The eigenvalues of this linearized system are found by solving the det(λ −)=0. Thus, for the coefficient matrix Φ, we have

This gives

The eigenvalues of the two by two submatrix are the most interesting. Consider the determinant

This has the complex roots (here the symbol is defined to be ).

Hence, we can get complex roots if

or and . Since an eigenvalue equation of the form

has roots

we have complex roots if (a−d)2−4b c<0. In particular, if a=d=α and b=β and c=β, then we have −4β 2<0 and the roots are complex. We can have complex roots for other choices of a, b, c and which implies conditions on the partials which is something we can explore, but we will not do that here. For our purposes, we will examine closely what happens when and The general solution to the system with the complex eigenvalues we have been discussing is this. The eigenvalues are and the complex conjugate pair α±β where and . The eigenvectors here are simple

We can then solve for and to find

Letting , we find

where δ is defined as tan−1(b/a). The full solution is then

where A, R, , β and δ determine a given model.

Here, we have =q 1 and =(1−p 0) . Hence, we roughly know at the time of the infection

Taking a ratio, we find

Hence, and

Finally, recall we have and ; thus, the oscillatory solutions for upregulation and free virus are

It seems unreasonable to us that the phase shift δ should be a constant; i.e. independent of . After all, the reasoning above is approximate and we should not think of this as actually fixed. So we will assume that all critical parameters here are proportional to . Our rough calculation showed us , so it seems reasonable that R is proportional to in general. Therefore, we now assume

for a new parameters r 1, r 2, r 3 and r 4. This leads to our estimate of the dependence of the upregulation and free virus on the initial dose : we have

A health model

Roughly speaking, if the total number of cells is , the number of healthy cells can be approximated by

We know

and so we are looking at deviations from the base values =0, =q 1 and =(1−p 0) . It follows we have

As discussed earlier, we have initially, =0, =p 2 p 0 and =p 2 p 0 . So we have

Thus, we have

Now collect all the terms involving and set that coefficient to Λ for convenience. Making this replacement, we have

This leads to the simplification

Now we have to compute these integrated transient values. We label them as for the transient integration; for the transient integration; and for the transient integration. We then have

The integration is straightforward

however, the integration are more complicated.

The calculation

To evaluate this term, we use integration by parts. Now we have to compute the integrated transient values required to find the health estimate. We have labeled the integrations as for the transient integration; and for the transient integration. We then have

Integration details

First, let’s calculate . To evaluate this term, we use integration by parts. We find

and so we can find as follows:

We can rewrite this is a much better form using our assumptions. First, rewrite as

Now we find . Another standard integration by parts gives

and so we can find also:

Model results

We know α, β, δ and R are really dependent of . For convenience of exposition, we drop the superscript in our calculations below

Finally, let’s define two new parameters, θ 1 and θ 2 as and . Using the above, we can rewrite (t) as

Using a standard reference triangle for the phase angle θ 2, we see and . We can then rewrite (t) again as

and using standard trigonometric identities, we then have

Recall, as showing in as shown in “Integration details”, we know

Another standard integration by parts similar to what was done in “Integration details” allows us to find . We note the same comment on the dependence of R, α, β and δ on holds still. Now using these values and the terms Q 1 and Q 2, we we can rewrite (t) as follows:

Now using the simplifications we obtained for α and β in terms of r 2 and r 3, we can rewrite this complicated expression as

Next, using the phase shift θ 2, we have

This then leads to our final form

Building the health model

Recall the health model is

Let and . Then we have

Now plug what we have found for our integrations. We have

Then we can rewrite as

Now put the together. We find

Let’s simplify some more using another phase shift. Define the phase angle ; then, we can rewrite the health like this.

This can be recast as

Next, we can combine the ratio into the new parameter and rewrite as to give

Finally, let . Then, we have the last form of the health estimate:

We could also assume the terms and are proportional to . We would model this by implying and . We then find

These parameters depend in complex ways on the initial virus dose and it is very difficult to tease out the details. The only data we have to help fit this model is the survival data, so the next step of seeing how this model of health gives rise to the experimental survival data requires additional analysis.

Results and discussions

We are going to graph Eq. 4 for various values of the parameters. First, we set the parameters as showin in Fig. 2.

Fig. 2

Set the parameters

Then, we define some auxiliary functions in Fig. 3.

Fig. 3

Define needed functions

Then, we generate a plot of health versus time for a range of choices of , here represented by S0 as shown in Fig. 4. For each plot, we find the minimum value and store it in the variable Min. We can then plot Min versus initial viral dose S0 and see if the plot looks like a traditional WNV survival curve. We then place all of this code into a MatLab function and use it in the usual way to generate the plots. For these choice of parameters, we can plot the minimal health curve which is shown in Fig. 5. Another choice of parameters, leads to an even better minimum health curve which captures the essence of the WNV survival plot. The parameter choices for this run are shown in Fig. 6. For these parameters values, all of the health plots versus time can be seen in Fig. 7 and show much oscillation. The corresponding minimal health curve is then shown in Fig. 8. We can make this look more like the real survival data by scaling the minimal health values and plotting them as a percentage. When we do this, we find a plot that is a bit easier to compare to the real data. It is not perfect, of course, but it is very interesting that we are capturing the essential quality of the data using our theoretical model. The percentage minimal values are shown in Fig. 9. From these experiments, it is clear what is happening. The model

Fig. 4

Generating survival plots

Fig. 5

Minimum heath versus viral dose

Fig. 6

Set the parameters for another run

Fig. 7

Heath plots versus time for many viral doses

Fig. 8

Minimum heath versus viral dose

Fig. 9

Minimal heath percentage versus viral dose

can be written in terms of decay and push - pull terms as follows:

Thus, we have (t) always decreases unless the push - pull terms counteract that decay. Hence, what is important is the term

can oscillate as viral load increases. To do this, it is important for the two terms in (t) to be out of phase. Hence, roughly speaking cos(r 3 t−r 4 −θ 2+θ 3) must be sometimes negative when cos(r 4 −θ 2−θ 3) is positive. This allows for an increase in health of approximately ξ s 1 where ξ is the difference between the two terms. This is possible when the two cos arguments are out of phase by about π radians. Note it is also important that the exponential term allows growth. The interaction dynamics are determined by

and we have argued that the appropriate algebraic signs for this coefficient matrix are

We can have complex eigenvalues and hence oscillating behavior if the signs were

but then the real part of the eigenvalues would be negative and we would have to model the exponential term as . The induced oscillations would then be damped and the potential for the WNV survival curve would vanish.

Collateral damage

We should also see oscillations in the collateral damage. Recall the collateral damage population is given by

We can then substitute for (t), (t) and (t) and obtain

Now, collect terms as we did in our earlier simplifications. We rewrite as

We can also introduce an additional phase shift, ϕ, as follows. It will be different from the phase shift

as here we only use the H 1 partials: . We rewrite as

We can then use the the usual cos laws of addition and subtraction of angles to repackage this as

Now define and rewrite as

Since collateral damage is initially zero, we have as our final form

Previously, we used the simplification

This needs to be modified to

Hence, although we can use some of parameter choices from our health simulations, it is difficult to compare completely as S 2, Λ 1 and ϕ are different. Our final collateral damage function is then

We can easily run a quick simulation to see if our prediction that the collateral damage will have oscillations is true or not. We use the function Collateral() for this which is listed in Fig. 10. The parameter values we use are as similar as possible to the ones we used in generated the survival curves, although the values of s 2, Λ 1 and ϕ are necessarily new choices. We ran the simulation with these parameter values and then plotted both the maximum and minimum collateral values versus the viral dose in Fig. 11. Note that there is variation in the collateral damage due to the nonlinear interactions between the upregulation, , and the free virus, .

Fig. 10

Collateral damage calculations

Fig. 11

Collateral damage versus viral dose

Conclusions

The presence of a form of self damage in our WNV infection model therefore appears to be a consequence of the nonlinear interactions in the , and model:

Note, if the two populations and coincide, one signal is unnecessary - say u - and this model reduces to a two dimensional model

and the chance of oscillation between the cellular population groups is lost. Hence, we can note some consequences and predictions due to our model.

The crucial assumption here is that the viral infections effect on the host divides into two parts. For a WNV infection, these two cell populations are the dividing and nondividing infected cells, and , respectively. We can envision other infectious agents or triggers that give rise to such a split response which then in principle could engender a similar collateral damage response which we interpret as an autoimmune reaction. So there is hope that this approach could perhaps give us insight into more general autoimmune responses. Note that Fig. 11 shows there is collateral damage that oscillates due to the infectious agent which here is WNV. It is clear that other triggering events, another virus or bacteria or even an environment toxin, could give rise to this behavior as well.

Specific to the WNV model, we assume that , , , , which then says the coefficient matrix of the linearized upregulation and free virus model has the form

This algebraic sign pattern itself can give rise to complex eigenvalues for the linearized nonlinear interaction model and we have not explored this more general problem. We have noted in our discussion in Section “Results and discussions” that if we did not have , we could still have oscillatory behavior but it would be damped and therefore it would not explain the data we see in the survival experiments. Here, we have posited specific relations that give rise to clearcut oscillations. We have assumed and which gives rise to the characteristic coefficient matrix

The assumptions above then give rise to a general model of how the minimal health changes with varying initial virus dose, which appears to allow us to approximate the data we can measure.

The exact replication of the data found in the biological situation is unlikely to occur. Indeed, a standard dose response curve of survival does not usually repeat exactly from experiment to experiment, except in the form of it, even when using genetically identical animals. The ragged dose response of mortality seen in WNV infection is similarly subject to biological variability. Clearly, with no virus there will be 100% survival, and with a large amount of virus there will be 100% death, as expected in a standard dose response curve. In between these two doses, however, the response to infection is subject to probability, which affects the outcome (survival or death). Thus, if the experiment were undertaken several times it would show the ragged form on each occasion, but not exactly the same percentage survival at each dose used. This implies that small biological differences at the starting point of infection, albeit in genetically identical mice, may subtend a large range of endpoints. It is of interest to note that bypassing the early initiation of the adaptive immune response, by inoculating virus intracranially, that the standard graded dose response occurs with WNV. This is because the replication of the virus overtakes the animal before an effective immune response can be generated, emphasising the role of the immune system in generating this ragged survival curve [12].

We have shown that we can build a reasonable model of how WNV infects a host’s cell in such a way that damage to the host can decrease, even though the inoculating viral dose increases.

Methods

Our model is a macro one and we believe it provides insight as to how we can model more general auto-immune reactions. We propose that for an infectious agent or trigger to cause oscillations in health it is required that the trigger causes alterations in two distinct cell populations. Then, if the nonlinear interactions between these two populations satisfies the conditions for damped oscillatory response we have mentioned here, we should see oscillations in the host health. We consider this work essentially a theoretical model and as we explore what we can do to extend the results to more general auto-immune settings, we hope that we can generate greater mechanistic insight into auto-immune disease.