A critical virus production rate for blow-up suppression in a haptotaxis model for oncolytic virotherapy

Youshan Tao,Michael Winkler

DOI: https://doi.org/10.1016/j.na.2020.111870

2020-01-01

Abstract:In a two-dimensional smoothly bounded domain, we consider the oncolytic virotherapy model given by ut=Δu−∇⋅(u∇v)−ρuz,vt=−(u+w)v,wt=DwΔw−w+uz,zt=DzΔz−z−uz+βw,with parameters Dw>0,Dz>0,β>0 and ρ≥0.

What problem does this paper attempt to address?

A Critical Virus Production Rate for Efficiency of Oncolytic Virotherapy

Youshan Tao,Michael Winkler

DOI: https://doi.org/10.1017/s0956792520000133

2021-01-01

European Journal of Applied Mathematics

Abstract:In a planar smoothly bounded domain$\Omega$, we consider the model for oncolytic virotherapy given by$$\left\{ \begin{array}{l} u_t = \Delta u - \nabla \cdot (u\nabla v) - uz, \\[1mm] v_t = - (u+w)v, \\[1mm] w_t = d_w \Delta w - w + uz, \\[1mm] z_t = d_z \Delta z - z - uz + \beta w, \end{array} \right.$$with positive parameters$ D_w $,$ D_z $and$\beta$. It is firstly shown that whenever$\beta \lt 1$, for any choice of$M \gt 0$, one can find initial data such that the solution of an associated no-flux initial-boundary value problem, well known to exist globally actually for any choice of$\beta \gt 0$, satisfies$$u\ge M \qquad \mbox{in } \Omega\times (0,\infty).$$If$\beta \gt 1$, however, then for arbitrary initial data the corresponding is seen to have the property that$$\liminf_{t\to\infty} \inf_{x\in\Omega} u(x,t)\le \frac{1}{\beta-1}.$$This may be interpreted as indicating that$\beta$plays the role of a critical virus replication rate with regard to efficiency of the considered virotherapy, with corresponding threshold value given by$\beta = 1$.
Asymptotic stability of spatial homogeneity in a haptotaxis model for oncolytic virotherapy

Youshan Tao,Michael Winkler

DOI: https://doi.org/10.1017/prm.2020.97

2022-01-01

Abstract:This study considers a model for oncolytic virotherapy, as given by the reaction-diffusion-taxis system {u(t) - Delta u - del(u del v) - rho uz, v(t) = -(u + w)v, w(t) = D-w Delta w - w + uz, z(t) = D-z Delta z - z - uz + beta w, in a smoothly bounded domain Omega subset of R-2, with parameters D-w > 0, D-z > 0, beta > 0 and rho >= 0. Previous analysis has asserted that for all reasonably regular initial data, an associated no-flux type initial-boundary value problem admits a global classical solution, and that this solution is bounded if beta < 1, whereas whenever beta > 1 and (1/vertical bar Omega vertical bar) integral(Omega) u(., 0) > 1/(beta-1), infinite-time blow-up occurs at least in the particular case when rho = 0. In order to provide an appropriate complement to this, the current study reveals that for any rho >= 0 and arbitrary beta > 0, at each prescribexl level gamma is an element of (0,1/(beta - 1)+) one can identify an L-infinity-neighbourhood of the homogeneous distribution (u, v, w, z) (gamma, 0, 0, 0) within which all initial data lead to globally bounded solutions that stabilize towards the constant equilibrium (u(infinity), 0 ,0, 0, 0) with some u(infinity) > 0.
Global Boundedness and Asymptotic Behavior in a Double Haptotaxis Model for Oncolytic Virotherapy

Xueyan Tao,Shulin Zhou

DOI: https://doi.org/10.1016/j.jmaa.2024.128522

IF: 1.417

2024-01-01

Journal of Mathematical Analysis and Applications

Abstract:This work studies a haptotactic cross -diffusion system modeling oncolytic virotherapy {u(t) = Delta u - V ( u V v ) - puz, v(t) = - (u + w) v - kappa v, w(t) = D-w Delta w-xi((sic))(w)(w v ) - w + puz, z(t) = D-z Delta z - z - puz + beta w, in a smooth bounded domain Omega subset of R-2 , with parameters D-w > 0, D-z > 0, xi(w) >= 0, p >= 0, kappa > 0 and beta >= 0. This system describes the interaction among uninfected and infected cancer cells, extracellular matrix and oncolytic viruses. It is rigorously proved that an associated no -flux initial -boundary value problem has a unique global classical solution which is bounded in ( L infinity (Omega)) 4 . Moreover, it is shown that this bounded solution can approach a spatially constant equilibrium in the large time limit under the additional assumption that 0 <= beta < 1. (c) 2024 Elsevier Inc. All rights are reserved, including those for text and data mining, AI training, and similar technologies.
Global Classical Solutions to a Doubly Haptotactic Cross-Diffusion System Modeling Oncolytic Virotherapy

Youshan Tao,Michael Winkler

DOI: https://doi.org/10.1016/j.jde.2019.10.046

IF: 2.615

2020-01-01

Journal of Differential Equations

Abstract:This work studies a haptotaxis system proposed as a model for oncolytic virotherapy, accounting for interaction between uninfected cancer cells, infected cancer cells, extracellular matrix (ECM) and oncolytic virus. In addition to random movement, both uninfected and infected tumor cells migrate haptotactically toward higher ECM densities; moreover, besides degrading the non-diffusible ECM upon contact the two cancer cell populations are subject to an infection-induced transition mechanism driven by virus particles which are released by infected cancer cells, and which assault the uninfected part of the tumor.The main results assert global classical solvability in an associated initial-boundary value problem posed in one- or two-dimensional domains with any given suitably regular initial data. This is achieved by discovering a quasi-Lyapunov functional structure that allows to appropriately cope with the presence of nonlinear zero-order interaction terms which apparently form the most significant additional mathematical challenge of the considered system in comparison to previously studied haptotaxis models.
A parabolic-hyperbolic free boundary problem modelling tumor treatment with virus

Youshan Tao,Hui Zhang

DOI: https://doi.org/10.1142/S0218202507001838

2011-01-01

Abstract:We consider a procedure for cancer therapy which consists of injecting replication-competent viruses into the tumor. The viruses infect tumor cells, replicate inside them, and eventually cause their death. As infected cells die, the viruses inside them are released and then proceed to infect adjacent tumor cells. However, a major factor influencing the efficacy of virus agents is the immune response that may limit the replication and spread of the replication-competent virus. The competition between tumor cells, a replication-competent virus and an immune response is modelled as a free boundary problem for a nonlinear system of partial differential equations, where the free boundary is the surface of the tumor. In this model, the immune response equation is a semilinear parabolic equation, including a chemotaxis term which is used to describe the movement of the immune response induced by gradients of the infected cell density. Under the assumption that the chemotactic sensitivity coefficient is small compared with the diffusion coefficient of the immune response, we prove the global existence and uniqueness of the solution of this free boundary problem. For large chemotactic coefficient, the global existence is still open.
Simulation of a Model of Tumors with Virus-therapy

Youshan Tao,Qian Guo

DOI: https://doi.org/10.1007/978-3-7643-7719-9_42

2007-01-01

Abstract:We consider a procedure for cancer therapy which consists of injecting replication-competent viruses into the tumor. The viruses infect tumor cells, replicate inside them, and eventually cause their death. As infected cells die, the viruses inside them are released and then proceed to infect adjacent tumor cells. However, a major factor influencing the efficacy of virus agents is the immune response that may limit the replication and spread of the replication-competent virus. The immune response is cytokine-mediated. The expression of viruses in tumor cells sensitize cells to lysis by the TNF (tumor necrosis factor) cytokine. The competition between tumor cells, a replication-competent virus and an immune response is modelled as a free boundary problem for a nonlinear system of partial differential equations, where the free boundary is the surface of the tumor. In this model, the immune response equation is a non-standard parabolic equation due to the chemotaxis (spatial gradients of diffusible chemicals) of the immune response. The purpose of this paper is to give the numerical methods for solving this kind of free boundary problems. Several simulation results are also given.
Critical Mass for Infinite-Time Blow-Up in a Haptotaxis System with Nonlinear Zero-Order Interaction

Youshan Tao,Michael Winkler

DOI: https://doi.org/10.3934/dcds.2020216

2021-01-01

Discrete and Continuous Dynamical Systems

Abstract:We consider the haptotaxis system {u(t = )Delta u - del . (u del u), u(t )= -(u+w)u, w(t) = D-w Delta w - w + uz, z(t) = D-z Delta z - z - uz + beta w, which arises as a simplified version of a recently proposed model for oncolytic virotherapy. When posed under no-flux boundary conditions in a smoothly bounded domain Omega subset of R-2, with positive parameters D-w, D-z and beta, and along with initial conditions involving suitably regular data, this system is known to admit global classical solutions. It is shown that with respect to infinite-time blow-up, this system exhibits a critical mass phenomenon related to the quantity m(e) := 1/(beta-1)(+) In fact, it is seen that each solution fulfilling 1/vertical bar Omega vertical bar f(Omega) u(., 0) > m(c) must be unbounded, and this is complemented by a boundedness result which inter alia asserts that for any choice of m < m(c) one can find a nontrivial set of solutions, particularly containing spatially heterogeneous solutions, each of which is bounded though satisfying 1/vertical bar Omega vertical bar f(Omega) u(.,0) = m.
Dampening effects on global boundedness and asymptotic behavior in an oncolytic virotherapy model

Xueyan Tao,Shulin Zhou

DOI: https://doi.org/10.1016/j.jde.2021.11.003

IF: 2.615

2022-01-01

Journal of Differential Equations

Abstract:This work studies a haptotactic cross-diffusion system modeling oncolytic virotherapy in two-dimensional domains, accounting for interaction among uninfected and infected cancer cells, extracellular matrix and oncolytic viruses. It is rigorously proved that an associated no-flux initial-boundary value problem has a unique global classical solution which is uniformly bounded under suitable assumptions. Moreover, it is shown that a bounded solution can approach spatially constant equilibrium in the large time limit. (c) 2021 Elsevier Inc. All rights reserved.
A general non-local delay model on oncolytic virus therapy

Zizi Wang,Qian Zhang,Yong Luo

DOI: https://doi.org/10.1016/j.apm.2021.09.045

IF: 5.336

2022-02-01

Applied Mathematical Modelling

Abstract:The oncolytic virus is regarded as a novel, powerful, and biologically safe method of cancer treatment. A general delay differential system was driven by the age-dependent model better to understand the interaction between tumor cells and viruses. General continuous functions F(x,y) and G(x) depict the tumor proliferation rate and virus infection rate. The critical threshold value R0 was calculated that plays a determinant role in whether virus therapy occurs. The non-local delay term ∫t−τtβG(x(θ))v(θ)e−α(t−θ)dθ makes our model hard to analyze when using the traditional eigenvalue method. The method combining implicit function theorem and comparison theorem is used to overcome this problem. Furthermore, we support the fact that virotherapy can lead to tumor remission by using the fluctuation method. Lastly, Bayesian information criterion was adopted to select a better model when fitting the experimental data.

mechanics,engineering, multidisciplinary,mathematics, interdisciplinary applications
Modeling oncolytic virus therapy with distributed delay and non-local diffusion

Zizi Wang

2024-02-21

Abstract:In the field of modeling the dynamics of oncolytic viruses, researchers often face the challenge of using specialized mathematical terms to explain uncertain biological phenomena. This paper introduces a basic framework for an oncolytic virus dynamics model with a general growth rate $\mathcal{F}$ and a general nonlinear incidence term $\mathcal{G}$. The construction and derivation of the model explain in detail the generation process and practical significance of the distributed time delays and non-local infection terms. The paper provides the existence and uniqueness of solutions to the model, as well as the existence of a global attractor. Furthermore, through two auxiliary linear partial differential equations, the threshold parameters $\sigma_1$ are determined for sustained tumor growth and $\lambda_1$ for successful viral invasion of tumor cells to analyze the global dynamic behavior of the model. Finally, we illustrate and analyze our abstract theoretical results through a specific example.

Mathematics
Blow-up Prevention by Logistic Damping in a Chemotaxis-May-Nowak Model for Virus Infection

Zhang, Qingshan

DOI: https://doi.org/10.1007/s00025-024-02183-7

2024-05-09

Results in Mathematics

Abstract:In this paper, we study the no-flux boundary initial-boundary problem for a three-component reaction-diffusion system originating from the classical May-Nowak model for viral infection

mathematics, applied
Mathematical Models Of Tumor Treatment With Virus

Youshan Tao,Dashun Wang

2008-01-01

Abstract:We consider several mathematical models of tumor treatment involved in using oncolytic virus. These models are free boundary problems of nonlinear systems of partial differential equations, where the surface of a tumor is the free boundary. The first model deals with the competition between tumor cells, a replication-competent virus and an immune response. The second model describes combined therapies against cancer using virus and inhibitor. Finally, we consider a model of cancer radiovirotherapy which requires not only injection of replication-competent virus but also administration of radiodide.
Mathematical Analysis of a Model of a Replication-competent Virus in Tumor Cells

Tao Youshan

DOI: https://doi.org/10.3321/j.issn:1003-3998.2006.02.004

2006-01-01

Acta Mathematica Scientia

Abstract:The author considers a procedure for cancer therapy which consists of injecting replication-competent viruses into the tumor. The viruses infect tumor cells, replicate insider them, and eventually cause their death. As infected cells die, the viruses inside them are released and then proceed to infect adjacent tumor cells. This process is modelled as a free boundary problem for a nonlinear system of hyperbolic differential equations, where the free boundary is the surface of the tumor. The unknowns are the densities of uninfected cells, infected cells, necrotic cells and the free virus particles, and the velocity of cells within the tumor as well as the free boundary r = R(t). The aim of this paper is to explore the conditions under which the tumor can be made to shrink to zero.
Analysis of a model of a virus that replicates selectively in tumor cells

Avner Friedman,Youshan Tao

DOI: https://doi.org/10.1007/s00285-003-0199-5

2003-01-01

Journal of Mathematical Biology

Abstract:. We consider a procedure for cancer therapy which consists of injecting replication-competent viruses into the tumor. The viruses infect tumor cells, replicate inside them, and eventually cause their death. As infected cells die, the viruses inside them are released and then proceed to infect adjacent tumor cells. This process is modelled as a free boundary problem for a nonlinear system of hyperbolic-parabolic differential equations, where the free boundary is the surface of the tumor. The unknowns are the densities of uninfected cells, infected cells, necrotic cells and free virus particles, and the velocity of cells within the tumor as well as the free boundary r = R ( t ). The purpose of this paper is to establish a rigorous mathematical analysis of the model, and to explore the reduction of the tumor size that can be achieved by this therapy.
Oscillations in a Spatial Oncolytic Virus Model

Arwa Abdulla Baabdulla,Thomas Hillen

DOI: https://doi.org/10.1007/s11538-024-01322-z

2024-06-20

Bulletin of Mathematical Biology

Abstract: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.

mathematical & computational biology,biology
Modeling the Virus-Induced Tumor-Specific Immune Response with Delay in Tumor Virotherapy.

Qian Li,Yanni Xiao

DOI: https://doi.org/10.1016/j.cnsns.2021.106196

IF: 4.186

2021-01-01

Communications in Nonlinear Science and Numerical Simulation

Abstract:It is urgently required to design novel cancer therapies which lead to permanent cancer eradication or cancer control. Oncolytic virotherapy is a promising cancer treatment strategy using genetically engineered viruses which can selectively infect, replicate in and kill tumor cells without harming normal cells. Due to the variable interactions between tumor cells and oncolytic viruses (OVs) as well as the accordingly immune response, the impact of tumor virotherapy on tumor control is prominent. We propose a novel mathematical modeling framework based on delay differential equation to study tumor virotherapy with mediated antitumor immunity by OVs, which incorporates complex tumor–virus–immune system interactions. We initially study the existence and local stability of equilibria, and theoretically and numerically investigate the local Hopf bifurcation from the positive equilibrium by considering the time delay as a bifurcation parameter. Further, we analyze the direction of Hopf bifurcation and the stability of bifurcating periodic solution. Main results show that the time delay can induce Hopf bifurcation and result in periodic oscillations, indicating that the delayed tumor-specific CTL response induced by OVs leads to complex dynamics and may significantly influence the development process of tumor growth. Our findings also provide insight into important aspects of the virotherapy, including the dependence of the efficacy on key factors. We find that enhancing the induction of antitumor immunity by OVs can reduce the complexity of dynamics by strengthening the stability and show considerable effects on combating tumor. Further, proper choice of OVs with relative strong ability to lyse tumor cells is beneficial to ultimately control the development process of tumor growth.
Questing for an optimal, universal viral agent for oncolytic virotherapy

L. R. Paiva,M. L. Martins,S. C. Ferreira

DOI: https://doi.org/10.1103/physreve.84.041918

IF: 2.707

2011-10-17

Physical Review E

Abstract:One of the most promising strategies to treat cancer is attacking it with viruses designed to exploit specific altered pathways. Here, the effects of oncolytic virotherapy on tumors having compact, papillary, and disconnected morphologies are investigated through computer simulations of a multiscale model coupling macroscopic reaction-diffusion equations for the nutrients with microscopic stochastic rules for the actions of individual cells and viruses. The interaction among viruses and tumor cells involves cell infection, intracellular virus replication, and the release of new viruses in the tissue after cell lysis. The evolution over time of both the viral load and cancer cell population, as well as the probabilities for tumor eradication, were evaluated for a range of multiplicities of infection, viral entries, and burst sizes. It was found that in immunosuppressed hosts, the antitumor efficacy of a virus is primarily determined by its entry efficiency, its replicative capacity within the tumor, and its ability to spread over the tissue. However, the optimal traits for oncolytic viruses depend critically on the tumor growth dynamics and do not necessarily include rapid replication, cytolysis, or spreading, currently assumed as necessary conditions for a successful therapeutic outcome. Our findings have potential implications on the design of new vectors for the viral therapy of cancer.

physics, fluids & plasmas, mathematical
Backward Hopf Bifurcation In A Mathematical Model For Oncolytic Virotherapy With The Infection Delay And Innate Immune Effects

Yuxiao Guo,Ben Niu,Jianjun Paul Tian

DOI: https://doi.org/10.1080/17513758.2019.1667443

2019-01-01

Journal of Biological Dynamics

Abstract:In this paper, we consider a system of delay differential equations that models the oncolytic virotherapy on solid tumours with the delay of viral infection in the presence of the innate immune response. We conduct qualitative and numerical analysis, and provide possible medical implications for our results. The system has four equilibrium solutions. Fixed point analysis indicates that increasing the burst size and infection rate of the viruses has positive contribution to the therapy. However, increasing the immune killing infection rate, the immune stimulation rate, or the immune killing virus rate may lead the treatment failed. The viral infection time delay induces backward Hopf bifurcations, which means that the therapy may fail before time delay increases passing through a Hopf bifurcation. The parameter analysis also shows how saddle-node and Hopf bifurcations occur as viral burst size and other parameters vary, which yields further insights into the dynamics of the virotherapy.
Treating cancerous cells with a continuous release of virus particles

Ahmed Msmali,Mark Nelson,Fahad Al Saadi

DOI: https://doi.org/10.21914/anziamj.v63.17108

2022-12-07

Abstract:We investigate a model for the treatment of a tumour through the application of a virus. In the original model it was assumed that the virus particles are released only at one time. Such a treatment strategy cannot eliminate a tumour, as the tumour-free steady-state solution is unstable except for pathological circumstances in which the tumour does not grow and/or the virus does not die. We extend the model by allowing the tumour to be treated by a continuous release of virus particles. We show that the scaled delivery rate has two threshold values: below the lower threshold the system evolves to a stable periodic solution; above the higher threshold the tumour is eradicated. References C. E. Engeland, J. P. W. Heidbuechel, R. P. Araujo, and A. L. Jenner. Improving immunovirotherapies: The intersection of mathematical modelling and experiments. ImmunoInformatics 6 (2022), p. 100011. doi: 10.1016/j.immuno.2022.100011 A. L. Jenner, A. C. F. Coster, P. S. Kim, and F. Frascoli. Treating cancerous cells with viruses. Lett. Biomath. 5.2 (2018), S117–S136. doi: 10.1080/23737867.2018.1440977 A. L. Jenner, F. Frascoli, C.-O. Yun, and P. S. Kim. Optimising hydrogel release profiles for viro-immunotherapy using oncolytic adenovirus expressing IL-12 and GM-CSF with immature dendritic cells. Appl. Sci. 10.8 (2020). doi: 10.3390/app10082872 J. P. Tian. The replicability of oncolytic virus: Defining conditions in tumor virotherapy. Math. Bio. Eng. 8.3 (2011), pp. 841–860. doi: 10.3934/mbe.2011.8.841

mathematics, applied

A critical virus production rate for blow-up suppression in a haptotaxis model for oncolytic virotherapy

A Critical Virus Production Rate for Efficiency of Oncolytic Virotherapy

Asymptotic stability of spatial homogeneity in a haptotaxis model for oncolytic virotherapy

Global Boundedness and Asymptotic Behavior in a Double Haptotaxis Model for Oncolytic Virotherapy

Global Classical Solutions to a Doubly Haptotactic Cross-Diffusion System Modeling Oncolytic Virotherapy

A parabolic-hyperbolic free boundary problem modelling tumor treatment with virus

Simulation of a Model of Tumors with Virus-therapy

Critical Mass for Infinite-Time Blow-Up in a Haptotaxis System with Nonlinear Zero-Order Interaction

Dampening effects on global boundedness and asymptotic behavior in an oncolytic virotherapy model

A general non-local delay model on oncolytic virus therapy

Modeling oncolytic virus therapy with distributed delay and non-local diffusion

Blow-up Prevention by Logistic Damping in a Chemotaxis-May-Nowak Model for Virus Infection

Mathematical Models Of Tumor Treatment With Virus

Mathematical Analysis of a Model of a Replication-competent Virus in Tumor Cells

Analysis of a model of a virus that replicates selectively in tumor cells

Oscillations in a Spatial Oncolytic Virus Model

Modeling the Virus-Induced Tumor-Specific Immune Response with Delay in Tumor Virotherapy.

Questing for an optimal, universal viral agent for oncolytic virotherapy

Backward Hopf Bifurcation In A Mathematical Model For Oncolytic Virotherapy With The Infection Delay And Innate Immune Effects

Treating cancerous cells with a continuous release of virus particles