Ahmadreza Ghaffarizadeh

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.

Progress Towards Computational 3-D Multicellular Systems Biology

Tumors cannot be understood in isolation from their microenvironment. Tumor and stromal cells change phenotype based upon biochemical and biophysical inputs from their surroundings, even as they interact with and remodel the microenvironment. Cancer should be investigated as an adaptive, multicellular system in a dynamical microenvironment. Computational modeling offers the potential to detangle this complex system, but the modeling platform must ideally account for tumor heterogeneity, substrate and signaling factor biotransport, cell and tissue biophysics, tissue and vascular remodeling, microvascular and interstitial flow, and links between all these sub-systems. Such a platform should leverage high-throughput experimental data, while using open data standards for reproducibility. In this chapter, we review advances by our groups in these key areas, particularly in advanced models of tissue mechanics and interstitial flow, open source simulation software, high-throughput phenotypic screening, and multicellular data standards. In the future, we expect a transformation of computational cancer biology from individual groups modeling isolated parts of cancer, to coalitions of groups combining compatible tools to simulate the 3-D multicellular systems biology of cancer tissues.

BioFVM: an efficient, parallelized diffusive transport solver for 3-D biological simulations

Motivation: Computational models of multicellular systems require solving systems of PDEs for release, uptake, decay and diffusion of multiple substrates in 3D, particularly when incorporating the impact of drugs, growth substrates and signaling factors on cell receptors and subcellular systems biology. Results: We introduce BioFVM, a diffusive transport solver tailored to biological problems. BioFVM can simulate release and uptake of many substrates by cell and bulk sources, diffusion and decay in large 3D domains. It has been parallelized with OpenMP, allowing efficient simulations on desktop workstations or single supercomputer nodes. The code is stable even for large time steps, with linear computational cost scalings. Solutions are first-order accurate in time and second-order accurate in space. The code can be run by itself or as part of a larger simulator. Availability and implementation: BioFVM is written in C ++ with parallelization in OpenMP. It is maintained and available for download at http://BioFVM.MathCancer.org and http://BioFVM.sf.net under the Apache License (v2.0). Contact: paul.macklin@usc.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Agent-based simulation of large tumors in 3-D microenvironments

Multicellular simulations of tumor growth in complex 3-D tissues, where data come from high content in vitro and bioengineered experiments, have gained significant attention by the cancer modeling community in recent years. Agent-based models are often selected for these problems because they can directly model and track cells’ states and their interactions with the microenvironment. We describe PhysiCell, a specific agent-based model that includes cell motion, cell cycling, and cell volume changes. The model has been performance tested on systems of 105 cells on desktop computers, and is expected to scale to 106 or more cells on single super-computer compute nodes. We plan an open source release of the software in early 2016 at PhysiCell.MathCancer.org

MultiCellDS: a community-developed standard for curating microenvironment-dependent multicellular data

Exchanging and understanding scientific data and their context represents a significant barrier to advancing research, especially with respect to information siloing. Maintaining information provenance and providing data curation and quality control help overcome common concerns and barriers to the effective sharing of scientific data. To address these problems in and the unique challenges of multicellular systems, we assembled a panel composed of investigators from several disciplines to create the MultiCellular Data Standard (MultiCellDS) with a use-case driven development process. The standard includes (1) digital cell lines, which are analogous to traditional biological cell lines, to record metadata, cellular microenvironment, and cellular phenotype variables of a biological cell line, (2) digital snapshots to consistently record simulation, experimental, and clinical data for multicellular systems, and (3) collections that can logically group digital cell lines and snapshots. We have created a MultiCellular DataBase (MultiCellDB) to store digital snapshots and the 200+ digital cell lines we have generated. MultiCellDS, by having a fixed standard, enables discoverability, extensibility, maintainability, searchability, and sustainability of data, creating biological applicability and clinical utility that permits us to identify upcoming challenges to uplift biology and strategies and therapies for improving human health.

High-throughput cancer hypothesis testing with an integrated PhysiCell-EMEWS workflow

Cancer is a complex, multiscale dynamical system, with interactions between tumor cells and non-cancerous systems. Therapies act on this cancer-host system, sometimes with unexpected results. Systematic investigation of mechanistic models could help identify the factors driving a treatment’s success or failure, but exploring mechanistic models over high-dimensional parameter spaces is computationally challenging. In this paper, we introduce a high throughput computing (HTC) framework that integrates a mechanistic 3-D multicellular simulator (PhysiCell) with an extreme-scale model exploration platform (EMEWS) to investigate high-dimensional parameter spaces. We show early results in adapting PhysiCell-EMEWS to 3-D cancer immunotherapy and show insights on therapeutic failure. We describe a PhysiCell-EMEWS work-flow for high-throughput cancer hypothesis testing, where thousands of mechanistic simulations are compared against data-driven error metrics to perform hypothesis optimization. We close by discussing novel applications to synthetic multicellular systems for cancer therapy.

MultiCellDS: a standard and a community for sharing multicellular data

Cell biology is increasingly focused on cellular heterogeneity and multicellular systems. To make the fullest use of experimental, clinical, and computational efforts, we need standardized data formats, community-curated “public data libraries”, and tools to combine and analyze shared data. To address these needs, our multidisciplinary community created MultiCellDS (MultiCellular Data Standard): an extensible standard, a library of digital cell lines and tissue snapshots, and support software. With the help of experimentalists, clinicians, modelers, and data and library scientists, we can grow this seed into a community-owned ecosystem of shared data and tools, to the benefit of basic science, engineering, and human health.

Estimating cell cycle model parameters using systems identification

A current challenge in data-driven mathematical modeling of cancer is identifying biologically-relevant parameters of mathematical models from sparse and often noisy experimental data of mixed types. We describe a cell cycle model and outline how to use the Optimization Toolbox in Matlab to estimate its timescale parameters, given flow cytometry and cell viability (synthetic) data, and illustrate the technique with simulated data. This technique can be similarly applied to a variety of cell cycle models, particularly as more laboratories begin to use high-content, quantitative cell screening and imaging platforms. An advanced version of this work (CellPD: cell line phenotype digitizer) will be released as open source in early 2016 at MultiCellDS.org.

Simulating multi-substrate diffusive transport in 3-D tissues with BioFVM

To simulate the spatiotemporal distribution of chemical compounds, we present BioFVM, an open-source reaction-diffusion equation solver using finite volume methods with motivation for biological applications. With various numerical solvers, we can simulate the interaction of dozens of compounds, including growth substrates, drugs, and signaling compounds in 3-D tissues, with cells by treating them as various source/sink terms. BioFVM has linear computational cost scalings and demonstrates first-order accuracy in time and second-order accuracy in space. Beyond simulating the transport of drugs and growth substrates in tissues, the ability to simulate dozens of compounds should make 3-D simulations of multicellular secretomics feasible.