Building multiscale models with PhysiBoSS, an agent-based modeling tool.

Multiscale models provide a unique tool for analyzing complex processes that study events occurring at different scales across space and time. In the context of biological systems, such models can simulate mechanisms happening at the intracellular level such as signaling, and at the extracellular level where cells communicate and coordinate with other cells. These models aim to understand the impact of genetic or environmental deregulation observed in complex diseases, describe the interplay between a pathological tissue and the immune system, and suggest strategies to revert the diseased phenotypes. The construction of these multiscale models remains a very complex task, including the choice of the components to consider, the level of details of the processes to simulate, or the fitting of the parameters to the data. One additional difficulty is the expert knowledge needed to program these models in languages such as C++ or Python, which may discourage the participation of non-experts. Simplifying this process through structured description formalisms-coupled with a graphical interface-is crucial in making modeling more accessible to the broader scientific community, as well as streamlining the process for advanced users. This article introduces three examples of multiscale models which rely on the framework PhysiBoSS, an add-on of PhysiCell that includes intracellular descriptions as continuous time Boolean models to the agent-based approach. The article demonstrates how to construct these models more easily, relying on PhysiCell Studio, the PhysiCell Graphical User Interface. A step-by-step tutorial is provided as Supplementary Material and all models are provided at https://physiboss.github.io/tutorial/.

Building multiscale models with PhysiBoSS, an agent-based modeling tool.

Multiscale models provide a unique tool for studying complex processes that study events occurring at different scales across space and time. In the context of biological systems, such models can simulate mechanisms happening at the intracellular level such as signaling, and at the extracellular level where cells communicate and coordinate with other cells. They aim to understand the impact of genetic or environmental deregulation observed in complex diseases, describe the interplay between a pathological tissue and the immune system, and suggest strategies to revert the diseased phenotypes. The construction of these multiscale models remains a very complex task, including the choice of the components to consider, the level of details of the processes to simulate, or the fitting of the parameters to the data. One additional difficulty is the expert knowledge needed to program these models in languages such as C++ or Python, which may discourage the participation of non-experts. Simplifying this process through structured description formalisms - coupled with a graphical interface - is crucial in making modeling more accessible to the broader scientific community, as well as streamlining the process for advanced users. This article introduces three examples of multiscale models which rely on the framework PhysiBoSS, an add-on of PhysiCell that includes intracellular descriptions as continuous time Boolean models to the agent-based approach. The article demonstrates how to easily construct such models, relying on PhysiCell Studio, the PhysiCell Graphical User Interface. A step-by-step tutorial is provided as a Supplementary Material and all models are provided at: https://physiboss.github.io/tutorial/.

PhysiCell Studio: a graphical tool to make agent-based modeling more accessible.

Defining a multicellular model can be challenging. There may be hundreds of parameters that specify the attributes and behaviors of objects. In the best case, the model will be defined using some format specification - a markup language - that will provide easy model sharing (and a minimal step toward reproducibility). PhysiCell is an open-source, physics-based multicellular simulation framework with an active and growing user community. It uses XML to define a model and, traditionally, users needed to manually edit the XML to modify the model. PhysiCell Studio is a tool to make this task easier. It provides a GUI that allows editing the XML model definition, including the creation and deletion of fundamental objects: cell types and substrates in the microenvironment. It also lets users build their model by defining initial conditions and biological rules, run simulations, and view results interactively. PhysiCell Studio has evolved over multiple workshops and academic courses in recent years, which has led to many improvements. There is both a desktop and cloud version. Its design and development has benefited from an active undergraduate and graduate research program. Like PhysiCell, the Studio is open-source software and contributions from the community are encouraged.

Merging Metabolic Modeling and Imaging for Screening Therapeutic Targets in Colorectal Cancer.

Cancer-associated fibroblasts (CAFs) play a key role in metabolic reprogramming and are well-established contributors to drug resistance in colorectal cancer (CRC). To exploit this metabolic crosstalk, we integrated a systems biology approach that identified key metabolic targets in a data-driven method and validated them experimentally. This process involved a novel machine learning-based method to computationally screen, in a high-throughput manner, the effects of enzyme perturbations predicted by a computational model of CRC metabolism. This approach reveals the network-wide effects of metabolic perturbations. Our results highlighted hexokinase (HK) as the crucial target, which subsequently became our focus for experimental validation using patient-derived tumor organoids (PDTOs). Through metabolic imaging and viability assays, we found that PDTOs cultured in CAF-conditioned media exhibited increased sensitivity to HK inhibition, confirming the model predictions. Our approach emphasizes the critical role of integrating computational and experimental techniques in exploring and exploiting CRC-CAF crosstalk.

Drug-loaded nanoparticles for cancer therapy: a high-throughput multicellular agent-based modeling study.

Interactions between biological systems and engineered nanomaterials have become an important area of study due to the application of nanomaterials in medicine. In particular, the application of nanomaterials for cancer diagnosis or treatment presents a challenging opportunity due to the complex biology of this disease spanning multiple time and spatial scales. A system-level analysis would benefit from mathematical modeling and computational simulation to explore the interactions between anticancer drug-loaded nanoparticles (NPs), cells, and tissues, and the associated parameters driving this system and a patient's overall response. Although a number of models have explored these interactions in the past, few have focused on simulating individual cell-NP interactions. This study develops a multicellular agent-based model of cancer nanotherapy that simulates NP internalization, drug release within the cell cytoplasm, "inheritance" of NPs by daughter cells at cell division, cell pharmacodynamic response to the intracellular drug, and overall drug effect on tumor dynamics. A large-scale parallel computational framework is used to investigate the impact of pharmacokinetic design parameters (NP internalization rate, NP decay rate, anticancer drug release rate) and therapeutic strategies (NP doses and injection frequency) on the tumor dynamics. In particular, through the exploration of NP "inheritance" at cell division, the results indicate that cancer treatment may be improved when NPs are inherited at cell division for cytotoxic chemotherapy. Moreover, smaller dosage of cytostatic chemotherapy may also improve inhibition of tumor growth when cell division is not completely inhibited. This work suggests that slow delivery by "heritable" NPs can drive new dimensions of nanotherapy design for more sustained therapeutic response.

Toward mechanistic medical digital twins: some use cases in immunology.

A fundamental challenge for personalized medicine is to capture enough of the complexity of an individual patient to determine an optimal way to keep them healthy or restore their health. This will require personalized computational models of sufficient resolution and with enough mechanistic information to provide actionable information to the clinician. Such personalized models are increasingly referred to as medical digital twins. Digital twin technology for health applications is still in its infancy, and extensive research and development is required. This article focuses on several projects in different stages of development that can lead to specific-and practical-medical digital twins or digital twin modeling platforms. It emerged from a two-day forum on problems related to medical digital twins, particularly those involving an immune system component. Open access video recordings of the forum discussions are available.

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.

Forum on immune digital twins: a meeting report.

Medical digital twins are computational models of human biology relevant to a given medical condition, which are tailored to an individual patient, thereby predicting the course of disease and individualized treatments, an important goal of personalized medicine. The immune system, which has a central role in many diseases, is highly heterogeneous between individuals, and thus poses a major challenge for this technology. In February 2023, an international group of experts convened for two days to discuss these challenges related to immune digital twins. The group consisted of clinicians, immunologists, biologists, and mathematical modelers, representative of the interdisciplinary nature of medical digital twin development. A video recording of the entire event is available. This paper presents a synopsis of the discussions, brief descriptions of ongoing digital twin projects at different stages of progress. It also proposes a 5-year action plan for further developing this technology. The main recommendations are to identify and pursue a small number of promising use cases, to develop stimulation-specific assays of immune function in a clinical setting, and to develop a database of existing computational immune models, as well as advanced modeling technology and infrastructure.

An agent-based modeling approach for lung fibrosis in response to COVID-19.

The severity of the COVID-19 pandemic has created an emerging need to investigate the long-term effects of infection on patients. Many individuals are at risk of suffering pulmonary fibrosis due to the pathogenesis of lung injury and impairment in the healing mechanism. Fibroblasts are the central mediators of extracellular matrix (ECM) deposition during tissue regeneration, regulated by anti-inflammatory cytokines including transforming growth factor beta (TGF-ß). The TGF-ß-dependent accumulation of fibroblasts at the damaged site and excess fibrillar collagen deposition lead to fibrosis. We developed an open-source, multiscale tissue simulator to investigate the role of TGF-ß sources in the progression of lung fibrosis after SARS-CoV-2 exposure, intracellular viral replication, infection of epithelial cells, and host immune response. Using the model, we predicted the dynamics of fibroblasts, TGF-ß, and collagen deposition for 15 days post-infection in virtual lung tissue. Our results showed variation in collagen area fractions between 2% and 40% depending on the spatial behavior of the sources (stationary or mobile), the rate of activation of TGF-ß, and the duration of TGF-ß sources. We identified M2 macrophages as primary contributors to higher collagen area fraction. Our simulation results also predicted fibrotic outcomes even with lower collagen area fraction when spatially-localized latent TGF-ß sources were active for longer times. We validated our model by comparing simulated dynamics for TGF-ß, collagen area fraction, and macrophage cell population with independent experimental data from mouse models. Our results showed that partial removal of TGF-ß sources changed the fibrotic patterns; in the presence of persistent TGF-ß sources, partial removal of TGF-ß from the ECM significantly increased collagen area fraction due to maintenance of chemotactic gradients driving fibroblast movement. The computational findings are consistent with independent experimental and clinical observations of collagen area fractions and cell population dynamics not used in developing the model. These critical insights into the activity of TGF-ß sources may find applications in the current clinical trials targeting TGF-ß for the resolution of lung fibrosis.

PhysiCell Studio: a graphical tool to make agent-based modeling more accessible.

Defining a multicellular model can be challenging. There may be hundreds of parameters that specify the attributes and behaviors of objects. Hopefully the model will be defined using some format specification, e.g., a markup language, that will provide easy model sharing (and a minimal step toward reproducibility). PhysiCell is an open source, physics-based multicellular simulation framework with an active and growing user community. It uses XML to define a model and, traditionally, users needed to manually edit the XML to modify the model. PhysiCell Studio is a tool to make this task easier. It provides a graphical user interface that allows editing the XML model definition, including the creation and deletion of fundamental objects, e.g., cell types and substrates in the microenvironment. It also lets users build their model by defining initial conditions and biological rules, run simulations, and view results interactively. PhysiCell Studio has evolved over multiple workshops and academic courses in recent years which has led to many improvements. Its design and development has benefited from an active undergraduate and graduate research program. Like PhysiCell, the Studio is open source software and contributions from the community are encouraged.

Digitize your Biology! Modeling multicellular systems through interpretable cell behavior.

Cells are fundamental units of life, constantly interacting and evolving as dynamical systems. While recent spatial multi-omics can quantitate individual cells' characteristics and regulatory programs, forecasting their evolution ultimately requires mathematical modeling. We develop a conceptual framework-a cell behavior hypothesis grammar-that uses natural language statements (cell rules) to create mathematical models. This allows us to systematically integrate biological knowledge and multi-omics data to make them computable. We can then perform virtual "thought experiments" that challenge and extend our understanding of multicellular systems, and ultimately generate new testable hypotheses. In this paper, we motivate and describe the grammar, provide a reference implementation, and demonstrate its potential through a series of examples in tumor biology and immunotherapy. Altogether, this approach provides a bridge between biological, clinical, and systems biology researchers for mathematical modeling of biological systems at scale, allowing the community to extrapolate from single-cell characterization to emergent multicellular behavior.

An agent-based modeling approach for lung fibrosis in response to COVID-19.

The severity of the COVID-19 pandemic has created an emerging need to investigate the long-term effects of infection on patients. Many individuals are at risk of suffering pulmonary fibrosis due to the pathogenesis of lung injury and impairment in the healing mechanism. Fibroblasts are the central mediators of extracellular matrix (ECM) deposition during tissue regeneration, regulated by anti-inflammatory cytokines including transforming growth factor beta (TGF-ß). The TGF-ß-dependent accumulation of fibroblasts at the damaged site and excess fibrillar collagen deposition lead to fibrosis. We developed an open-source, multiscale tissue simulator to investigate the role of TGF-ß sources in the progression of lung fibrosis after SARS-CoV-2 exposure, intracellular viral replication, infection of epithelial cells, and host immune response. Using the model, we predicted the dynamics of fibroblasts, TGF-ß, and collagen deposition for 15 days post-infection in virtual lung tissue. Our results showed variation in collagen area fractions between 2% and 40% depending on the spatial behavior of the sources (stationary or mobile), the rate of activation of TGF-ß, and the duration of TGF-ß sources. We identified M2 macrophages as primary contributors to higher collagen area fraction. Our simulation results also predicted fibrotic outcomes even with lower collagen area fraction when spatially-localized latent TGF-ß sources were active for longer times. We validated our model by comparing simulated dynamics for TGF-ß, collagen area fraction, and macrophage cell population with independent experimental data from mouse models. Our results showed that partial removal of TGF-ß sources changed the fibrotic patterns; in the presence of persistent TGF-ß sources, partial removal of TGF-ß from the ECM significantly increased collagen area fraction due to maintenance of chemotactic gradients driving fibroblast movement. The computational findings are consistent with independent experimental and clinical observations of collagen area fractions and cell population dynamics not used in developing the model. These critical insights into the activity of TGF-ß sources may find applications in the current clinical trials targeting TGF-ß for the resolution of lung fibrosis. Author summary: COVID-19 survivors are at risk of lung fibrosis as a long-term effect. Lung fibrosis is the excess deposition of tissue materials in the lung that hinder gas exchange and can collapse the whole organ. We identified TGF-ß as a critical regulator of fibrosis. We built a model to investigate the mechanisms of TGF-ß sources in the process of fibrosis. Our results showed spatial behavior of sources (stationary or mobile) and their activity (activation rate of TGF-ß, longer activation of sources) could lead to lung fibrosis. Current clinical trials for fibrosis that target TGF-ß need to consider TGF-ß sources' spatial properties and activity to develop better treatment strategies.

Exploring approaches for predictive cancer patient digital twins: Opportunities for collaboration and innovation.

We are rapidly approaching a future in which cancer patient digital twins will reach their potential to predict cancer prevention, diagnosis, and treatment in individual patients. This will be realized based on advances in high performance computing, computational modeling, and an expanding repertoire of observational data across multiple scales and modalities. In 2020, the US National Cancer Institute, and the US Department of Energy, through a trans-disciplinary research community at the intersection of advanced computing and cancer research, initiated team science collaborative projects to explore the development and implementation of predictive Cancer Patient Digital Twins. Several diverse pilot projects were launched to provide key insights into important features of this emerging landscape and to determine the requirements for the development and adoption of cancer patient digital twins. Projects included exploring approaches to using a large cohort of digital twins to perform deep phenotyping and plan treatments at the individual level, prototyping self-learning digital twin platforms, using adaptive digital twin approaches to monitor treatment response and resistance, developing methods to integrate and fuse data and observations across multiple scales, and personalizing treatment based on cancer type. Collectively these efforts have yielded increased insights into the opportunities and challenges facing cancer patient digital twin approaches and helped define a path forward. Given the rapidly growing interest in patient digital twins, this manuscript provides a valuable early progress report of several CPDT pilot projects commenced in common, their overall aims, early progress, lessons learned and future directions that will increasingly involve the broader research community.

Agent-based computational modeling of glioblastoma predicts that stromal density is central to oncolytic virus efficacy.

Oncolytic viruses (OVs) are emerging cancer immunotherapy. Despite notable successes in the treatment of some tumors, OV therapy for central nervous system cancers has failed to show efficacy. We used an ex vivo tumor model developed from human glioblastoma tissue to evaluate the infiltration of herpes simplex OV rQNestin (oHSV-1) into glioblastoma tumors. We next leveraged our data to develop a computational, model of glioblastoma dynamics that accounts for cellular interactions within the tumor. Using our computational model, we found that low stromal density was highly predictive of oHSV-1 therapeutic success, suggesting that the efficacy of oHSV-1 in glioblastoma may be determined by stromal-to-tumor cell regional density. We validated these findings in heterogenous patient samples from brain metastatic adenocarcinoma. Our integrated modeling strategy can be applied to suggest mechanisms of therapeutic responses for central nervous system cancers and to facilitate the successful translation of OVs into the clinic.

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.

Elucidating tumor-stromal metabolic crosstalk in colorectal cancer through integration of constraint-based models and LC-MS metabolomics.

Colorectal cancer (CRC) is a major cause of morbidity and mortality in the United States. Tumor-stromal metabolic crosstalk in the tumor microenvironment promotes CRC development and progression, but exactly how stromal cells, in particular cancer-associated fibroblasts (CAFs), affect the metabolism of tumor cells remains unknown. Here we take a data-driven approach to investigate the metabolic interactions between CRC cells and CAFs, integrating constraint-based modeling and metabolomic profiling. Using metabolomics data, we perform unsteady-state parsimonious flux balance analysis to infer flux distributions for central carbon metabolism in CRC cells treated with or without CAF-conditioned media. We find that CAFs reprogram CRC metabolism through stimulation of glycolysis, the oxidative arm of the pentose phosphate pathway (PPP), and glutaminolysis, as well as inhibition of the tricarboxylic acid cycle. To identify potential therapeutic targets, we simulate enzyme knockouts and find that CAF-treated CRC cells are especially sensitive to inhibitions of hexokinase and glucose-6-phosphate, the rate limiting steps of glycolysis and oxidative PPP. Our work gives mechanistic insights into the metabolic interactions between CRC cells and CAFs and provides a framework for testing hypotheses towards CRC-targeted therapies.

Maraviroc inhibits SARS-CoV-2 multiplication and s-protein mediated cell fusion in cell culture.

In an effort to identify therapeutic intervention strategies for the treatment of COVID-19, we have investigated a selection of FDA-approved small molecules and biologics that are commonly used to treat other human diseases. A investigation into 18 small molecules and 3 biologics was conducted in cell culture and the impact of treatment on viral titer was quantified by plaque assay. The investigation identified 4 FDA-approved small molecules, Maraviroc, FTY720 (Fingolimod), Atorvastatin and Nitazoxanide that were able to inhibit SARS-CoV-2 infection. Confocal microscopy with over expressed S-protein demonstrated that Maraviroc reduced the extent of S-protein mediated cell fusion as observed by fewer multinucleate cells in the context of drug-treatment. Mathematical modeling of drug-dependent viral multiplication dynamics revealed that prolonged drug treatment will exert an exponential decrease in viral load in a multicellular/tissue environment. Taken together, the data demonstrate that Maraviroc, Fingolimod, Atorvastatin and Nitazoxanide inhibit SARS-CoV-2 in cell culture.

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.