Studying the capability of different cancer hallmarks to initiate tumor growth using a cellular automaton simulation. Application in a cancer stem cell context

tugHall: a simulator of cancer-cell evolution based on the hallmarks of cancer and tumor-related genes

Summary

The flood of recent cancer genomic data requires a coherent model that can sort out the findings to systematically explain clonal evolution and the resultant intra-tumor heterogeneity (ITH). Here, we present a new mathematical model designed to computationally simulate the evolution of cancer cells. The model connects the well-known hallmarks of cancer with the specific mutational states of tumor-related genes. The cell behavior phenotypes are stochastically determined, and the hallmarks probabilistically interfere with the phenotypic probabilities. In turn, the hallmark variables depend on the mutational states of tumor-related genes. Thus, our software can deepen our understanding of cancer-cell evolution and generation of ITH.

Availability and implementation

The open-source code is available in the repository https://github.com/nagornovys/Cancer_cell_evolution.

Contact

mamkato@ncc.go.jp

Supplementary information

Supplementary data are available at Bioinformatics online.

A Lock Free Approach To Parallelize The Cellular Potts Model: Application To Ductal Carcinoma In Situ

In the field of computational biology, in order to simulate multiscale biological systems, the Cellular Potts Model (CPM) has been used, which determines the actions that simulated cells can perform by determining a hamiltonian of energy that takes into account the influence that neighboring cells exert, under a wide range of parameters. There are some proposals in the literature that parallelize the CPM; in all cases, either lock-based techniques or other techniques that require large amounts of information to be disseminated among parallel tasks are used to preserve data coherence. In both cases, computational performance is limited. This work proposes an alternative approach for the parallelization of the model that uses transactional memory to maintain the coherence of the information. A Java implementation has been applied to the simulation of the ductal adenocarcinoma of breast in situ (DCIS). Times and speedups of the simulated execution of the model on the cluster of our university are analyzed. The results show a good speedup.

Cellular Automata Modeling of Stem-Cell-Driven Development of Tissue in the Nervous System

Mathematical and computational modeling enables biologists to integrate data from observations and experiments into a theoretical framework. In this review, we describe how developmental processes associated with stem-cell-driven growth of tissue in both the embryonic and adult nervous system can be modeled using cellular automata (CA). A cellular automaton is defined by its discrete nature in time, space, and state. The discrete space is represented by a uniform grid or lattice containing agents that interact with other agents within their local neighborhood. This possibility of local interactions of agents makes the cellular automata approach particularly well suited for studying through modeling how complex patterns at the tissue level emerge from fundamental developmental processes (such as proliferation, migration, differentiation, and death) at the single-cell level. As part of this review, we provide a primer for how to define biologically inspired rules governing these processes so that they can be implemented into a CA model. We then demonstrate the power of the CA approach by presenting simulations (in the form of figures and movies) based on building models of three developmental systems: the formation of the enteric nervous system through invasion by neural crest cells; the growth of normal and tumorous neurospheres induced by proliferation of adult neural stem/progenitor cells; and the neural fate specification through lateral inhibition of embryonic stem cells in the neurogenic region of Drosophila.

m-CALP - Yet another way of generating handwritten data through evolution for pattern recognition

Recently a cellular automata learning and prediction (CALP) model was proposed that used evolution to generate handwritten data for pattern recognition. However, the proposed method by CALP was mechanical and not intuitive, as it did not consider the degree of image deformation that takes place during evolution. The problem with this approach is that it evolved handwritten shapes to a completely different form. We propose a new data generation through evolution method called m-CALP, that evolves handwritten shapes using an objective function that also considers the degree of image deformation as it evolves. We replicated the exact same experimental setup of CALP - with the same 5 handwritten data sets, and classifiers with the exact same settings, and then used it to record the performance of our evolved data. We show that data evolved using our method maintains the interesting properties of a pattern. CALP was shown to have performed better than the state-of-the-art synthetic over sampling methods such as SMOTE and BORDERLINE-SMOTE, and m-CALP performs better than CALP in most cases, in the same environment. Data generation using evolution is a promising field, as it allows researchers to develop trainers and classifier for even small-sized training data. Therefore, we also show that evolving small-sized training data to generate more data using m-CALP performs better than both CALP and larger sized training data.

Cell adhesion heterogeneity reinforces tumour cell dissemination: novel insights from a mathematical model

BACKGROUND

Cancer cell invasion, dissemination, and metastasis have been linked to an epithelial-mesenchymal transition (EMT) of individual tumour cells. During EMT, adhesion molecules like E-cadherin are downregulated and the decrease of cell-cell adhesion allows tumour cells to dissociate from the primary tumour mass. This complex process depends on intracellular cues that are subject to genetic and epigenetic variability, as well as extrinsic cues from the local environment resulting in a spatial heterogeneity in the adhesive phenotype of individual tumour cells. Here, we use a novel mathematical model to study how adhesion heterogeneity, influenced by intrinsic and extrinsic factors, affects the dissemination of tumour cells from an epithelial cell population. The model is a multiscale cellular automaton that couples intracellular adhesion receptor regulation with cell-cell adhesion.

RESULTS

Simulations of our mathematical model indicate profound effects of adhesion heterogeneity on tumour cell dissemination. In particular, we show that a large variation of intracellular adhesion receptor concentrations in a cell population reinforces cell dissemination, regardless of extrinsic cues mediated through the local cell density. However, additional control of adhesion receptor concentration through the local cell density, which can be assumed in healthy cells, weakens the effect. Furthermore, we provide evidence that adhesion heterogeneity can explain the remarkable differences in adhesion receptor concentrations of epithelial and mesenchymal phenotypes observed during EMT and might drive early dissemination of tumour cells.

CONCLUSIONS

Our results suggest that adhesion heterogeneity may be a universal trigger to reinforce cell dissemination in epithelial cell populations. This effect can be at least partially compensated by a control of adhesion receptor regulation through neighbouring cells. Accordingly, our findings explain how both an increase in intra-tumour adhesion heterogeneity and the loss of control through the local environment can promote tumour cell dissemination.

REVIEWERS

This article was reviewed by Hanspeter Herzel, Thomas Dandekar and Marek Kimmel.

A Six-Stage Workflow for Robust Application of Systems Pharmacology

Quantitative and systems pharmacology (QSP) is increasingly being applied in pharmaceutical research and development. One factor critical to the ultimate success of QSP is the establishment of commonly accepted language, technical criteria, and workflows. We propose an integrated workflow that bridges conceptual objectives with underlying technical detail to support the execution, communication, and evaluation of QSP projects.

In silico study of the impact of cancer stem cell dynamics and radiobiological hypoxia on tumour response to hyperfractionated radiotherapy

OBJECTIVES

Advanced head and neck carcinomas (HNCs) are aggressive tumours, mainly due to hypoxia and a cancer stem cell (CSC) subpopulation. The aim of this study was to simulate tumour growth and behaviour during radiotherapy of three HNC groups (governed by different growth kinetics, hypoxia levels and CSC division pattern) to determine correlation between resistance factors and responses to hyperfractionated radiotherapy.

METHODS

An in silico HNC model was developed based on biologically realistic input parameters. During radiotherapy simulation, three parameters were studied: growth kinetics, hypoxia and probability of CSC symmetrical division. Both independent and combined effects on tumour response to hyperfractionated radiotherapy were assessed.

RESULTS

Oxic and very mildly hypoxic HNCs were revealed to be controlled by hyperfractionated radiotherapy, irrespective of growth kinetics and CSC division pattern. Moderately hypoxic tumours had different responses to radiotherapy: while slowly proliferating HNCs were still controllable, tumours with higher cell turnover were more resistant. In rapidly proliferating tumours, the number of fractions needed for tumour control increased exponentially with the probability of CSC symmetrical division, whereas in moderately growing HNC, this behaviour was linear. Severely hypoxic tumours could not be controlled by radiotherapy alone. Tumours with CSCs in a severely hypoxic niche required adjuvant therapies to be eradicated.

CONCLUSIONS

Growth kinetics strongly influence tumour responses to treatment. Slowly growing tumours showed linear dependence between dose and hypoxia/CSC, whereas rapidly growing tumours followed exponential behaviour.

Analysis of behaviour transitions in tumour growth using a cellular automaton simulation

The authors used computational biology as an approach for analysing the emergent dynamics of tumour growth at cellular level. They applied cellular automata for modelling the behaviour of cells when the main cancer cell hallmarks are present. Their model is oriented to mimic the development of multicellular spheroids of tumour cells. In their modelling, cells have a genome associated with the different cancer hallmarks, indicating if those are acquired as a consequence of mutations. The presence of the cancer hallmarks defines cell states and cell mitotic behaviours. These hallmarks are associated with a series of parameters, and depending on their values and the activation of the hallmarks in each of the cells, the system can evolve to different dynamics. With the simulation tool the authors performed an analysis of the first phases of cancer growth, using different and alternative strategies: firstly, studying the evolution of cancer cells and hallmarks in different representative situations regarding initial conditions and parameters, analysing the relative importance of the hallmarks for tumour progression; secondly, being the focus of this work, inspecting the behaviour transitions when the cancer cells are killed with a given probability during the cellular system progression.

Treatment Analysis in a Cancer Stem Cell Context Using a Tumor Growth Model Based on Cellular Automata

Cancer can be viewed as an emergent behavior in terms of complex system theory and artificial life, Cellular Automata (CA) being the tool most used for studying and characterizing the emergent behavior. Different approaches with CA models were used to model cancer growth. The use of the abstract model of acquired cancer hallmarks permits the direct modeling at cellular level, where a cellular automaton defines the mitotic and apoptotic behavior of cells, and allows for an analysis of different dynamics of the cellular system depending on the presence of the different hallmarks. A CA model based on the presence of hallmarks in the cells, which includes a simulation of the behavior of Cancer Stem Cells (CSC) and their implications for the resultant growth behavior of the multicellular system, was employed. This modeling of cancer growth, in the avascular phase, was employed to analyze the effect of cancer treatments in a cancer stem cell context. The model clearly explains why, after treatment against non-stem cancer cells, the regrowth capability of CSCs generates a faster regrowth of tumor behavior, and also shows that a continuous low-intensity treatment does not favor CSC proliferation and differentiation, thereby allowing an unproblematic control of future tumor regrowth. The analysis performed indicates that, contrary to the current attempts at CSC control, trying to make CSC proliferation more difficult is an important point to consider, especially in the immediate period after a standard treatment for controlling non-stem cancer cell proliferation.

What Tumor Dynamics Modeling Can Teach us About Exploiting the Stem-Cell View for Better Cancer Treatment

The cancer stem cell hypothesis is that in human solid cancers, only a small proportion of the cells, the cancer stem cells (CSCs), are self-renewing; the vast majority of the cancer cells are unable to sustain tumor growth indefinitely on their own. In recent years, discoveries have led to the concentration, if not isolation, of putative CSCs. The evidence has mounted that CSCs do exist and are important. This knowledge may promote better understanding of treatment resistance, create opportunities to test agents against CSCs, and open up promise for a fresh approach to cancer treatment. The first clinical trials of new anti-CSC agents are completed, and many others follow. Excitement is mounting that this knowledge will lead to major improvements, even breakthroughs, in treating cancer. However, exploitation of this phenomenon may be more successful if informed by insights into the population dynamics of tumor development. We revive some ideas in tumor dynamics modeling to extract some guidance in designing anti-CSC treatment regimens and the clinical trials that test them.