Cancer-induced immunosuppression can enable effectiveness of immunotherapy through bistability generation: A mathematical and computational examination

Order-of-Mutation Effects on Cancer Progression: Models for Myeloproliferative Neoplasm

In some patients with myeloproliferative neoplasms (MPN), two genetic mutations are often found: JAK2 V617F and one in the TET2 gene. Whether one mutation is present influences how the other subsequent mutation will affect the regulation of gene expression. In other words, when a patient carries both mutations, the order of when they first arose has been shown to influence disease progression and prognosis. We propose a nonlinear ordinary differential equation, the Moran process, and Markov chain models to explain the non-additive and non-commutative mutation effects on recent clinical observations of gene expression patterns, proportions of cells with different mutations, and ages at diagnosis of MPN. Combined, these observations are used to shape our modeling framework. Our key proposal is that bistability in gene expression provides a natural explanation for many observed order-of-mutation effects. We also propose potential experimental measurements that can be used to confirm or refute predictions of our models.

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.

Order-of-mutation effects on cancer progression: models for myeloproliferative neoplasm

In some patients with myeloproliferative neoplasms (MPN), two genetic mutations are often found, JAK2 V617F and one in the TET2 gene. Whether or not one mutation is present will influence how the other subsequent mutation affects the regulation of gene expression. When both mutations are present, the order of their occurrence has been shown to influence disease progression and prognosis. We propose a nonlinear ordinary differential equation (ODE), Moran process, and Markov chain models to explain the non-additive and non-commutative mutation effects on recent clinical observations of gene expression patterns, proportions of cells with different mutations, and ages at diagnosis of MPN. These observations consistently shape our modeling framework. Our key proposal is that bistability in gene expression provides a natural explanation for many observed order-of-mutation effects. We also propose potential experimental measurements that can be used to confirm or refute predictions of our models.

Order-of-mutation effects on cancer progression: models for myeloproliferative neoplasm

In some patients with myeloproliferative neoplasms (MPN), two genetic mutations are often found, JAK2 V617F and one in the TET2 gene. Whether or not one mutation is present will influence how the other subsequent mutation affects the regulation of gene expression. When both mutations are present, the order of their occurrence has been shown to influence disease progression and prognosis. We propose a nonlinear ordinary differential equation (ODE), Moran process, and Markov chain models to explain the non-additive and non-commutative mutation effects on recent clinical observations of gene expression patterns, proportions of cells with different mutations, and ages at diagnosis of MPN. These observations consistently shape our modeling framework. Our key proposal is that bistability in gene expression provides a natural explanation for many observed order-of-mutation effects. We also propose potential experimental measurements that can be used to confirm or refute predictions of our models.

Modelling the Tumour Microenvironment, but What Exactly Do We Mean by "Model"?

The Oxford English Dictionary includes 17 definitions for the word "model" as a noun and another 11 as a verb. Therefore, context is necessary to understand the meaning of the word model. For instance, "model railways" refer to replicas of railways and trains at a smaller scale and a "model student" refers to an exemplary individual. In some cases, a specific context, like cancer research, may not be sufficient to provide one specific meaning for model. Even if the context is narrowed, specifically, to research related to the tumour microenvironment, "model" can be understood in a wide variety of ways, from an animal model to a mathematical expression. This paper presents a review of different "models" of the tumour microenvironment, as grouped by different definitions of the word into four categories: model organisms, in vitro models, mathematical models and computational models. Then, the frequencies of different meanings of the word "model" related to the tumour microenvironment are measured from numbers of entries in the MEDLINE database of the United States National Library of Medicine at the National Institutes of Health. The frequencies of the main components of the microenvironment and the organ-related cancers modelled are also assessed quantitatively with specific keywords. Whilst animal models, particularly xenografts and mouse models, are the most commonly used "models", the number of these entries has been slowly decreasing. Mathematical models, as well as prognostic and risk models, follow in frequency, and these have been growing in use.

Effects of immunosuppressants on T-cell dynamics: Understanding from a generic coarse-grained immune network model

Long-term immunosuppressive therapy is a drug regimen often used to lower aggressive immune responses in various chronic inflammatory diseases. However, such long-term therapy leading to immune suppression may trigger other adverse reactions in the immune system. The rising concern regarding the optimal dose and duration of such treatment has motivated us to understand non-classical immunomodulatory responses induced by various immunosuppressive steroid and secosteroid drugs such as glucocorticoid and vitamin D supplements. The immunomodulatory actions of such immunosuppressants (that govern the adaptive immune response) are often mediated through their characteristic control over CD4+ T-cells involving pro- and antiinflammatory T-cells. Several early studies attempted to decode temporal and dose-dependent behaviors of such pro- and anti-inflammatory T-cells using the chemical dynamics approach. We first summarize these early works. Then, we develop a minimal coarse-grained kinetic network model to capture the commonality in their immunomodulatory functions. This generic model successfully reproduces the characteristic dynamical features, including the clinical latency period in long-term T-cell dynamics. The temporal behavior of T-cells is found to be sensitive to specific rate parameters and doses of immunosuppressants. The steady-state analysis reflects the transition from an early classified weakly regulated (autoimmune-prone) immune state to a strongly regulated state (immunocompromised state), separated by an intervening state of moderate/balanced regulation. An optimal dose and duration are essential in rescuing balanced immune regulation. This review elucidates how developing a simple generic coarse-grained immune network model may provide immense information that helps diagnose inefficacy in adaptive immune function before and after administering immunosuppressants such as glucocorticoid or vitamin D.

Bistability in a model of tumor-immune system interactions with an oncolytic viral therapy

A mathematical model of tumor-immune system interactions with an oncolytic virus therapy for which the immune system plays a twofold role against cancer cells is derived. The immune cells can kill cancer cells but can also eliminate viruses from the therapy. In addition, immune cells can either be stimulated to proliferate or be impaired to reduce their growth by tumor cells. It is shown that if the tumor killing rate by immune cells is above a critical value, the tumor can be eradicated for all sizes, where the critical killing rate depends on whether the immune system is immunosuppressive or proliferative. For a reduced tumor killing rate with an immunosuppressive immune system, that bistability exists in a large parameter space follows from our numerical bifurcation study. Depending on the tumor size, the tumor can either be eradicated or be reduced to a size less than its carrying capacity. However, reducing the viral killing rate by immune cells always increases the effectiveness of the viral therapy. This reduction may be achieved by manipulating certain genes of viruses via genetic engineering or by chemical modification of viral coat proteins to avoid detection by the immune cells.