A numerical model of EMT6/Ro spheroid dynamics under irradiation: calibration and estimation of the underlying irradiation-induced cell survival probability.

We present extensions to our quasi-2D cellular automata spheroid model that add a cellular kinetics module together with an irradiation and repair module. Significantly, our approach is not based on the Linear Quadratic (LQ) model, instead, we propose a simple two-parameter, algorithmic model which captures the essential biological features of irradiation-induced cell death, repair and associated cell cycle delays. This approach allows us to estimate directly the underlying irradiation-induced cell survival probability. We present the calibration of this extended model both with and without the application of single irradiation doses to the commonly studied (in vitro) EMT6/Ro (mammary carcinoma) cell line. A comparison of the estimated underlying cell survival probability with the in vitro survival probability data confirms the expected differences in the measures.

A quantitative cellular automaton model of in vitro multicellular spheroid tumour growth.

We report numerical results from a 2D cellular automaton (CA) model describing the dynamics of the in vitro cultivated multicellular spheroid obtained from EMT6/Ro (mammary carcinoma) cell line. Significantly, the CA model relaxes the often assumed one-to-one correspondence between cells and CA sites so as to correctly model the peripheral mitotic boundary region, and to enable the study of necrosis in large avascular tumours. By full calibration and scaling to available experimental data, the model produces with good accuracy experimentally comparable data on a range of bulk tumour kinetics and necrosis measures. Our main finding is that the metabolic production of H(+) ions is not sufficient to cause central necrosis prior to the sub-viable nutrient-deficient stage of tumour development being reached. Thus, the model suggests that an additional process is required to explain the experimentally observable onset of necrosis prior to the non-viable nutrient-deficient point being reached.

Preface.

This volume was inspired by the topics presented at the international conference "Micro and Macro Systems in Life Sciences" which was held on Jun 8-12, 2015 in Bedlewo, Poland. System biology is an approach which tries to understand how micro systems, at the molecular and cellular levels, affect macro systems such as organs, tissue and populations. Thus it is not surprising that a major theme of this volume evolves around cancer and its treatment. Articles on this topic include models for tumor induced angiogenesis, without and with delays, metastatic niche of the bone marrow, drug resistance and metronomic chemotherapy, and virotherapy of glioma. Methods range from dynamical systems to optimal control. Another well represented topic of this volume is mathematical modeling in epidemiology. Mathematical approaches to modeling and control of more specific diseases like malaria, Ebola or human papillomavirus are discussed as well as a more general approaches to the SEIR, and even more general class of models in epidemiology, by using the tools of optimal control and optimization. The volume also brings up challenges in mathematical modeling of other diseases such as tuberculosis. Partial differential equations combined with numerical approaches are becoming important tools in modeling not only tumor growth and treatment, but also other diseases, such as fibrosis of the liver, and atherosclerosis and its associated blood flow dynamics, and our volume presents a state of the art approach on these topics. Understanding mathematics behind the cell motion, appearance of the special patterns in various cell populations, and age structured mutations are among topics addressed inour volume. A spatio-temporal models of synthetic genetic oscillators brings the analysis to the gene level which is the focus of much of current biological research. Mathematics can help biologists to explain the collective behavior of bacterial, a topic that is also presented here. Finally some more across the discipline topics are being addresses, which can appear as a challenge in studying problems in systems biology on all, macro, meso and micro levels. They include numerical approaches to stochastic wave equation arising in modeling Brownian motion, discrete velocity models, many particle approximations as well as very important aspect on the connection between discrete measurement and the construction of the models for various phenomena, particularly the one involving delays. With the variety of biological topics and their mathematical approaches we very much hope that the reader of the Mathematical Biosciences and Engineering will find this volume interesting and inspirational for their own research.

A matter of timing: identifying significant multi-dose radiotherapy improvements by numerical simulation and genetic algorithm search.

Multi-dose radiotherapy protocols (fraction dose and timing) currently used in the clinic are the product of human selection based on habit, received wisdom, physician experience and intra-day patient timetabling. However, due to combinatorial considerations, the potential treatment protocol space for a given total dose or treatment length is enormous, even for relatively coarse search; well beyond the capacity of traditional in-vitro methods. In constrast, high fidelity numerical simulation of tumor development is well suited to the challenge. Building on our previous single-dose numerical simulation model of EMT6/Ro spheroids, a multi-dose irradiation response module is added and calibrated to the effective dose arising from 18 independent multi-dose treatment programs available in the experimental literature. With the developed model a constrained, non-linear, search for better performing cadidate protocols is conducted within the vicinity of two benchmarks by genetic algorithm (GA) techniques. After evaluating less than 0.01% of the potential benchmark protocol space, candidate protocols were identified by the GA which conferred an average of 9.4% (max benefit 16.5%) and 7.1% (13.3%) improvement (reduction) on tumour cell count compared to the two benchmarks, respectively. Noticing that a convergent phenomenon of the top performing protocols was their temporal synchronicity, a further series of numerical experiments was conducted with periodic time-gap protocols (10 h to 23 h), leading to the discovery that the performance of the GA search candidates could be replicated by 17-18 h periodic candidates. Further dynamic irradiation-response cell-phase analysis revealed that such periodicity cohered with latent EMT6/Ro cell-phase temporal patterning. Taken together, this study provides powerful evidence towards the hypothesis that even simple inter-fraction timing variations for a given fractional dose program may present a facile, and highly cost-effecitive means of significantly improving clinical efficacy.

Gompertz model with delays and treatment: mathematical analysis.

In this paper we study the delayed Gompertz model, as a typical model of tumor growth, with a term describing external interference that can reflect a treatment, e.g. chemotherapy. We mainly consider two types of delayed models, the one with the delay introduced in the per capita growth rate (we call it the single delayed model) and the other with the delay introduced in the net growth rate (the double delayed model). We focus on stability and possible stability switches with increasing delay for the positive steady state. Moreover, we study a Hopf bifurcation, including stability of arising periodic solutions for a constant treatment. The analytical results are extended by numerical simulations for a pharmacokinetic treatment function.

A simple model of carcinogenic mutations with time delay and diffusion.

In the paper we consider a system of delay differential equations (DDEs) of Lotka-Volterra type with diffusion reflecting mutations from normal to malignant cells. The model essentially follows the idea of Ahangar and Lin (2003) where mutations in three different environmental conditions, namely favorable, competitive and unfavorable, were considered. We focus on the unfavorable conditions that can result from a given treatment, e.g. chemotherapy. Included delay stands for the interactions between benign and other cells. We compare the dynamics of ODEs system, the system with delay and the system with delay and diffusion. We mainly focus on the dynamics when a positive steady state exists. The system which is globally stable in the case without the delay and diffusion is destabilized by increasing delay, and therefore the underlying kinetic dynamics becomes oscillatory due to a Hopf bifurcation for appropriate values of the delay. This suggests the occurrence of spatially non-homogeneous periodic solutions for the system with the delay and diffusion.

Model of tumour angiogenesis -- analysis of stability with respect to delays.

In the paper we consider the model of tumour angiogenesis process proposed by Bodnar and Fory (2009). The model combines ideas of Hahnfeldt et al. (1999) and Agur et al. (2004) describing the dynamics of tumour, angiogenic proteins and effective vessels density. Presented analysis is focused on the dependance of the model dynamics on delays introduced to the system. These delays reflect time lags in the proliferation/death term and the vessel formation/regression response to stimuli. It occurs that the dynamics strongly depends on the model parameters and the behaviour independent of the delays magnitude as well as multiple stability switches with increasing delay can be obtained.