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 10⁵-10⁶ 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. This paper introduces PhysiCell: an open source, agent-based modeling framework for 3-D multicellular simulations. It includes a standard library of sub-models for cell fluid and solid volume changes, cycle progression, apoptosis, necrosis, mechanics, and motility. PhysiCell is directly coupled to a biotransport solver to simulate many diffusing substrates and cell-secreted signals. Each cell can dynamically update its phenotype based on its microenvironmental conditions. Users can customize or replace the included sub-models. PhysiCell runs on a variety of platforms (Linux, OSX, and Windows) with few software dependencies. Its computational cost scales linearly in the number of cells. It is feasible to simulate 500,000 cells on quad-core desktop workstations, and millions of cells on single HPC compute nodes. We demonstrate PhysiCell by simulating the impact of necrotic core biomechanics, 3-D geometry, and stochasticity on hanging drop tumor spheroids (HDS) and ductal carcinoma in situ (DCIS) of the breast. We demonstrate contact- and chemokine-based interactions among multiple cell types with examples in synthetic multicellular bioengineering, cancer heterogeneity, and cancer immunology. We developed PhysiCell to help the scientific community tackle multicellular systems biology problems involving many interacting cells in multi-substrate microenvironments. PhysiCell is also an independent, cross-platform codebase for replicating results from other simulators.

Colorectal cancer (CRC) progression is a complex process that is not well understood. We describe an in vitro organ-on-chip model that emulates in vivo tissue structure and the tumor microenvironment (TME) to better understand intravasation, an early step in metastasis. The CRC-on-chip incorporates fluid flow and peristalsis-like cyclic stretching and consists of endothelial and epithelial compartments, separated by a porous membrane. On-chip imaging and effluent analyses are used to interrogate CRC progression and the resulting cellular heterogeneity. Mass spectrometry-based metabolite profiles are indicative of a CRC disease state. Tumor cells intravasate from the epithelial channel to the endothelial channel, revealing differences in invasion between aggressive and non-aggressive tumor cells. Tuning the TME by peristalsis-like mechanical forces, the epithelial:endothelial interface, and the addition of fibroblasts influences the invasive capabilities of tumor cells. The CRC-on-chip is a tunable human-relevant model system and a valuable tool to study early invasive events in cancer.

Targeted agents have improved the efficacy of chemotherapy for cancer patients, however, there remains a lack of understanding of how these therapies affect the unsuspecting bystanders of the stromal microenvironment. Cetuximab, a monoclonal antibody therapy targeting the epidermal growth factor receptor (EGFR), is given in combination with chemotherapy as the standard of care for a subset of metastatic colorectal cancer patients. The overall response to this treatment is underwhelming and, while genetic mutations that confer resistance have been identified, it is still not known why this drug is ineffective for some patients. We discovered that cancer-associated fibroblasts (CAFs), a major cellular subset of the tumor stroma, can provide a source of cancer cell resistance. Specifically, we observed that upon treatment with cetuximab, CAFs increased their secretion of EGF, which was sufficient to render neighboring cancer cells resistant to cetuximab treatment through sustained mitogen-activated protein kinases (MAPK) signaling. Furthermore, we show the cetuximab-induced EGF secretion to be specific to CAFs and not to cancer cells or normal fibroblasts. Altogether, this work emphasizes the importance of the tumor microenvironment and considering the potential unintended consequences of therapeutically targeting cancer-driving proteins on non-tumorigenic cell types.

3D cell culture models have been developed to better mimic the physiological environments that exist in human diseases. As such, these models are advantageous over traditional 2D cultures for screening drug compounds. However, the practicalities of transitioning from 2D to 3D drug treatment studies pose challenges with respect to analysis methods. Patient-derived tumor organoids (PDTOs) possess unique features given their heterogeneity in size, shape, and growth patterns. A detailed assessment of the length scale at which PDTOs should be evaluated (i.e., individual cell or organoid-level analysis) has not been done to our knowledge. Therefore, using dynamic confocal live cell imaging and data analysis methods we examined tumor cell growth rates and drug response behaviors in colorectal cancer (CRC) PDTOs. High-resolution imaging of H2B-GFP-labeled organoids with DRAQ7 vital dye permitted tracking of cellular changes, such as cell birth and death events, in individual organoids. From these same images, we measured morphological features of the 3D objects, including volume, sphericity, and ellipticity. Sphericity and ellipticity were used to evaluate intra- and interpatient tumor organoid heterogeneity. We found a strong correlation between organoid live cell number and volume. Linear growth rate calculations based on volume or live cell counts were used to determine differential responses to therapeutic interventions. We showed that this approach can detect different types of drug effects (cytotoxic vs cytostatic) in PDTO cultures. Overall, our imaging-based quantification workflow results in multiple parameters that can provide patient- and drug-specific information for screening applications.

Tumor progression results from a complex interplay between cellular heterogeneity, treatment response, microenvironment and heterocellular interactions. Existing approaches to characterize this interplay suffer from an inability to distinguish between multiple cell types, often lack environmental context and are unable to perform multiplex phenotypic profiling of cell populations. Here we present a high-throughput platform for characterizing, with single-cell resolution, the dynamic phenotypic responses (i.e. morphology changes, proliferation, apoptosis) of heterogeneous cell populations both during standard growth and in response to multiple, co-occurring selective pressures. The speed of this platform enables a thorough investigation of the impacts of diverse selective pressures including genetic alterations, therapeutic interventions, heterocellular components and microenvironmental factors. The platform has been applied to both 2D and 3D culture systems and readily distinguishes between (1) cytotoxic versus cytostatic cellular responses; and (2) changes in morphological features over time and in response to perturbation. These important features can directly influence tumor evolution and clinical outcome. Our image-based approach provides a deeper insight into the cellular dynamics and heterogeneity of tumors (or other complex systems), with reduced reagents and time, offering advantages over traditional biological assays.

Three-quarters of compounds that enter clinical trials fail to make it to market due to safety or efficacy concerns. This statistic strongly suggests a need for better screening methods that result in improved translatability of compounds during the preclinical testing period. Patient-derived organoids have been touted as a promising 3D preclinical model system to impact the drug discovery pipeline, particularly in oncology. However, assessing drug efficacy in such models poses its own set of challenges, and traditional cell viability readouts fail to leverage some of the advantages that the organoid systems provide. Consequently, phenotypically evaluating complex 3D cell culture models remains difficult due to intra- and inter-patient organoid size differences, cellular heterogeneities, and temporal response dynamics. Here, we present an image-based high-content assay that provides object level information on 3D patient-derived tumor organoids without the need for vital dyes. Leveraging computer vision, we segment and define organoids as independent regions of interest and obtain morphometric and textural information per organoid. By acquiring brightfield images at different timepoints in a robust, non-destructive manner, we can track the dynamic response of individual organoids to various drugs. Furthermore, to simplify the analysis of the resulting large, complex data files, we developed a web-based data visualization tool, the Organoizer, that is available for public use. Our work demonstrates the feasibility and utility of using imaging, computer vision and machine learning to determine the vital status of individual patient-derived organoids without relying upon vital dyes, thus taking advantage of the characteristics offered by this preclinical model system.