1
|
L Rocha H, Aguilar B, Getz M, Shmulevich I, Macklin P. A multiscale model of immune surveillance in micrometastases gives insights on cancer patient digital twins. NPJ Syst Biol Appl 2024; 10:144. [PMID: 39627216 PMCID: PMC11614875 DOI: 10.1038/s41540-024-00472-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2024] [Accepted: 11/15/2024] [Indexed: 12/06/2024] Open
Abstract
Metastasis is the leading cause of death in patients with cancer, driving considerable scientific and clinical interest in immunosurveillance of micrometastases. We investigated this process by creating a multiscale mathematical model to study the interactions between the immune system and the progression of micrometastases in general epithelial tissue. We analyzed the parameter space of the model using high-throughput computing resources to generate over 100,000 virtual patient trajectories. We demonstrated that the model could recapitulate a wide variety of virtual patient trajectories, including uncontrolled growth, partial response, and complete immune response to tumor growth. We classified the virtual patients and identified key patient parameters with the greatest effect on the simulated immunosurveillance. We highlight the lessons derived from this analysis and their impact on the nascent field of cancer patient digital twins (CPDTs). While CPDTs could enable clinicians to systematically dissect the complexity of cancer in each individual patient and inform treatment choices, our work shows that key challenges remain before we can reach this vision. In particular, we show that there remain considerable uncertainties in immune responses, unreliable patient stratification, and unpredictable personalized treatment. Nonetheless, we also show that in spite of these challenges, patient-specific models suggest strategies to increase control of clinically undetectable micrometastases even without complete parameter certainty.
Collapse
Affiliation(s)
- Heber L Rocha
- Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA
| | | | - Michael Getz
- Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA
| | | | - Paul Macklin
- Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA.
| |
Collapse
|
2
|
Luque LM, Carlevaro CM, Rodriguez-Lomba E, Lomba E. In silico study of heterogeneous tumour-derived organoid response to CAR T-cell therapy. Sci Rep 2024; 14:12307. [PMID: 38811838 PMCID: PMC11137006 DOI: 10.1038/s41598-024-63125-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Accepted: 05/24/2024] [Indexed: 05/31/2024] Open
Abstract
Chimeric antigen receptor (CAR) T-cell therapy is a promising immunotherapy for treating cancers. This method consists in modifying the patients' T-cells to directly target antigen-presenting cancer cells. One of the barriers to the development of this type of therapies, is target antigen heterogeneity. It is thought that intratumour heterogeneity is one of the leading determinants of therapeutic resistance and treatment failure. While understanding antigen heterogeneity is important for effective therapeutics, a good therapy strategy could enhance the therapy efficiency. In this work we introduce an agent-based model (ABM), built upon a previous ABM, to rationalise the outcomes of different CAR T-cells therapies strategies over heterogeneous tumour-derived organoids. We found that one dose of CAR T-cell therapy should be expected to reduce the tumour size as well as its growth rate, however it may not be enough to completely eliminate it. Moreover, the amount of free CAR T-cells (i.e. CAR T-cells that did not kill any cancer cell) increases as we increase the dosage, and so does the risk of side effects. We tested different strategies to enhance smaller dosages, such as enhancing the CAR T-cells long-term persistence and multiple dosing. For both approaches an appropriate dosimetry strategy is necessary to produce "effective yet safe" therapeutic results. Moreover, an interesting emergent phenomenon results from the simulations, namely the formation of a shield-like structure of cells with low antigen expression. This shield turns out to protect cells with high antigen expression. Finally we tested a multi-antigen recognition therapy to overcome antigen escape and heterogeneity. Our studies suggest that larger dosages can completely eliminate the organoid, however the multi-antigen recognition increases the risk of side effects. Therefore, an appropriate small dosages dosimetry strategy is necessary to improve the outcomes. Based on our results, it is clear that a proper therapeutic strategy could enhance the therapies outcomes. In that direction, our computational approach provides a framework to model treatment combinations in different scenarios and to explore the characteristics of successful and unsuccessful treatments.
Collapse
Affiliation(s)
- Luciana Melina Luque
- Centre for Regenerative Medicine, University of Edinburgh, Edinburgh, EH16 4UU, UK.
| | - Carlos Manuel Carlevaro
- Instituto de Física de Líquidos y Sistemas Biológicos, Consejo Nacional de Investigaciones Científicas y Técnicas, 1900, La Plata, Argentina
- Departamento de Ingeniería Mecánica, Universidad Tecnológica Nacional, Facultad Regional La Plata, 1900, La Plata, Argentina
| | | | - Enrique Lomba
- Instituto de Química Física Blas Cabrera, Consejo Superior de Investigaciones Científicas, 28006, Madrid, Spain
| |
Collapse
|
3
|
Metzcar J, Jutzeler CR, Macklin P, Köhn-Luque A, Brüningk SC. A review of mechanistic learning in mathematical oncology. Front Immunol 2024; 15:1363144. [PMID: 38533513 PMCID: PMC10963621 DOI: 10.3389/fimmu.2024.1363144] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Accepted: 02/20/2024] [Indexed: 03/28/2024] Open
Abstract
Mechanistic learning refers to the synergistic combination of mechanistic mathematical modeling and data-driven machine or deep learning. This emerging field finds increasing applications in (mathematical) oncology. This review aims to capture the current state of the field and provides a perspective on how mechanistic learning may progress in the oncology domain. We highlight the synergistic potential of mechanistic learning and point out similarities and differences between purely data-driven and mechanistic approaches concerning model complexity, data requirements, outputs generated, and interpretability of the algorithms and their results. Four categories of mechanistic learning (sequential, parallel, extrinsic, intrinsic) of mechanistic learning are presented with specific examples. We discuss a range of techniques including physics-informed neural networks, surrogate model learning, and digital twins. Example applications address complex problems predominantly from the domain of oncology research such as longitudinal tumor response predictions or time-to-event modeling. As the field of mechanistic learning advances, we aim for this review and proposed categorization framework to foster additional collaboration between the data- and knowledge-driven modeling fields. Further collaboration will help address difficult issues in oncology such as limited data availability, requirements of model transparency, and complex input data which are embraced in a mechanistic learning framework.
Collapse
Affiliation(s)
- John Metzcar
- Intelligent Systems Engineering, Luddy School of Informatics, Computing, and Engineering, Bloomington, IN, United States
- Informatics, Luddy School of Informatics, Computing, and Engineering, Bloomington, IN, United States
| | - Catherine R. Jutzeler
- Department of Health Sciences and Technology (D-HEST), Eidgenössische Technische Hochschule Zürich (ETH), Zürich, Switzerland
- Swiss Institute of Bioinformatics (SIB), Lausanne, Switzerland
| | - Paul Macklin
- Intelligent Systems Engineering, Luddy School of Informatics, Computing, and Engineering, Bloomington, IN, United States
| | - Alvaro Köhn-Luque
- Oslo Centre for Biostatistics and Epidemiology, Faculty of Medicine, University of Oslo, Oslo, Norway
- Oslo Centre for Biostatistics and Epidemiology, Research Support Services, Oslo University Hospital, Oslo, Norway
| | - Sarah C. Brüningk
- Department of Health Sciences and Technology (D-HEST), Eidgenössische Technische Hochschule Zürich (ETH), Zürich, Switzerland
- Swiss Institute of Bioinformatics (SIB), Lausanne, Switzerland
| |
Collapse
|
4
|
Coulier A, Singh P, Sturrock M, Hellander A. Systematic comparison of modeling fidelity levels and parameter inference settings applied to negative feedback gene regulation. PLoS Comput Biol 2022; 18:e1010683. [PMID: 36520957 PMCID: PMC9799300 DOI: 10.1371/journal.pcbi.1010683] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2022] [Revised: 12/29/2022] [Accepted: 10/25/2022] [Indexed: 12/23/2022] Open
Abstract
Quantitative stochastic models of gene regulatory networks are important tools for studying cellular regulation. Such models can be formulated at many different levels of fidelity. A practical challenge is to determine what model fidelity to use in order to get accurate and representative results. The choice is important, because models of successively higher fidelity come at a rapidly increasing computational cost. In some situations, the level of detail is clearly motivated by the question under study. In many situations however, many model options could qualitatively agree with available data, depending on the amount of data and the nature of the observations. Here, an important distinction is whether we are interested in inferring the true (but unknown) physical parameters of the model or if it is sufficient to be able to capture and explain available data. The situation becomes complicated from a computational perspective because inference needs to be approximate. Most often it is based on likelihood-free Approximate Bayesian Computation (ABC) and here determining which summary statistics to use, as well as how much data is needed to reach the desired level of accuracy, are difficult tasks. Ultimately, all of these aspects-the model fidelity, the available data, and the numerical choices for inference-interplay in a complex manner. In this paper we develop a computational pipeline designed to systematically evaluate inference accuracy for a wide range of true known parameters. We then use it to explore inference settings for negative feedback gene regulation. In particular, we compare a detailed spatial stochastic model, a coarse-grained compartment-based multiscale model, and the standard well-mixed model, across several data-scenarios and for multiple numerical options for parameter inference. Practically speaking, this pipeline can be used as a preliminary step to guide modelers prior to gathering experimental data. By training Gaussian processes to approximate the distance function values, we are able to substantially reduce the computational cost of running the pipeline.
Collapse
Affiliation(s)
- Adrien Coulier
- Department of Information Technology, Uppsala University, Uppsala, Sweden
| | - Prashant Singh
- Science for Life Laboratory, Department of Information Technology, Uppsala University, Uppsala, Sweden
| | - Marc Sturrock
- Department of Physiology, Royal College of Surgeons in Ireland, Dublin, Ireland
| | - Andreas Hellander
- Department of Information Technology, Uppsala University, Uppsala, Sweden
- * E-mail:
| |
Collapse
|
5
|
Shojaee P, Mornata F, Deutsch A, Locati M, Hatzikirou H. The impact of tumor associated macrophages on tumor biology under the lens of mathematical modelling: A review. Front Immunol 2022; 13:1050067. [PMID: 36439180 PMCID: PMC9685623 DOI: 10.3389/fimmu.2022.1050067] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2022] [Accepted: 10/18/2022] [Indexed: 09/10/2023] Open
Abstract
In this article, we review the role of mathematical modelling to elucidate the impact of tumor-associated macrophages (TAMs) in tumor progression and therapy design. We first outline the biology of TAMs, and its current application in tumor therapies, and their experimental methods that provide insights into tumor cell-macrophage interactions. We then focus on the mechanistic mathematical models describing the role of macrophages as drug carriers, the impact of macrophage polarized activation on tumor growth, and the role of tumor microenvironment (TME) parameters on the tumor-macrophage interactions. This review aims to identify the synergies between biological and mathematical approaches that allow us to translate knowledge on fundamental TAMs biology in addressing current clinical challenges.
Collapse
Affiliation(s)
- Pejman Shojaee
- Centre for Information Services and High Performance Computing, Technische Universität (TU) Dresden, Dresden, Germany
| | - Federica Mornata
- Leukocyte Biology Lab, IRCCS Humanitas Research Hospital, Rozzano, Italy
| | - Andreas Deutsch
- Centre for Information Services and High Performance Computing, Technische Universität (TU) Dresden, Dresden, Germany
| | - Massimo Locati
- Leukocyte Biology Lab, IRCCS Humanitas Research Hospital, Rozzano, Italy
- Department of Medical Biotechnologies and Translational Medicine, Universitàdegli Studi di Milano, Milan, Italy
| | - Haralampos Hatzikirou
- Centre for Information Services and High Performance Computing, Technische Universität (TU) Dresden, Dresden, Germany
- Mathematics Department, Khalifa University, Abu Dhabi, United Arab Emirates
- Healthcare Engineering Innovation Centre (HEIC), Khalifa University, Abu Dhabi, United Arab Emirates
| |
Collapse
|
6
|
Madamanchi A, Thomas M, Magana A, Heiland R, Macklin P. Supporting Computational Apprenticeship Through Educational and Software Infrastructure: A Case Study in a Mathematical Oncology Research Lab. PRIMUS 2022; 32:446-467. [PMID: 35197716 PMCID: PMC8863170 DOI: 10.1080/10511970.2021.1881849] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
There is growing awareness of the need for mathematics and computing to quantitatively understand the complex dynamics and feedbacks in the life sciences. Although several institutions and research groups are conducting pioneering multidisciplinary research, communication and education across fields remain a bottleneck. The opportunity is ripe for using education research-supported mechanisms of cross-disciplinary training at the intersection of mathematics, computation, and biology. This case study uses the computational apprenticeship theoretical framework to describe the efforts of a computational biology lab to rapidly prototype, test, and refine a mentorship infrastructure for undergraduate research experiences. We describe the challenges, benefits, and lessons learned, as well as the utility of the computational apprenticeship framework in supporting computational/math students learning and contributing to biology, and biologists in learning computational methods. We also explore implications for undergraduate classroom instruction and cross-disciplinary scientific communication.
Collapse
|
7
|
Getz M, Wang Y, An G, Asthana M, Becker A, Cockrell C, Collier N, Craig M, Davis CL, Faeder JR, Ford Versypt AN, Mapder T, Gianlupi JF, Glazier JA, Hamis S, Heiland R, Hillen T, Hou D, Islam MA, Jenner AL, Kurtoglu F, Larkin CI, Liu B, Macfarlane F, Maygrundter P, Morel PA, Narayanan A, Ozik J, Pienaar E, Rangamani P, Saglam AS, Shoemaker JE, Smith AM, Weaver JJA, Macklin P. Iterative community-driven development of a SARS-CoV-2 tissue simulator. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2021:2020.04.02.019075. [PMID: 32511322 PMCID: PMC7239052 DOI: 10.1101/2020.04.02.019075] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The 2019 novel coronavirus, SARS-CoV-2, is a pathogen of critical significance to international public health. Knowledge of the interplay between molecular-scale virus-receptor interactions, single-cell viral replication, intracellular-scale viral transport, and emergent tissue-scale viral propagation is limited. Moreover, little is known about immune system-virus-tissue interactions and how these can result in low-level (asymptomatic) infections in some cases and acute respiratory distress syndrome (ARDS) in others, particularly with respect to presentation in different age groups or pre-existing inflammatory risk factors. Given the nonlinear interactions within and among each of these processes, multiscale simulation models can shed light on the emergent dynamics that lead to divergent outcomes, identify actionable "choke points" for pharmacologic interventions, screen potential therapies, and identify potential biomarkers that differentiate patient outcomes. Given the complexity of the problem and the acute need for an actionable model to guide therapy discovery and optimization, we introduce and iteratively refine a prototype of a multiscale model of SARS-CoV-2 dynamics in lung tissue. The first prototype model was built and shared internationally as open source code and an online interactive model in under 12 hours, and community domain expertise is driving regular refinements. In a sustained community effort, this consortium is integrating data and expertise across virology, immunology, mathematical biology, quantitative systems physiology, cloud and high performance computing, and other domains to accelerate our response to this critical threat to international health. More broadly, this effort is creating a reusable, modular framework for studying viral replication and immune response in tissues, which can also potentially be adapted to related problems in immunology and immunotherapy.
Collapse
|
8
|
Rocha HL, Godet I, Kurtoglu F, Metzcar J, Konstantinopoulos K, Bhoyar S, Gilkes DM, Macklin P. A persistent invasive phenotype in post-hypoxic tumor cells is revealed by fate mapping and computational modeling. iScience 2021; 24:102935. [PMID: 34568781 PMCID: PMC8449249 DOI: 10.1016/j.isci.2021.102935] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Revised: 04/23/2021] [Accepted: 07/29/2021] [Indexed: 12/03/2022] Open
Abstract
Hypoxia is a critical factor in solid tumors that has been associated with cancer progression and aggressiveness. We recently developed a hypoxia fate mapping system to trace post-hypoxic cells within a tumor for the first time. This approach uses an oxygen-dependent fluorescent switch and allowed us to measure key biological features such as oxygen distribution, cell proliferation, and migration. We developed a computational model to investigate the motility and phenotypic persistence of hypoxic and post-hypoxic cells during tumor progression. The cellular behavior was defined by phenotypic persistence time, cell movement bias, and the fraction of cells that respond to an enhanced migratory stimulus. This work combined advanced cell tracking and imaging techniques with mathematical modeling, to reveal that a persistent invasive migratory phenotype that develops under hypoxia is required for cellular escape into the surrounding tissue, promoting the formation of invasive structures (“plumes”) that expand toward the oxygenated tumor regions. A fluorescent fate mapping system allows tracking of hypoxic and post-hypoxic cells Computational modeling predicts the formation of post-hypoxic invasive plumes Simulations show post-hypoxic cells must maintain persistant migration to form plumes Tracking cells exposed to intratumoral hypoxia confirms persistent migration
Collapse
Affiliation(s)
- Heber L Rocha
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN 47408, USA
| | - Inês Godet
- Department of Oncology, The Sidney Kimmel Comprehensive Cancer Center, The Johns Hopkins University School of Medicine, Baltimore, MD 21231, USA.,Department of Chemical and Biomolecular Engineering, The Johns Hopkins University, Baltimore, MD 21218, USA
| | - Furkan Kurtoglu
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN 47408, USA
| | - John Metzcar
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN 47408, USA.,Department of Informatics, Indiana University, Bloomington, IN 47408, USA
| | - Kali Konstantinopoulos
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN 47408, USA
| | - Soumitra Bhoyar
- Department of Oncology, The Sidney Kimmel Comprehensive Cancer Center, The Johns Hopkins University School of Medicine, Baltimore, MD 21231, USA.,Department of Chemical and Biomolecular Engineering, The Johns Hopkins University, Baltimore, MD 21218, USA
| | - Daniele M Gilkes
- Department of Oncology, The Sidney Kimmel Comprehensive Cancer Center, The Johns Hopkins University School of Medicine, Baltimore, MD 21231, USA.,Department of Chemical and Biomolecular Engineering, The Johns Hopkins University, Baltimore, MD 21218, USA.,Cellular and Molecular Medicine Program, The Johns Hopkins University School of Medicine, Baltimore, MD 21231, USA
| | - Paul Macklin
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN 47408, USA
| |
Collapse
|
9
|
Bustamante DJ, Basile EJ, Hildreth BM, Browning NW, Jensen SA, Moldovan L, Petrache HI, Moldovan NI. Biofabrication of spheroids fusion-based tumor models: computational simulation of glucose effects. Biofabrication 2021; 13. [PMID: 33498017 DOI: 10.1088/1758-5090/abe025] [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: 06/30/2020] [Accepted: 01/26/2021] [Indexed: 01/04/2023]
Abstract
In vitro tumor models consisting of cell spheroids are increasingly used for mechanistic studies and pharmacological testing. However, unless vascularized, the availability of nutrients such as glucose to deeper layers of multicellular aggregates is limited. In addition, recent developments in cells-only biofabrication (e.g. 'scaffold-free bioprinting'), allow the creation of more complex spheroid-based structures, further exposing the cells to nutrient deprivation within these constructs. To explore the impact of glucose availability on such tumor-like structures, we used the CompuCell3D (CC3D) platform for modeling of tumor spheroids. By monitoring the types of cells, fusing pairs geometry and the distance between spheroids centers of mass, we made novel heuristic observations on how binary- and multi-spheroid fusions are impacted by glucose availability. At limiting glucose concentrations mimicking hypoglycemia we noted an abrupt collapse of the tumor spheroids, unexpectedly amplified by the contact with normal cell spheroids. At higher glucose concentrations, we found an increased intermixing of cancerous cells, strong anti-phase oscillations between proliferating and quiescent tumor cells and a structural instability of fusing tumor spheroids, leading to their re-fragmentation. In a model of tumor microenvironment composed of normal cell spheroids fusing around a tumoral one, the competition for glucose lead to either the tumor's disappearance, or to its steady expansion. Moreover, the invasion of this microenvironment by individual tumor cells was also strongly depended on the available glucose. In conclusion, we demonstrate the value of computational simulations for anticipating the properties of biofabricated tumor models, and in generating testable hypotheses regarding the relationship between cancer, nutrition and diabetes.
Collapse
Affiliation(s)
- David J Bustamante
- IUPUI BME, 799 W. Michigan Street, Indianapolis, Indiana, 46202-5195, UNITED STATES
| | - Elijah J Basile
- IUPUI, 301 University Boulevard, Indianapolis, Indiana, 46202-5146, UNITED STATES
| | - Brady M Hildreth
- IUPUI, 301 University Boulevard, Indianapolis, Indiana, 46202-5146, UNITED STATES
| | - Nathan W Browning
- IUPUI, 301 University Boulevard, Indianapolis, Indiana, 46202-5146, UNITED STATES
| | - S Alexander Jensen
- IUPUI, 301 University Boulevard Suite, Indianapolis, Indiana, 46202-5146, UNITED STATES
| | - Leni Moldovan
- Surgery, Indiana University School of Medicine, 1481 W. 10th St., Room D-2008, Indianapolis, Indiana, 46202-5114, UNITED STATES
| | - Horia I Petrache
- Department of Physics, Indiana University - Purdue University at Indianapolis, 402 N. Blackford Street, LD 154, Indianapolis, Indiana, 46202, UNITED STATES
| | - Nicanor I Moldovan
- VA Medical Center, Indiana University Purdue University at Indianapolis, 1481 W. 10th St., Room C6128, Indianapolis, Indiana, 46202-5143, UNITED STATES
| |
Collapse
|
10
|
Perrillat-Mercerot A, Guillevin C, Miranville A, Guillevin R. Using mathematics in MRI data management for glioma assesment. J Neuroradiol 2019; 48:282-290. [PMID: 31811826 DOI: 10.1016/j.neurad.2019.11.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2019] [Revised: 10/25/2019] [Accepted: 11/26/2019] [Indexed: 12/01/2022]
Abstract
Our aim is to review the mathematical tools usefulness in MR data management for glioma diagnosis and treatment optimization. MRI does not give access to organs variations in hours or days. However a lot of multiparametric data are generated. Mathematics could help to override this paradox, the aim of this article is to show how. We first make a review on mathematical modelling using equations. Afterwards we present statistical analysis. We provide detailed examples in both sections. We finally conclude, giving some clues on in silico models.
Collapse
Affiliation(s)
- A Perrillat-Mercerot
- UMR CNRS 7348, SP2MI, équipe DACTIM-MIS, laboratoire de mathématiques et applications, université de Poitiers, boulevard Marie-et-Pierre-Curie, Téléport 2, 86962 Chasseneuil Futuroscope cedex, France.
| | - C Guillevin
- UMR CNRS 7348, SP2MI, équipe DACTIM-MIS, laboratoire de mathématiques et applications, université de Poitiers, boulevard Marie-et-Pierre-Curie, Téléport 2, 86962 Chasseneuil Futuroscope cedex, France; CHU de Poitiers, 2, rue de la Milétrie, 86021 Poitiers, France.
| | - A Miranville
- UMR CNRS 7348, SP2MI, équipe DACTIM-MIS, laboratoire de mathématiques et applications, université de Poitiers, boulevard Marie-et-Pierre-Curie, Téléport 2, 86962 Chasseneuil Futuroscope cedex, France.
| | - R Guillevin
- UMR CNRS 7348, SP2MI, équipe DACTIM-MIS, laboratoire de mathématiques et applications, université de Poitiers, boulevard Marie-et-Pierre-Curie, Téléport 2, 86962 Chasseneuil Futuroscope cedex, France; CHU de Poitiers, 2, rue de la Milétrie, 86021 Poitiers, France.
| |
Collapse
|
11
|
Macklin P. Key challenges facing data-driven multicellular systems biology. Gigascience 2019; 8:giz127. [PMID: 31648301 PMCID: PMC6812467 DOI: 10.1093/gigascience/giz127] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2018] [Revised: 09/27/2019] [Accepted: 09/30/2019] [Indexed: 12/17/2022] Open
Abstract
Increasingly sophisticated experiments, coupled with large-scale computational models, have the potential to systematically test biological hypotheses to drive our understanding of multicellular systems. In this short review, we explore key challenges that must be overcome to achieve robust, repeatable data-driven multicellular systems biology. If these challenges can be solved, we can grow beyond the current state of isolated tools and datasets to a community-driven ecosystem of interoperable data, software utilities, and computational modeling platforms. Progress is within our grasp, but it will take community (and financial) commitment.
Collapse
Affiliation(s)
- Paul Macklin
- Department of Intelligent Systems Engineering, Indiana University, 700 N Woodlawn Ave, Bloomington, IN 47408, USA
| |
Collapse
|
12
|
Ozik J, Collier N, Heiland R, An G, Macklin P. Learning-accelerated discovery of immune-tumour interactions. MOLECULAR SYSTEMS DESIGN & ENGINEERING 2019; 4:747-760. [PMID: 31497314 PMCID: PMC6690424 DOI: 10.1039/c9me00036d] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2019] [Accepted: 04/18/2019] [Indexed: 05/04/2023]
Abstract
We present an integrated framework for enabling dynamic exploration of design spaces for cancer immunotherapies with detailed dynamical simulation models on high-performance computing resources. Our framework combines PhysiCell, an open source agent-based simulation platform for cancer and other multicellular systems, and EMEWS, an open source platform for extreme-scale model exploration. We build an agent-based model of immunosurveillance against heterogeneous tumours, which includes spatial dynamics of stochastic tumour-immune contact interactions. We implement active learning and genetic algorithms using high-performance computing workflows to adaptively sample the model parameter space and iteratively discover optimal cancer regression regions within biological and clinical constraints.
Collapse
Affiliation(s)
- Jonathan Ozik
- Decision and Infrastructure Sciences , Argonne National Laboratory , 9700 S. Cass Ave , Lemont , IL 60439 , USA .
- Consortium for Advanced Science and Engineering , University of Chicago , 5801 S. Ellis Ave. , Chicago , IL 60637 , USA
| | - Nicholson Collier
- Decision and Infrastructure Sciences , Argonne National Laboratory , 9700 S. Cass Ave , Lemont , IL 60439 , USA .
- Consortium for Advanced Science and Engineering , University of Chicago , 5801 S. Ellis Ave. , Chicago , IL 60637 , USA
| | - Randy Heiland
- Intelligent Systems Engineering , Indiana University , 700 N. Woodlawn Avenue Bloomington , IN 47408 , USA
| | - Gary An
- The University of Vermont Medical Center , 111 Colchester Avenue , Burlington , VT 05401 , USA
| | - Paul Macklin
- Intelligent Systems Engineering , Indiana University , 700 N. Woodlawn Avenue Bloomington , IN 47408 , USA
| |
Collapse
|
13
|
Metzcar J, Wang Y, Heiland R, Macklin P. A Review of Cell-Based Computational Modeling in Cancer Biology. JCO Clin Cancer Inform 2019; 3:1-13. [PMID: 30715927 PMCID: PMC6584763 DOI: 10.1200/cci.18.00069] [Citation(s) in RCA: 205] [Impact Index Per Article: 34.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/28/2018] [Indexed: 12/14/2022] Open
Abstract
Cancer biology involves complex, dynamic interactions between cancer cells and their tissue microenvironments. Single-cell effects are critical drivers of clinical progression. Chemical and mechanical communication between tumor and stromal cells can co-opt normal physiologic processes to promote growth and invasion. Cancer cell heterogeneity increases cancer's ability to test strategies to adapt to microenvironmental stresses. Hypoxia and treatment can select for cancer stem cells and drive invasion and resistance. Cell-based computational models (also known as discrete models, agent-based models, or individual-based models) simulate individual cells as they interact in virtual tissues, which allows us to explore how single-cell behaviors lead to the dynamics we observe and work to control in cancer systems. In this review, we introduce the broad range of techniques available for cell-based computational modeling. The approaches can range from highly detailed models of just a few cells and their morphologies to millions of simpler cells in three-dimensional tissues. Modeling individual cells allows us to directly translate biologic observations into simulation rules. In many cases, individual cell agents include molecular-scale models. Most models also simulate the transport of oxygen, drugs, and growth factors, which allow us to link cancer development to microenvironmental conditions. We illustrate these methods with examples drawn from cancer hypoxia, angiogenesis, invasion, stem cells, and immunosurveillance. An ecosystem of interoperable cell-based simulation tools is emerging at a time when cloud computing resources make software easier to access and supercomputing resources make large-scale simulation studies possible. As the field develops, we anticipate that high-throughput simulation studies will allow us to rapidly explore the space of biologic possibilities, prescreen new therapeutic strategies, and even re-engineer tumor and stromal cells to bring cancer systems under control.
Collapse
|
14
|
Schumacher LJ. Neural crest migration with continuous cell states. J Theor Biol 2019; 481:84-90. [PMID: 30707976 DOI: 10.1016/j.jtbi.2019.01.029] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2018] [Revised: 01/24/2019] [Accepted: 01/28/2019] [Indexed: 01/09/2023]
Abstract
Models of cranial neural crest cell migration in cell-induced (or self-generated) gradients have included a division of labour into leader and follower migratory states, which undergo chemotaxis and contact guidance, respectively. Despite validated utility of these models through experimental perturbation of migration in the chick embryo and gene expression analysis showing relevant heterogeneity at the single cell level, an often raised question has been whether the discrete cell states are necessary, or if a continuum of cell behaviours offers a functionally equivalent description. Here we argue that this picture is supported by recent single-cell transcriptome data. Motivated by this, we implement two versions of a continuous-state model: (1) signal choice and (2) signal combination. We find that the cell population migrates further than in the discrete-state model and than in experimental observations. We further show that the signal combination model, but not the signal choice model, can be successfully adjusted to experimentally plausible regimes by reducing the chemoattractant consumption parameter. Thus we show an equivalently plausible, experimentally motivated, model of neural crest cell migration.
Collapse
Affiliation(s)
- Linus J Schumacher
- MCR Centre for Regenerative Medicine, University of Edinburgh, United Kingdom.
| |
Collapse
|
15
|
Ghaffarizadeh A, Heiland R, Friedman SH, Mumenthaler SM, Macklin P. PhysiCell: An open source physics-based cell simulator for 3-D multicellular systems. PLoS Comput Biol 2018; 14:e1005991. [PMID: 29474446 PMCID: PMC5841829 DOI: 10.1371/journal.pcbi.1005991] [Citation(s) in RCA: 233] [Impact Index Per Article: 33.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2016] [Revised: 03/07/2018] [Accepted: 01/19/2018] [Indexed: 02/07/2023] Open
Abstract
Many multicellular systems problems can only be understood by studying how cells move, grow, divide, interact, and die. Tissue-scale dynamics emerge from systems of many interacting cells as they respond to and influence their microenvironment. The ideal "virtual laboratory" for such multicellular systems simulates both the biochemical microenvironment (the "stage") and many mechanically and biochemically interacting cells (the "players" upon the stage). PhysiCell-physics-based multicellular simulator-is an open source agent-based simulator that provides both the stage and the players for studying many interacting cells in dynamic tissue microenvironments. It builds upon a multi-substrate biotransport solver to link cell phenotype to multiple diffusing substrates and signaling factors. It includes biologically-driven sub-models for cell cycling, apoptosis, necrosis, solid and fluid volume changes, mechanics, and motility "out of the box." The C++ code has minimal dependencies, making it simple to maintain and deploy across platforms. PhysiCell has been parallelized with OpenMP, and its performance scales linearly with the number of cells. Simulations up to 105-106 cells are feasible on quad-core desktop workstations; larger simulations are attainable on single HPC compute nodes. We demonstrate PhysiCell by simulating the impact of necrotic core biomechanics, 3-D geometry, and stochasticity on the dynamics of hanging drop tumor spheroids and ductal carcinoma in situ (DCIS) of the breast. We demonstrate stochastic motility, chemical and contact-based interaction of multiple cell types, and the extensibility of PhysiCell with examples in synthetic multicellular systems (a "cellular cargo delivery" system, with application to anti-cancer treatments), cancer heterogeneity, and cancer immunology. PhysiCell is a powerful multicellular systems simulator that will be continually improved with new capabilities and performance improvements. It also represents a significant independent code base for replicating results from other simulation platforms. The PhysiCell source code, examples, documentation, and support are available under the BSD license at http://PhysiCell.MathCancer.org and http://PhysiCell.sf.net.
Collapse
Affiliation(s)
- Ahmadreza Ghaffarizadeh
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, California, United States of America
| | - Randy Heiland
- Intelligent Systems Engineering, Indiana University, Bloomington, Indiana, United States of America
| | - Samuel H. Friedman
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, California, United States of America
- Opto-Knowledge Systems, Inc., Torrance, California, United States of America
| | - Shannon M. Mumenthaler
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, California, United States of America
| | - Paul Macklin
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, California, United States of America
- Intelligent Systems Engineering, Indiana University, Bloomington, Indiana, United States of America
| |
Collapse
|