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

Multiscale modelling of vascular tumour growth in 3D: the roles of domain size & boundary condition

We investigate a three-dimensional multiscale model of vascular tumour growth, which couples blood flow, angiogenesis, vascular remodelling, nutrient/growth factor transport, movement of, and interactions between, normal and tumour cells, and nutrient-dependent cell cycle dynamics within each cell. In particular, we determine how the domain size, aspect ratio and initial vascular network influence the tumour's growth dynamics and its long-time composition. We establish whether it is possible to extrapolate simulation results obtained for small domains to larger ones, by constructing a large simulation domain from a number of identical subdomains, each subsystem initially comprising two parallel parent vessels, with associated cells and diffusible substances. We find that the subsystem is not representative of the full domain and conclude that, for this initial vessel geometry, interactions between adjacent subsystems contribute to the overall growth dynamics. We then show that extrapolation of results from a small subdomain to a larger domain can only be made if the subdomain is sufficiently large and is initialised with a sufficiently complex vascular network. Motivated by these results, we perform simulations to investigate the tumour's response to therapy and show that the probability of tumour elimination in a larger domain can be extrapolated from simulation results on a smaller domain. Finally, we demonstrate how our model may be combined with experimental data, to predict the spatio-temporal evolution of a vascular tumour

Systems biologists seek fuller integration of systems biology approaches in new cancer research programs

Systems biology takes an interdisciplinary approach to the systematic study of complex interactions in biological systems. This approach seeks to decipher the emergent behaviors of complex systems rather than focusing only on their constituent properties. As an increasing number of examples illustrate the value of systems biology approaches to understand the initiation, progression, and treatment of cancer, systems biologists from across Europe and the United States hope for changes in the way their field is currently perceived among cancer researchers. In a recent EU-US workshop, supported by the European Commission, the German Federal Ministry for Education and Research, and the National Cancer Institute of the NIH, the participants discussed the strengths, weaknesses, hurdles, and opportunities in cancer systems biology

Colloquium: Mechanical formalisms for tissue dynamics

The understanding of morphogenesis in living organisms has been renewed by tremendous progressin experimental techniques that provide access to cell-scale, quantitative information both on theshapes of cells within tissues and on the genes being expressed. This information suggests that ourunderstanding of the respective contributions of gene expression and mechanics, and of their crucialentanglement, will soon leap forward. Biomechanics increasingly benefits from models, which assistthe design and interpretation of experiments, point out the main ingredients and assumptions, andultimately lead to predictions. The newly accessible local information thus calls for a reflectionon how to select suitable classes of mechanical models. We review both mechanical ingredientssuggested by the current knowledge of tissue behaviour, and modelling methods that can helpgenerate a rheological diagram or a constitutive equation. We distinguish cell scale ("intra-cell")and tissue scale ("inter-cell") contributions. We recall the mathematical framework developpedfor continuum materials and explain how to transform a constitutive equation into a set of partialdifferential equations amenable to numerical resolution. We show that when plastic behaviour isrelevant, the dissipation function formalism appears appropriate to generate constitutive equations;its variational nature facilitates numerical implementation, and we discuss adaptations needed in thecase of large deformations. The present article gathers theoretical methods that can readily enhancethe significance of the data to be extracted from recent or future high throughput biomechanicalexperiments.Comment: 33 pages, 20 figures. This version (26 Sept. 2015) contains a few corrections to the published version, all in Appendix D.2 devoted to large deformation

Agent based modelling helps in understanding the rules by which fibroblasts support keratinocyte colony formation

Background: Autologous keratincoytes are routinely expanded using irradiated mouse fibroblasts and bovine serum for clinical use. With growing concerns about the safety of these xenobiotic materials, it is desirable to culture keratinocytes in media without animal derived products. An improved understanding of epithelial/mesenchymal interactions could assist in this. Methodology/Principal Findings: A keratincyte/fibroblast o-culture model was developed by extending an agent-based keratinocyte colony formation model to include the response of keratinocytes to both fibroblasts and serum. The model was validated by comparison of the in virtuo and in vitro multicellular behaviour of keratinocytes and fibroblasts in single and co-culture in Greens medium. To test the robustness of the model, several properties of the fibroblasts were changed to investigate their influence on the multicellular morphogenesis of keratinocyes and fibroblasts. The model was then used to generate hypotheses to explore the interactions of both proliferative and growth arrested fibroblasts with keratinocytes. The key predictions arising from the model which were confirmed by in vitro experiments were that 1) the ratio of fibroblasts to keratinocytes would critically influence keratinocyte colony expansion, 2) this ratio needed to be optimum at the beginning of the co-culture, 3) proliferative fibroblasts would be more effective than irradiated cells in expanding keratinocytes and 4) in the presence of an adequate number of fibroblasts, keratinocyte expansion would be independent of serum. Conclusions: A closely associated computational and biological approach is a powerful tool for understanding complex biological systems such as the interactions between keratinocytes and fibroblasts. The key outcome of this study is the finding that the early addition of a critical ratio of proliferative fibroblasts can give rapid keratinocyte expansion without the use of irradiated mouse fibroblasts and bovine serum

Growth Based Morphogenesis of Vertebrate Limb Bud

Many genes and their regulatory relationships are involved in developmental phenomena. However, by chemical information alone, we cannot fully understand changing organ morphologies through tissue growth because deformation and growth of the organ are essentially mechanical processes. Here, we develop a mathematical model to describe the change of organ morphologies through cell proliferation. Our basic idea is that the proper specification of localized volume source (e.g., cell proliferation) is able to guide organ morphogenesis, and that the specification is given by chemical gradients. We call this idea “growth-based morphogenesis.” We find that this morphogenetic mechanism works if the tissue is elastic for small deformation and plastic for large deformation. To illustrate our concept, we study the development of vertebrate limb buds, in which a limb bud protrudes from a flat lateral plate and extends distally in a self-organized manner. We show how the proportion of limb bud shape depends on different parameters and also show the conditions needed for normal morphogenesis, which can explain abnormal morphology of some mutants. We believe that the ideas shown in the present paper are useful for the morphogenesis of other organs

Fragment-Based Learning of Visual Object Categories in Non-Human Primates

When we perceive a visual object, we implicitly or explicitly associate it with an object category we know. Recent research has shown that the visual system can use local, informative image fragments of a given object, rather than the whole object, to classify it into a familiar category. We have previously reported, using human psychophysical studies, that when subjects learn new object categories using whole objects, they incidentally learn informative fragments, even when not required to do so. However, the neuronal mechanisms by which we acquire and use informative fragments, as well as category knowledge itself, have remained unclear. Here we describe the methods by which we adapted the relevant human psychophysical methods to awake, behaving monkeys and replicated key previous psychophysical results. This establishes awake, behaving monkeys as a useful system for future neurophysiological studies not only of informative fragments in particular, but also of object categorization and category learning in general

The relationship between stimulus intensity and response amplitude for the photopic negative response of the flash electroretinogram

The aim of this study was to investigate the relationship between stimulus intensity and response amplitude for the photopic negative response (PhNR) of the flash ERG. Specific aims were (i) to determine whether a generalized Naka-Rushton function provided a good fit to the intensity-response data and (ii) to determine the variability of the parameters of the best-fitting Naka-Rushton models. Electroretinograms were recorded in 18 participants, on two occasions, using both DTL fibre and skin active electrodes, in response to Ganzfeld red stimuli (Lee filter "terry red") ranging in stimulus strength from -1.30 to 0.53 log cd.s.m(-2) (0.28-2.11 log phot td.s) presented over a steady blue background (Schott glass filter BG28; 3.9 log scot td). PhNR amplitude was measured from b-wave peak and from pre-stimulus baseline. The Naka-Rushton function was fitted to all intensity-response data, and parameters, 'n', 'Vmax' and 'K' were obtained. Coefficients of variation (CoV), and inter-ocular and inter-session limits of agreement (LoA) were calculated for both Naka-Rushton parameters. A generalized Naka-Rushton function was found to provide a good fit to the intensity-response data, except at the highest stimulus intensity, where a reduction in amplitude occurred in many individuals. The 'Vmax' parameter was less variable than 'K' for all intensity-response data. Variability was lower for DTL than skin electrodes, and for peak-to-trough PhNR measurements, compared to baseline-to-trough. This study has demonstrated for the first time that the Naka-Rushton model provides a useful means of quantifying the intensity-response relationship of the PhNR

A Sub-Cellular Viscoelastic Model for Cell Population Mechanics

Understanding the biomechanical properties and the effect of biomechanical force on epithelial cells is key to understanding how epithelial cells form uniquely shaped structures in two or three-dimensional space. Nevertheless, with the limitations and challenges posed by biological experiments at this scale, it becomes advantageous to use mathematical and ‘in silico’ (computational) models as an alternate solution. This paper introduces a single-cell-based model representing the cross section of a typical tissue. Each cell in this model is an individual unit containing several sub-cellular elements, such as the elastic plasma membrane, enclosed viscoelastic elements that play the role of cytoskeleton, and the viscoelastic elements of the cell nucleus. The cell membrane is divided into segments where each segment (or point) incorporates the cell's interaction and communication with other cells and its environment. The model is capable of simulating how cells cooperate and contribute to the overall structure and function of a particular tissue; it mimics many aspects of cellular behavior such as cell growth, division, apoptosis and polarization. The model allows for investigation of the biomechanical properties of cells, cell-cell interactions, effect of environment on cellular clusters, and how individual cells work together and contribute to the structure and function of a particular tissue. To evaluate the current approach in modeling different topologies of growing tissues in distinct biochemical conditions of the surrounding media, we model several key cellular phenomena, namely monolayer cell culture, effects of adhesion intensity, growth of epithelial cell through interaction with extra-cellular matrix (ECM), effects of a gap in the ECM, tensegrity and tissue morphogenesis and formation of hollow epithelial acini. The proposed computational model enables one to isolate the effects of biomechanical properties of individual cells and the communication between cells and their microenvironment while simultaneously allowing for the formation of clusters or sheets of cells that act together as one complex tissue

A Multicellular Model of Intestinal Crypt Buckling and Fission

Crypt fission is an in vivo tissue deformation process that is involved in both intestinal homeostasis and colorectal tumourigenesis. Despite its importance, the mechanics underlying crypt fission are currently poorly understood. Recent experimental development of organoids, organ-like buds cultured from crypt stem cells in vitro, has shown promise in shedding light on crypt fission. Drawing inspiration from observations of organoid growth and fission in vivo, we develop a computational model of a deformable epithelial tissue layer. Results from in silico experiments show the stiffness of cells and the proportions of cell subpopulations affect the nature of deformation in the epithelial layer. In particular, we find that increasing the proportion of stiffer cells in the layer increases the likelihood of crypt fission occurring. This is in agreement with and helps explain recent experimental work