When Seeing Isn't Believing: How Math Can Guide Our Interpretation of Measurements and Experiments

A multiscale model of immune surveillance in micrometastases gives insights on cancer patient digital twins

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.

In silico study of heterogeneous tumour-derived organoid response to CAR T-cell therapy

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.

A review of mechanistic learning in mathematical oncology

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.

Systematic comparison of modeling fidelity levels and parameter inference settings applied to negative feedback gene regulation

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.

The impact of tumor associated macrophages on tumor biology under the lens of mathematical modelling: A review

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.

Supporting Computational Apprenticeship Through Educational and Software Infrastructure: A Case Study in a Mathematical Oncology Research Lab

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.

Iterative community-driven development of a SARS-CoV-2 tissue simulator

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.

A persistent invasive phenotype in post-hypoxic tumor cells is revealed by fate mapping and computational modeling

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

Biofabrication of spheroids fusion-based tumor models: computational simulation of glucose effects

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.

Using mathematics in MRI data management for glioma assesment

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.

Key challenges facing data-driven multicellular systems biology

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.

Learning-accelerated discovery of immune-tumour interactions

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.

A Review of Cell-Based Computational Modeling in Cancer Biology

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.

Neural crest migration with continuous cell states

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.

PhysiCell: An open source physics-based cell simulator for 3-D multicellular systems

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.