Modelling and analysis of dynamics of viral infection of cells and of interferon resistance

Interferons are active biomolecules, which help fight viral infections by spreading from infected to uninfected cells and activate effector molecules, which confer resistance from the virus on cells. We propose a new model of dynamics of viral infection, including endocytosis, cell death, production of interferon and development of resistance. The novel element is a specific biologically justified mechanism of interferon action, which results in dynamics different from other infection models. The model reflects conditions prevailing in liquid cultures (ideal mixing), and the absence of cells or virus influx from outside. The basic model is a nonlinear system of five ordinary differential equations. For this variant, it is possible to characterise global behaviour, using a conservation law. Analytic results are supplemented by computational studies. The second variant of the model includes age-of-infection structure of infected cells, which is described by a transport-type partial differential equation for infected cells. The conclusions are: (i) If virus mortality is included, the virus becomes eventually extinct and subpopulations of uninfected and resistant cells are established. (ii) If virus mortality is not included, the dynamics may lead to extinction of uninfected cells. (iii) Switching off the interferon defense results in a decrease of the sum total of uninfected and resistant cells. (iv) Infection-age structure of infected cells may result in stabilisation or destabilisation of the system, depending on detailed assumptions. Our work seems to constitute the first comprehensive mathematical analysis of the cell-virus-interferon system based on biologically plausible hypotheses.

[1]  J. Hale Functional Differential Equations , 1971 .

[2]  O. Diekmann Mathematical Epidemiology of Infectious Diseases , 1996 .

[3]  S. Baron,et al.  IFN-beta mediates coordinate expression of antigen-processing genes in RSV-infected pulmonary epithelial cells. , 2001, American journal of physiology. Lung cellular and molecular physiology.

[4]  Allan R. Brasier,et al.  Retinoic Acid-Inducible Gene I Mediates Early Antiviral Response and Toll-Like Receptor 3 Expression in Respiratory Syncytial Virus-Infected Airway Epithelial Cells , 2006, Journal of Virology.

[5]  G. Stépán Retarded dynamical systems : stability and characteristic functions , 1989 .

[6]  Alan S. Perelson,et al.  Viral and Latent Reservoir Persistence in HIV-1–Infected Patients on Therapy , 2006, PLoS Comput. Biol..

[7]  R. Garofalo,et al.  Ribavirin Treatment Up-Regulates Antiviral Gene Expression via the Interferon-Stimulated Response Element in Respiratory Syncytial Virus-Infected Epithelial Cells , 2003, Journal of Virology.

[8]  John Yin,et al.  Quantifying Viral Propagation in Vitro: Toward a Method for Characterization of Complex Phenotypes , 2001, Biotechnology progress.

[9]  O. Diekmann,et al.  The Dynamics of Physiologically Structured Populations , 1986 .

[10]  James B. Rawlings,et al.  Dynamics of viral infections: incorporating both the intracellular and extracellular levels , 2005, Comput. Chem. Eng..

[11]  Horst R. Thieme,et al.  Mathematics in Population Biology , 2003 .

[12]  Urszula Foryś,et al.  BIOLOGICAL DELAY SYSTEMS AND THE MIKHAILOV CRITERION OF STABILITY , 2004 .

[13]  John Yin,et al.  Arrested spread of vesicular stomatitis virus infections in vitro depends on interferon-mediated antiviral activity. , 2005, Biotechnology and bioengineering.

[14]  A. Perelson,et al.  Viral dynamics and response differences in HCV‐infected African American and white patients treated with IFN and ribavirin , 2003, Hepatology.

[15]  S. Lunel,et al.  Delay Equations. Functional-, Complex-, and Nonlinear Analysis , 1995 .

[16]  E. Coddington,et al.  Theory of Ordinary Differential Equations , 1955 .