**DOI:**10.1128/JVI.03395-12

## ABSTRACT

Theiler murine encephalomyelitis virus (TMEV) infection of a mouse's central nervous system is biphasic: first the virus infects motor neurons (acute phase), and this is followed by a chronic phase in which the virus infects glial cells (primarily microglia and macrophages [Mϕ]) of the spinal cord white matter, leading to inflammation and demyelination. As such, TMEV-induced demyelinating disease in mice provides a highly relevant experimental animal model for multiple sclerosis. Mathematical models have proven valuable in understanding the *in vivo* dynamics of persistent virus infections, such as HIV-1, hepatitis B virus, and hepatitis C virus infections. However, viral dynamic modeling has not been used for understanding TMEV infection. We constructed the first mathematical model of TMEV-host kinetics during acute and early chronic infections in mice and fit measured viral kinetic data with the model. The data fitting allowed us to estimate several unknown parameters, including the following: the rate of infection of neurons, 0.5 × 10^{−8} to 5.6 × 10^{−8} day^{−1}; the percent reduction of the infection rate due to the presence of virus-specific antibodies, which reaches 98.5 to 99.9% after day 15 postinfection (p.i.); the half-life of infected neurons, 0.1 to 1.2 days; and a cytokine-enhanced macrophage source rate of 25 to 350 Mϕ/day into the spinal cord starting at 10.9 to 12.9 days p.i. The model presented here is a first step toward building a comprehensive model for TMEV-induced demyelinating disease. Moreover, the model can serve as an important tool in understanding TMEV infectious mechanisms and may prove useful in evaluating antivirals and/or therapeutic modalities to prevent or inhibit demyelination.

## INTRODUCTION

One of the few available experimental animal models of virus-induced demyelination is Theiler murine encephalomyelitis virus (TMEV) infection in mice, which has been recognized as an experimental analog of multiple sclerosis (MS) (1, 2). TMEV belongs to the genus Cardiovirus in the family Picornaviridae (3, 4). It is a highly cytolytic nonenveloped (or naked) virus which consists of a spherical protein shell that encapsidates a single positive-strand RNA genome of about 8,100 nucleotides (3, 4). TMEV strains have been divided into two groups: high neurovirulence and Theiler original (TO) (5), or low neurovirulence. The high-neurovirulence group includes virus strains that cause a rapidly fatal encephalitis in mice, whereas members of the TO group have much lower virulence and establish chronic (or persistent) viral infection, inflammation, and demyelinating disease.

Persistent infection with virus strains of the TO group is characterized by a biphasic central nervous system (CNS) disease (6). During the first phase (termed here the acute phase), the virus infects sensory and motor neurons and causes an acute but mild encephalomyelitis that lasts for 1 to 2 weeks. The acute phase is followed by a second phase (termed here the chronic phase), during which the virus infects glial cells, primarily microglia and macrophages (Mϕ), of the spinal cord white matter (7, 8). Unlike viral infections with agents such as human immunodeficiency virus type 1 (HIV-1), hepatitis B virus (HBV), and hepatitis C virus (HCV), the unique feature of TMEV infection is that the dominant cell type of infection changes during the transition from the acute to the chronic phase of infection. The temporal dynamics of TMEV RNA replication in the CNS of susceptible and resistant strains of mice has been examined by quantitative RT-PCR and correlated with host immune responses (9). During the acute phase of infection in both susceptible and resistant mice, levels of viral replication peak and then decrease at approximately 5 days postinfection (p.i.) in parallel with the appearance of virus-specific antibodies (9) and CD8^{+} T cells (10–12). However, after about 2 weeks p.i., viral RNA numbers begin to increase again only in the spinal cords of susceptible mice. High-viral-genome equivalents in spinal cords are observed only in susceptible strains of mice, such as SJL/J mice developing demyelinating disease (9, 13).

During the acute phase, TMEV epitopes are presented to Th1 CD4^{+} lymphocytes by antigen-presenting cells in peripheral lymphoid organs, which causes them to divide and expand to effector CD4^{+} T cells, which in turn secrete cytokines and chemokines (14). As a result, monocytes that are recruited into the CNS differentiate into Mϕ, some of which become susceptible to TMEV infection, while others are activated and become resistant to TMEV infection (because they produce type I IFN [15, 16]) and are lost due to death and/or emigration from the spinal cord (reviewed in reference 17). A restriction of TMEV infectious virus production was observed in Mϕ isolated from mouse CNS (18) and also in murine-infected Mϕ in culture. Interestingly, this restriction is likely to reflect a block later in the virus life cycle rather than in TMEV RNA replication, translation, and polyprotein processing (19, 20). Thus, infection of Mϕ is characterized by unaffected levels of TMEV RNA replication and protein synthesis but with markedly reduced levels of infectious virus. Because TMEV is a cytolytic virus and leads to apoptosis, persistence of the infection requires continuous virus spread in the face of host immune responses (e.g., virus-specific antibodies and CD8^{+} T cells) and markedly reduced production of free infectious particles. Thus, infection via phagocytosis of infected apoptotic debris (or blebs) containing infectious viral RNA by Mϕ and other glial cells such as oligodendrocytes has been hypothesized, as the presence of virus-neutralizing antibodies (36) does not seem to prevent virus spread during the chronic phase (reviewed in reference 17).

Mathematical models have proven valuable in understanding the *in vivo* dynamics of HIV-1 (21–23), HBV (24–26), and HCV (27, 28), viruses that all cause persistent infections, as well as the dynamics of influenza A (29–31) that causes an acute infection. Mathematical modeling is also improving our understanding of intracellular viral genome dynamics (32–34), T-cell dynamics, and the quantitative events that underlie the immune response to pathogens (21). By comparing mathematical models with the data obtained from patients being treated with antiviral drugs, it has been possible to determine many quantitative features of these infections. The most dramatic finding has been that even though HIV-1-, HBV-, and HCV-induced diseases occur on a time scale of one or more decades, there are also very rapid dynamic processes that occur in hours to days, as well as slower processes that occur in weeks to months. Dynamic modeling and parameter estimation techniques have uncovered these important features of HIV-1, HBV, and HCV infections, with a subsequent impact on treatment protocols in patients. A mathematical model of chronic herpes simplex virus infection, which resides in neurons that innervate the mucosal tissue of the genital tract in infected humans, was also developed and suggested a control point of viral infectivity that could serve as a vulnerable target for therapeutic intervention (35).

In this paper, we describe a mathematical model of TMEV infection and estimate unknown viral and host parameters. The results provide insight into early (i.e., during the first 30 days p.i.) TMEV-host dynamics and may serve as an important starting point in modeling long-term (beyond 30 days of infection) demyelinating disease and possibly inform future therapeutic approaches for MS.

## MATERIALS AND METHODS

Model description.Our model of TMEV infection is based on the infection scheme shown in Fig. 1. We assume that during acute phase of infection TMEV primarily infects motor neurons (36), and equations 1 to 3 apply, while during the chronic (or persistent) phase TMEV mainly infects Mϕ via phagocytosis of infectious debris, and equations 4 to 6 apply. The infectious debris units are generated upon death of infected neurons and Mϕ.

Our model of the acute phase of infection considers uninfected motor neurons (*T*), infected motor neurons (*I*), and free virions (*V*). Uninfected motor neurons survive for essentially the life of the host, and hence a death term is not included in equation 1. Motor neurons, when infected by TMEV at rate β_{n}, release virions at rate *p _{n}* per cell and die by cytolysis at rate δ

_{n}. Free virions are cleared from the CNS at rate

*c*per virion. The factor η(

*t*), where 0 ≤ η(

*t*) ≤ 1, describes the reduction of the infection rate due to virus-specific antibodies (equations 7 and 8), which are detected as early as day 5 p.i. (9).

Due to acute TMEV infection, monocytes are recruited into the CNS and differentiate into susceptible Mϕ at rate *s*(*t*) (equation 9) at time *t _{M}* p.i. (9, 37), triggering the transition from the acute to the persistent phase of infection. Our model of the chronic phase infection includes susceptible Mϕ, infected Mϕ (

*I*

_{M}_{ϕ}), and infectious debris (

*I*). Susceptible Mϕ are lost at rate δ

_{D}_{Mϕ}due to activation of Mϕ, degradation/death and/or emigration out of the spinal cord. Infected Mϕ are lost at rate δ

_{IMϕ}, due to apoptosis and/or emigration out of spinal cord. The viral RNA (vRNA) present in apoptotic remnants (or blebs) of infected cells allow infection of susceptible Mϕ at rate β via phagocytosis (based on the hypothetical mechanism described in Fig. 2 in reference 17). Infectious debris units (

*I*) are produced upon death/apoptosis of infected Mϕ or neurons and are cleared at rate

_{D}*c*per unit. The numbers of

_{ID}*I*units produced per infected neuron (

_{D}*I*) and per infected Mϕ (

*I*

_{Mϕ}) upon their death are α

_{I}and α

_{IMϕ}, respectively. Since TMEV production of infectious viral particles is markedly reduced in infected Mϕ, due to a block in viral assembly, we did not consider virus release from infected Mϕ.

Experimental data.We fit our model to our previously published data on BeAn strain TMEV RNA levels measured in susceptible SJL/J mouse spinal cords (Fig. 2) and serum virus-specific antibody titers (Fig. 3) (9). In some infected mice, typical acute/chronic infection symptoms were not observed, as the viral RNA (vRNA) level was probably too low to generate symptoms. This may be due to the low infection titers used in that study. In general, with low infection titers or low viral loads, 10 to 15% of mice either do not develop persistent infection or have delayed kinetics of persistent infection by several months (9, 38). As such, we excluded the data points corresponding to <700 copy equivalents/μg of total spinal cord RNA for simulation and model fits (i.e., the cutoff observed in Fig. 5 in reference 9). The term “copy equivalents” is used because viral genome copies are measured by real-time RT-PCR for short viral genome sequences, which are assumed to represent the full-length viral genome (∼8,100 nucleotides). Since a mouse is sacrificed to measure its vRNA level, there are no longitudinal kinetic data for an individual mouse. Because there is variation between mice, the data are noisy, and thus we used the mean of log_{10} vRNA copy equivalents/mg of spinal cord RNA at each measurement time. To quantify the variation, we calculated the standard deviation (SD) of log_{10} vRNA copy equivalents. Some measurement time points included four vRNA measurements. When this was the case, the SD of log_{10} vRNA copy equivalents was generally higher than 0.8. When there was only one vRNA measurement at a given time point, or two vRNA measurements at a given time point that differ by less than 0.8, we used 0.8 as the default SD. The means and SDs of log_{10} vRNA copy equivalents are shown in Fig. 2.

Time-dependent immune responses.We assume that the measured serum virus-specific antibody titer (Fig. 3) reflects the total virus-specific antibody in the spinal cord (36). In order to estimate *h*(*t*) in equations 1 and 2, i.e., the effect of virus-specific antibody in reducing TMEV infection in the spinal cord, we fit the virus-specific antibody levels shown in Fig. 3 by
*t* that is <5, *A _{b}*(

*t*) = 0 (since virus-specific antibody is only detected after day 5), κ is the maximum antibody expansion rate, and

*A*

_{b}^{max}is the maximum antibody level in the spinal cord. We used Mathematica 8.0 to find the values of κ and

*A*

_{b}^{max}that generated the best fit curve of the virus-specific antibody data (Fig. 3).

We assume that the main effect of antibodies is to block the infection of susceptible cells, with blocking efficacy determined as
*K _{A}* is the association constant between virus-specific antibody and free virions in the spinal cord, as previously used (39).

The influx of monocytes into the CNS results from virus-specific Th1 and Th17 T lymphocyte responses leading to cytokine and chemokine recruitment of monocytes into the CNS, which then differentiate into susceptible Mϕ (40). Thus, the susceptible Mϕ source rate in the model is
*t*_{M} is time of immune-mediated Mϕ influx, *s*_{0} is baseline source rate, and *s*_{+} is source rate after *t*_{M}. Since before *t _{M}* there is no immune-mediated influx of susceptible Mϕ, we set

*s*

_{0}at 0. The time of immune-mediated Mϕ influx,

*t*

_{M}, is about 2 weeks p.i. (9, 37) and was estimated here while fitting the model to the experimental data.

Estimation of fixed parameters and fitting procedures.The number of motor neurons in the spinal cord before infection is approximately 30,000 (41). The spinal cord contains ∼100 μg of total RNA. The first vRNA data, collected on day 1 p.i. (9), were 10^{3} copy equivalents/μg or 10^{5} copies in the whole spinal cord, so we started the simulations with a value for *V* of 10^{5} copies on day 1 p.i., assuming there are no infected neurons, susceptible Mϕ, infected Mϕ, or infectious debris at the beginning of a simulation. An infected neuron contains about 10^{5} copy equivalents of viral RNA (42), and each infected Mϕ also contains about 10^{5} copy equivalents of viral RNA (18, 37). For simplicity, we assumed here that this level is reached immediately once a neuron or a macrophage is infected, since vRNA reaches its intracellular maximum levels (10^{5}) within 4 to 8 h p.i. *in vitro* (43). The loss rate of uninfected Mϕ in mouse spinal cord is unknown. Here we adopted a loss rate constant of 0.1/day, measured in the irises of rat eyes (44), in our calculations. We also assumed that infected Mϕ die faster than uninfected ones, i.e., δ_{IMϕ} > δ_{Mϕ}, as estimated *in vitro* (45). Trottier et al. (37) found that TMEV-infected cell lines such as BHK-21 (baby hamster kidney cells) and M1-D (macrophages) produce 3.5 × 10^{5} and 2.6 × 10^{5} copy equivalents of vRNA/cell, respectively. Assuming that neurons produce slightly higher numbers of viral genomes than BHK-21 or M1-D cells (e.g., because neurons are larger cells than BHK-21), we chose a value for *p _{n}* of 5 × 10

^{5}δ

_{n}. We also confirmed that changing 5 × 10

^{5}to 1 × 10

^{5}copy equivalents/cell leads to a negligible change in the results. Since the system is finite, if the number of cells of a particular type or infectious debris units in a mouse spinal cord is less than 1, it may generate unrealistic results. To avoid such an artifact, we treated the change of the number as 0 if the change was less than a cutoff value of Δ

*t*/10 days, where Δ

*t*is the time step in the simulations and 10 days is the approximate time of the acute phase or transition phase in our simulations. In other words, if the time step was 1 day, then if the number of cells changed by less than 0.1 during that time step, the change was assumed to be negligible. Further, if the rate of change was constant, this would imply that the change in number was less than 1 over the full 10-day acute phase. This cutoff value differs from 1, allowing accumulation of cell and debris numbers during the small time steps in the simulation.

Since total vRNA was measured in the spinal cord, the total vRNA computed by the model is
^{5} vRNA copies (42), an infected Mϕ contains 10^{5} vRNA copies (18, 37), and the average number of vRNA copies a unit of infectious debris contains is defined as *n _{D}*. Flaccid paralysis was not observed in the infected mice (9). Although the exact association between flaccid paralysis symptoms and the magnitude of damaged/lost sensory and motor neurons is not known, we considered only simulation runs where the total loss of motor neurons due to acute TMEV infection in SJL mice did not exceed 50%.

Monte Carlo filtering.To estimate a reasonable range of the unknown model parameters, we randomly chose many sets of parameters and calculated *V*(*t*) numerically. We claim a set of parameters is good if the predicted *V*(*t*) generated a small enough weighted root mean squared error when compared with the average measured values (*V̂*) i.e.,
_{10} *V*(*t _{i}*) is the predicted values, log

_{10}

*V̂*(

*t*) is the average of experimental values, σ(

_{i}*t*) is the SD of the data at

_{i}*t*, and

_{i}*n*is the number of data points to fit. We used Latin hypercube sampling; i.e., the logarithmic values of parameters were first chosen randomly from an assigned interval, e.g., [0, 5]. We performed a large number (>10

^{4}) of simulations to collect histograms of parameter values that give good fits of the data in equation 11. If the histogram of a given parameter is bell-shaped, the range of the parameter we report (Table 1) is given by the half-peak positions in the histogram (Fig. 4). There are three parameters for which we assign knowledge-based intervals. The interval [−1, 1] was chosen for log

_{10}δ

_{IMϕ}, because we expect that δ

_{IMϕ}is slightly larger than the death rate constant of uninfected Mϕ (δ

_{Mϕ}) of 0.1/day. The interval [0, 3] is chosen for log

_{10}α

_{I}and log

_{10}α

_{IMϕ}, where α

_{I}and α

_{IMϕ}are the number of debris units created upon the death of an infected neuron and Mϕ, respectively. The number of vRNA copies in a debris unit (

*n*) is

_{D}*f*× 10

^{5}/α

_{IMϕ}, where

*f*, chosen from among [1, 0.5, 0.1, 0.01], is the average fraction of the ∼10

^{5}vRNA copies in an infected Mϕ that remain intact upon loss of the cell. The logarithmic values of the other parameters, β

_{n},

*c*, δ

_{n},

*K*,

_{A}*t*

_{M},

*s*

_{0}, β,

*c*, and

_{ID}*s*

_{+}, were chosen randomly from assigned large intervals in the first round as described in Table 1. To optimize the choice of parameter intervals in the next rounds of simulations, we adjusted the parameter intervals in every round, aiming to determine the ranges of parameters more precisely. For instance, if the peak of the histogram fell out of the chosen parameter interval, we moved the interval to cover the peak; if the half peak range occupied only a small part of the interval, we shrank the interval. Because there are many unknown parameters, we estimated their values in two steps. The viral dynamics during the acute phase (i.e.,

*t*<

*t*) is controlled by parameters β

_{M}_{n},

*p*

_{n},

*c*, δ

_{n}, and

*K*

_{A}(Fig. 1). With the chosen parameter intervals, we performed 1.7 × 10

^{6}simulations of the model up to 11 days p.i. in order to collect 2 × 10

^{5}sets of the five acute phase parameters mentioned above that generated good fits with the expression

*J*< 0.55. Here we used a lower cutoff than in equation 11, because the data during acute phase are less noisy than those during the chronic phase. The estimation of these five unknown parameters concluded the first step of model fit with the data. In the second step, we used these 2 × 10

^{5}sets of acute-phase parameters and varied the other model parameters (β,

*c*, and

_{ID}*s*

_{+}) randomly to look for parameter sets that fit the overall data with good quality. We ran 7.0 × 10

^{4}simulations in order to collect 6,000 sets of parameters that met the fitting criteria (equation 11). As a result, we produced a histogram of 6,000 values for each parameter (Fig. 4).

## RESULTS

The model described by equations 1 to 9 and Fig. 1 characterize the dynamic transition of TMEV infection from acute to chronic phases in the mouse spinal cord. We fit the model to published data (9) of vRNA copy equivalents in mouse spinal cords (Fig. 2) and serum virus-specific antibody titers (Fig. 3). Using a Monte Carlo filtering approach, as described in Materials and Methods, we obtained a range of parameters that generated good fits of the model to the data; an example of a good fit is shown in Fig. 2, and the histograms of representative parameters that can generate fits of good quality are shown in Fig. 4. The parameter ranges obtained are summarized in Table 1. Importantly, some key parameters fall within a narrow range, as described below.

Rate of infection of neurons by free TMEV particles.During the first 5 days of infection, before virus-specific antibodies were detected, viral RNA increased steadily due to an increase in infection of neurons. The rate constant that describes neuron infection, β_{n}, was found to be constrained to a limited range, 0.5 × 10^{−8} to 5.6 × 10^{−8} day^{−1}. A higher or lower value of β_{n} fails to fit the data. To get an intuitive understanding of this rate constant, consider what happens on day 1 p.i., when there are approximately 10^{5} virions in the spinal cord, which contains approximately 30,000 motor neurons. Using this value of β_{n}, an infection event occurs every 0.14 to 1.6 h on average. Thus, the model predicts that the initial infection events are slow but then occur more rapidly as the viral load increases.

Effectiveness of TMEV-specific antibodies in blocking infection.Between days 6 and 11 p.i., the average total vRNA level decreased (Fig. 2), which is explained in the model to be due to blocking of new virus infection of motor neurons by virus-specific antibodies and the TMEV-induced cytolytic loss of already-infected neurons. To estimate the virus-specific antibody kinetic parameters during TMEV infection, we obtained the best fit of equation 7 to the virus-specific level shown in Fig. 3. We estimate that the antibody levels increased after day 5 p.i. at a rate (κ) of 1.2/day and reached a maximum titer of *A _{b}*

^{max}= 10

^{4.2}after day 15 p.i. From the Monte Carlo filtering,

*K*is in the range 0.4 × 10

_{A}^{−2}to 5.0 × 10

^{−2}/antibody/spinal cord volume. From equation 8, this range of

*K*corresponds to a value of η of 98.5 to 99.9% reduction of the infection rate after day 15 p.i.

_{A}Death rate and virion production rate of infected neurons.The death rate constant of infected neurons (δ_{n}) was found to be in the range 0.6 day^{−1} to 7.1 day^{−1}, which corresponds to a half-life of 0.1 to 1.2 days. From this quantity, we deduce that the infected neurons produce virions at a rate (*p _{n}*) of 0.3 × 10

^{6}to 3.5 × 10

^{6}virions/cell/day.

Timing of Mϕ source activation.The TMEV-immune-mediated activation of the Mϕ source rate (*t _{M}*) is estimated to occur between 10.9 and 12.9 days p.i. This is in agreement with the observation that the vRNA level was lowest at day 11 p.i. and became significantly higher at day 13. The increased Mϕ source rate provides more target cells, and hence vRNA would be expected to increase.

Activated source of susceptible Mϕ.The total viral RNA reaches a steady state after about 15 days p.i. (Fig. 2). The simulation results suggest that the source rate of susceptible Mϕ (*s*_{+}) is in the range 25 to 350 macrophages/day. Assuming that the infection of susceptible Mϕ is much faster than death/loss of Mϕ, we expect that almost all susceptible Mϕ become infected; hence, the source rate of infected Mϕ is also *s*_{+}. Since an infected Mϕ dies at rate δ_{IMϕ}, the steady-state number of infected Mϕ is *s*_{+}/δ_{IMϕ}. At steady state, there are ∼10^{7} vRNA molecules in the spinal cord, which corresponds to an *s*_{+}/δ_{IMϕ} of ∼10^{2} infected Mϕ assuming each infected Mϕ contains 10^{5} copy equivalents of vRNA (18, 37). If *s*_{+}/δ_{IMϕ} is 10^{2} and δ_{IMϕ} is 0.1 to 10 day^{−1}, the source of Mϕ needs to be 10 to 1,000 day^{−1}, in agreement with the result for *s*_{+} of 25 to 350 day^{−1}.

Parameters without narrowed-range estimates.Some parameters are not estimated well by our model and simulations. The fraction of vRNA that remains in infectious debris, *f*, is chosen from [1, 0.5, 0.1, 0.01], but all four values have similar likelihood based on fitting the model to the data. We find large ranges (over 4 log units) of possible values for the infection rate constant β of Mϕ by infectious debris via phagocytosis and for the clearance rate constant *c _{ID}* of infectious debris. We also find that significant correlation exists between these two rate constants (Fig. 5). Obviously, a higher clearance rate of infectious debris can be compensated for by a higher rate at which infectious debris infects Mϕ; therefore, the values for one can span a large range as long as the other value changes accordingly.

## DISCUSSION

TMEV infection in the mouse has been extensively studied for more than 3 decades. In order to provide new insights into the virus-host dynamics during the acute and early chronic phases, we developed a mathematical model of TMEV infection and estimated unknown viral and host parameters. These parameter estimates include the neuron infection rate (β_{n}), which is 0.5 × 10^{−8} to 5.6 × 10^{−8} day^{−1}, the reduction of infection by virus-specific antibodies (η), which reaches 98.5 to 99.9% after day 15 p.i., the half-life of infected neurons, which is 0.1 to 1.2 days, and the immune-mediated Mϕ source rate, which is 25 to 350 macrophages/day starting at 10.9 to 12.9 days p.i.

In contrast to a chronic infection by a noncytolytic RNA virus (such as HCV), in which the host cells may live for days or weeks (27), persistence of cytolytic RNA viruses (in which host cells die rapidly, i.e., on a time scale of hours) requires dynamic cell-to-cell infection in the presence of host immunity in order to maintain the infection (46). Here we predict that virus-specific antibodies efficiently block cell-free infection by more than 98.4%. This model prediction may be supported by a prior study that indicated the presence of virus-neutralizing antibodies (36), but further studies are needed to verify that neutralizing-antibody kinetics approximate the measured virus-specific antibodies modeled here. However, we cannot rule out the possibility that virus-specific antibodies also enhance the clearance of virus, as has been observed with antibodies to HBV (33). Indeed, model simulations cannot distinguish between enhanced virus clearance and the blocking of infection by virus-specific antibodies (data not shown). Together with a significant reduction of infectious virus production by infected Mϕ compared to that by infected motor neurons, TMEV hypothetically needs to have an alternative route of infection in the presence of adaptive immune responses, such as phagocytosis of infectious debris, which is modeled here.

There is no available estimate of the influx rate of Mϕ into the spinal cord (or their basal number) before infection or during the acute phase (i.e., before a *t _{M}* of ∼11 days p.i.). If their basal influx number is >2 Mϕ/day and they are susceptible to infection and productive, as are immune-recruited Mϕ that enter the spinal cord after

*t*, then our model simulations predict that the transient viral decline observed during the transition from acute to chronic phases will disappear (data not shown). It is plausible that most susceptible Mϕ have already died during the acute phase and that without an enhanced influx rate of new Mϕ, their number is very low at the time of transition to the chronic phase. Thus, for simplicity, we assumed here that no Mϕ susceptible to TMEV infection were present (and that their influx rate [

_{M}*s*

_{0}] is 0/day) during the acute phase. Interestingly, Rossi et al. (8) selectively depleted Mϕ in infected SJL mice by intravenously administering mannosylated liposomes containing dichloromethylene diphosphonate, which are phagocytosed by Mϕ and then kill them, at various times from days 7 to 19 p.i. when macrophage infiltration of the CNS occurs and then sacrificing the mice on day 21 p.i. No viral RNA was detected by hybridization in 70% of treated mice (

*n*= 16), but viral RNA was detected in 100% of untreated mice (

*n*= 13), suggesting that Mϕ infiltration into the CNS early in the chronic phase plays an important role in viral persistence.

Acute TMEV infection is initiated by intracerebral inoculation of mice. The infection is characterized by a transition from an acute neural stage in the gray matter (brain and spinal cord) to a chronic macrophage and glial infection in the white matter (spinal cord) (13). After 30 to 60 days, high viral RNA and proinflammatory cytokine mRNA levels were observed only in susceptible SJL mice that developed demyelinating disease (9). Therefore, we modeled the CNS site, i.e., the spinal cord, where TMEV infection takes place both acutely and then chronically.

We assumed in the model that motor neurons are the primary site of acute infection, due to a prior study in which we showed that more virus is present in anterior horn (motor) neurons than in sensory neurons (36). In addition, although we modeled viral persistence in Mϕ, which is considered the main site of persistence, TMEV has been found also to replicate in other cells, such as oligodendrocytes (47–49) and astrocytes (47). However, since the dynamics of TMEV infection in these cells and the type of virus-induced death (apoptosis versus necrosis) are still under investigation, the role of these cells in TMEV persistence is unknown. It is anticipated that future studies will shed light on the role of sensory neurons, oligodendrocytes, and astrocytes, and the model here will be modified accordingly.

During the acute phase, the majority of CNS-infiltrating CD8^{+} T cells in susceptible mice are TMEV capsid specific (10, 12, 50). Interestingly, Lyman et al. (12) showed that at day 8 p.i., ≤1% of splenic CD8^{+} T cells are reactive to any of the TMEV-capsid epitopes in susceptible SJL/J mice, compared to 12% in resistant C57BL/6 mice. The lower level of peripheral virus-specific CD8^{+} T cells in susceptible SJL/J mice resulted in >3-fold more infiltrating virus-specific CD8^{+} T cells in the CNS in resistant C57BL/6 mice than in susceptible mice. As CD8^{+} T cells can kill TMEV-infected neurons, our estimate of the short infected-neuron half-life (0.1 to 1.2 days) may not be due solely to TMEV cytopathic effects but may also include effects of infiltrating virus-specific CD8^{+} T cells. The limited amount of data obtained to date on the temporal relationship between virus-specific CD8^{+} T cells and viral RNA copies during the acute and chronic phases of infection precluded incorporating the role of CD8^{+} T cells into the current mathematical model. In future studies, it will be necessary to fully characterize and model both TMEV RNA and virus-specific CD8^{+} T cell kinetics during the transition from the acute to the chronic phase in both susceptible and resistant mice.

In summary, we have developed a mathematical model for TMEV during the acute and early chronic phases of infection in a susceptible strain of mice. The model reproduces the virus-host dynamics during the transition from the acute to the chronic phase of infection and allowed us to estimate, for the first time, dynamic parameters of TMEV infection. Since the mouse model of TMEV infection provides an experimental analog of MS in humans, which is a chronic immune-mediated demyelinating disease of the CNS, the next step will be to incorporate later virus-host events (beyond ∼30 days) to create one comprehensive model of the complete TMEV infection process, including the clinical expression of demyelinating disease. The model developed here (and the future comprehensive one) can serve as an important tool in understanding TMEV kinetics, the antiviral immune response(s) (9, 51), and demyelinating disease during therapeutic modalities (38, 52–58), and it may prove useful in designing and evaluating new drugs for use against MS, should a viral cause be established.

## ACKNOWLEDGMENTS

This work was performed under the auspices of the U.S. Department of Energy and was supported by NIH grants R56/R01-AI078881, P20-GM103452, AI028433, NS065945, and OD011095, as well as NSF grant PHY11-25915.

## FOOTNOTES

- Received 18 December 2012.
- Accepted 23 January 2013.
- Accepted manuscript posted online 30 January 2013.

- Copyright © 2013, American Society for Microbiology. All Rights Reserved.