Mathematical Modelling for the Role of CD4+T Cells in Tumor-Immune Interactions

A Mathematical Model of TCR-T Cell Therapy for Cervical Cancer

Engineered T cell receptor (TCR)-expressing T (TCR-T) cells are intended to drive strong anti-tumor responses upon recognition of the specific cancer antigen, resulting in rapid expansion in the number of TCR-T cells and enhanced cytotoxic functions, causing cancer cell death. However, although TCR-T cell therapy against cancers has shown promising results, it remains difficult to predict which patients will benefit from such therapy. We develop a mathematical model to identify mechanisms associated with an insufficient response in a mouse cancer model. We consider a dynamical system that follows the population of cancer cells, effector TCR-T cells, regulatory T cells (Tregs), and "non-cancer-killing" TCR-T cells. We demonstrate that the majority of TCR-T cells within the tumor are "non-cancer-killing" TCR-T cells, such as exhausted cells, which contribute little or no direct cytotoxicity in the tumor microenvironment (TME). We also establish two important factors influencing tumor regression: the reversal of the immunosuppressive TME following depletion of Tregs, and the increased number of effector TCR-T cells with antitumor activity. Using mathematical modeling, we show that certain parameters, such as increasing the cytotoxicity of effector TCR-T cells and modifying the number of TCR-T cells, play important roles in determining outcomes.

Mathematical model for IL-2-based cancer immunotherapy

A basic mathematical model for IL-2-based cancer immunotherapy is proposed and studied. Our analysis shows that the outcome of therapy is mainly determined by three parameters, the relative death rate of CD4+ T cells, the relative death rate of CD8+ T cells, and the dose of IL-2 treatment. Minimal equilibrium tumor size can be reached with a large dose of IL-2 in the case that CD4+ T cells die out. However, in cases where CD4+ and CD8+ T cells persist, the final tumor size is independent of the IL-2 dose and is given by the relative death rate of CD4+ T cells. Two groups of in silico clinical trials show some short-term behaviors of IL-2 treatment. IL-2 administration can slow the proliferation of CD4+ T cells, while high doses for a short period of time over several days transiently increase the population of CD8+ T cells during treatment before it recedes to its equilibrium. IL-2 administration for a short period of time over many days suppresses the tumor population for a longer time before approaching its steady-state levels. This implies that intermittent administration of IL-2 may be a good strategy for controlling tumor size.

The effect of chemotaxis on T-cell regulatory dynamics

Autoimmune diseases, such as Multiple Sclerosis, are often modelled through the dynamics of T-cell interactions. However, the spatial aspect of such diseases, and how dynamics may result in spatially heterogeneous outcomes, is often overlooked. We consider the effects of diffusion and chemotaxis on T-cell regulatory dynamics using a three-species model of effector and regulator T-cell populations, along with a chemical signalling agent. While diffusion alone cannot lead to instability and spatial patterning, the inclusion of chemotaxis can result in multiple forms of instability that produce highly complicated spatiotemporal behaviour. The parameter regimes in which different instabilities occur are determined through linear stability analysis and the full dynamics is explored through numerical simulation.

Aging, Inflammation, and Comorbidity in Cancers-A General In Silico Study Exemplified by Myeloproliferative Malignancies

(1) Background: We consider dormant, pre-cancerous states prevented from developing into cancer by the immune system. Inflammatory morbidity may compromise the immune system and cause the pre-cancer to escape into an actual cancerous development. The immune deficiency described is general, but the results may vary across specific cancers due to different variances (2) Methods: We formulate a general conceptual model to perform rigorous in silico consequence analysis. Relevant existing data for myeloproliferative malignancies from the literature are used to calibrate the in silico computations. (3) Results and conclusions: The hypothesis suggests a common physiological origin for many clinical and epidemiological observations in relation to cancers in general. Examples are the observed age-dependent prevalence for hematopoietic cancers, a general mechanism-based explanation for why the risk of cancer increases with age, and how somatic mutations in general, and specifically seen in screenings of citizens, sometimes are non-increased or even decrease when followed over time. The conceptual model is used to characterize different groups of citizens and patients, describing different treatment responses and development scenarios.

Analytical Models of Intra- and Extratumoral Cell Interactions at Avascular Stage of Growth in the Presence of Targeted Chemotherapy

In this study, we propose a set of nonlinear differential equations to model the dynamic growth of avascular stage tumors, considering nutrient supply from underlying tissue, innate immune response, contact inhibition of cell migration, and interactions with a chemotherapeutic agent. The model has been validated against available experimental data from the literature for tumor growth. We assume that the size of the modeled tumor is already detectable, and it represents all clinically observed existent cell populations; initial conditions are selected accordingly. Numerical results indicate that the tumor size and regression significantly depend on the strength of the host immune system. The effect of chemotherapy is investigated, not only within the malignancy, but also in terms of the responding immune cells and healthy tissue in the vicinity of a tumor.

Evaluation of innate and adaptive immune system interactions in the tumor microenvironment via a 3D continuum model

Immune cells in the tumor microenvironment (TME) are known to affect tumor growth, vascularization, and extracellular matrix (ECM) deposition. Marked interest in system-scale analysis of immune species interactions within the TME has encouraged progress in modeling tumor-immune interactions in silico. Due to the computational cost of simulating these intricate interactions, models have typically been constrained to representing a limited number of immune species. To expand the capability for system-scale analysis, this study develops a three-dimensional continuum mixture model of tumor-immune interactions to simulate multiple immune species in the TME. Building upon a recent distributed computing implementation that enables efficient solution of such mixture models, major immune species including monocytes, macrophages, natural killer cells, dendritic cells, neutrophils, myeloid-derived suppressor cells (MDSC), cytotoxic, helper, regulatory T-cells, and effector and regulatory B-cells and their interactions are represented in this novel implementation. Immune species extravasate from blood vasculature, undergo chemotaxis toward regions of high chemokine concentration, and influence the TME in proportion to locally defined levels of stimulation. The immune species contribute to the production of angiogenic and tumor growth factors, promotion of myofibroblast deposition of ECM, upregulation of angiogenesis, and elimination of living and dead tumor species. The results show that this modeling approach offers the capability for quantitative insight into the modulation of tumor growth by diverse immune-tumor interactions and immune-driven TME effects. In particular, MDSC-mediated effects on tumor-associated immune species' activation levels, volume fraction, and influence on the TME are explored. Longer term, linking of the model parameters to particular patient tumor information could simulate cancer-specific immune responses and move toward a more comprehensive evaluation of immunotherapeutic strategies.

Can the Kuznetsov Model Replicate and Predict Cancer Growth in Humans?

Several mathematical models to predict tumor growth over time have been developed in the last decades. A central aspect of such models is the interaction of tumor cells with immune effector cells. The Kuznetsov model (Kuznetsov et al. in Bull Math Biol 56(2):295–321, 1994) is the most prominent of these models and has been used as a basis for many other related models and theoretical studies. However, none of these models have been validated with large-scale real-world data of human patients treated with cancer immunotherapy. In addition, parameter estimation of these models remains a major bottleneck on the way to model-based and data-driven medical treatment. In this study, we quantitatively fit Kuznetsov’s model to a large dataset of 1472 patients, of which 210 patients have more than six data points, by estimating the model parameters of each patient individually. We also conduct a global practical identifiability analysis for the estimated parameters. We thus demonstrate that several combinations of parameter values could lead to accurate data fitting. This opens the potential for global parameter estimation of the model, in which the values of all or some parameters are fixed for all patients. Furthermore, by omitting the last two or three data points, we show that the model can be extrapolated and predict future tumor dynamics. This paves the way for a more clinically relevant application of mathematical tumor modeling, in which the treatment strategy could be adjusted in advance according to the model’s future predictions.

A Mathematical Study of the Role of tBregs in Breast Cancer

A model for the mathematical study of immune response to breast cancer is proposed and studied, both analytically and numerically. It is a simplification of a complex one, recently introduced by two of the present authors. It serves for a compact study of the dynamical role in cancer promotion of a relatively recently described subgroup of regulatory B cells, which are evoked by the tumour.

Dynamics of an Antitumour Model with Pulsed Radioimmunotherapy

In this paper, an antitumour model for characterising radiotherapy and immunotherapy processes at different fixed times is proposed. The global attractiveness of the positive periodic solution for each corresponding subsystem is proved with the integral inequality technique. Then, based on the differentiability of the solutions with respect to the initial values, the eigenvalues of the Jacobian matrix at a fixed point corresponding to the tumour-free periodic solution are determined, resulting in a sufficient condition for local stability. The solutions to the ordinary differential equations are compared, the threshold condition for the global attractiveness of the tumour-free periodic solution is provided in terms of an indicator $R_0$ , and the permanence of a system with at least one tumour-present periodic solution is investigated. Furthermore, the effects of the death rate, effector cell injection dosage, therapeutic period, and effector cell activation rate on indicator $R_0$ are determined through numerical simulations, and the results indicate that radioimmunotherapy is more effective than either radiotherapy or immunotherapy alone.

Eosinophils and melanoma: Implications for immunotherapy

New therapies such as immune checkpoint blockers (ICB) have offered extended survival to patients affected by advanced melanoma. However, ICBs have demonstrated debilitating side effects on the joints, liver, lungs, skin, and gut. Several biomarkers have been identified for their role in predicting which patients better tolerate ICBs. Still, these biomarkers are limited by immunologic and genetic heterogeneity and the complexity of translation into clinical practice. Recent observational studies have suggested eosinophil counts, and serum levels of eosinophil cationic protein are significantly associated with prolonged survival in advanced-stage melanoma. It is likely that eosinophils thereby modulate treatment response through mechanisms yet to be explored. Here, we review the functionality of eosinophils, their oncogenic role in melanoma and discuss how these mechanisms may influence patient response to ICBs and their implications in clinical practice.

Mathematical modeling of tumor-immune system interactions: the effect of rituximab on breast cancer immune response

tBregs are a newly discovered subcategory of B regulatory cells, which are generated by breast cancer, resulting in the increase of Tregs and therefore in the death of NK cells. In this study, we use a mathematical and computational approach to investigate the complex interactions between the aforementioned cells as well as CD8⁺ T cells, CD4⁺ T cells and B cells. Furthermore, we use data fitting to prove that the functional response regarding the lysis of breast cancer cells by NK cells has a ratio-dependent form. Additionally, we include in our model the concentration of rituximab - a monoclonal antibody that has been suggested as a potential breast cancer therapy - and test its effect, when the standard, as well as experimental dosages, are administered.

Simulation Model for Hashimoto Autoimmune Thyroiditis Disease

Hashimoto thyroiditis (HT) is a pathology that often causes a gradual thyroid insufficiency in affected patients due to the autoimmune destruction of this gland. The cellular immune response mediated by T helper lymphocytes TH1 and TH17 can induce the HT disease. In this pathologic condition, there is an imbalance between the TH17 and Treg lymphocytes as well as a gut microbiota dysfunction. The objective of this work was to describe the interactions of the cell subpopulations that participate in HT. To achieve this goal, we generated a mathematical model that allowed the simulation of different scenarios for the dynamic interaction between thyroid cells, the immune system, and the gut microbiota. We used a hypothetical-deductive design of mathematical modeling based on a system of ordinary differential equations, where the state variables are the TH1, TH17, and Treg lymphocytes, the thyrocytes, and the bacteria from gut microbiota. This work generated a compartmental model of the cellular immune response occurring in the thyroid gland. It was observed that TH1 and TH17 lymphocytes could increase the immune cells' activity, as well as activate effector cells directly and trigger the apoptosis and inflammation processes of healthy thyrocytes indirectly. Likewise, the model showed that a reduction in Treg lymphocytes could increase the activity of TH17 lymphocytes when an imbalance of the gut microbiota composition occurred. The numerical results highlight the TH1, TH17, and bacterial balance of the gut microbiota activities as important factors for the development of HT disease.

Lyapunov function and global asymptotic stability for a new multiscale viral dynamics model incorporating the immune system response: Implemented upon HCV

In this paper, a new mathematical model is formulated that describes the interaction between uninfected cells, infected cells, viruses, intracellular viral RNA, Cytotoxic T-lymphocytes (CTLs), and antibodies. Hence, the model contains certain biological relations that are thought to be key factors driving this interaction which allow us to obtain precise logical conclusions. Therefore, it improves our perception, that would otherwise not be possible, to comprehend the pathogenesis, to interpret clinical data, to control treatment, and to suggest new relations. This model can be used to study viral dynamics in patients for a wide range of infectious diseases like HIV, HPV, HBV, HCV, and Covid-19. Though, analysis of a new multiscale HCV model incorporating the immune system response is considered in detail, the analysis and results can be applied for all other viruses. The model utilizes a transformed multiscale model in the form of ordinary differential equations (ODE) and incorporates into it the interaction of the immune system. The role of CTLs and the role of antibody responses are investigated. The positivity of the solutions is proven, the basic reproduction number is obtained, and the equilibrium points are specified. The stability at the equilibrium points is analyzed based on the Lyapunov invariance principle. By using appropriate Lyapunov functions, the uninfected equilibrium point is proven to be globally asymptotically stable when the reproduction number is less than one and unstable otherwise. Global stability of the infected equilibrium points is considered, and it has been found that each equilibrium point has a specific domain of stability. Stability regions could be overlapped and a bistable equilibria could be found, which means the coexistence of two stable equilibrium points. Hence, the solution converges to one of them depending on the initial conditions.

Potential of Immunotherapies in Treating Hematological Cancer-Infection Comorbidities-A Mathematical Modelling Approach

BACKGROUND

The immune system attacks threats like an emerging cancer or infections like COVID-19 but it also plays a role in dealing with autoimmune disease, e.g., inflammatory bowel diseases, and aging. Malignant cells may tend to be eradicated, to appraoch a dormant state or escape the immune system resulting in uncontrolled growth leading to cancer progression. If the immune system is busy fighting a cancer, a severe infection on top of it may compromise the immunoediting and the comorbidity may be too taxing for the immune system to control.

METHOD

A novel mechanism based computational model coupling a cancer-infection development to the adaptive immune system is presented and analyzed. The model maps the outcome to the underlying physiological mechanisms and agree with numerous evidence based medical observations.

RESULTS AND CONCLUSIONS

Progression of a cancer and the effect of treatments depend on the cancer size, the level of infection, and on the efficiency of the adaptive immune system. The model exhibits bi-stability, i.e., virtual patient trajectories gravitate towards one of two stable steady states: a dormant state or a full-blown cancer-infection disease state. An infectious threshold curve exists and if infection exceed this separatrix for sufficiently long time the cancer escapes. Thus, early treatment is vital for remission and severe infections may instigate cancer progression. CAR T-cell Immunotherapy may sufficiently control cancer progression back into a dormant state but the therapy significantly gains efficiency in combination with antibiotics or immunomodulation.

Dynamical Analysis of a Multiscale Model of Hepatitis C Virus Infection Using a Transformed ODEs Model

A mathematically identical ordinary differential equations (ODEs) model was derived from a multiscale partial differential equations (PDEs) model of hepatitis c virus infection, which helps to overcome the limitations of the PDE model in clinical data analysis. We have discussed about basic properties of the system and found the basic reproduction number of the system. A condition for the local stability of the uninfected and the infected steady states is presented. The local stability analysis of the model shows that the system is asymptotically stable at the disease-free equilibrium point when the basic reproduction number is less than one. When the basic reproduction number is greater than one endemic equilibrium point exists, and the local stability analysis proves that this point is asymptotically stable. Numerical sensitivity analysis based on model parameters is performed and therefore the result describes the influence of each parameter on the basic reproduction number.

Using Padé Approximant Method to Solve the Mathematical Model of Tumor-Immune Interactions

A mathematical model, in the form of a system of nonlinear ordinary differential equations, that describes the interaction between tumor cells and effective immune cells is proposed. An exact solution cannot be found to this system like many other nonlinear systems. Yet, approximate analytical solution is explored. This solution should have a large interval of convergence to be acceptable because the interaction can take many days to reach its steady state. Power series method is used to obtain a series solution. In this process, some auxiliary variables are used to transform the system of equations to polynomial form. However, this solution has a small radius of convergence, therefore, Padé approximant method is used to extend the domain of convergence. Hence, the obtained approximate analytical solution is valid over a large interval and has a remarkable accuracy when compared with numerical solution.