1
|
Kolokotroni E, Abler D, Ghosh A, Tzamali E, Grogan J, Georgiadi E, Büchler P, Radhakrishnan R, Byrne H, Sakkalis V, Nikiforaki K, Karatzanis I, McFarlane NJB, Kaba D, Dong F, Bohle RM, Meese E, Graf N, Stamatakos G. A Multidisciplinary Hyper-Modeling Scheme in Personalized In Silico Oncology: Coupling Cell Kinetics with Metabolism, Signaling Networks, and Biomechanics as Plug-In Component Models of a Cancer Digital Twin. J Pers Med 2024; 14:475. [PMID: 38793058 PMCID: PMC11122096 DOI: 10.3390/jpm14050475] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2024] [Revised: 04/11/2024] [Accepted: 04/17/2024] [Indexed: 05/26/2024] Open
Abstract
The massive amount of human biological, imaging, and clinical data produced by multiple and diverse sources necessitates integrative modeling approaches able to summarize all this information into answers to specific clinical questions. In this paper, we present a hypermodeling scheme able to combine models of diverse cancer aspects regardless of their underlying method or scale. Describing tissue-scale cancer cell proliferation, biomechanical tumor growth, nutrient transport, genomic-scale aberrant cancer cell metabolism, and cell-signaling pathways that regulate the cellular response to therapy, the hypermodel integrates mutation, miRNA expression, imaging, and clinical data. The constituting hypomodels, as well as their orchestration and links, are described. Two specific cancer types, Wilms tumor (nephroblastoma) and non-small cell lung cancer, are addressed as proof-of-concept study cases. Personalized simulations of the actual anatomy of a patient have been conducted. The hypermodel has also been applied to predict tumor control after radiotherapy and the relationship between tumor proliferative activity and response to neoadjuvant chemotherapy. Our innovative hypermodel holds promise as a digital twin-based clinical decision support system and as the core of future in silico trial platforms, although additional retrospective adaptation and validation are necessary.
Collapse
Affiliation(s)
- Eleni Kolokotroni
- In Silico Oncology and In Silico Medicine Group, Institute of Communication and Computer Systems, School of Electrical and Computer Engineering, National Technical University of Athens, 157 80 Zografos, Greece;
| | - Daniel Abler
- Department of Oncology, Geneva University Hospitals and University of Geneva, 1205 Geneva, Switzerland;
- Department of Oncology, Lausanne University Hospital and University of Lausanne, 1011 Lausanne, Switzerland
| | - Alokendra Ghosh
- Department of Chemical and Biomolecular Engineering, Department of Bioengineering, University of Pennsylvania, Philadelphia, PA 19104, USA; (A.G.); (R.R.)
| | - Eleftheria Tzamali
- Institute of Computer Science, Foundation for Research and Technology—Hellas, 70013 Heraklion, Greece; (E.T.); (V.S.); (K.N.); (I.K.)
| | - James Grogan
- Irish Centre for High End Computing, University of Galway, H91 TK33 Galway, Ireland;
| | - Eleni Georgiadi
- In Silico Oncology and In Silico Medicine Group, Institute of Communication and Computer Systems, School of Electrical and Computer Engineering, National Technical University of Athens, 157 80 Zografos, Greece;
- Biomedical Engineering Department, University of West Attica, 12243 Egaleo, Greece
| | | | - Ravi Radhakrishnan
- Department of Chemical and Biomolecular Engineering, Department of Bioengineering, University of Pennsylvania, Philadelphia, PA 19104, USA; (A.G.); (R.R.)
| | - Helen Byrne
- Mathematical Institute, University of Oxford, Oxford OX1 2JD, UK;
| | - Vangelis Sakkalis
- Institute of Computer Science, Foundation for Research and Technology—Hellas, 70013 Heraklion, Greece; (E.T.); (V.S.); (K.N.); (I.K.)
| | - Katerina Nikiforaki
- Institute of Computer Science, Foundation for Research and Technology—Hellas, 70013 Heraklion, Greece; (E.T.); (V.S.); (K.N.); (I.K.)
| | - Ioannis Karatzanis
- Institute of Computer Science, Foundation for Research and Technology—Hellas, 70013 Heraklion, Greece; (E.T.); (V.S.); (K.N.); (I.K.)
| | | | - Djibril Kaba
- Department of Computer Science and Technology, University of Bedfordshire, Luton LU1 3JU, UK;
| | - Feng Dong
- Department of Computer & Information Sciences, University of Strathclyde, Glasgow G1 1XH, UK;
| | - Rainer M. Bohle
- Department of Pathology, Saarland University, 66421 Homburg, Germany;
| | - Eckart Meese
- Department of Human Genetics, Saarland University, 66421 Homburg, Germany;
| | - Norbert Graf
- Department of Paediatric Oncology and Haematology, Saarland University, 66421 Homburg, Germany;
| | - Georgios Stamatakos
- In Silico Oncology and In Silico Medicine Group, Institute of Communication and Computer Systems, School of Electrical and Computer Engineering, National Technical University of Athens, 157 80 Zografos, Greece;
| |
Collapse
|
2
|
Hamis S, Somervuo P, Ågren JA, Tadele DS, Kesseli J, Scott JG, Nykter M, Gerlee P, Finkelshtein D, Ovaskainen O. Spatial cumulant models enable spatially informed treatment strategies and analysis of local interactions in cancer systems. J Math Biol 2023; 86:68. [PMID: 37017776 PMCID: PMC10076412 DOI: 10.1007/s00285-023-01903-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Revised: 01/13/2023] [Accepted: 03/09/2023] [Indexed: 04/06/2023]
Abstract
Theoretical and applied cancer studies that use individual-based models (IBMs) have been limited by the lack of a mathematical formulation that enables rigorous analysis of these models. However, spatial cumulant models (SCMs), which have arisen from theoretical ecology, describe population dynamics generated by a specific family of IBMs, namely spatio-temporal point processes (STPPs). SCMs are spatially resolved population models formulated by a system of differential equations that approximate the dynamics of two STPP-generated summary statistics: first-order spatial cumulants (densities), and second-order spatial cumulants (spatial covariances). We exemplify how SCMs can be used in mathematical oncology by modelling theoretical cancer cell populations comprising interacting growth factor-producing and non-producing cells. To formulate model equations, we use computational tools that enable the generation of STPPs, SCMs and mean-field population models (MFPMs) from user-defined model descriptions (Cornell et al. Nat Commun 10:4716, 2019). To calculate and compare STPP, SCM and MFPM-generated summary statistics, we develop an application-agnostic computational pipeline. Our results demonstrate that SCMs can capture STPP-generated population density dynamics, even when MFPMs fail to do so. From both MFPM and SCM equations, we derive treatment-induced death rates required to achieve non-growing cell populations. When testing these treatment strategies in STPP-generated cell populations, our results demonstrate that SCM-informed strategies outperform MFPM-informed strategies in terms of inhibiting population growths. We thus demonstrate that SCMs provide a new framework in which to study cell-cell interactions, and can be used to describe and perturb STPP-generated cell population dynamics. We, therefore, argue that SCMs can be used to increase IBMs' applicability in cancer research.
Collapse
Affiliation(s)
- Sara Hamis
- Tampere Institute for Advanced Study, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland.
- Department of Biological and Environmental Science, University of Jyväskylä, Jyväskylä, Finland.
| | - Panu Somervuo
- Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland
| | - J Arvid Ågren
- Department of Evolutionary Biology, Uppsala University, Uppsala, Sweden
- Department of Translational Hematology and Oncology Research, Cleveland Clinic, Cleveland, OH, USA
| | - Dagim Shiferaw Tadele
- Department of Translational Hematology and Oncology Research, Cleveland Clinic, Cleveland, OH, USA
- Department for Medical Genetics, Oslo University Hospital, Ullevål, Oslo, Norway
| | - Juha Kesseli
- Prostate Cancer Research Center, Faculty of Medicine and Health Technology, Tampere University and Tays Cancer Centre, Tampere, Finland
| | - Jacob G Scott
- Department of Translational Hematology and Oncology Research, Cleveland Clinic, Cleveland, OH, USA
- Case Western Reserve School of Medicine, Cleveland, OH, USA
- Department of Radiation Oncology, Cleveland Clinic, Cleveland, OH, USA
| | - Matti Nykter
- Prostate Cancer Research Center, Faculty of Medicine and Health Technology, Tampere University and Tays Cancer Centre, Tampere, Finland
- Foundation for the Finnish Cancer Institute, Helsinki, Finland
| | - Philip Gerlee
- Mathematical Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Mathematical Sciences, University of Gothenburg, Gothenburg, Sweden
| | - Dmitri Finkelshtein
- Department of Mathematics, Faculty of Science and Engineering, Swansea University, Swansea, UK
| | - Otso Ovaskainen
- Department of Biological and Environmental Science, University of Jyväskylä, Jyväskylä, Finland
- Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland
- Department of Biology, Centre for Biodiversity Dynamics, Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
3
|
Mohammadrezaei D, Moghimi N, Vandvajdi S, Powathil G, Hamis S, Kohandel M. Predicting and elucidating the post-printing behavior of 3D printed cancer cells in hydrogel structures by integrating in-vitro and in-silico experiments. Sci Rep 2023; 13:1211. [PMID: 36681762 PMCID: PMC9867702 DOI: 10.1038/s41598-023-28286-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2022] [Accepted: 01/16/2023] [Indexed: 01/22/2023] Open
Abstract
A key feature distinguishing 3D bioprinting from other 3D cell culture techniques is its precise control over created structures. This property allows for the high-resolution fabrication of biomimetic structures with controlled structural and mechanical properties such as porosity, permeability, and stiffness. However, analyzing post-printing cellular dynamics and optimizing their functions within the 3D fabricated environment is only possible through trial and error and replicating several experiments. This issue motivated the development of a cellular automata model for the first time to simulate post-printing cell behaviour within the 3D bioprinted construct. To improve our model, we bioprinted a 3D construct using MDA-MB-231 cell-laden hydrogel and evaluated cellular functions, including viability and proliferation in 11 days. The results showed that our model successfully simulated the 3D bioprinted structure and captured in-vitro observations. We demonstrated that in-silico model could predict and elucidate post-printing biological functions for different initial cell numbers in bioink and different bioink formulations with gelatine and alginate, without replicating several costly and time-consuming in-vitro measurements. We believe such a computational framework will substantially impact 3D bioprinting's future application. We hope this study inspires researchers to further realize how an in-silico model might be utilized to advance in-vitro 3D bioprinting research.
Collapse
Affiliation(s)
- Dorsa Mohammadrezaei
- Department of Applied Mathematics, University of Waterloo, 200 University Ave West, Waterloo, ON, N2L 3G1, Canada.
| | - Nafiseh Moghimi
- Department of Applied Mathematics, University of Waterloo, 200 University Ave West, Waterloo, ON, N2L 3G1, Canada
| | - Shadi Vandvajdi
- Department of Applied Mathematics, University of Waterloo, 200 University Ave West, Waterloo, ON, N2L 3G1, Canada
| | - Gibin Powathil
- Department of Mathematics, Faculty of Science and Engineering, Swansea University, Swansea, UK
| | - Sara Hamis
- School of Mathematics and Statistics, University of St Andrews, St Andrews, UK
| | - Mohammad Kohandel
- Department of Applied Mathematics, University of Waterloo, 200 University Ave West, Waterloo, ON, N2L 3G1, Canada
| |
Collapse
|
4
|
Butner JD, Dogra P, Chung C, Pasqualini R, Arap W, Lowengrub J, Cristini V, Wang Z. Mathematical modeling of cancer immunotherapy for personalized clinical translation. NATURE COMPUTATIONAL SCIENCE 2022; 2:785-796. [PMID: 38126024 PMCID: PMC10732566 DOI: 10.1038/s43588-022-00377-z] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 11/14/2022] [Indexed: 12/23/2023]
Abstract
Encouraging advances are being made in cancer immunotherapy modeling, especially in the key areas of developing personalized treatment strategies based on individual patient parameters, predicting treatment outcomes and optimizing immunotherapy synergy when used in combination with other treatment approaches. Here we present a focused review of the most recent mathematical modeling work on cancer immunotherapy with a focus on clinical translatability. It can be seen that this field is transitioning from pure basic science to applications that can make impactful differences in patients' lives. We discuss how researchers are integrating experimental and clinical data to fully inform models so that they can be applied for clinical predictions, and present the challenges that remain to be overcome if widespread clinical adaptation is to be realized. Lastly, we discuss the most promising future applications and areas that are expected to be the focus of extensive upcoming modeling studies.
Collapse
Affiliation(s)
- Joseph D. Butner
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
| | - Prashant Dogra
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
| | - Caroline Chung
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Renata Pasqualini
- Rutgers Cancer Institute of New Jersey, Newark, NJ, USA
- Department of Radiation Oncology, Division of Cancer Biology, Rutgers New Jersey Medical School, Newark, NJ, USA
| | - Wadih Arap
- Rutgers Cancer Institute of New Jersey, Newark, NJ, USA
- Department of Medicine, Division of Hematology/Oncology, Rutgers New Jersey Medical School, Newark, NJ, USA
| | - John Lowengrub
- Department of Mathematics, University of California at Irvine, Irvine, CA, USA
| | - Vittorio Cristini
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
- Neal Cancer Center, Houston Methodist Research Institute, Houston, TX, USA
- Department of Imaging Physics, University of Texas MD Anderson Cancer Center, Houston, TX, USA
- Physiology, Biophysics, and Systems Biology Program, Graduate School of Medical Sciences, Weill Cornell Medicine, New York, NY, USA
| | - Zhihui Wang
- Mathematics in Medicine Program, Houston Methodist Research Institute, Houston, TX, USA
- Neal Cancer Center, Houston Methodist Research Institute, Houston, TX, USA
- Department of Physiology and Biophysics, Weill Cornell Medicine, New York, NY, USA
| |
Collapse
|
5
|
Deb D, Zhu S, LeBlanc MJ, Danino T. Assessing chemotherapy dosing strategies in a spatial cell culture model. Front Oncol 2022; 12:980770. [PMID: 36505801 PMCID: PMC9729937 DOI: 10.3389/fonc.2022.980770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2022] [Accepted: 11/07/2022] [Indexed: 11/25/2022] Open
Abstract
Predicting patient responses to chemotherapy regimens is a major challenge in cancer treatment. Experimental model systems coupled with quantitative mathematical models to calculate optimal dose and frequency of drugs can enable improved chemotherapy regimens. Here we developed a simple approach to track two-dimensional cell colonies composed of chemo-sensitive and resistant cell populations via fluorescence microscopy and coupled this to computational model predictions. Specifically, we first developed multiple 4T1 breast cancer cell lines resistant to varying concentrations of doxorubicin, and demonstrated how heterogeneous populations expand in a two-dimensional colony. We subjected cell populations to varied dose and frequency of chemotherapy and measured colony growth. We then built a mathematical model to describe the dynamics of both chemosensitive and chemoresistant populations, where we determined which number of doses can produce the smallest tumor size based on parameters in the system. Finally, using an in vitro model we demonstrated multiple doses can decrease overall colony growth as compared to a single dose at the same total dose. In the future, this system can be adapted to optimize dosing strategies in the setting of heterogeneous cell types or patient derived cells with varied chemoresistance.
Collapse
Affiliation(s)
- Dhruba Deb
- Department of Biomedical Engineering, Columbia University, New York, NY, United States
| | - Shu Zhu
- Department of Biomedical Engineering, Columbia University, New York, NY, United States
| | - Michael J. LeBlanc
- Department of Biomedical Engineering, Columbia University, New York, NY, United States
| | - Tal Danino
- Department of Biomedical Engineering, Columbia University, New York, NY, United States,Data Science Institute, Columbia University, New York, NY, United States,Herbert Irving Comprehensive Cancer Center, Columbia University, New York, NY, United States,*Correspondence: Tal Danino,
| |
Collapse
|
6
|
Wu C, Lorenzo G, Hormuth DA, Lima EABF, Slavkova KP, DiCarlo JC, Virostko J, Phillips CM, Patt D, Chung C, Yankeelov TE. Integrating mechanism-based modeling with biomedical imaging to build practical digital twins for clinical oncology. BIOPHYSICS REVIEWS 2022; 3:021304. [PMID: 35602761 PMCID: PMC9119003 DOI: 10.1063/5.0086789] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Accepted: 04/29/2022] [Indexed: 12/11/2022]
Abstract
Digital twins employ mathematical and computational models to virtually represent a physical object (e.g., planes and human organs), predict the behavior of the object, and enable decision-making to optimize the future behavior of the object. While digital twins have been widely used in engineering for decades, their applications to oncology are only just emerging. Due to advances in experimental techniques quantitatively characterizing cancer, as well as advances in the mathematical and computational sciences, the notion of building and applying digital twins to understand tumor dynamics and personalize the care of cancer patients has been increasingly appreciated. In this review, we present the opportunities and challenges of applying digital twins in clinical oncology, with a particular focus on integrating medical imaging with mechanism-based, tissue-scale mathematical modeling. Specifically, we first introduce the general digital twin framework and then illustrate existing applications of image-guided digital twins in healthcare. Next, we detail both the imaging and modeling techniques that provide practical opportunities to build patient-specific digital twins for oncology. We then describe the current challenges and limitations in developing image-guided, mechanism-based digital twins for oncology along with potential solutions. We conclude by outlining five fundamental questions that can serve as a roadmap when designing and building a practical digital twin for oncology and attempt to provide answers for a specific application to brain cancer. We hope that this contribution provides motivation for the imaging science, oncology, and computational communities to develop practical digital twin technologies to improve the care of patients battling cancer.
Collapse
Affiliation(s)
- Chengyue Wu
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, USA
| | | | | | | | - Kalina P. Slavkova
- Department of Physics, The University of Texas at Austin, Austin, Texas 78712, USA
| | | | | | - Caleb M. Phillips
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Debra Patt
- Texas Oncology, Austin, Texas 78731, USA
| | - Caroline Chung
- Department of Radiation Oncology, MD Anderson Cancer Center, University of Texas, Houston, Texas 77030, USA
| | | |
Collapse
|
7
|
Demetriades M, Zivanovic M, Hadjicharalambous M, Ioannou E, Ljujic B, Vucicevic K, Ivosevic Z, Dagovic A, Milivojevic N, Kokkinos O, Bauer R, Vavourakis V. Interrogating and Quantifying In Vitro Cancer Drug Pharmacodynamics via Agent-Based and Bayesian Monte Carlo Modelling. Pharmaceutics 2022; 14:pharmaceutics14040749. [PMID: 35456583 PMCID: PMC9029523 DOI: 10.3390/pharmaceutics14040749] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2022] [Revised: 03/25/2022] [Accepted: 03/28/2022] [Indexed: 02/01/2023] Open
Abstract
The effectiveness of chemotherapy in cancer cell regression is often limited by drug resistance, toxicity, and neoplasia heterogeneity. However, due to the significant complexities entailed by the many cancer growth processes, predicting the impact of interference and symmetry-breaking mechanisms is a difficult problem. To quantify and understand more about cancer drug pharmacodynamics, we combine in vitro with in silico cancer models. The anti-proliferative action of selected cytostatics is interrogated on human colorectal and breast adenocarcinoma cells, while an agent-based computational model is employed to reproduce experiments and shed light on the main therapeutic mechanisms of each chemotherapeutic agent. Multiple drug administration scenarios on each cancer cell line are simulated by varying the drug concentration, while a Bayesian-based method for model parameter optimisation is employed. Our proposed procedure of combining in vitro cancer drug screening with an in silico agent-based model successfully reproduces the impact of chemotherapeutic drugs in cancer growth behaviour, while the mechanisms of action of each drug are characterised through model-derived probabilities of cell apoptosis and division. We suggest that our approach could form the basis for the prospective generation of experimentally-derived and model-optimised pharmacological variables towards personalised cancer therapy.
Collapse
Affiliation(s)
- Marios Demetriades
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
| | - Marko Zivanovic
- Department of Science, Institute for Information Technologies Kragujevac, University of Kragujevac, 34000 Kragujevac, Serbia; (M.Z.); (N.M.)
| | - Myrianthi Hadjicharalambous
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
| | - Eleftherios Ioannou
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
| | - Biljana Ljujic
- Faculty of Medical Sciences, Human Genetics, University of Kragujevac, 34000 Kragujevac, Serbia; (B.L.); (Z.I.)
| | - Ksenija Vucicevic
- Department for Pharmaceutical Technologies, Faculty of Medical Sciences, University of Kragujevac, 34000 Kragujevac, Serbia;
| | - Zeljko Ivosevic
- Faculty of Medical Sciences, Human Genetics, University of Kragujevac, 34000 Kragujevac, Serbia; (B.L.); (Z.I.)
| | - Aleksandar Dagovic
- Oncology and Radiotherapy Centre, Faculty of Medical Sciences, University of Kragujevac, 34000 Kragujevac, Serbia;
| | - Nevena Milivojevic
- Department of Science, Institute for Information Technologies Kragujevac, University of Kragujevac, 34000 Kragujevac, Serbia; (M.Z.); (N.M.)
| | - Odysseas Kokkinos
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
| | - Roman Bauer
- Department of Computer Science, University of Surrey, Guilford GU2 7XH, UK
- Correspondence: (R.B.); or (V.V.)
| | - Vasileios Vavourakis
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 2109, Cyprus; (M.D.); (M.H.); (E.I.); (O.K.)
- Department of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, UK
- Correspondence: (R.B.); or (V.V.)
| |
Collapse
|
8
|
de los Reyes AA, Kim Y. Optimal regulation of tumour-associated neutrophils in cancer progression. ROYAL SOCIETY OPEN SCIENCE 2022; 9:210705. [PMID: 35127110 PMCID: PMC8808100 DOI: 10.1098/rsos.210705] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Accepted: 11/19/2021] [Indexed: 06/14/2023]
Abstract
In a tumour microenvironment, tumour-associated neutrophils could display two opposing differential phenotypes: anti-tumour (N1) and pro-tumour (N2) effector cells. Converting N2 to N1 neutrophils provides innovative therapies for cancer treatment. In this study, a mathematical model for N1-N2 dynamics describing the cancer survival and immune inhibition in response to TGF-β and IFN-β is considered. The effects of exogenous intervention of TGF-β inhibitor and IFN-β are examined in order to enhance N1 recruitment to combat tumour progression. Our approach employs optimal control theory to determine drug infusion protocols that could minimize tumour volume with least administration cost possible. Four optimal control scenarios corresponding to different therapeutic strategies are explored, namely, TGF-β inhibitor control only, IFN-β control only, concomitant TGF-β inhibitor and IFN-β controls, and alternating TGF-β inhibitor and IFN-β controls. For each scheme, different initial conditions are varied to depict different pathophysiological condition of a cancer patient, leading to adaptive treatment schedule. TGF-β inhibitor and IFN-β drug dosages, total drug amount, infusion times and relative cost of drug administrations are obtained under various circumstances. The control strategies achieved could guide in designing individualized therapeutic protocols.
Collapse
Affiliation(s)
- Aurelio A. de los Reyes
- Biomedical Mathematics Group, Pioneer Research Center for Mathematical and Computational Sciences, Institute for Basic Science, Daejeon 34126, Republic of Korea
- Institute of Mathematics, University of the Philippines Diliman, Quezon City 1101, Philippines
| | - Yangjin Kim
- Department of Mathematics, Konkuk University, Seoul 05029, Republic of Korea
- Mathematical Biosciences Institute, Columbus, OH 43210, USA
| |
Collapse
|
9
|
Italia M, Dercole F, Lucchetti R. Optimal chemotherapy counteracts cancer adaptive resistance in a cell-based, spatially-extended, evolutionary model. Phys Biol 2022; 19. [PMID: 35100568 DOI: 10.1088/1478-3975/ac509c] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2021] [Accepted: 01/31/2022] [Indexed: 11/12/2022]
Abstract
Most aggressive cancers are incurable due to their fast evolution of drug resistance. We model cancer growth and adaptive response in a simplified cell-based (CB) setting, assuming a genetic resistance to two chemotherapeutic drugs. We show that optimal administration protocols can steer cells resistance and turned it into a weakness for the disease. Our work extends the population-based (PB) model proposed by Orlando et al. (Physical Biology, 2012), in which a homogeneous population of cancer cells evolves according to a fitness landscape. The landscape models three types of trade-offs, differing on whether the cells are more, less, or equal effective when generalizing resistance to two drugs as opposed to specializing to a single one. The CB framework allows us to include genetic heterogeneity, spatial competition, and drugs diffusion, as well as realistic administration protocols. By calibrating our model on Orlando et al.'s assumptions, we show that dynamical protocols that alternate the two drugs minimize the cancer size at the end of (or at mid-points during) treatment. These results significantly differ from those obtained with the homogeneous model---suggesting static protocols under the pro-generalizing and neutral allocation trade-offs---highlighting the important role of spatial and genetic heterogeneities. Our work is the first attempt to search for optimal treatments in a CB setting, a step forward toward realistic clinical applications.
Collapse
Affiliation(s)
- Matteo Italia
- Electronics, Information and Bioengineering, Politecnico di Milano, Via Ponzio 34/5, Milano, 20133, ITALY
| | - Fabio Dercole
- Electronics, Information and Bioengineering, Politecnico di Milano, Via Ponzio 34/5, Milano, 20133, ITALY
| | - Roberto Lucchetti
- Mathematics, Politecnico di Milano, Via Edoardo Bonardi, 9, Milano, 20133, ITALY
| |
Collapse
|
10
|
Ioannidis JP. Pre-registration of mathematical models. Math Biosci 2022; 345:108782. [DOI: 10.1016/j.mbs.2022.108782] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2021] [Revised: 01/13/2022] [Accepted: 01/13/2022] [Indexed: 11/28/2022]
|
11
|
Quantifying ERK activity in response to inhibition of the BRAFV600E-MEK-ERK cascade using mathematical modelling. Br J Cancer 2021; 125:1552-1560. [PMID: 34621046 PMCID: PMC8608797 DOI: 10.1038/s41416-021-01565-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2021] [Revised: 09/06/2021] [Accepted: 09/21/2021] [Indexed: 01/07/2023] Open
Abstract
BACKGROUND Simultaneous inhibition of multiple components of the BRAF-MEK-ERK cascade (vertical inhibition) has become a standard of care for treating BRAF-mutant melanoma. However, the molecular mechanism of how vertical inhibition synergistically suppresses intracellular ERK activity, and consequently cell proliferation, are yet to be fully elucidated. METHODS We develop a mechanistic mathematical model that describes how the mutant BRAF inhibitor, dabrafenib, and the MEK inhibitor, trametinib, affect BRAFV600E-MEK-ERK signalling. The model is based on a system of chemical reactions that describes cascade signalling dynamics. Using mass action kinetics, the chemical reactions are re-expressed as ordinary differential equations that are parameterised by in vitro data and solved numerically to obtain the temporal evolution of cascade component concentrations. RESULTS The model provides a quantitative method to compute how dabrafenib and trametinib can be used in combination to synergistically inhibit ERK activity in BRAFV600E-mutant melanoma cells. The model elucidates molecular mechanisms of vertical inhibition of the BRAFV600E-MEK-ERK cascade and delineates how elevated BRAF concentrations generate drug resistance to dabrafenib and trametinib. The computational simulations further suggest that elevated ATP levels could be a factor in drug resistance to dabrafenib. CONCLUSIONS The model can be used to systematically motivate which dabrafenib-trametinib dose combinations, for treating BRAFV600E-mutated melanoma, warrant experimental investigation.
Collapse
|
12
|
Targeting Cellular DNA Damage Responses in Cancer: An In Vitro-Calibrated Agent-Based Model Simulating Monolayer and Spheroid Treatment Responses to ATR-Inhibiting Drugs. Bull Math Biol 2021; 83:103. [PMID: 34459993 PMCID: PMC8405495 DOI: 10.1007/s11538-021-00935-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2020] [Accepted: 08/10/2021] [Indexed: 11/26/2022]
Abstract
We combine a systems pharmacology approach with an agent-based modelling approach to simulate LoVo cells subjected to AZD6738, an ATR (ataxia–telangiectasia-mutated and rad3-related kinase) inhibiting anti-cancer drug that can hinder tumour proliferation by targeting cellular DNA damage responses. The agent-based model used in this study is governed by a set of empirically observable rules. By adjusting only the rules when moving between monolayer and multi-cellular tumour spheroid simulations, whilst keeping the fundamental mathematical model and parameters intact, the agent-based model is first parameterised by monolayer in vitro data and is thereafter used to simulate treatment responses in in vitro tumour spheroids subjected to dynamic drug delivery. Spheroid simulations are subsequently compared to in vivo data from xenografts in mice. The spheroid simulations are able to capture the dynamics of in vivo tumour growth and regression for approximately 8 days post-tumour injection. Translating quantitative information between in vitro and in vivo research remains a scientifically and financially challenging step in preclinical drug development processes. However, well-developed in silico tools can be used to facilitate this in vitro to in vivo translation, and in this article, we exemplify how data-driven, agent-based models can be used to bridge the gap between in vitro and in vivo research. We further highlight how agent-based models, that are currently underutilised in pharmaceutical contexts, can be used in preclinical drug development.
Collapse
|
13
|
Scott M, Żychaluk K, Bearon RN. A mathematical framework for modelling 3D cell motility: applications to glioblastoma cell migration. MATHEMATICAL MEDICINE AND BIOLOGY-A JOURNAL OF THE IMA 2021; 38:333-354. [PMID: 34189581 DOI: 10.1093/imammb/dqab009] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Revised: 06/09/2021] [Accepted: 06/10/2021] [Indexed: 11/14/2022]
Abstract
The collection of 3D cell tracking data from live images of micro-tissues is a recent innovation made possible due to advances in imaging techniques. As such there is increased interest in studying cell motility in 3D in vitro model systems but a lack of rigorous methodology for analysing the resulting data sets. One such instance of the use of these in vitro models is in the study of cancerous tumours. Growing multicellular tumour spheroids in vitro allows for modelling of the tumour microenvironment and the study of tumour cell behaviours, such as migration, which improves understanding of these cells and in turn could potentially improve cancer treatments. In this paper, we present a workflow for the rigorous analysis of 3D cell tracking data, based on the persistent random walk model, but adaptable to other biologically informed mathematical models. We use statistical measures to assess the fit of the model to the motility data and to estimate model parameters and provide confidence intervals for those parameters, to allow for parametrization of the model taking correlation in the data into account. We use in silico simulations to validate the workflow in 3D before testing our method on cell tracking data taken from in vitro experiments on glioblastoma tumour cells, a brain cancer with a very poor prognosis. The presented approach is intended to be accessible to both modellers and experimentalists alike in that it provides tools for uncovering features of the data set that may suggest amendments to future experiments or modelling attempts.
Collapse
Affiliation(s)
- M Scott
- Department of Mathematical Sciences, University of Liverpool, Liverpool L69 7ZL, UK
| | - K Żychaluk
- Department of Mathematical Sciences, University of Liverpool, Liverpool L69 7ZL, UK
| | - R N Bearon
- Department of Mathematical Sciences, University of Liverpool, Liverpool L69 7ZL, UK
| |
Collapse
|
14
|
Nardini JT, Stolz BJ, Flores KB, Harrington HA, Byrne HM. Topological data analysis distinguishes parameter regimes in the Anderson-Chaplain model of angiogenesis. PLoS Comput Biol 2021; 17:e1009094. [PMID: 34181657 PMCID: PMC8270459 DOI: 10.1371/journal.pcbi.1009094] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Revised: 07/09/2021] [Accepted: 05/18/2021] [Indexed: 12/27/2022] Open
Abstract
Angiogenesis is the process by which blood vessels form from pre-existing vessels. It plays a key role in many biological processes, including embryonic development and wound healing, and contributes to many diseases including cancer and rheumatoid arthritis. The structure of the resulting vessel networks determines their ability to deliver nutrients and remove waste products from biological tissues. Here we simulate the Anderson-Chaplain model of angiogenesis at different parameter values and quantify the vessel architectures of the resulting synthetic data. Specifically, we propose a topological data analysis (TDA) pipeline for systematic analysis of the model. TDA is a vibrant and relatively new field of computational mathematics for studying the shape of data. We compute topological and standard descriptors of model simulations generated by different parameter values. We show that TDA of model simulation data stratifies parameter space into regions with similar vessel morphology. The methodologies proposed here are widely applicable to other synthetic and experimental data including wound healing, development, and plant biology.
Collapse
Affiliation(s)
- John T. Nardini
- Department of Mathematics, North Carolina State University, Raleigh, North Carolina, United States of America
| | | | - Kevin B. Flores
- Department of Mathematics, North Carolina State University, Raleigh, North Carolina, United States of America
| | - Heather A. Harrington
- Mathematical Institute, University of Oxford, Oxford, United Kingdom
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, United Kingdom
| | - Helen M. Byrne
- Mathematical Institute, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
15
|
Sosa-Marrero C, de Crevoisier R, Hernandez A, Fontaine P, Rioux-Leclercq N, Mathieu R, Fautrel A, Paris F, Acosta O. Towards a Reduced In Silico Model Predicting Biochemical Recurrence After Radiotherapy in Prostate Cancer. IEEE Trans Biomed Eng 2021; 68:2718-2729. [PMID: 33460366 DOI: 10.1109/tbme.2021.3052345] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
OBJECTIVE Purposes of this work were i) to develop an in silico model of tumor response to radiotherapy, ii) to perform an exhaustive sensitivity analysis in order to iii) propose a simplified version and iv) to predict biochemical recurrence with both the comprehensive and the reduced model. METHODS A multiscale computational model of tumor response to radiotherapy was developed. It integrated the following radiobiological mechanisms: oxygenation, including hypoxic death; division of tumor cells; VEGF diffusion driving angiogenesis; division of healthy cells and oxygen-dependent response to irradiation, considering, cycle arrest and mitotic catastrophe. A thorough sensitivity analysis using the Morris screening method was performed on 21 prostate computational tissues. Tumor control probability (TCP) curves of the comprehensive model and 15 reduced versions were compared. Logistic regression was performed to predict biochemical recurrence after radiotherapy on 76 localized prostate cancer patients using an output of the comprehensive and the reduced models. RESULTS No significant difference was found between the TCP curves of the comprehensive and a simplified version which only considered oxygenation, division of tumor cells and their response to irradiation. Biochemical recurrence predictions using the comprehensive and the reduced models improved those made from pre-treatment imaging parameters (AUC = 0.81 ± 0.02 and 0.82 ± 0.02 vs. 0.75 ± 0.03, respectively). CONCLUSION A reduced model of tumor response to radiotherapy able to predict biochemical recurrence in prostate cancer was obtained. SIGNIFICANCE This reduced model may be used in the future to optimize personalized fractionation schedules.
Collapse
|
16
|
Hamis S, Kohandel M, Dubois LJ, Yaromina A, Lambin P, Powathil GG. Combining hypoxia-activated prodrugs and radiotherapy in silico: Impact of treatment scheduling and the intra-tumoural oxygen landscape. PLoS Comput Biol 2020; 16:e1008041. [PMID: 32745136 PMCID: PMC7425994 DOI: 10.1371/journal.pcbi.1008041] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2019] [Revised: 08/13/2020] [Accepted: 06/11/2020] [Indexed: 12/30/2022] Open
Abstract
Hypoxia-activated prodrugs (HAPs) present a conceptually elegant approach to not only overcome, but better yet, exploit intra-tumoural hypoxia. Despite being successful in vitro and in vivo, HAPs are yet to achieve successful results in clinical settings. It has been hypothesised that this lack of clinical success can, in part, be explained by the insufficiently stringent clinical screening selection of determining which tumours are suitable for HAP treatments. Taking a mathematical modelling approach, we investigate how tumour properties and HAP-radiation scheduling influence treatment outcomes in simulated tumours. The following key results are demonstrated in silico: (i) HAP and ionising radiation (IR) monotherapies may attack tumours in dissimilar, and complementary, ways. (ii) HAP-IR scheduling may impact treatment efficacy. (iii) HAPs may function as IR treatment intensifiers. (iv) The spatio-temporal intra-tumoural oxygen landscape may impact HAP efficacy. Our in silico framework is based on an on-lattice, hybrid, multiscale cellular automaton spanning three spatial dimensions. The mathematical model for tumour spheroid growth is parameterised by multicellular tumour spheroid (MCTS) data.
Collapse
Affiliation(s)
- Sara Hamis
- School of Mathematics and Statistics, University of St Andrews, St Andrews, Scotland
- Department of Mathematics, College of Science, Swansea University, Swansea, Wales, United Kingdom
| | - Mohammad Kohandel
- Department of Applied Mathematics, University of Waterloo, Waterloo, Canada
| | - Ludwig J. Dubois
- The M-Lab, Department of Precision Medicine, GROW - School for Oncology and Developmental Biology, Maastricht University, Maastricht, The Netherlands
| | - Ala Yaromina
- The M-Lab, Department of Precision Medicine, GROW - School for Oncology and Developmental Biology, Maastricht University, Maastricht, The Netherlands
| | - Philippe Lambin
- The M-Lab, Department of Precision Medicine, GROW - School for Oncology and Developmental Biology, Maastricht University, Maastricht, The Netherlands
| | - Gibin G. Powathil
- Department of Mathematics, College of Science, Swansea University, Swansea, Wales, United Kingdom
| |
Collapse
|
17
|
|
18
|
A general mathematical framework for understanding the behavior of heterogeneous stem cell regeneration. J Theor Biol 2020; 492:110196. [PMID: 32067937 DOI: 10.1016/j.jtbi.2020.110196] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2019] [Revised: 12/28/2019] [Accepted: 02/11/2020] [Indexed: 11/21/2022]
Abstract
Stem cell heterogeneity is essential for homeostasis in tissue development. This paper establishes a general mathematical framework to model the dynamics of stem cell regeneration with cell heterogeneity and random transitions of epigenetic states. The framework generalizes the classical G0 cell cycle model and incorporates the epigenetic states of individual cells represented by a continuous multidimensional variable. In the model, the kinetic rates of cell behaviors, including proliferation, differentiation, and apoptosis, are dependent on their epigenetic states, and the random transitions of epigenetic states between cell cycles are represented by an inheritance probability function that describes the conditional probability of cell state changes. Moreover, the model can be extended to include genotypic changes and describe the process of gene mutation-induced tumor development. The proposed mathematical framework provides a generalized formula that helps us to understand various dynamic processes of stem cell regeneration, including tissue development, degeneration, and abnormal growth.
Collapse
|
19
|
Heudobler D, Lüke F, Vogelhuber M, Klobuch S, Pukrop T, Herr W, Gerner C, Pantziarka P, Ghibelli L, Reichle A. Anakoinosis: Correcting Aberrant Homeostasis of Cancer Tissue-Going Beyond Apoptosis Induction. Front Oncol 2019; 9:1408. [PMID: 31921665 PMCID: PMC6934003 DOI: 10.3389/fonc.2019.01408] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2019] [Accepted: 11/28/2019] [Indexed: 12/16/2022] Open
Abstract
The current approach to systemic therapy for metastatic cancer is aimed predominantly at inducing apoptosis of cancer cells by blocking tumor-promoting signaling pathways or by eradicating cell compartments within the tumor. In contrast, a systems view of therapy primarily considers the communication protocols that exist at multiple levels within the tumor complex, and the role of key regulators of such systems. Such regulators may have far-reaching influence on tumor response to therapy and therefore patient survival. This implies that neoplasia may be considered as a cell non-autonomous disease. The multi-scale activity ranges from intra-tumor cell compartments, to the tumor, to the tumor-harboring organ to the organism. In contrast to molecularly targeted therapies, a systems approach that identifies the complex communications networks driving tumor growth offers the prospect of disrupting or "normalizing" such aberrant communicative behaviors and therefore attenuating tumor growth. Communicative reprogramming, a treatment strategy referred to as anakoinosis, requires novel therapeutic instruments, so-called master modifiers to deliver concerted tumor growth-attenuating action. The diversity of biological outcomes following pro-anakoinotic tumor therapy, such as differentiation, trans-differentiation, control of tumor-associated inflammation, etc. demonstrates that long-term tumor control may occur in multiple forms, inducing even continuous complete remission. Accordingly, pro-anakoinotic therapies dramatically extend the repertoire for achieving tumor control and may activate apoptosis pathways for controlling resistant metastatic tumor disease and hematologic neoplasia.
Collapse
Affiliation(s)
- Daniel Heudobler
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Florian Lüke
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Martin Vogelhuber
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Sebastian Klobuch
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Tobias Pukrop
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Wolfgang Herr
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Christopher Gerner
- Institut for Analytical Chemistry, Faculty Chemistry, University Vienna, Vienna, Austria
| | - Pan Pantziarka
- The George Pantziarka TP53 Trust, London, United Kingdom
- Anticancer Fund, Brussels, Belgium
| | - Lina Ghibelli
- Department Biology, Università di Roma Tor Vergata, Rome, Italy
| | - Albrecht Reichle
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| |
Collapse
|
20
|
Park DS, Robertson-Tessi M, Luddy KA, Maini PK, Bonsall MB, Gatenby RA, Anderson ARA. The Goldilocks Window of Personalized Chemotherapy: Getting the Immune Response Just Right. Cancer Res 2019; 79:5302-5315. [PMID: 31387920 DOI: 10.1158/0008-5472.can-18-3712] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2018] [Revised: 05/20/2019] [Accepted: 08/01/2019] [Indexed: 11/16/2022]
Abstract
The immune system is a robust and often untapped accomplice of many standard cancer therapies. A majority of tumors exist in a state of immune tolerance where the patient's immune system has become insensitive to the cancer cells. Because of its lymphodepleting effects, chemotherapy has the potential to break this tolerance. To investigate this, we created a mathematical modeling framework of tumor-immune dynamics. Our results suggest that optimal chemotherapy scheduling must balance two opposing objectives: maximizing tumor reduction while preserving patient immune function. Successful treatment requires therapy to operate in a "Goldilocks Window" where patient immune health is not overly compromised. By keeping therapy "just right," we show that the synergistic effects of immune activation and chemotherapy can maximize tumor reduction and control. SIGNIFICANCE: To maximize the synergy between chemotherapy and antitumor immune response, lymphodepleting therapy must be balanced in a "Goldilocks Window" of optimal dosing.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/79/20/5302/F1.large.jpg.
Collapse
Affiliation(s)
- Derek S Park
- Department of Zoology, University of Oxford, Oxford, United Kingdom. .,Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, Florida
| | - Mark Robertson-Tessi
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, Florida
| | - Kimberly A Luddy
- Comparative Immunology Group, School of Biochemistry and Immunology, Trinity College Dublin, Dublin, Ireland.,Department of Cancer Physiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, Florida
| | - Philip K Maini
- Mathematical Institute, University of Oxford, Oxford, Oxfordshire, United Kingdom
| | | | - Robert A Gatenby
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, Florida.,Department of Radiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, Florida
| | - Alexander R A Anderson
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, Florida.
| |
Collapse
|