Emergent Properties of Tumor Microenvironment in a Real-life Model of Multicell Tumor Spheroids

Multicellular tumor spheroids are an important {\it in vitro} model of the pre-vascular phase of solid tumors, for sizes well below the diagnostic limit: therefore a biophysical model of spheroids has the ability to shed light on the internal workings and organization of tumors at a critical phase of their development. To this end, we have developed a computer program that integrates the behavior of individual cells and their interactions with other cells and the surrounding environment. It is based on a quantitative description of metabolism, growth, proliferation and death of single tumor cells, and on equations that model biochemical and mechanical cell-cell and cell-environment interactions. The program reproduces existing experimental data on spheroids, and yields unique views of their microenvironment. Simulations show complex internal flows and motions of nutrients, metabolites and cells, that are otherwise unobservable with current experimental techniques, and give novel clues on tumor development and strong hints for future therapies.Comment: 20 pages, 10 figures. Accepted for publication in PLOS One. The published version contains links to a supplementary text and three video file

Bridging the gap between the micro- and the macro-world of tumors

At present it is still quite difficult to match the vast knowledge on the behavior of individual tumor cells with macroscopic measurements on clinical tumors. On the modeling side, we already know how to deal with many molecular pathways and cellular events, using systems of differential equations and other modeling tools, and ideally, we should be able to extend such a mathematical description up to the level of large tumor masses. An extended model should thus help us forecast the behavior of large tumors from our basic knowledge of microscopic processes. Unfortunately, the complexity of these processes makes it very difficult -- probably impossible -- to develop comprehensive analytical models. We try to bridge the gap with a simulation program which is based on basic biochemical and biophysical processes -- thereby building an effective computational model -- and in this paper we describe its structure, endeavoring to make the description sufficiently detailed and yet understandable.Comment: 24 pages, 10 figures. Accepted for publication in AIP Advances, in the special issue on the physics of cance

Implementing vertex dynamics models of cell populations in biology within a consistent computational framework

The dynamic behaviour of epithelial cell sheets plays a central role during development, growth, disease and wound healing. These processes occur as a result of cell adhesion, migration, division, differentiation and death, and involve multiple processes acting at the cellular and molecular level. Computational models offer a useful means by which to investigate and test hypotheses about these processes, and have played a key role in the study of cell–cell interactions. However, the necessarily complex nature of such models means that it is difficult to make accurate comparison between different models, since it is often impossible to distinguish between differences in behaviour that are due to the underlying model assumptions, and those due to differences in the in silico implementation of the model. In this work, an approach is described for the implementation of vertex dynamics models, a discrete approach that represents each cell by a polygon (or polyhedron) whose vertices may move in response to forces. The implementation is undertaken in a consistent manner within a single open source computational framework, Chaste, which comprises fully tested, industrial-grade software that has been developed using an agile approach. This framework allows one to easily change assumptions regarding force generation and cell rearrangement processes within these models. The versatility and generality of this framework is illustrated using a number of biological examples. In each case we provide full details of all technical aspects of our model implementations, and in some cases provide extensions to make the models more generally applicable

The Epitheliome: agent-based modelling of the social behaviour of cells

We have developed a new computational modelling paradigm for predicting the emergent behaviour resulting from the interaction of cells in epithelial tissue. As proof-of-concept, an agent-based model, in which there is a one-to-one correspondence between biological cells and software agents, has been coupled to a simple physical model. Behaviour of the computational model is compared with the growth characteristics of epithelial cells in monolayer culture, using growth media with low and physiological calcium concentrations. Results show a qualitative fit between the growth characteristics produced by the simulation and the in vitro cell models

Development and validation of computational models of cellular interaction

In this paper we take the view that computational models of biological systems should satisfy two conditions – they should be able to predict function at a systems biology level, and robust techniques of validation against biological models must be available. A modelling paradigm for developing a predictive computational model of cellular interaction is described, and methods of providing robust validation against biological models are explored, followed by a consideration of software issues

Towards a digital model of zebrafish embryogenesis. Integration of Cell Tracking and Gene Expression Quantification

We elaborate on a general framework composed of a set of computational tools to accurately quantificate cellular position and gene expression levels throughout early zebrafish embryogenesis captured over a time-lapse series of in vivo 3D images. Our modeling strategy involves nuclei detection, cell geometries extraction, automatic gene levels quantification and cell tracking to reconstruct cell trajectories and lineage tree which describe the animal development. Each cell in the embryo is then precisely described at each given time t by a vector composed of the cell 3D spatial coordinates (x; y; z) along with its gene expression level g. This comprehensive description of the embryo development is used to assess the general connection between genetic expression and cell movement. We also investigate genetic expression propagation between a cell and its progeny in the lineage tree. More to the point, this paper focuses on the evolution of the expression pattern of transcriptional factor goosecoid (gsc) through the gastrulation process between 6 and 9 hours post fertilization (hpf

Colorectal Cancer Through Simulation and Experiment

Colorectal cancer has continued to generate a huge amount of research interest over several decades, forming a canonical example of tumourigenesis since its use in Fearon and Vogelstein’s linear model of genetic mutation. Over time, the field has witnessed a transition from solely experimental work to the inclusion of mathematical biology and computer-based modelling. The fusion of these disciplines has the potential to provide valuable insights into oncologic processes, but also presents the challenge of uniting many diverse perspectives. Furthermore, the cancer cell phenotype defined by the ‘Hallmarks of Cancer’ has been extended in recent times and provides an excellent basis for future research. We present a timely summary of the literature relating to colorectal cancer, addressing the traditional experimental findings, summarising the key mathematical and computational approaches, and emphasising the role of the Hallmarks in current and future developments. We conclude with a discussion of interdisciplinary work, outlining areas of experimental interest which would benefit from the insight that mathematical and computational modelling can provide