key: cord-1046026-uyim084w
authors: Iwami, Shingo; Takeuchi, Yasuhiro; Iwamoto, Kentaro; Naruo, Yoshimi; Yasukawa, Masahiro
title: A mathematical design of vector vaccine against autoimmune disease
date: 2009-02-07
journal: Journal of Theoretical Biology
DOI: 10.1016/j.jtbi.2008.09.038
sha: f0ea24dbad758b4ef547b8e2c54bd4bb721caf01
doc_id: 1046026
cord_uid: uyim084w

Abstract Viruses have been implicated in the initiation, progression, and exacerbation of several human autoimmune diseases. Evidence also exists that viruses can protect against autoimmune disease. Several proposed mechanisms explain the viral effects. One mechanism is “molecular mimicry” which represents a shared immunologic epitope with a microbe and the host. We consider, using a simple mathematical model, whether and how a viral infection with molecular mimicry can be beneficial or detrimental for autoimmune disease. Furthermore, we consider the possibility of development of a vector therapeutic vaccine that can relieve autoimmune disease symptoms. Our findings demonstrate that vaccine therapy success necessitates (i) appropriate immune response function, (ii) appropriate affinities with self and non-self antigen, and (iii) a replicative vector vaccine. Moreover, the model shows that the viral infection can cause autoimmune relapses.

The concept of autoimmunity was first predicted by Nobel Laureate Paul Ehrlich at the start of the twentieth century: he described it as ''horror autotoxicus'' (Janewa et al., 2004) . His experiments led him to conclude that the immune system is normally focused on responding to foreign materials; it has an inherent tendency to avoid attacking self tissues. Nevertheless, when this process goes wrong, the immune system can attack self tissues, resulting in autoimmune disease (Bell and Bird, 2005) . Autoimmune diseases, such as rheumatoid arthritis and multiple sclerosis (MS), can create life-long disability and increased mortality.

Even in 2008, we do not completely understand the primary initiators or causes of many of these autoimmune diseases. No single theory or mechanism can adequately explain all features or pathogeneses of autoimmune diseases. The clinically wide spectrum of autoimmune diseases is best considered as the mosaic of autoimmunity. The many factors involved are genetic, hormonal, immunological, and environmental (Deodhar, 1992) . In particular, virus infections have long been associated with autoimmune diseases, whether MS, diabetes, or myocarditis (Fujinami, 2004; Fujinami et al., 2006; Horwitz and Sarvetnick, 1999; von Herrath and Oldstone, 1996) . Several mechanisms have been proposed to explain virus triggers of autoimmune diseases. One mechanism is ''molecular mimicry'' (Fujinami, 2001; Fujinami et al., 2006; Libbey and Fujinami, 2002) , which represents a shared immunologic epitope with a microbe and the host. For example, molecular mimicry has been assigned a presumptive role in the pathogenesis of several human diseases, including insulin-dependent diabetes mellitus type-1 (IDDM), ankylosing spondylitis, Guillain-Barre syndrome, primary biliary cirrhosis, and MS (Fujinami, 2001; Horwitz and Sarvetnick, 1999; Janewa et al., 2004) . The immune response to the virus cross-reacts with self because of molecular mimicry, which engenders autoimmunity as follows (Fujinami et al., 2006) ; virus-infected-antigen presenting cells (APCs) present viral peptides in the contexts of MHC class I and II, respectively, to naive CD8 þ T cells and CD4 þ T cells. Activation of T cells engenders IFN-g production, which further activates APCs, leading to production of IL-12, a potent Tcell-differentiating cytokine. Effector CD4 þ T cells release proinflammatory cytokines such as IFN-g and IL-2, thereby stimulating T cells to differentiate into effector T cells. Activated T cells can also secrete IFN-g and TNF, which can engender macrophage activation. The activated macrophages in turn release TNF, nitric oxide, and reactive oxygen intermediates, which can kill infected cells and uninfected cells. Dead and dying cells are then phagocytosed by macrophages and dendritic cells (DCs), which can present self antigens to autoreactive CD4 þ T cells. Similarly, effector CD8 þ T cells can kill infected cells via perforin and granzyme granules. Cell debris is taken up by APCs, which can present self antigens to autoreactive CD8 þ T cells. The generation of such cells can engender autoimmune responses with enhanced inflammation if not modulated by regulatory T cells releasing IL-10 and/or TGF-g. Consequently, patients can develop autoimmune disease through virus-induced autoimmunity.

Although some viruses can modulate the development of autoimmune disease as discussed above, interestingly, some experimental evidence exists for experimental allergic encephalomyelitis (EAE) (Barnett et al., 1996; Fujinami et al., 2006) , which is an experimental model of MS, that viruses can protect against autoimmune disease. Possible mechanisms of protecting against autoimmune disease are considered as ''altered peptide ligand'', which activates regulatory cells that modulate the disease (Barnett et al., 1996; Fujinami et al., 2006) , and ''activationinduced cell death'' (AICD) which engenders anergy or unresponsiveness of T cells (Fujinami, 2001) . These imply that viruses having molecular mimicry with self proteins are useful to vaccinate against autoimmune disease. Using molecular biology and DNA manipulation methods, it has also been possible to express mimic proteins in adequate live vectors (Arnon and Ben-Yedidia, 2003) and thereby design transgenic vector vaccines (Janewa et al., 2004; Roitt et al., 1998) against autoimmune disease. The development of vaccines has been an important contribution of autoimmune disease therapy and public health.

Herein, we construct a simple mathematical model based on the autoimmune disease model proposed in Iwami et al. (2007a) and Iwami (2007) . We consider a viral infection that can induce cross-reactive immune responses with self antigen caused by molecular mimicry. Our model suggests that the viral infection can induce various symptoms of autoimmune disease such as relapse. Furthermore, we propose that a form of immune response function determines whether a viral infection can be beneficial or detrimental. Using the model, we consider the possibility of development of a vector vaccine that can relieve autoimmune disease symptoms.

The breaking of tolerance or unresponsiveness to self antigens, involving the activation of autoreactive lymphocytes, is a critical event in the pathogenesis of autoimmune disease (von Herrath and Oldstone, 1996) . The molecular mimicry theory has become an important paradigm to explain the triggering of autoaggressive T lymphocytes (Anderton, 2006) . Viruses and microbial agents might possess protein structures or shapes that mimic normal host self proteins. An immune response elicited against the pathogen will eliminate it and will cross-react with one or more self antigens that share determinants with the agent (von Herrath and Oldstone, 1996) . The cross-reactive immune response can break a tolerance for self antigens and might engender autoimmune disease. We consider a viral infection that can induce cross-reactive immune responses with a self antigen caused by molecular mimicry (see Fig. 1 ). To explore effects of the viral infection and dynamical behavior of the vicious cycle of autoimmunity, we propose the following mathematical model based on the autoimmune disease model proposed in Iwami et al. (2007a) and Iwami (2007) ;

Variables T, D, I, and V, respectively, signify the number of target (uninfected) cells, damaged cells (which implies a concentration of self antigen), cross-reactive immune cells, and viral agents with molecular mimicry. The immune responses eliminate target cells and viral agents at a rate ofb 1 and b 2 , respectively. We assume a ''Malthusian growth rate'' k in viral agents, which decay at a rate u as considered in Nowak et al. (1991) . The parameters a and g represent the decay rate of damaged and cross-reactive immune cells, respectively. The function gðTÞ is the ''target cell growth function''. In Iwami et al. (2007a) , we consider two target cell growth functions g 1 ðTÞ ¼ l À mT, which means a constant growth and g 2 ðTÞ ¼ l À mT þ pTð1 À T=LÞ which means a logistic growth. A form of this function considerably affects the dynamics of autoimmune disease (Iwami et al., 2007a; Iwami, 2007) . However, to examine the effect of viral infection with molecular mimicry specifically, we consider that target cells are always constant T ¼ T Ã (we omit target cell dynamics). This assumption is justified by some homeostatic ability of organs, which exactly balances the growth of target cells with the damage. A modeling approach of this kind is used in Nowak and May (1994) . On the other hand, the assumption is not biologically unnatural if the target cells are abundant in the organ. We leave the inclusion of target cell dynamics as a subject for future work. We can model the vicious cycle in the following form:

where b 1 ¼b 1 T Ã . Therein, b 1 and b 2 , respectively, represent a degree of affinity to self and viral antigens (non-self antigen). A property of molecular mimicry can be described using these parameters. The function f ðD; VÞ is the ''immune response function'' (see Iwami et al., 2007a) . In general, the proliferation ability is saturated for a sufficiently large amount of antigens (Borghans and de Boer, 1995; Borghans et al., 1998; De Boer and Perelson, 1995) . We consider that the proliferation of immune response is dependent on the total number of self and viral antigens as

The forms of f 1 and f 2 can be regarded, respectively, as functional responses of Holling types II and III. Parameters m and h, respectively, signify the maximum proliferation rate and the efficiency of the proliferation.

The mechanism of immune proliferation remains unclear. However, using mathematical and experimental models, many theoretical immunologists study the proliferation function to obtain better knowledge related to immune response (Borghans et al., 1999; Wodarz and Jansen, 2001; Wodarz and Thomsen, 2005) . The function might depend on the immune cell type, genetic and physiological factors, and so on. As described herein, we investigate two immune response functions f 1 and f 2 which are not biologically unnatural. The difference of the immune response function strongly affects dynamic behaviors of model (1), as discussed below.

We consider in the context that immune response function is convex form f 1 . Therefore, our mathematical model is the following:

We must consider two different situations for model (2): (i) k À uo0 and (ii) k À u40. If (i) holds, then lim t!1 VðtÞ ¼ 0, which implies that the virus cannot maintain its replication in the host.

Then the system has two possible equilibria:

Equilibria E h and E a , respectively, represent the ''healthy state'' and the ''autoimmune disease state''. If (ii) holds, then the virus can persist in the host and we have one more equilibrium,

which represents the ''complication state'' (i.e., the patient develops autoimmune disease with the viral infection). Furthermore, the number of viral agents can explode ðlim t!1 VðtÞ ¼ 1Þ under (ii) and the dynamics of model (2) converges to some steady state

which represents the ''infection state'' (see Appendix C). The exact expressions for equilibria are referred to Appendix A. Fig. 2 portrays the existence and stability conditions of these equilibria in model (2). For detailed mathematical analysis of the equilibria, see Appendix A.

3.1.1. Symptoms of autoimmune disease in model (2) We consider that the number of damaged cells, such as D a or D c , represents the level of autoimmune disease progression. A strong immune affinity with self antigen (large b 1 ) engenders the development of autoimmune disease E a (see Fig. 2 (c)) when the virus cannot establish the persistence infection ðk À uo0Þ. Because dD a =db 1 40, an increase of the affinity with self antigen deteriorates autoimmune disease (which means that autoimmune disease is worsening; so is the patient's condition). In other words, if we can appropriately reduce the immune affinity with self antigen using some drugs or therapies, we can delay autoimmune disease progression (a marked decrease of the affinity with self antigen ðb 1 oahg=mÞ can engender immune tolerance E h ).

On the other hand, when the virus can maintain its replications ðk À u40Þ, strong immune affinities with self and viral antigen (large b 1 and b 2 ) also engender the development of autoimmune disease without viral replication (see Fig. 2 (d)). However, if the immune affinity with self antigen is weak, the patient develops autoimmune disease with viral infection E c . Because dD a =db 1 40 and dD c =db 1 40, a decrease of the affinity with self antigen also reduces autoimmune disease progression. The immune tolerance cannot occur under viral persistence (the viral infection initiates autoimmune disease) even if we can reduce the immune affinity with the self antigen. A decrease of the affinity with viral antigen engenders deterioration of autoimmune disease progression because dD c =db 2 o0 whenever the patient is in a complication state. Furthermore, a decrease of the affinity with viral antigen engenders unlimited replications of viral agents and a deterioration E 1 (note that D 1 4D a ).

We consider in the context that immune response function is sigmoid form f 2 . Therefore, our mathematical model is the following:

We must also consider two different situations for model (3):

The virus cannot maintain its replication in the host under (i). Then the system has three possible equilibria:

0Þ, but E À a is not biologically appropriate because the equilibrium is always unstable even if it exists (see Appendix B). The equilibria E h and E þ a also, respectively, represent the ''healthy state'' and the ''autoimmune disease state''. On the other hand, if the virus can persist in the host ((ii)), then we have one more equilibrium, If k À u40, then model (3) represents various dynamical behaviors (see Fig. 4 , which corresponds to Fig. 3(d) ). The parameter region (I) represents that E þ a is stable; the region (V) represents that E Ã c is stable. In addition, E þ a and E Ã c are stable simultaneously (bistability) in region (III). Therefore, the orbit converges to E þ a or E Ã c according to its initial values. In parameter region (IV), we can observe periodic behavior by Hopf bifurcation of E Ã c (all orbits except E Ã c converge to the limit cycle). Furthermore, the limit cycle and E þ a are stable simultaneously in region (II); the orbit also converges to the limit cycle or E þ a according to its initial values (see Fig. 5 ). However, if we choose a parameter set near Gðb 2 Þ in (II), we can numerically confirm that almost all orbits converge to E þ a because the amplitude of the periodic orbit increases as a parameter set approaches Gðb 2 Þ, the orbit crosses the stable manifold of E þ a , the periodic orbit vanishes and the orbit converges to E þ a . This phenomenon is also observed and particularly explained in Iwami (2007) . Therefore, we can classify the dynamical behavior of model (3) under k À u40. For detailed mathematical analysis of the equilibria, see Appendix B.

We also consider that the number of damaged cells, such as D þ a or D Ã c , represents a level of autoimmune disease progression. A strong immune affinity with self antigen ðb 1 Þ tends to develop into autoimmune disease E þ a according to initial values (see Fig. 3 (c)) when the virus cannot establish persistent infection ðk À uo0Þ. The healthy state E h is always stable, which implies that the immune tolerance can occur after development of autoimmune disease if we can remove the damaged cells and the immune responses (simulations not shown). Furthermore, because dD þ a =db 1 40, an appropriate decrease of the immune affinity with self antigen can delay autoimmune disease progression (a considerable decrease of the affinity with self antigen ðb 1 o2ahg=mÞ can engender the immune tolerance E h ).

On the other hand, when the virus can maintain its replications ðk À u40Þ, strong immune affinities with self and viral antigen (large b 1 and b 2 ) also engender the development of autoimmune disease without viral infection ( Fig. 4(I) ). However, if the immune affinity with self antigen is weak, the patient develops autoimmune disease with a viral infection (Fig. 4(V) ). Because dD þ a =db 1 40 and dD Ã c =db 1 40, a decrease of the affinity with self antigen also reduces autoimmune disease progression in regions (I) and (V). A decrease of the affinity with viral antigen engenders

Existence and stability conditions of the equilibria in model (2): panels (a) and (b), respectively, present the existence conditions of the equilibria with k À uo0 and k À u40. Horizontal and vertical stripes, respectively, signify the existence regions for Ea and Ec. Equilibrium E h always exists. Panels (c) and (d), respectively, present the stability conditions of the equilibria with k À uo0 and k À u40. The horizontal and vertical stripes, respectively, signify the stability regions for Ea and Ec . Shadows and asterisks, respectively, denote the stability and attractive regions for E h (or E1). Function Hðb 2 Þ is referred from Appendix A.

deterioration of autoimmune disease progression because dD Ã c =db 2 o0 in region (V) (a marked decrease of affinity with viral antigen engenders unlimited viral replications). Interestingly, because E þ a and E Ã c are stable simultaneously in region (III), the symptoms of autoimmune disease depend on the patients' states. Furthermore, in regions (II) and (IV), the relapse of autoimmune disease, which is a common symptom of autoimmune disease, can occur. Actually, the symptoms of autoimmune disease also depend on the patients' states in (II) (patients with affinities near Gðb 2 Þ in (II) tend not to represent relapse symptoms, as discussed above). Consequently, the viral infection prevents immune tolerance and engenders the relapse pattern of autoimmune disease.

Virus infection can initiate or accelerate autoimmune disease via epitope spreading (Libbey and Fujinami, 2002; Miller et al., 1995) and molecular mimicry, thereby engendering the development of an inflammatory region with activated APCs and possible presentation of a self antigen (Fujinami et al., 2006; von Herrath and Oldstone, 1996; von Herrath et al., 2003) . However, several interesting experimental examples exist for prevention of autoimmune disease caused by viral infections. Possible mechanisms of prevention caused by viral infections are immunosuppression, chemokine gradients, apoptosis of autoaggressive lymphocytes, Fig. 3 . Existence and stability conditions of the equilibria in model (3). Panels (a) and (b), respectively, represent the existence conditions of equilibria with k À uo0 and k À u40. The horizontal and vertical stripes, respectively, mark the existence regions for E þ a and E Ã c . Equilibrium E h always exists. Panels (c) and (d), respectively, represent the stability conditions of equilibria with k À uo0 and k À u40. Horizontal and vertical stripes, respectively, signify the stability regions for E þ a and E Ã c . Shadows and asterisks, respectively, denote the stability and attractive regions for E h (or E1). Functions Fðb 2 Þ and Gðb 2 Þ are referred from Appendix B. and so on (Fujinami et al., 2006) . Furthermore, experimental evidence indicates that a viral infection with molecular mimicry can provide protection from EAE (Fujinami, 2001; Fujinami et al., 2006) . Reasons for the protection are said to be altered peptide ligand, AICD, and so on, which implies that viruses having molecular mimicry with self proteins are useful for vaccination against autoimmune disease.

We consider whether infection by a virus having molecular mimicry with a self antigen that can replicate itself ðk À u40Þ is beneficial or detrimental for autoimmune disease. First we assume that autoimmune disease has already developed before the patient is infected with the virus having molecular mimicry with self proteins ðVð0Þ ¼ 0Þ. In model (2), we assume that E a is stable ðb 1 4agh=mÞ with Vð0Þ ¼ 0, which corresponds to Fig. 2(c) , irrespective of the sign of k À u (we remark that the infection with k À u40 persists and initiates autoimmune disease even if b 1 oagh=m). Furthermore, to avoid the obvious result, we assume that the immune affinity with viral antigen is high ðb 2 4gðk À uÞ=mÞ. Actually, if the affinity with viral antigen is low ðb 2 ogðk À uÞ=mÞ, then the viral population explodes (see Appendix C) and the number of damaged cells ðD 1 Þ represents its maximum value, which implies that the viral infection detrimental affects the patients. We also exclude the high affinity with self antigen ðb 1 4Hðb 2 ÞÞ because the viral infection cannot affect the disease progression. Consequently, we assume the following conditions: k À u40; agh m ob 1 oHðb 2 Þ; b 2 4 gðk À uÞ m .

We evaluate the effect of viral infection as follows. Let Vð0Þ40; then the stable equilibrium changes from E a to E c (the patients get the viral infection). We then have the following relations:

b 1 oHðb 2 Þ ) D a oD c ; I a oI c .

They imply that the viral infection always imparts a detrimental effect on patients and accelerates autoimmune disease on the region because the viral infection increases damaged cells ðD a oD c Þ and cross-reactive immune responses ðI a oI c Þ (see Fig. 6 ).

On the other hand, in model (3), we assume that E þ a is stable ðb 1 42agh=mÞ with Vð0Þ ¼ 0, which fundamentally corresponds to Fig. 3(c) , irrespective of a sign of k À u. Furthermore, let the immune affinity with viral antigen be high ðb 2 4gðk À uÞ=mÞ. We also exclude high affinity with self antigen ðb 1 4Gðb 2 ÞÞ because of a neutral effect of the infection. To avoid a bad prognosis such as a relapse of autoimmune disease, we also exclude the possibility of the relapse caused by viral infection (b 2 42gðk À uÞ=m and Fðb 2 Þob 1 oGðb 2 Þ; (II)(IV) in Fig. 4) . Consequently, we assume the gðk À uÞ m .

Let Vð0Þ40 and then the convergence equilibrium changes from E þ a to E Ã c (the patients get the viral infection) or it does not change because the region (III) in Fig. 4 represents bistability of E þ a and E Ã c . When the viral infection changes the convergence equilibrium, we evaluate the effect of viral infection as shown below.

The expressions presented above imply that the viral infection can give patients a beneficial effect and relieve autoimmune disease symptoms for the former case because the viral infection engenders a decrease of damaged cells and cross-reactive immune responses (Fig. 6) .

Vaccines are, by definition, prophylactic, but in recent years we saw an interest in developing therapeutic vaccines, in infectious diseases (for diseases such as AIDS, tuberculosis, and peptic ulcer), in cancer (a variety of approaches to combat different kinds of cancer), and in autoimmune diseases (a definite success in developing a drug/vaccine against MS and hopes for myasthenia gravis, lupus and diabetes) (Arnon and Ben-Yedidia, 2003) . Using molecular biology and DNA manipulation methods, it is possible to produce a therapeutic vaccine against autoimmune disease.

We consider the possibility of development of a vector vaccine having molecular mimicry with self proteins. After emergence of autoimmunity, the vector vaccine can reproduce itself, induce a cross-reactive immune response, and be removed by the immune response in patients as a similar mechanism for viral infection with molecular mimicry (see Fig. 7 ). Therefore, we can consider that the vaccine can be described similarly as the viral infection in terms of the mathematical model. Fujinami (2001) and Fujinami et al. (2006) , explained that a possible mechanism of protection from EAE by the viral infection might be suppression of autoreactive immune cells caused by regulatory cells, AICD, and so on. However, our results demonstrate that a viral infection in model (3) can be effective for autoimmune disease without involving these suppressive effects. This shows that, when the immune response function is f 2 , the vector vaccine can reduce the cross-reactive immune response ðI þ a 4I Ã c Þ and relieve symptoms of autoimmune disease ðD þ a 4D Ã c Þ, which implies that the vector vaccine can be effective by itself even if we consider no additional suppressive abilities in our immune system. Therefore, to make the vector vaccine effective, we at least require that the immune proliferation be a sigmoid function such as f 2 (we explain the immune response function in Discussion). The affinities with self and non-self antigen are in the beneficial region, as shown in Fig. 6 . Actually, our immune system has high affinity with non-self antigen and low affinity with self antigen because, in the process of differentiation thymic, lymphocytes undergo positive and negative selection. Positive selection generates a functional T cell repertoire restricted to self MHC expressed on the epithelial cells of the thymic cortex. Negative selection eliminates T cells that are aggressively reactive to self antigen (Goldrath and Bevan, 1999; Janewa et al., 2004) . Consequently, the beneficial region, which has high affinity with non-self antigen and low affinity with self antigen, is biologically realistic. Furthermore, the vector vaccine must replicate effectively in autoimmune disease patients, which means that the orbit converges to E Ã c . We remark that the vaccine might be unable to replicate (the orbit converges to E þ a ) according to patients' states because of bistability. Therefore, the vector vaccine might have to be used with other immunosuppressive drugs to switch the patient state from E þ a to E Ã c . Although many restrictions discussed above exist for the success of therapeutic vaccine, we can theoretically design the transgenic vector vaccine (because of these stringent restrictions, it is also true that the transgenic vector vaccine might only be effective to certain specific patient states). Moreover, if we consider the additional effects for suppression of autoreactive immune cells, these restrictions can be relaxed. However, we leave the inclusion of additional effects as a subject for future work.

Increasing evidence exists that infectious agents play an important role in autoimmune disease (Janewa et al., 2004; Roitt et al., 1998) . Viruses are an important factor that can precipitate autoimmune disease by various mechanisms (Fujinami, 2004; Pender, 2004; von Herrath et al., 2003) . For example, viruses that stimulate the production of IL-12, such as herpes simplex virus, human herpesvirus, influenza virus, and coronavirus, have been isolated from or have been associated with exacerbation of MS (Fujinami, 2001; Libbey and Fujinami, 2002; Monteiro et al., 1998; Whitton and Fujinami, 1999) . On the other hand, viruses can abrogate an ongoing autoimmune reaction by inducing apoptosis of autoreactive cells, by secreting various cytokines, or by immune suppression (Fujinami, 2001) . Nevertheless, it is difficult to obtain direct evidence for virus-induced initiation of autoimmune disease or protection from autoimmune disease because we are all infected by multiple viruses (Fujinami et al., 2006) .

Using the simple mathematical model, we analyzed whether and how a viral infection having molecular mimicry with self proteins can impart a detrimental or beneficial effect to autoimmune disease patients (see Fig. 6 ). Furthermore, we consider the possibility of development of a therapeutic vaccine against autoimmune disease. Our findings suggest that the success of therapeutic vaccine necessitates the following: (i) patients have an appropriate immune response function, such as f 2 ; (ii) affinities with self and non-self antigen are in Fig. 7 . A transgenic vector vaccine having molecular mimicry with self proteins; after emergence of autoimmunity, the vaccine can induce a cross-reactive immune response as a similar mechanism for viral infection with molecular mimicry. the beneficial region shown in Fig. 6; (iii) the vector vaccine replicates in vivo. Although these restrictions can be relaxed if we consider additional immunological effects such as immune suppression, memory response, and exhaustion of effector cells, we found the distinct possibility of designing a therapeutic vaccine.

In Iwami et al. (2007a) and Iwami (2007) , we investigated the influence of functional form of immune proliferation on autoimmune disease symptoms. We demonstrated that sigmoid function such as f 2 can induce a bistable structure and periodic behavior. As described in this paper, we omit the target cell dynamics and add the viral replication cycle compared with the model considered in Iwami et al. (2007a) . The difference of immune response functions also strongly affects the autoimmune disease symptoms. If the virus cannot establish a persistent infection ðk À uo0Þ, then f 2 induces bistability (Fig. 3(c) ) but f 1 does not, then (Fig. 2(c) ). On the other hand, if the virus can replicate itself ðk À u40Þ, then f 2 induces bistability and a periodic orbit (Fig. 3(d) ), but f 1 does not, then (Fig. 2(d) ). Consequently, the sigmoid function f 2 , which is biologically more reasonable (because APCs only slightly induce immune cells when only a few antigens exist, Iwami et al., 2007a) than the convex function f 1 represents various symptoms such as relapse in autoimmune disease. That slight inducement implies that the various symptoms of autoimmune disease might be related with the function. However, the results might not be sufficiently robust to the form of immune proliferation function. We need to know an appropriate proliferation function that has biological relevance. For example, APCs are known to induce immune cells only slightly when many antigens exist (high zone tolerance) (Janewa et al., 2004) . For that reason, we might have to use a bell-shaped proliferation function to consider high zone tolerance instead of the proliferation as an increasing function of antigen load (De Boer et al., 1993 . Consequently, a more complete understanding of immune proliferation must be the foundation for the development of a therapeutic vaccine.

Viruses trigger autoimmune disease, but they are also likely to be important for reactivation of autoimmunity (viruses can behave as reactivators of autoimmune disease) (Horwitz and Sarvetnick, 1999) . Some clinical data show that viral infections trigger MS relapse (Andersen et al., 1993) . It has been suggested that a determinant spreading (self-antigen diversification) is a relapse mechanism (Lehmann et al., 1993) . Furthermore, in Borghans et al. (1998) , they showed that T cell regulatory circuitry induces autoimmune relapse using a simple mathematical model. Another possible mechanism of the relapse is considered as crossreactive immune responses through a process of molecular mimicry. In our model, relapses can also occur under the viral infection with molecular mimicry (in the absence of determinant spreading). The immune response function f 2 can induce a limit cycle that corresponds to relapse of autoimmune disease only when the virus can replicate itself (see Fig. 4 ). This implies that the viral infection engenders the relapse. Therefore, our model supports that the cross-reactive immune response is also a relapse mechanism.

This study highlights the immune response functions and molecular mimicry with the self antigen to investigate the possibility of vector vaccine development. Although our model might over-simplify complex interactions in autoimmune disease, its simplicity illustrates the general and qualitative properties of viral infection with molecular mimicry. This prediction should be verified using actual experiments. Our model is a starting point, but must include more immunological factors such as immunosuppression and apoptosis to merit further theoretical and mathematical investigation.

We analyze the existence and stability conditions of the equilibria in (2). The model has three possible equilibria:

It might be readily apparent that E h always exists. The existence conditions of E a and E c are as follows:

It is noteworthy that E c cannot exist in R 3 þ if k À uo0. In Fig. 2(a) and (b), we summarize these conditions in the b 1 2b 2 plane.

Hereinafter we explain the stability of these equilibria in detail. The Jacobian matrix of (2) at E h is 

That is, if the above conditions hold, then E h is locally asymptotically stable (LAS); otherwise E h is unstable. We can show that E h is globally asymptotically stable (GAS) using similar methods to those described in Iwami et al. (2007b) . The Jacobian matrix of (2) at E a is The characteristic equation of JðE a Þ is 

We remark that the above second condition holds with k À uo0 if the first is satisfied, which implies that E a is always LAS if it exists. Let

It is clear that Hðb 2 Þ4agh=m when k À u40. Consequently, E a is LAS if b 1 4Hðb 2 Þ. We can show that E a is GAS under k À uo0 using similar methods to those explained in Iwami et al. (2007b) .

The Jacobian matrix of (2) at E c is

Àb 2 V c 0 2 6 6 6 4 3 7 7 7 5 .

The characteristic equation of JðE c Þ is

Therefore, from the Routh-Hurwitz criterion, the stability of E c is determined by the sign of a 1 a 2 À a 3 . Because the following relation holds:

Direct but tedious calculations yield a 1 a 2 À a 3 40 b 2 À gðk À uÞ=m 2 1 À 1 mb 2 g 2 ðk À uÞ 2 gðk À uÞ þ aða þ gÞ b 2 À gðk À uÞ=m  gðk À uÞ gðk À uÞ þ aða þ gÞ þ b 2 b 2 À gðk À uÞ=m aða þ gÞ gðk À uÞ þ aða þ gÞ .

It is noteworthy that gðk À uÞ gðk À uÞ þ aða þ gÞ þ b 2 b 2 À gðk À uÞ=m a ða þ gÞ gðk À uÞ þ aða þ gÞ

That relation implies that a 1 a 2 À a 3 40 ( b 1 o agh m b 2 b 2 À gðk À uÞ=m . Therefore, we can conclude that E c is always LAS whenever it exists. In Fig. 2(c) and (d), we present these conditions in the b 1 2b 2 plane.

We also analyze the existence and stability conditions of the equilibria in (3). The model has four possible equilibria as follows:

Clearly, E h always exists. The existence condition of E AE a and E Ã c are as follows:

It might be readily apparent that E Ã c cannot exist in R 3 þ if k À uo0.

In Fig. 3(a) and (b), we present these conditions in the b 1 2b 2 plane. Hereinafter, we explain the stability of these equilibria in detail. From a direct calculation, the eigenvalues of JðE h Þ are Àa, Àg, and k À u, where J is the Jacobian matrix of (3). Those results imply that if k À uo0, then E h is always LAS; otherwise, E h is unstable.

The Jacobian matrix of (2) at E AE a is The characteristic equation of JðE AE a Þ is

As a result of tedious but straightforward calculations, we can show that ag À 2mb 1 h 2 D À a =ðh 2 þ D À a Þ 2 o0 and ag À 2mb 1 h 2 D þ a = ðh 2 þ D þ a Þ 2 40. That relation implies that E À a is always unstable if it exists. Consequently, the stability of E þ a is determined by a sign of k À u À b 2 I þ a , which is an eigenvalue of JðE þ a Þ. If k À uo0, then E þ a is LAS. On the other hand, if k À u40, then we have

Therefore, if ðk À uÞ=b 2 À m=2go0, then E þ a is LAS. Although, if ðk À uÞ=b 2 À m=2g40, then the following relations hold: 

Viral infections trigger multiple sclerosis relapses: a prospective seroepidemiological study

Avoiding autoimmune disease-T cells know their limits

Old and new vaccine approaches

Virus encoding an encephalitogenic peptide protects mice from experimental allergic encephalomyelitis

A minimal model for T-cell vaccination

T cell vaccination in experimental autoimmune encephalomyelitis: a mathematical model

Competition for antigen site during T cell proliferation: a mathematical interpretation of in vivo data

Towards a general function describing T cell proliferation

Immune network behavior-I. From stationary state to limit cycle oscillations

A new bell-shaped function for idiotypic interactions based on cross-linking

Autoimmune diseases: overview and current concepts of pathogenesis

Viruses and autoimmune disease-two sides of the same coin?

Molecular mimicry, bystander activation, or viral persistence: infections and autoimmune disease

Selecting and maintaining a diverse T-cell repertoire

Viruses, host responses, and autoimmunity

Potential mechanisms of relapse in autoimmune disease

Dynamical properties of autoimmune disease models: tolerance, flare-up, dormancy

Avian-human influenza epidemic model

Immunobiology: The Immune System in Health and Disease

Determinant spreading and the dynamics of the autoimmune T-cell repertoire

Are virus infections triggers for autoimmune disease?

Blockade of CD28/B7-1 interaction prevents epitope spreading and clinical relapses of murine EAE

Role of interleukin-12 in primary influenza virus infection

Superinfection and the evolution of parasite virulence

Antigenic diversity thresholds and the development of AIDS

Infection of autoreactive B lymphocytes with EBV, causing chronic autoimmune diseases

Virus-induced autoimmune disease

Microorganisms and autoimmunity: making the barren field fertile?

Viruses as triggers of autoimmunity: facts and fantasies

The role of T cell help for anti-viral CTL responses

Does programmed CTL proliferation optimize virus control?

The authors would like to thank anonymous referees for very helpful suggestions and comments which improved the quality of this paper and study. SI was supported by Research Fellowships of the Japan Society for the Promotion of Science for Young Scientists.

Here we assume that b 2 4gðk À uÞ=m; otherwise the above condition does not hold. Therefore, if the following conditions hold, gðk À uÞ m ob 2 o 2gðk À uÞ m :The Jacobian matrix of (2) at E Ã c isTherefore, from the Routh-Hurwitz criterion, the stability of E Ã c is determined by the sign of a 1 a 2 À a 3 . Because the following relation holds:we can obtain a 1 a 2 À a 3 ¼ agða þ gÞ þ 2g 2 fmb 2 À gðk À uÞgðk À uÞ mb 2 À gðk À uÞ þ aða þ gÞ a 2b 1 fmb 2 À gðk À uÞg ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi gðk À uÞ p ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi mb 2 À gðk À uÞ p mhb 2 2 . Therefore, if the following conditions hold, b 1 o ab 2 gh½mab 2 ða þ gÞ þ 2gðk À uÞfmb 2 À gðk À uÞg 2faða þ gÞ þ gðk À uÞg ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi gðk À uÞ p fmb 2 À gðk À uÞg 3=2 , then E Ã c is LAS. We define the following functions.Fðb 2 Þ ¼ ab 2 gh½mab 2 ða þ gÞ þ 2gðk À uÞfmb 2 À gðk À uÞg 2faða þ gÞ þ gðk À uÞg ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi gðk À uÞ p fmb 2 À gðk À uÞg 3=2 ,Because we have In Fig. 3(c) and (d), we present these conditions in the b 1 2b 2 plane.

We show that VðtÞ in (1) blows up for sufficiently large t if b 2 ogðk À uÞ=m and k À u40. Because maxff 1 ; f 2 gpm, we have the following inequality:I 0 ¼ f ðD; VÞ À gIpm À gI.Using a comparison theorem for ordinary differential equations, we can obtain IðtÞp Ið0Þ À m g e Àgt þ m g . Therefore, we can evaluateBecause b 2 ogðk À uÞ=m and k À u40, there exist some T such as k À u À b 2 ðIð0Þ À m=gÞe Àgt À mb 2 =g40 for any t4T, which implies that lim t!1 VðtÞ ¼ 1. Consequently, we can show that lim t!1 IðtÞ ¼ m=g and lim t!1 DðtÞ ¼ mb 1 =ag. All orbits in (1) converge to E 1 if b 2 ogðk À uÞ=m and k À u40.