Lattice-Gas Cellular Automaton Models for Biology: From Fluids to Cells

Lattice-gas cellular automaton (LGCA) and lattice Boltzmann (LB) models are promising models for studying emergent behaviour of transport and interaction processes in biological systems. In this chapter, we will emphasise the use of LGCA/LB models and the derivation and analysis of LGCA models ranging from the classical example dynamics of fluid flow to clotting phenomena in cerebral aneurysms and the invasion of tumour cell

Lattice-gas cellular automata for the analysis of cancer invasion

Cancer cells display characteristic traits acquired in a step-wise manner during carcinogenesis. Some of these traits are autonomous growth, induction of angiogenesis, invasion and metastasis. In this thesis, the focus is on one of the latest stages of tumor progression, tumor invasion. Tumor invasion emerges from the combined effect of tumor cell-cell and cell-microenvironment interactions, which can be studied with the help of mathematical analysis. Cellular automata (CA) can be viewed as simple models of self-organizing complex systems in which collective behavior can emerge out of an ensemble of many interacting "simple" components. In particular, we focus on an important class of CA, the so-called lattice-gas cellular automata (LGCA). In contrast to traditional CA, LGCA provide a straightforward and intuitive implementation of particle transport and interactions. Additionally, the structure of LGCA facilitates the mathematical analysis of their behavior. Here, the principal tools of mathematical analysis of LGCA are the mean-field approximation and the corresponding Lattice Boltzmann equation. The main objective of this thesis is to investigate important aspects of tumor invasion, under the microscope of mathematical modeling and analysis: Impact of the tumor environment: We introduce a LGCA as a microscopic model of tumor cell migration together with a mathematical description of different tumor environments. We study the impact of the various tumor environments (such as extracellular matrix) on tumor cell migration by estimating the tumor cell dispersion speed for a given environment. Effect of tumor cell proliferation and migration: We study the effect of tumor cell proliferation and migration on the tumor’s invasive behavior by developing a simplified LGCA model of tumor growth. In particular, we derive the corresponding macroscopic dynamics and we calculate the tumor’s invasion speed in terms of tumor cell proliferation and migration rates. Moreover, we calculate the width of the invasive zone, where the majority of mitotic activity is concentrated, and it is found to be proportional to the invasion speed. Mechanisms of tumor invasion emergence: We investigate the mechanisms for the emergence of tumor invasion in the course of cancer progression. We conclude that the response of a microscopic intracellular mechanism (migration/proliferation dichotomy) to oxygen shortage, i.e. hypoxia, maybe responsible for the transition from a benign (proliferative) to a malignant (invasive) tumor. Computing in vivo tumor invasion: Finally, we propose an evolutionary algorithm that estimates the parameters of a tumor growth LGCA model based on time-series of patient medical data (in particular Magnetic Resonance and Diffusion Tensor Imaging data). These parameters may allow to reproduce clinically relevant tumor growth scenarios for a specific patient, providing a prediction of the tumor growth at a later time stage.Krebszellen zeigen charakteristische Merkmale, die sie in einem schrittweisen Vorgang während der Karzinogenese erworben haben. Einige dieser Merkmale sind autonomes Wachstum, die Induktion von Angiogenese, Invasion und Metastasis. Der Schwerpunkt dieser Arbeit liegt auf der Tumorinvasion, einer der letzten Phasen der Tumorprogression. Die Tumorinvasion ensteht aus der kombinierten Wirkung von den Wechselwirkungen Tumorzelle-Zelle und Zelle-Mikroumgebung, die mit die Hilfe von mathematischer Analyse untersucht werden können. Zelluläre Automaten (CA) können als einfache Modelle von selbst-organisierenden komplexen Systemen betrachtet werden, in denen kollektives Verhalten aus einer Kombination von vielen interagierenden "einfachen" Komponenten entstehen kann. Insbesondere konzentrieren wir uns auf eine wichtige CA-Klasse, die sogenannten Zelluläre Gitter-Gas Automaten (LGCA). Im Gegensatz zu traditionellen CA bieten LGCA eine einfache und intuitive Umsetzung der Teilchen und Wechselwirkungen. Zusätzlich erleichtert die Struktur der LGCA die mathematische Analyse ihres Verhaltens. Die wichtigsten Werkzeuge der mathematischen Analyse der LGCA sind hier die Mean-field Approximation und die entsprechende Lattice - Boltzmann - Gleichung. Das wichtigste Ziel dieser Arbeit ist es, wichtige Aspekte der Tumorinvasion unter dem Mikroskop der mathematischen Modellierung und Analyse zu erforschen: Auswirkungen der Tumorumgebung: Wir stellen einen LGCA als mikroskopisches Modell der Tumorzellen-Migration in Verbindung mit einer mathematischen Beschreibung der verschiedenen Tumorumgebungen vor. Wir untersuchen die Auswirkungen der verschiedenen Tumorumgebungen (z. B. extrazellulären Matrix) auf die Migration von Tumorzellen dürch Schätzung der Tumorzellen-Dispersionsgeschwindigkeit in einem gegebenen Umfeld. Wirkung von Tumor-Zellenproliferation und Migration: Wir untersuchen die Wirkung von Tumorzellenproliferation und Migration auf das invasive Verhalten der Tumorzellen durch die Entwicklung eines vereinfachten LGCA Tumorwachstumsmodells. Wir leiten die entsprechende makroskopische Dynamik und berechnen die Tumorinvasionsgeschwindigkeit im Hinblick auf die Tumorzellenproliferation- und Migrationswerte. Darüber hinaus berechnen wir die Breite der invasiven Zone, wo die Mehrheit der mitotischer Aktivität konzentriert ist, und es wird festgestellt, dass diese proportional zu den Invasionsgeschwindigkeit ist. Mechanismen der Tumorinvasion Entstehung: Wir untersuchen Mechanismen, die für die Entstehung von Tumorinvasion im Verlauf des Krebs zuständig sind. Wir kommen zu dem Schluss, dass die Reaktion eines mikroskopischen intrazellulären Mechanismus (Migration/Proliferation Dichotomie) zu Sauerstoffmangel, d.h. Hypoxie, möglicheweise für den Übergang von einem gutartigen (proliferative) zu einer bösartigen (invasive) Tumor verantwortlich ist. Berechnung der in-vivo Tumorinvasion: Schließlich schlagen wir einen evolutionären Algorithmus vor, der die Parameter eines LGCA Modells von Tumorwachstum auf der Grundlage von medizinischen Daten des Patienten für mehrere Zeitpunkte (insbesondere die Magnet-Resonanz-und Diffusion Tensor Imaging Daten) ermöglicht. Diese Parameter erlauben Szenarien für einen klinisch relevanten Tumorwachstum für einen bestimmten Patienten zu reproduzieren, die eine Vorhersage des Tumorwachstums zu einem späteren Zeitpunkt möglich machen

Travelling wave analysis of a mathematical model of glioblastoma growth

In this paper we analyse a previously proposed cell-based model of glioblastoma (brain tumour) growth, which is based on the assumption that the cancer cells switch phenotypes between a proliferative and motile state (Gerlee and Nelander, PLoS Comp. Bio., 8(6) 2012). The dynamics of this model can be described by a system of partial differential equations, which exhibits travelling wave solutions whose wave speed depends crucially on the rates of phenotypic switching. We show that under certain conditions on the model parameters, a closed form expression of the wave speed can be obtained, and using singular perturbation methods we also derive an approximate expression of the wave front shape. These new analytical results agree with simulations of the cell-based model, and importantly show that the inverse relationship between wave front steepness and speed observed for the Fisher equation no longer holds when phenotypic switching is considered.Comment: Corrected error in the equation for the Jacobia

Cell-Based Modeling

A cell-based model is a simulation model that predicts collective behavior of cell-clusters from the behavior and interactions of individual cells. The inputs to a cell-based model are cell behaviors as observed in experiments or deriving from single cell models, including the cellular responses to cues from the micro-environment. The cell behaviors are encoded in a set of biologically plausible rules that the simulated cells will follow. The outputs of a cell-based model are the patterns and behaviors that follow indirectly from the cell behaviors and the cellular interactions. Cell-based models resemble agent-based models, but typically contain more biophysically-detailed descriptions of the individual cells

Mesoscopic Methods in Engineering and Science

(First paragraph) Matter, conceptually classified into fluids and solids, can be completely described by the microscopic physics of its constituent atoms or molecules. However, for most engineering applications a macroscopic or continuum description has usually been sufficient, because of the large disparity between the spatial and temporal scales relevant to these applications and the scales of the underlying molecular dynamics. In this case, the microscopic physics merely determines material properties such as the viscosity of a fluid or the elastic constants of a solid. These material properties cannot be derived within the macroscopic framework, but the qualitative nature of the macroscopic dynamics is usually insensitive to the details of the underlying microscopic interactions

Simulating tissue mechanics with Agent Based Models: concepts and perspectives

International audienceIn this paper we present an overview of agent based models that are used to simulate mechanical and physiological phenomena in cells and tissues, and we discuss underlying concepts, limitations and future perspectives of these models. As the interest in cell and tissue mechanics increase, agent based models are becoming more common the modeling community. We overview the physical aspects, complexity, shortcomings and capabilities of the major agent based model categories: lattice-based models (cellular automata, lattice gas cellular automata, cellular Potts models), off-lattice models (center based models, deformable cell models, vertex models), and hybrid discrete-continuum models. In this way, we hope to assist future researchers in choosing a model for the phenomenon they want to model and understand. The article also contains some novel results

Evolutionary games on graphs

Game theory is one of the key paradigms behind many scientific disciplines from biology to behavioral sciences to economics. In its evolutionary form and especially when the interacting agents are linked in a specific social network the underlying solution concepts and methods are very similar to those applied in non-equilibrium statistical physics. This review gives a tutorial-type overview of the field for physicists. The first three sections introduce the necessary background in classical and evolutionary game theory from the basic definitions to the most important results. The fourth section surveys the topological complications implied by non-mean-field-type social network structures in general. The last three sections discuss in detail the dynamic behavior of three prominent classes of models: the Prisoner's Dilemma, the Rock-Scissors-Paper game, and Competing Associations. The major theme of the review is in what sense and how the graph structure of interactions can modify and enrich the picture of long term behavioral patterns emerging in evolutionary games.Comment: Review, final version, 133 pages, 65 figure

Mean field mutation dynamics and the continuous Luria-Delbr\"uck distribution

The Luria-Delbr\"uck mutation model has a long history and has been mathematically formulated in several different ways. Here we tackle the problem in the case of a continuous distribution using some mathematical tools from nonlinear statistical physics. Starting from the classical formulations we derive the corresponding differential models and show that under a suitable mean field scaling they correspond to generalized Fokker-Planck equations for the mutants distribution whose solutions are given by the corresponding Luria-Delbr\"uck distribution. Numerical results confirming the theoretical analysis are also presented

ZIH-Info

- Erweiterung des VPN-Dienstes - Aktivierung der Diskquoten - Kommerzielle Software-Lizenzen - HRSK-Wartungsarbeiten - Uni-Tag am 21. Mai 2011 - ZIH-Kolloquium - Neue ZIH-Publikationen - Veranstaltunge