skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on December 1, 2026

Title: Mathematical modeling of multicellular tumor spheroids quantifies inter-patient and intra-tumor heterogeneity
Abstract In the study of brain tumors, patient-derived three-dimensional sphere cultures provide an important tool for studying emerging treatments. The growth of such spheroids depends on the combined effects of proliferation and migration of cells, but it is challenging to make accurate distinctions between increase in cell number versus the radial movement of cells. To address this, we formulate a novel model in the form of a system of two partial differential equations (PDEs) incorporating both migration and growth terms, and show that it more accurately fits our data compared to simpler PDE models. We show that traveling-wave speeds are strongly associated with population heterogeneity. Having fitted the model to our dataset we show that a subset of the cell lines are best described by a “Go-or-Grow”-type model, which constitutes a special case of our model. Finally, we investigate whether our fitted model parameters are correlated with patient age and survival.  more » « less
Award ID(s):
2320244
PAR ID:
10616603
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Nature Partner Journals
Date Published:
Journal Name:
npj Systems Biology and Applications
Volume:
11
Issue:
1
ISSN:
2056-7189
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Colorectal cancer, a significant cause of cancer-related mortality, often exhibits drug resistance, highlighting the need for improved tumor models to advance personalized drug testing and precision therapy. We generated organoids from primary colorectal cancer cells cultured through the conditional reprogramming technique, establishing a framework to perform short-term drug testing studies on patient-derived cells. To model interactions with stromal cells in the tumor microenvironment, we combined cancer cell organoids with carcinoma-associated fibroblasts, a cell type implicated in disease progression and drug resistance. Our organotypic models revealed that carcinoma-associated fibroblasts promote cancer cell proliferation and stemness primarily through hepatocyte growth factor–MET paracrine signaling and activation of cyclin-dependent kinases. Disrupting these tumor–stromal interactions reduced organoid size while limiting oncogenic signals and cancer stemness. Leveraging this tumor model, we identified effective drug combinations targeting colorectal cancer cells and their tumorigenic activities. Our study highlights a path to incorporate patient-derived cells and tumor–stromal interactions into a drug testing workflow that could identify effective therapies for individual patients. 
    more » « less
  2. The CDK4/6 inhibitor palbociclib blocks cell cycle progression in Estrogen receptor–positive, human epidermal growth factor 2 receptor–negative (ER+/HER2−) breast tumor cells. Despite the drug’s success in improving patient outcomes, a small percentage of tumor cells continues to divide in the presence of palbociclib—a phenomenon we refer to as fractional resistance. It is critical to understand the cellular mechanisms underlying fractional resistance because the precise percentage of resistant cells in patient tissue is a strong predictor of clinical outcomes. Here, we hypothesize that fractional resistance arises from cell-to-cell differences in core cell cycle regulators that allow a subset of cells to escape CDK4/6 inhibitor therapy. We used multiplex, single-cell imaging to identify fractionally resistant cells in both cultured and primary breast tumor samples resected from patients. Resistant cells showed premature accumulation of multiple G1 regulators including E2F1, retinoblastoma protein, and CDK2, as well as enhanced sensitivity to pharmacological inhibition of CDK2 activity. Using trajectory inference approaches, we show how plasticity among cell cycle regulators gives rise to alternate cell cycle “paths” that allow individual tumor cells to escape palbociclib treatment. Understanding drivers of cell cycle plasticity, and how to eliminate resistant cell cycle paths, could lead to improved cancer therapies targeting fractionally resistant cells to improve patient outcomes. 
    more » « less
  3. Wodarz, Dominik (Ed.)
    Patient-derived tumor organoids (PDTOs) are novel cellular models that maintain the genetic, phenotypic and structural features of patient tumor tissue and are useful for studying tumorigenesis and drug response. When integrated with advanced 3D imaging and analysis techniques, PDTOs can be used to establish physiologically relevant high-throughput and high-content drug screening platforms that support the development of patient-specific treatment strategies. However, in order to effectively leverage high-throughput PDTO observations for clinical predictions, it is critical to establish a quantitative understanding of the basic properties and variability of organoid growth dynamics. In this work, we introduced an innovative workflow for analyzing and understanding PDTO growth dynamics, by integrating a high-throughput imaging deep learning platform with mathematical modeling, incorporating flexible growth laws and variable dormancy times. We applied the workflow to colon cancer organoids and demonstrated that organoid growth is well-described by the Gompertz model of growth. Our analysis showed significant intrapatient heterogeneity in PDTO growth dynamics, with the initial exponential growth rate of an organoid following a lognormal distribution within each dataset. The level of intrapatient heterogeneity varied between patients, as did organoid growth rates and dormancy times of single seeded cells. Our work contributes to an emerging understanding of the basic growth characteristics of PDTOs, and it highlights the heterogeneity in organoid growth both within and between patients. These results pave the way for further modeling efforts aimed at predicting treatment response dynamics and drug resistance timing. 
    more » « less
  4. Abstract Recent years have seen a tremendous growth of interest in understanding the role that the adaptive immune system could play in interdicting tumor progression. In this context, it has been shown that the density of adaptive immune cells inside a solid tumor serves as a favorable prognostic marker across different types of cancer. The exact mechanisms underlying the degree of immune cell infiltration is largely unknown. Here, we quantify the temporal dynamics of the density profile of activated immune cells around a solid tumor spheroid. We propose a computational model incorporating immune cells with active, persistent movement and a proliferation rate that depends on the presence of cancer cells, and show that the model able to reproduce semi-quantitatively the experimentally measured infiltration profile. Studying the density distribution of immune cells inside a solid tumor can help us better understand immune trafficking in the tumor micro-environment, hopefully leading towards novel immunotherapeutic strategies. 
    more » « less
  5. Abstract Collagen fibers in the 3D tumor microenvironment (TME) exhibit complex alignment landscapes that are critical in directing cell migration through a process called contact guidance. Previous in vitro work studying this phenomenon has focused on quantifying cell responses in uniformly aligned environments. However, the TME also features short‐range gradients in fiber alignment that result from cell‐induced traction forces. Although the influence of graded biophysical taxis cues is well established, cell responses to physiological alignment gradients remain largely unexplored. In this work, fiber alignment gradients in biopsy samples are characterized and recreated using a new microfluidic biofabrication technique to achieve tunable sub‐millimeter to millimeter scale gradients. This study represents the first successful engineering of continuous alignment gradients in soft, natural biomaterials. Migration experiments on graded alignment show that human umbilical vein endothelial cells (HUVECs) exhibit increased directionality, persistence, and speed compared to uniform and unaligned fiber architectures. Similarly, patterned MDA‐MB‐231 breast cancer cell aggregates exhibit biased migration toward increasing fiber alignment, suggesting a role for alignment gradients as a taxis cue. This user‐friendly approach, requiring no specialized equipment, is anticipated to offer new insights into the biophysical cues that cells interpret as they traverse the extracellular matrix (ECM), with broad applicability in healthy and diseased tissue environments. 
    more » « less