Abstract
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