A mathematical model for the immune-mediated theory of metastasis

A Role of Effector CD 8 + T Cells Against Circulating Tumor Cells Cloaked with Platelets: Insights from a Mathematical Model

Cancer metastasis accounts for a majority of cancer-related deaths worldwide. Metastasis occurs when the primary tumor sheds cells into the blood and lymphatic circulation, thereby becoming circulating tumor cells (CTCs) that transverse through the circulatory system, extravasate the circulation and establish a secondary distant tumor. Accumulating evidence suggests that circulating effector CD 8 + T cells are able to recognize and attack arrested or extravasating CTCs, but this important antitumoral effect remains largely undefined. Recent studies highlighted the supporting role of activated platelets in CTCs's extravasation from the bloodstream, contributing to metastatic progression. In this work, a simple mathematical model describes how the primary tumor, CTCs, activated platelets and effector CD 8 + T cells participate in metastasis. The stability analysis reveals that for early dissemination of CTCs, effector CD 8 + T cells can present or keep secondary metastatic tumor burden at low equilibrium state. In contrast, for late dissemination of CTCs, effector CD 8 + T cells are unlikely to inhibit secondary tumor growth. Moreover, global sensitivity analysis demonstrates that the rate of the primary tumor growth, intravascular CTC proliferation, as well as the CD 8 + T cell proliferation, strongly affects the number of the secondary tumor cells. Additionally, model simulations indicate that an increase in CTC proliferation greatly contributes to tumor metastasis. Our simulations further illustrate that the higher the number of activated platelets on CTCs, the higher the probability of secondary tumor establishment. Intriguingly, from a mathematical immunology perspective, our simulations indicate that if the rate of effector CD 8 + T cell proliferation is high, then the secondary tumor formation can be considerably delayed, providing a window for adjuvant tumor control strategies. Collectively, our results suggest that the earlier the effector CD 8 + T cell response is enhanced the higher is the probability of preventing or delaying secondary tumor metastases.

Metastasis Models: Thermodynamics and Complexity

The thermodynamic formalism of nonequilibrium systems together with the theory of complex systems and systems biology offer an appropriate theoretical framework to explain the complexity observed at the macroscopic level in physiological phenomena. In turn, they allow the establishment of an appropriate conceptual and operational framework to address the study of phenomena such as the emergence and evolution of cancer.This chapter is organized as follows: In Subheading 1, an integrated vision of these disciplines is offered for the characterization of the emergence and evolution of cancer, seen as a nonlinear dynamic system, temporally and spatially self-organized out of thermodynamic equilibrium. The development of the various mathematical models and different techniques and approaches used in the characterization of cancer metastasis is presented in Subheading 2. Subheading 3 is devoted to the time course of cancer metastasis, with particular emphasis on the epithelial-mesenchymal transition (EMT henceforth) as well as chronotherapeutic treatments. In Subheading 4, models of the spatial evolution of cancer metastasis are presented. Finally, in Subheading 5, some conclusions and remarks are presented.

Mathematical Modeling Support for Lung Cancer Therapy-A Short Review

The paper presents a review of models that can be used to describe dynamics of lung cancer growth and its response to treatment at both cell population and intracellular processes levels. To address the latter, models of signaling pathways associated with cellular responses to treatment are overviewed. First, treatment options for lung cancer are discussed, and main signaling pathways and regulatory networks are briefly reviewed. Then, approaches used to model specific therapies are discussed. Following that, models of intracellular processes that are crucial in responses to therapies are presented. The paper is concluded with a discussion of the applicability of the presented approaches in the context of lung cancer.

Intermetastatic and Intrametastatic Heterogeneity Shapes Adaptive Therapy Cycling Dynamics

Adaptive therapies that alternate between drug applications and drug-free vacations can exploit competition between sensitive and resistant cells to maximize the time to progression. However, optimal dosing schedules depend on the properties of metastases, which are often not directly measurable in clinical practice. Here, we proposed a framework for estimating features of metastases through tumor response dynamics during the first adaptive therapy treatment cycle. Longitudinal prostate-specific antigen (PSA) levels in 16 patients with metastatic castration-resistant prostate cancer undergoing adaptive androgen deprivation treatment were analyzed to investigate relationships between cycle dynamics and clinical variables such as Gleason score, the change in the number of metastases over a cycle, and the total number of cycles over the course of treatment. The first cycle of adaptive therapy, which consists of a response period (applying therapy until 50% PSA reduction), and a regrowth period (removing treatment until reaching initial PSA levels), delineated several features of the computational metastatic system: larger metastases had longer cycles; a higher proportion of drug-resistant cells slowed the cycles; and a faster cell turnover rate sped up drug response time and slowed regrowth time. The number of metastases did not affect cycle times, as response dynamics were dominated by the largest tumors rather than the aggregate. In addition, systems with higher intermetastasis heterogeneity responded better to continuous therapy and correlated with dynamics from patients with high or low Gleason scores. Conversely, systems with higher intrametastasis heterogeneity responded better to adaptive therapy and correlated with dynamics from patients with intermediate Gleason scores.

SIGNIFICANCE

Multiscale mathematical modeling combined with biomarker dynamics during adaptive therapy helps identify underlying features of metastatic cancer to inform treatment decisions.

Patient-Specific Mathematical Model of the Clear Cell Renal Cell Carcinoma Microenvironment

The interactions between cells and molecules in the tumor microenvironment can give insight into the initiation and progression of tumors and their optimal treatment options. In this paper, we developed an ordinary differential equation (ODE) mathematical model of the interaction network of key players in the clear cell renal cell carcinoma (ccRCC) microenvironment. We then performed a global gradient-based sensitivity analysis to investigate the effects of the most sensitive parameters of the model on the number of cancer cells. The results indicate that parameters related to IL-6 have high a impact on cancer cell growth, such that decreasing the level of IL-6 can remarkably slow the tumor's growth.

Modeling cancer immunoediting in tumor microenvironment with system characterization through the ising-model Hamiltonian

Background and objective

Cancer Immunoediting (CI) describes the cellular-level interaction between tumor cells and the Immune System (IS) that takes place in the Tumor Micro-Environment (TME). CI is a highly dynamic and complex process comprising three distinct phases (Elimination, Equilibrium and Escape) wherein the IS can both protect against cancer development as well as, over time, promote the appearance of tumors with reduced immunogenicity. Herein we present an agent-based model for the simulation of CI in the TME, with the objective of promoting the understanding of this process.

Methods

Our model includes agents for tumor cells and for elements of the IS. The actions of these agents are governed by probabilistic rules, and agent recruitment (including cancer growth) is modeled via logistic functions. The system is formalized as an analogue of the Ising model from statistical mechanics to facilitate its analysis. The model was implemented in the Netlogo modeling environment and simulations were performed to verify, illustrate and characterize its operation.

Results

A main result from our simulations is the generation of emergent behavior in silico that is very difficult to observe directly in vivo or even in vitro. Our model is capable of generating the three phases of CI; it requires only a couple of control parameters and is robust to these. We demonstrate how our simulated system can be characterized through the Ising-model energy function, or Hamiltonian, which captures the “energy” involved in the interaction between agents and presents it in clear and distinct patterns for the different phases of CI.

Conclusions

The presented model is very flexible and robust, captures well the behaviors of the target system and can be easily extended to incorporate more variables such as those pertaining to different anti-cancer therapies. System characterization via the Ising-model Hamiltonian is a novel and powerful tool for a better understanding of CI and the development of more effective treatments. Since data of CI at the cellular level is very hard to procure, our hope is that tools such as this may be adopted to shed light on CI and related developing theories.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04731-w.

Investigating key cell types and molecules dynamics in PyMT mice model of breast cancer through a mathematical model

The most common kind of cancer among women is breast cancer. Understanding the tumor microenvironment and the interactions between individual cells and cytokines assists us in arriving at more effective treatments. Here, we develop a data-driven mathematical model to investigate the dynamics of key cell types and cytokines involved in breast cancer development. We use time-course gene expression profiles of a mouse model to estimate the relative abundance of cells and cytokines. We then employ a least-squares optimization method to evaluate the model's parameters based on the mice data. The resulting dynamics of the cells and cytokines obtained from the optimal set of parameters exhibit a decent agreement between the data and predictions. We perform a sensitivity analysis to identify the crucial parameters of the model and then perform a local bifurcation on them. The results reveal a strong connection between adipocytes, IL6, and the cancer population, suggesting them as potential targets for therapies.

Implications of immune-mediated metastatic growth on metastatic dormancy, blow-up, early detection, and treatment

Metastatic seeding of distant organs can occur in the very early stages of primary tumor development. Once seeded, these micrometastases may enter a dormant phase that can last decades. Curiously, the surgical removal of the primary tumor can stimulate the accelerated growth of distant metastases, a phenomenon known as metastatic blow-up. Recent clinical evidence has shown that the immune response can have strong tumor promoting effects. In this work, we investigate if the pro-tumor effects of the immune response can have a significant contribution to metastatic dormancy and metastatic blow-up. We develop an ordinary differential equation model of the immune-mediated theory of metastasis. We include both anti- and pro-tumor immune effects, in addition to the experimentally observed phenomenon of tumor-induced immune cell phenotypic plasticity. Using geometric singular perturbation analysis, we derive a rather simple model that captures the main processes and, at the same time, can be fully analyzed. Literature-derived parameter estimates are obtained, and model robustness is demonstrated through a time dependent sensitivity analysis. We determine conditions under which the parameterized model can successfully explain both metastatic dormancy and blow-up. The results confirm the significant active role of the immune system in the metastatic process. Numerical simulations suggest a novel measure to predict the occurrence of future metastatic blow-up in addition to new potential avenues for treatment of clinically undetectable micrometastases.

The Antitumor Efficacy of β-Elemene by Changing Tumor Inflammatory Environment and Tumor Microenvironment

Inflammatory mediators and inflammatory cells in the inflammatory microenvironment promote the transformation of normal cells to cancer cells in the early stage of cancer, promote the growth and development of cancer cells, and induce tumor immune escape. The monomeric active ingredient β-elemene is extracted from the traditional Chinese medicine Curcuma wenyujin and has been proven to have good anti-inflammatory and antitumor activities in clinical applications for more than 20 years in China. Recent studies have found that this traditional Chinese medicine plays a vital role in macrophage infiltration and M2 polarization, as well as in regulating immune disorders, and it even regulates the transcription factors NF-κB and STAT3 to alter inflammation, tumorigenesis, and development. In addition, β-elemene regulates not only different inflammatory factors (such as TNF-α, IFN, TGF-β, and IL-6/10) but also oxidative stress in vivo and in vitro. The excellent anti-inflammatory and antitumor effects of β-elemene and its ability to alter the inflammatory microenvironment of tumors have been gradually elaborated. Although the study of monomeric active ingredients in traditional Chinese medicines is insufficient in terms of quality and quantity, the pharmacological effects of more active ingredients of traditional Chinese medicines will be revealed after β-elemene.