Oscillations in a Spatial Oncolytic Virus Model.

Virotherapy treatment is a new and promising target therapy that selectively attacks cancer cells without harming normal cells. Mathematical models of oncolytic viruses have shown predator-prey like oscillatory patterns as result of an underlying Hopf bifurcation. In a spatial context, these oscillations can lead to different spatio-temporal phenomena such as hollow-ring patterns, target patterns, and dispersed patterns. In this paper we continue the systematic analysis of these spatial oscillations and discuss their relevance in the clinical context. We consider a bifurcation analysis of a spatially explicit reaction-diffusion model to find the above mentioned spatio-temporal virus infection patterns. The desired pattern for tumor eradication is the hollow ring pattern and we find exact conditions for its occurrence. Moreover, we derive the minimal speed of travelling invasion waves for the cancer and for the oncolytic virus. Our numerical simulations in 2-D reveal complex spatial interactions of the virus infection and a new phenomenon of a periodic peak splitting. An effect that we cannot explain with our current methods.

Homogenization of a reaction diffusion equation can explain influenza A virus load data.

We study the influence of spatial heterogeneity on the antiviral activity of mouse embryonic fibroblasts (MEF) infected with influenza A. MEF of type Ube1L-/- are composed of two distinct sub-populations, the strong type that sustains a strong viral infection and the weak type, sustaining a weak viral load. We present new data on the virus load infection of Ube1L-/-, which have been micro-printed in a checker board pattern of different sizes of the inner squares. Surprisingly, the total viral load at one day after inoculation significantly depends on the sizes of the inner squares. We explain this observation by using a reaction diffusion model and we show that mathematical homogenization can explain the observed inhomogeneities. If the individual patches are large, then the growth rate and the carrying capacity will be the arithmetic means of the patches. For finer and finer patches the average growth rate is still the arithmetic mean, however, the carrying capacity uses the harmonic mean. While fitting the PDE to the experimental data, we also predict that a discrepancy in virus load would be unobservable after only half a day. Furthermore, we predict the viral load in different inner squares that had not been measured in our experiment and the travelling distance the virions can reach after one day.