Noise-assisted interactions of tumor and immune cells

Curating models from BioModels: Developing a workflow for creating OMEX files

The reproducibility of computational biology models can be greatly facilitated by widely adopted standards and public repositories. We examined 50 models from the BioModels Database and attempted to validate the original curation and correct some of them if necessary. For each model, we reproduced these published results using Tellurium. Once reproduced we manually created a new set of files, with the model information stored by the Systems Biology Markup Language (SBML), and simulation instructions stored by the Simulation Experiment Description Markup Language (SED-ML), and everything included in an Open Modeling EXchange (OMEX) file, which could be used with a variety of simulators to reproduce the same results. On the one hand, the reproducibility procedure of 50 models developed a manual workflow that we would use to build an automatic platform to help users more easily curate and verify models in the future. On the other hand, these exercises allowed us to find the limitations and possible enhancement of the current curation and tooling to verify and curate models.

Curating models from BioModels: Developing a workflow for creating OMEX files

The reproducibility of computational biology models can be greatly facilitated by widely adopted standards and public repositories. We examined 50 models from the BioModels Database and attempted to validate the original curation and correct some of them if necessary. For each model, we reproduced these published results using Tellurium. Once reproduced we manually created a new set of files, with the model information stored by the Systems Biology Markup Language (SBML), and simulation instructions stored by the Simulation Experiment Description Markup Language (SED-ML), and everything included in an Open Modeling EXchange (OMEX) file, which could be used with a variety of simulators to reproduce the same results. On the one hand, the reproducibility procedure of 50 models developed a manual workflow that we would use to build an automatic platform to help users more easily curate and verify models in the future. On the other hand, these exercises allowed us to find the limitations and possible enhancement of the current curation and tooling to verify and curate models.

Survival induced by short-scale autocorrelations in birth-death-diffusion models with negative growth

In a birth, death, and diffusion process, the extinction-survival transition occurs when the average net growth rate is zero. For instance, in the presence of normally distributed time-varying stochastic growth rates with no autocorrelation, the transition indeed occurs at zero net growth rates. In contrast, if the growth rates are constant in time, a large enough variance in the growth rate will systematically ensure the survival of the global population even in a small system and, more importantly, even with a negative net growth rate. We here show that, surprisingly, for any intermediate temporal autocorrelation, any length of correlation, and any negative average growth rate, the same result holds. We test this argument on exponential and power law autocorrelation models and propose a simple condition for the growth rate variance at the transition.

Persistent Correlation in Cellular Noise Determines Longevity of Viral Infections

The slowly decaying viral dynamics, even after 2-3 weeks from diagnosis, is one of the characteristics of COVID-19 infection that is still unexplored in theoretical and experimental studies. This long-lived characteristic of viral infections in the framework of inherent variations or noise present at the cellular level is often overlooked. Therefore, in this work, we aim to understand the effect of these variations by proposing a stochastic non-Markovian model that not only captures the coupled dynamics between the immune cells and the virus but also enables the study of the effect of fluctuations. Numerical simulations of our model reveal that the long-range temporal correlations in fluctuations dictate the long-lived dynamics of a viral infection and, in turn, also affect the rates of immune response. Furthermore, predictions of our model system are in agreement with the experimental viral load data of COVID-19 patients from various countries.

A noble extended stochastic logistic model for cell proliferation with density-dependent parameters

Cell proliferation often experiences a density-dependent intrinsic proliferation rate (IPR) and negative feedback from growth-inhibiting molecules in culture media. The lack of flexible models with explanatory parameters fails to capture such a proliferation mechanism. We propose an extended logistic growth law with the density-dependent IPR and additional negative feedback. The extended parameters of the proposed model can be interpreted as density-dependent cell-cell cooperation and negative feedback on cell proliferation. Moreover, we incorporate further density regulation for flexibility in the model through environmental resistance on cells. The proposed growth law has similarities with the strong Allee model and harvesting phenomenon. We also develop the stochastic analog of the deterministic model by representing possible heterogeneity in growth-inhibiting molecules and environmental perturbation of the culture setup as correlated multiplicative and additive noises. The model provides a conditional maximum sustainable stable cell density (MSSCD) and a new fitness measure for proliferative cells. The proposed model shows superiority to the logistic law after fitting to real cell culture datasets. We illustrate both conditional MSSCD and the new cell fitness for a range of parameters. The cell density distributions reveal the chance of overproliferation, underproliferation, or decay for different parameter sets under the deterministic and stochastic setups.

Mitigation of rare events in multistable systems driven by correlated noise

We consider rare transitions induced by colored noise excitation in multistable systems. We show that undesirable transitions can be mitigated by a simple time-delay feedback control if the control parameters are judiciously chosen. We devise a parsimonious method for selecting the optimal control parameters, without requiring any Monte Carlo simulations of the system. This method relies on a new nonlinear Fokker-Planck equation whose stationary response distribution is approximated by a rapidly convergent iterative algorithm. In addition, our framework allows us to accurately predict, and subsequently suppress, the modal drift and tail inflation in the controlled stationary distribution. We demonstrate the efficacy of our method on two examples, including an optical laser model perturbed by multiplicative colored noise.

A nonlinear mathematical model of cell-mediated immune response for tumor phenotypic heterogeneity

Human cancers display intra-tumor heterogeneity in many phenotypic features, such as expression of cell surface receptors, growth, and angiogenic, proliferative, and immunogenic factors, which represent obstacles to a successful immune response. In this paper, we propose a nonlinear mathematical model of cancer immunosurveillance that takes into account some of these features based on cell-mediated immune responses. The model describes phenomena that are seen in vivo, such as tumor dormancy, robustness, immunoselection over tumor heterogeneity (also called "cancer immunoediting") and strong sensitivity to initial conditions in the composition of tumor microenvironment. The results framework has as common element the tumor as an attractor for abnormal cells. Bifurcation analysis give us as tumor attractors fixed-points, limit cycles and chaotic attractors, the latter emerging from period-doubling cascade displaying Feigenbaum's universality. Finally, we simulated both elimination and escape tumor scenarios by means of a stochastic version of the model according to the Doob-Gillespie algorithm.

Mathematical modeling of tumor-immune cell interactions

The anti-tumor activity of the immune system is increasingly recognized as critical for the mounting of a prolonged and effective response to cancer growth and invasion, and for preventing recurrence following resection or treatment. As the knowledge of tumor-immune cell interactions has advanced, experimental investigation has been complemented by mathematical modeling with the goal to quantify and predict these interactions. This succinct review offers an overview of recent tumor-immune continuum modeling approaches, highlighting spatial models. The focus is on work published in the past decade, incorporating one or more immune cell types and evaluating immune cell effects on tumor progression. Due to their relevance to cancer, the following immune cells and their combinations are described: macrophages, Cytotoxic T Lymphocytes, Natural Killer cells, dendritic cells, T regulatory cells, and CD4+ T helper cells. Although important insight has been gained from a mathematical modeling perspective, the development of models incorporating patient-specific data remains an important goal yet to be realized for potential clinical benefit.

Downregulated Tim-3 expression is responsible for the incidence and development of colorectal cancer

The present study aimed to investigate the role of T cell immunoglobulin domain and mucin-3 (Tim-3) in its gene and protein forms in colorectal cancer (CRC), and to verify the significance of Tim-3 expression in patients with CRC. A prospective analysis of 258 patients with CRC and 246 normal controls was conducted between December 2012 and June 2015. Intestinal samples were collected, including of CRC tissues, paracancerous tissues and normal colon mucosa tissues. Peripheral venous blood samples were also collected. Polymerase chain reaction (PCR) amplification, reverse transcription-quantitative PCR (RT-qPCR) and western blot analysis was performed for the detection and evaluation of Tim-3 gene and protein in various tissues. PCR analysis indicated that the T and G alleles of -882C/T and 4259T/G are associated with a significantly increased risk of CRC. Following the confirmation of Tim-3 expression in CRC tissues, RT-qPCR detection and western blot analysis revealed clear downregulation of Tim-3 mRNA and protein expression in the blood and tissue samples obtained from patients with CRC, as compared with in the corresponding control samples. Similar trends of decreased Tim-3 mRNA levels and protein expression were observed in CRC tissues compared with in the paracancerous and the normal colon mucosa tissues. In addition, the mRNA and protein expression levels in the paracancerous tissues were lower than those in the normal colon mucosa tissues. Furthermore, significantly lower Tim-3 mRNA levels were observed in patients with a tumor size >5 cm, a poor differentiation degree, higher tumor-node-metastasis stage (stage III-IV), and lymph node and distant metastasis. Collectively, genetic changes to Tim-3, expressed as polymorphisms in Tim-3, and decreased mRNA/protein expression may be partially responsible for the incidence and development of CRC.

Role of the interpretation of stochastic calculus in systems with cross-correlated Gaussian white noises

We derive the Fokker-Planck equation for multivariable Langevin equations with cross-correlated Gaussian white noises for an arbitrary interpretation of the stochastic differential equation. We formulate the conditions when the solution of the Fokker-Planck equation does not depend on which stochastic calculus is adopted. Further, we derive an equivalent multivariable Ito stochastic differential equation for each possible interpretation of the multivariable Langevin equation. To demonstrate the usefulness and significance of these general results, we consider the motion of Brownian particles. We study in detail the stability conditions for harmonic oscillators with two white noises, one of which is additive, random forcing, and the other, which accounts for fluctuations of either the damping or the spring coefficient, is multiplicative. We analyze the role of cross correlation in terms of the different noise interpretations and confirm the theoretical predictions via numerical simulations. We stress the interest of our results for numerical simulations of stochastic differential equations with an arbitrary interpretation of the stochastic integrals.