2
|
Phan TA, Nguyen HD, Tian JP. Deterministic and stochastic modeling for PDGF-driven gliomas reveals a classification of gliomas. J Math Biol 2021; 83:22. [PMID: 34345961 DOI: 10.1007/s00285-021-01647-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2020] [Revised: 04/10/2021] [Accepted: 07/18/2021] [Indexed: 11/25/2022]
Abstract
Motivated by our study of infiltrating dynamics of immune cells into tumors, we propose a stochastic model in terms of Ito stochastic differential equations to study how two parameters, the chemoattractant production rate and the chemotactic coefficient, influence immune cell migration and how these parameters distinguish two types of gliomas. We conduct a detailed analysis of the stochastic model and its deterministic counterpart. The deterministic model can differentiate two types of gliomas according to the range of the chemoattractant production rate as two equilibrium solutions, while the stochastic model also can differentiate two types of gliomas according to the ranges of the chemoattractant production rate and chemotactic coefficient with thresholds as one non-zero ergodic invariant measure and one weak persistent state when the noise intensities are small. When the noise intensities are large comparing with the chemotactic coefficient, there is only one type of glioma that corresponds to a non-zero ergodic invariant measure. Using our experimental data, numerical simulations are carried out to demonstrate properties of our models, and we give medical interpretations and implications for our analytical results and numerical simulations. This study also confirms some of our results about IDH gliomas.
Collapse
Affiliation(s)
- Tuan Anh Phan
- Department of Mathematical Sciences, New Mexico State University, Las Cruces, NM, 88001, USA.,Institute for Modeling Collaboration and Innovation, 875 Perimeter Drive, MS 1122, Moscow, ID, 83844, USA
| | - Hai Dang Nguyen
- Department of Mathematics, The University of Alabama, Tuscaloosa, AL, 35401, USA
| | - Jianjun Paul Tian
- Department of Mathematical Sciences, New Mexico State University, Las Cruces, NM, 88001, USA.
| |
Collapse
|
3
|
Senekal NS, Mahasa KJ, Eladdadi A, de Pillis L, Ouifki R. Natural Killer Cells Recruitment in Oncolytic Virotherapy: A Mathematical Model. Bull Math Biol 2021; 83:75. [PMID: 34008149 DOI: 10.1007/s11538-021-00903-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Accepted: 04/20/2021] [Indexed: 01/17/2023]
Abstract
In this paper, we investigate how natural killer (NK) cell recruitment to the tumor microenvironment (TME) affects oncolytic virotherapy. NK cells play a major role against viral infections. They are, however, known to induce early viral clearance of oncolytic viruses, which hinders the overall efficacy of oncolytic virotherapy. Here, we formulate and analyze a simple mathematical model of the dynamics of the tumor, OV and NK cells using currently available preclinical information. The aim of this study is to characterize conditions under which the synergistic balance between OV-induced NK responses and required viral cytopathicity may or may not result in a successful treatment. In this study, we found that NK cell recruitment to the TME must take place neither too early nor too late in the course of OV infection so that treatment will be successful. NK cell responses are most influential at either early (partly because of rapid response of NK cells to viral infections or antigens) or later (partly because of antitumoral ability of NK cells) stages of oncolytic virotherapy. The model also predicts that: (a) an NK cell response augments oncolytic virotherapy only if viral cytopathicity is weak; (b) the recruitment of NK cells modulates tumor growth; and (c) the depletion of activated NK cells within the TME enhances the probability of tumor escape in oncolytic virotherapy. Taken together, our model results demonstrate that OV infection is crucial, not just to cytoreduce tumor burden, but also to induce the stronger NK cell response necessary to achieve complete or at least partial tumor remission. Furthermore, our modeling framework supports combination therapies involving NK cells and OV which are currently used in oncolytic immunovirotherapy to treat several cancer types.
Collapse
Affiliation(s)
- Noma Susan Senekal
- Department of Mathematics and Computer Science, National University of Lesotho, Roma, Maseru, Lesotho.
| | - Khaphetsi Joseph Mahasa
- Department of Mathematics and Computer Science, National University of Lesotho, Roma, Maseru, Lesotho
| | | | | | - Rachid Ouifki
- Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria, South Africa
| |
Collapse
|
4
|
Al-Tuwairqi SM, Al-Johani NO, Simbawa EA. Modeling dynamics of cancer radiovirotherapy. J Theor Biol 2020; 506:110405. [PMID: 32738266 DOI: 10.1016/j.jtbi.2020.110405] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Revised: 12/21/2019] [Accepted: 07/09/2020] [Indexed: 10/23/2022]
Abstract
Advances in genetic engineering have paved the way for a new therapy for cancer, which is called virotherapy. This treatment uses genetically engineered viruses which selectively infect, replicate in, and destroy cancer cells without damaging normal cells. Furthermore, current research and clinical trials have indicated that these viruses can be delivered as single agents or in combination with other therapies. In this paper, we propose systems of ordinary differential equations for modeling the dynamics of aggressive tumor growth under radiovirotherapy treatment. We divide the treatment period into two phases; consequently, we present two mathematical models. First, we formulate the virotherapy model as Phase I of the treatment. Then we extend the model to include radiotherapy in combination with virotherapy as Phase II of the treatment. Comprehensive qualitative analyses of both models are conducted. Furthermore, numerical experiments are performed in order to support the analytical results. An analysis of the parameters is also carried out to investigate their effects on the outcome of the treatment. Overall, the analytical results reveal that radiovirotherapy is more effective than, and a good alternative to, virotherapy, as it is capable of eradicating tumors completely.
Collapse
Affiliation(s)
| | - Najwa O Al-Johani
- Mathematics department, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Eman A Simbawa
- Mathematics department, King Abdulaziz University, Jeddah, Saudi Arabia
| |
Collapse
|
5
|
Pooladvand P, Yun CO, Yoon AR, Kim PS, Frascoli F. The role of viral infectivity in oncolytic virotherapy outcomes: A mathematical study. Math Biosci 2020; 334:108520. [PMID: 33290764 DOI: 10.1016/j.mbs.2020.108520] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2020] [Revised: 10/15/2020] [Accepted: 12/01/2020] [Indexed: 10/22/2022]
Abstract
A model capturing the dynamics between virus and tumour cells in the context of oncolytic virotherapy is presented and analysed. The ability of the virus to be internalised by uninfected cells is described by an infectivity parameter, which is inferred from available experimental data. The parameter is also able to describe the effects of changes in the tumour environment that affect viral uptake from tumour cells. Results show that when a virus is inoculated inside a growing tumour, strategies for enhancing infectivity do not lead to a complete eradication of the tumour. Within typical times of experiments and treatments, we observe the onset of oscillations, which always prevent a full destruction of the tumour mass. These findings are in good agreement with available laboratory results. Further analysis shows why a fully successful therapy cannot exist for the proposed model and that care must be taken when designing and engineering viral vectors with enhanced features. In particular, bifurcation analysis reveals that creating longer lasting virus particles or using strategies for reducing infected cell lifespan can cause unexpected and unwanted surges in the overall tumour load over time. Our findings suggest that virotherapy alone seems unlikely to be effective in clinical settings unless adjuvant strategies are included.
Collapse
Affiliation(s)
- Pantea Pooladvand
- School of Mathematics and Statistics, The University of Sydney, Sydney, NSW 2006, Australia.
| | - Chae-Ok Yun
- Department of Bioengineering, Collage of Engineering, Hanyang University, Seoul, South Korea; Institute of Nano Science and Technology (INST), Hanyang University, Seoul, South Korea
| | - A-Rum Yoon
- Department of Bioengineering, Collage of Engineering, Hanyang University, Seoul, South Korea; Institute of Nano Science and Technology (INST), Hanyang University, Seoul, South Korea
| | - Peter S Kim
- School of Mathematics and Statistics, The University of Sydney, Sydney, NSW 2006, Australia
| | - Federico Frascoli
- Department of Mathematics, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, VIC 3122, Australia
| |
Collapse
|
6
|
Chelliah V, Lazarou G, Bhatnagar S, Gibbs JP, Nijsen M, Ray A, Stoll B, Thompson RA, Gulati A, Soukharev S, Yamada A, Weddell J, Sayama H, Oishi M, Wittemer-Rump S, Patel C, Niederalt C, Burghaus R, Scheerans C, Lippert J, Kabilan S, Kareva I, Belousova N, Rolfe A, Zutshi A, Chenel M, Venezia F, Fouliard S, Oberwittler H, Scholer-Dahirel A, Lelievre H, Bottino D, Collins SC, Nguyen HQ, Wang H, Yoneyama T, Zhu AZX, van der Graaf PH, Kierzek AM. Quantitative Systems Pharmacology Approaches for Immuno-Oncology: Adding Virtual Patients to the Development Paradigm. Clin Pharmacol Ther 2020; 109:605-618. [PMID: 32686076 PMCID: PMC7983940 DOI: 10.1002/cpt.1987] [Citation(s) in RCA: 44] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2020] [Accepted: 07/06/2020] [Indexed: 12/12/2022]
Abstract
Drug development in oncology commonly exploits the tools of molecular biology to gain therapeutic benefit through reprograming of cellular responses. In immuno‐oncology (IO) the aim is to direct the patient’s own immune system to fight cancer. After remarkable successes of antibodies targeting PD1/PD‐L1 and CTLA4 receptors in targeted patient populations, the focus of further development has shifted toward combination therapies. However, the current drug‐development approach of exploiting a vast number of possible combination targets and dosing regimens has proven to be challenging and is arguably inefficient. In particular, the unprecedented number of clinical trials testing different combinations may no longer be sustainable by the population of available patients. Further development in IO requires a step change in selection and validation of candidate therapies to decrease development attrition rate and limit the number of clinical trials. Quantitative systems pharmacology (QSP) proposes to tackle this challenge through mechanistic modeling and simulation. Compounds’ pharmacokinetics, target binding, and mechanisms of action as well as existing knowledge on the underlying tumor and immune system biology are described by quantitative, dynamic models aiming to predict clinical results for novel combinations. Here, we review the current QSP approaches, the legacy of mathematical models available to quantitative clinical pharmacologists describing interaction between tumor and immune system, and the recent development of IO QSP platform models. We argue that QSP and virtual patients can be integrated as a new tool in existing IO drug development approaches to increase the efficiency and effectiveness of the search for novel combination therapies.
Collapse
Affiliation(s)
| | | | | | | | | | - Avijit Ray
- Abbvie Inc., North Chicago, Illinois, USA
| | | | | | - Abhishek Gulati
- Astellas Pharma Global Development Inc./Astellas Pharma Inc., Northbrook, Illinois, USA.,Astellas Pharma Global Development Inc./Astellas Pharma Inc., Tokyo or Tsukuba-shi, Japan
| | - Serguei Soukharev
- Astellas Pharma Global Development Inc./Astellas Pharma Inc., Northbrook, Illinois, USA.,Astellas Pharma Global Development Inc./Astellas Pharma Inc., Tokyo or Tsukuba-shi, Japan
| | - Akihiro Yamada
- Astellas Pharma Global Development Inc./Astellas Pharma Inc., Northbrook, Illinois, USA.,Astellas Pharma Global Development Inc./Astellas Pharma Inc., Tokyo or Tsukuba-shi, Japan
| | - Jared Weddell
- Astellas Pharma Global Development Inc./Astellas Pharma Inc., Northbrook, Illinois, USA.,Astellas Pharma Global Development Inc./Astellas Pharma Inc., Tokyo or Tsukuba-shi, Japan
| | - Hiroyuki Sayama
- Astellas Pharma Global Development Inc./Astellas Pharma Inc., Northbrook, Illinois, USA.,Astellas Pharma Global Development Inc./Astellas Pharma Inc., Tokyo or Tsukuba-shi, Japan
| | - Masayo Oishi
- Astellas Pharma Global Development Inc./Astellas Pharma Inc., Northbrook, Illinois, USA.,Astellas Pharma Global Development Inc./Astellas Pharma Inc., Tokyo or Tsukuba-shi, Japan
| | | | | | | | | | | | | | | | - Irina Kareva
- EMD Serono, Merck KGaA, Billerica, Massachusetts, USA
| | | | - Alex Rolfe
- EMD Serono, Merck KGaA, Billerica, Massachusetts, USA
| | - Anup Zutshi
- EMD Serono, Merck KGaA, Billerica, Massachusetts, USA
| | | | | | | | | | | | | | - Dean Bottino
- Millennium Pharmaceuticals Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Ltd., Cambridge, Massachusetts, USA
| | - Sabrina C Collins
- Millennium Pharmaceuticals Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Ltd., Cambridge, Massachusetts, USA
| | - Hoa Q Nguyen
- Millennium Pharmaceuticals Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Ltd., Cambridge, Massachusetts, USA
| | - Haiqing Wang
- Millennium Pharmaceuticals Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Ltd., Cambridge, Massachusetts, USA
| | - Tomoki Yoneyama
- Millennium Pharmaceuticals Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Ltd., Cambridge, Massachusetts, USA
| | - Andy Z X Zhu
- Millennium Pharmaceuticals Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Ltd., Cambridge, Massachusetts, USA
| | | | | |
Collapse
|
9
|
Investigating Macrophages Plasticity Following Tumour-Immune Interactions During Oncolytic Therapies. Acta Biotheor 2019; 67:321-359. [PMID: 31410657 PMCID: PMC6825040 DOI: 10.1007/s10441-019-09357-9] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2018] [Accepted: 08/02/2019] [Indexed: 12/22/2022]
Abstract
Over the last few years, oncolytic virus therapy has been recognised as a promising approach in cancer treatment, due to the potential of these viruses to induce systemic anti-tumour immunity and selectively killing tumour cells. However, the effectiveness of these viruses depends significantly on their interactions with the host immune responses, both innate (e.g., macrophages, which accumulate in high numbers inside solid tumours) and adaptive (e.g., \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\hbox {CD8}^{+}$$\end{document}CD8+ T cells). In this article, we consider a mathematical approach to investigate the possible outcomes of the complex interactions between two extreme types of macrophages (M1 and M2 cells), effector \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\hbox {CD8}^{+}$$\end{document}CD8+ T cells and an oncolytic Vesicular Stomatitis Virus (VSV), on the growth/elimination of B16F10 melanoma. We discuss, in terms of VSV, \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\hbox {CD8}^{+}$$\end{document}CD8+ and macrophages levels, two different types of immune responses which could ensure tumour control and eventual elimination. We show that both innate and adaptive anti-tumour immune responses, as well as the oncolytic virus, could be very important in delaying tumour relapse and eventually eliminating the tumour. Overall this study supports the use mathematical modelling to increase our understanding of the complex immune interaction following oncolytic virotherapies. However, the complexity of the model combined with a lack of sufficient data for model parametrisation has an impact on the possibility of making quantitative predictions.
Collapse
|