May 2004
Volume 45, Issue 13
ARVO Annual Meeting Abstract  |   May 2004
A Mathematical Model of Acute Herpes Simplex Virus Corneal Infection
Author Affiliations & Notes
  • H.W. Thompson
    Ophthalmology, LSU Eye Center, New Orleans, LA
  • D. Stark
    Ophthalmology, LSU Eye Center, New Orleans, LA
  • J.M. Hill
    Ophthalmology, LSU Eye Center, New Orleans, LA
  • Footnotes
    Commercial Relationships  H.W. Thompson, None; D. Stark, None; J.M. Hill, None.
  • Footnotes
    Support  NEI06311 (JMH), EY02377 (LSU Eye Center Core Grant)
Investigative Ophthalmology & Visual Science May 2004, Vol.45, 1660. doi:
  • Views
  • Share
  • Tools
    • Alerts
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      H.W. Thompson, D. Stark, J.M. Hill; A Mathematical Model of Acute Herpes Simplex Virus Corneal Infection . Invest. Ophthalmol. Vis. Sci. 2004;45(13):1660.

      Download citation file:

      © ARVO (1962-2015); The Authors (2016-present)

  • Supplements

Abstract: : Purpose: To devise a mathematical model of corneal cell turnover mimicking corneal cell dynamics and to extend this model with further equations to model viral infection, proliferation, and clearance to reproduce observed dynamics of viral plaque forming unit (PFU) production. Methods: The mathematical model is a non–linear system of six differential equations with terms that describe the corneal epithelium’s response to wounding and viral infection. Three equations for cellular dynamics of the corneal epithelium without viral infection describe 1) the rate of cell proliferation p(t), 2) the rate of quiescent cell production q(t), and 3) the combined amount of growth factors c(t). The equations are 1) dp(t)/dt=a1p(t)c(t)–a2p(t)c(t)+s(c), where s(c) is the term for stem cell mitosis and it =a3c(t) if c(t)<kc0 and s(c)=a3kc0 if c(t)≥kc0; 2) dq(t)/dt=a2p(t)c(t)–a5q(t)c(t); and 3) dc(t)/dt=a7a9c(t)(q(t)+p(t))–a10c(t). At equilibrium: a1p0c0–a2p0c0+a3c0=0; a2p0c0–a5q0c0=0; and a7–a9c0(q0+p0)–a10c0=0. By using the equilibrium equations and estimates from the literature for p0, q0 and c0,we obtain a1, a2, a3, a5, a7, a9, and a10. The system of equations behaves as expected if cells are removed (corneal wounding); proliferation (healing) continues until a stable number of cells is reached (equilibrium reestablished). Three further equations are added to the system to model the dynamics of viral infection: 1) an equation for the rate of cells infected but not yet producing virus di(t)/dt=a4p(t)v(t)+a6q(t)v(t)–a11i(t); 2) an equation for cells actively producing virus dj(t)/dt=a11i(t)–a12j(t); and 3) an equation for virus (PFU) in the tears dv(t)/dt=a13j(t)–a14p(t)v(t)–a15q(t)v(t)–a16tv(t)–a17v(t). Results: Given a range of viral inoculation concentrations, the dynamics of PFU observed in experimental data are closely duplicated by the model predictions. Extensive graphic analysis of the model predictions for experimental observations will be presented. Conclusions: Numerical analytic solutions of non–linear equation systems can duplicate the behavior of the corneal epithelium. This model can make accurate predictions under various conditions of PFU load and differences in viral replication efficiency.

Keywords: herpes simplex virus • cornea: epithelium • computational modeling 

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.