1
|
Nakano T, Okaie Y, Kinugasa Y, Koujin T, Suda T, Hiraoka Y, Haraguchi T. Roles of Remote and Contact Forces in Epithelial Cell Structure Formation. Biophys J 2020; 118:1466-1478. [PMID: 32097624 PMCID: PMC7091513 DOI: 10.1016/j.bpj.2020.01.037] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2019] [Revised: 01/25/2020] [Accepted: 01/29/2020] [Indexed: 12/28/2022] Open
Abstract
Cancer cells collectively form a large-scale structure for their growth. In this article, we report that HeLa cells, epithelial-like human cervical cancer cells, aggressively migrate on Matrigel and form a large-scale structure in a cell-density-dependent manner. To explain the experimental results, we develop a simple model in which cells interact and migrate using the two fundamentally different types of force, remote and contact forces, and show how cells form a large-scale structure. We demonstrate that the simple model reproduces experimental observations, suggesting that the remote and contact forces considered in this work play a major role in large-scale structure formation of HeLa cells. This article provides important evidence that cancer cells form a large-scale structure and develops an understanding into the poorly understood mechanisms of their structure formation.
Collapse
Affiliation(s)
- Tadashi Nakano
- Institute for Datability Science, Osaka University, Suita, Japan.
| | - Yutaka Okaie
- Institute for Datability Science, Osaka University, Suita, Japan
| | - Yasuha Kinugasa
- Graduate School of Frontier Biosciences, Osaka University, Suita, Japan
| | - Takako Koujin
- Advanced ICT Research Institute Kobe, National Institute of Information and Communications Technology, Kobe, Japan
| | | | - Yasushi Hiraoka
- Graduate School of Frontier Biosciences, Osaka University, Suita, Japan; Advanced ICT Research Institute Kobe, National Institute of Information and Communications Technology, Kobe, Japan
| | - Tokuko Haraguchi
- Graduate School of Frontier Biosciences, Osaka University, Suita, Japan; Advanced ICT Research Institute Kobe, National Institute of Information and Communications Technology, Kobe, Japan
| |
Collapse
|
2
|
Oelz D, Khataee H, Czirok A, Neufeld Z. Polarization wave at the onset of collective cell migration. Phys Rev E 2019; 100:032403. [PMID: 31640045 PMCID: PMC6894614 DOI: 10.1103/physreve.100.032403] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2018] [Indexed: 01/08/2023]
Abstract
Collective cell migration underlies morphogenesis, tissue regeneration, and cancer progression. How the biomechanical coupling between epithelial cells triggers and coordinates the collective migration is an open question. Here, we develop a one-dimensional model for an epithelial monolayer which predicts that after the onset of migration at an open boundary, cells in the bulk of the epithelium are gradually recruited into outward-directed motility, exhibiting traveling-wave-like behavior. We find an exact formula for the speed of this motility wave proportional to the square root of the cells' contractility, which accounts for cortex tension and adhesion between adjacent cells.
Collapse
Affiliation(s)
- Dietmar Oelz
- School of Mathematics and Physics, The University of Queensland, St. Lucia, Brisbane, QLD 4072, Australia
| | - Hamid Khataee
- School of Mathematics and Physics, The University of Queensland, St. Lucia, Brisbane, QLD 4072, Australia
| | - Andras Czirok
- Department of Biological Physics, Eotvos University, Budapest, 1053, Hungary
- Department of Anatomy and Cell Biology, University of Kansas Medical Center, Kansas City, KS 66160, USA
| | - Zoltan Neufeld
- School of Mathematics and Physics, The University of Queensland, St. Lucia, Brisbane, QLD 4072, Australia
| |
Collapse
|
3
|
A stochastic model of corneal epithelium maintenance and recovery following perturbation. J Math Biol 2018; 78:1245-1276. [PMID: 30478759 DOI: 10.1007/s00285-018-1308-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2018] [Revised: 11/01/2018] [Indexed: 01/15/2023]
Abstract
Various biological studies suggest that the corneal epithelium is maintained by active stem cells located in the limbus, the so-called limbal epithelial stem cell hypothesis. While numerous mathematical models have been developed to describe corneal epithelium wound healing, only a few have explored the process of corneal epithelium homeostasis. In this paper we present a purposefully simple stochastic mathematical model based on a chemical master equation approach, with the aim of clarifying the main factors involved in the maintenance process. Model analysis provides a set of constraints on the numbers of stem cells, division rates, and the number of division cycles required to maintain a healthy corneal epithelium. In addition, our stochastic analysis reveals noise reduction as the epithelium approaches its homeostatic state, indicating robustness to noise. Finally, recovery is analysed in the context of perturbation scenarios.
Collapse
|
4
|
Gavagnin E, Ford MJ, Mort RL, Rogers T, Yates CA. The invasion speed of cell migration models with realistic cell cycle time distributions. J Theor Biol 2018; 481:91-99. [PMID: 30219568 DOI: 10.1016/j.jtbi.2018.09.010] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2018] [Revised: 09/07/2018] [Accepted: 09/10/2018] [Indexed: 01/02/2023]
Abstract
Cell proliferation is typically incorporated into stochastic mathematical models of cell migration by assuming that cell divisions occur after an exponentially distributed waiting time. Experimental observations, however, show that this assumption is often far from the real cell cycle time distribution (CCTD). Recent studies have suggested an alternative approach to modelling cell proliferation based on a multi-stage representation of the CCTD. In this paper we investigate the connection between the CCTD and the speed of the collective invasion. We first state a result for a general CCTD, which allows the computation of the invasion speed using the Laplace transform of the CCTD. We use this to deduce the range of speeds for the general case. We then focus on the more realistic case of multi-stage models, using both a stochastic agent-based model and a set of reaction-diffusion equations for the cells' average density. By studying the corresponding travelling wave solutions, we obtain an analytical expression for the speed of invasion for a general N-stage model with identical transition rates, in which case the resulting cell cycle times are Erlang distributed. We show that, for a general N-stage model, the Erlang distribution and the exponential distribution lead to the minimum and maximum invasion speed, respectively. This result allows us to determine the range of possible invasion speeds in terms of the average proliferation time for any multi-stage model.
Collapse
Affiliation(s)
- Enrico Gavagnin
- Department of Mathematical Sciences University of Bath, Claverton Down, Bath, BA2 7AY, UK.
| | - Matthew J Ford
- Centre for Research in Reproduction and Development McGill University, Montréal, H3G 1Y6, Québec
| | - Richard L Mort
- Division of Biomedical and Life Sciences Faculty of Health and Medicine Lancaster University, Bailrigg, Lancaster LA1 4YG, UK
| | - Tim Rogers
- Department of Mathematical Sciences University of Bath, Claverton Down, Bath, BA2 7AY, UK
| | - Christian A Yates
- Department of Mathematical Sciences University of Bath, Claverton Down, Bath, BA2 7AY, UK
| |
Collapse
|
5
|
Johnston ST, Baker RE, Simpson MJ. Filling the gaps: A robust description of adhesive birth-death-movement processes. Phys Rev E 2016; 93:042413. [PMID: 27176335 DOI: 10.1103/physreve.93.042413] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2015] [Indexed: 11/07/2022]
Abstract
Existing continuum descriptions of discrete adhesive birth-death-movement processes provide accurate predictions of the average discrete behavior for limited parameter regimes. Here we present an alternative continuum description in terms of the dynamics of groups of contiguous occupied and vacant lattice sites. Our method provides more accurate predictions, is valid in parameter regimes that could not be described by previous continuum descriptions, and provides information about the spatial clustering of occupied sites. Furthermore, we present a simple analytic approximation of the spatial clustering of occupied sites at late time, when the system reaches its steady-state configuration.
Collapse
Affiliation(s)
- Stuart T Johnston
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| | - Ruth E Baker
- Mathematical Institute, University of Oxford, Oxford, United Kingdom
| | - Matthew J Simpson
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| |
Collapse
|
6
|
Johnston ST, Simpson MJ, Baker RE. Modelling the movement of interacting cell populations: A moment dynamics approach. J Theor Biol 2015; 370:81-92. [DOI: 10.1016/j.jtbi.2015.01.025] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2014] [Revised: 01/16/2015] [Accepted: 01/20/2015] [Indexed: 01/26/2023]
|
7
|
Castiglione F, Pappalardo F, Bianca C, Russo G, Motta S. Modeling biology spanning different scales: an open challenge. BIOMED RESEARCH INTERNATIONAL 2014; 2014:902545. [PMID: 25143952 PMCID: PMC4124842 DOI: 10.1155/2014/902545] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/16/2014] [Accepted: 06/25/2014] [Indexed: 02/03/2023]
Abstract
It is coming nowadays more clear that in order to obtain a unified description of the different mechanisms governing the behavior and causality relations among the various parts of a living system, the development of comprehensive computational and mathematical models at different space and time scales is required. This is one of the most formidable challenges of modern biology characterized by the availability of huge amount of high throughput measurements. In this paper we draw attention to the importance of multiscale modeling in the framework of studies of biological systems in general and of the immune system in particular.
Collapse
Affiliation(s)
- Filippo Castiglione
- Institute for Applied Mathematics, National Research Council of Italy, Rome, Italy
| | | | - Carlo Bianca
- Theoretical Physics of Condensed Matter, Sorbonne Universities, UPMC Univ Paris 6, 75252 Paris Cedex 05, France
- UMR 7600 LPTMC, CNRS, 75252 Paris Cedex 05, France
| | - Giulia Russo
- Department of Pharmaceutical Sciences, University of Catania, Catania, Italy
| | - Santo Motta
- Department of Mathematics and Computer Science, University of Catania, 95125 Catania, Italy
| |
Collapse
|
8
|
Middleton AM, Fleck C, Grima R. A continuum approximation to an off-lattice individual-cell based model of cell migration and adhesion. J Theor Biol 2014; 359:220-32. [PMID: 24972155 DOI: 10.1016/j.jtbi.2014.06.011] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2013] [Revised: 06/05/2014] [Accepted: 06/10/2014] [Indexed: 11/24/2022]
Abstract
Cell-cell adhesion plays a key role in the collective migration of cells and in determining correlations in the relative cell positions and velocities. Recently, it was demonstrated that off-lattice individual cell based models (IBMs) can accurately capture the correlations observed experimentally in a migrating cell population. However, IBMs are often computationally expensive and difficult to analyse mathematically. Traditional continuum-based models, in contrast, are amenable to mathematical analysis and are computationally less demanding, but typically correspond to a mean-field approximation of cell migration and so ignore cell-cell correlations. In this work, we address this problem by using an off-lattice IBM to derive a continuum approximation which does take into account correlations. We furthermore show that a mean-field approximation of the off-lattice IBM leads to a single partial integro-differential equation of the same form as proposed by Sherratt and co-workers to model cell adhesion. The latter is found to be only effective at approximating the ensemble averaged cell number density when mechanical interactions between cells are weak. In contrast, the predictions of our novel continuum model for the time-evolution of the ensemble cell number density distribution and of the density-density correlation function are in close agreement with those obtained from the IBM for a wide range of mechanical interaction strengths. In particular, we observe 'front-like' propagation of cells in simulations using both our IBM and our continuum model, but not in the continuum model simulations obtained using the mean-field approximation.
Collapse
Affiliation(s)
- Alistair M Middleton
- University of Heidelberg, Im Neuenheimer Feld 267, 69120 Heidelberg, Germany; Center for Biological Systems Analysis, University of Freiburg, Habsburgerstr. 49, 79104 Freiburg, Germany
| | - Christian Fleck
- Laboratory for Systems and Synthetic Biology, Dreijenplein 10, 6703HB Wageningen, The Netherlands
| | - Ramon Grima
- SynthSys and School of Biological Sciences, University of Edinburgh, EH9 3JR Edinburgh, UK.
| |
Collapse
|
9
|
Meyer M, Schimansky-Geier L, Romanczuk P. Active Brownian agents with concentration-dependent chemotactic sensitivity. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:022711. [PMID: 25353513 DOI: 10.1103/physreve.89.022711] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/30/2013] [Indexed: 06/04/2023]
Abstract
We study a biologically motivated model of overdamped, autochemotactic Brownian agents with concentration-dependent chemotactic sensitivity. The agents in our model move stochastically and produce a chemical ligand at their current position. The ligand concentration obeys a reaction-diffusion equation and acts as a chemoattractant for the agents, which bias their motion towards higher concentrations of the dynamically altered chemical field. We explore the impact of concentration-dependent response to chemoattractant gradients on large-scale pattern formation, by deriving a coarse-grained macroscopic description of the individual-based model, and compare the conditions for emergence of inhomogeneous solutions for different variants of the chemotactic sensitivity. We focus primarily on the so-called receptor-law sensitivity, which models a nonlinear decrease of chemotactic sensitivity with increasing ligand concentration. Our results reveal qualitative differences between the receptor law, the constant chemotactic response, and the so-called log law, with respect to stability of the homogeneous solution, as well as the emergence of different patterns (labyrinthine structures, clusters, and bubbles) via spinodal decomposition or nucleation. We discuss two limiting cases, where the model can be reduced to the dynamics of single species: (I) the agent density governed by a density-dependent effective diffusion coefficient and (II) the ligand field with an effective bistable, time-dependent reaction rate. In the end, we turn to single clusters of agents, studying domain growth and determining mean characteristics of the stationary inhomogeneous state. Analytical results are confirmed and extended by large-scale GPU simulations of the individual based model.
Collapse
Affiliation(s)
- Marcel Meyer
- Department of Physics, Humboldt Universität zu Berlin, Newtonstraße 15, 12489 Berlin, Germany
| | - Lutz Schimansky-Geier
- Department of Physics, Humboldt Universität zu Berlin, Newtonstraße 15, 12489 Berlin, Germany
| | - Pawel Romanczuk
- Physikalisch-Technische Bundesanstalt, Abbestraße 2-12, 10587 Berlin, Germany
| |
Collapse
|
10
|
Effects of bursty protein production on the noisy oscillatory properties of downstream pathways. Sci Rep 2013; 3:2438. [PMID: 23942456 PMCID: PMC3743060 DOI: 10.1038/srep02438] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2013] [Accepted: 07/30/2013] [Indexed: 01/06/2023] Open
Abstract
Experiments show that proteins are translated in sharp bursts; similar bursty phenomena have been observed for protein import into compartments. Here we investigate the effect of burstiness in protein expression and import on the stochastic properties of downstream pathways. We consider two identical pathways with equal mean input rates, except in one pathway proteins are input one at a time and in the other proteins are input in bursts. Deterministically the dynamics of these two pathways are indistinguishable. However the stochastic behavior falls in three categories: (i) both pathways display or do not display noise-induced oscillations; (ii) the non-bursty input pathway displays noise-induced oscillations whereas the bursty one does not; (iii) the reverse of (ii). We derive necessary conditions for these three cases to classify systems involving autocatalysis, trimerization and genetic feedback loops. Our results suggest that single cell rhythms can be controlled by regulation of burstiness in protein production.
Collapse
|
11
|
Walpole J, Papin JA, Peirce SM. Multiscale computational models of complex biological systems. Annu Rev Biomed Eng 2013; 15:137-54. [PMID: 23642247 DOI: 10.1146/annurev-bioeng-071811-150104] [Citation(s) in RCA: 121] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Integration of data across spatial, temporal, and functional scales is a primary focus of biomedical engineering efforts. The advent of powerful computing platforms, coupled with quantitative data from high-throughput experimental methodologies, has allowed multiscale modeling to expand as a means to more comprehensively investigate biological phenomena in experimentally relevant ways. This review aims to highlight recently published multiscale models of biological systems, using their successes to propose the best practices for future model development. We demonstrate that coupling continuous and discrete systems best captures biological information across spatial scales by selecting modeling techniques that are suited to the task. Further, we suggest how to leverage these multiscale models to gain insight into biological systems using quantitative biomedical engineering methods to analyze data in nonintuitive ways. These topics are discussed with a focus on the future of the field, current challenges encountered, and opportunities yet to be realized.
Collapse
Affiliation(s)
- Joseph Walpole
- Department of Biomedical Engineering, University of Virginia, Charlottesville, VA 22908, USA
| | | | | |
Collapse
|
12
|
Holzhütter HG, Drasdo D, Preusser T, Lippert J, Henney AM. The virtual liver: a multidisciplinary, multilevel challenge for systems biology. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2012; 4:221-35. [PMID: 22246674 DOI: 10.1002/wsbm.1158] [Citation(s) in RCA: 81] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
The liver is the central metabolic organ in human physiology, with functions that are fundamentally important to the detoxification of xenobiotics (drugs), the maintenance of homeostasis of numerous blood metabolites, and the production of mediators of the acute phase response. Liver toxicity, whether actual or implied is the reason for the failure of a significant proportion of many promising novel medicines that consequently never reach the market, and diseases such as atherosclerosis, diabetes, and fatty liver diseases, that are a major burden on current health resources, are directly linked to functional and structural disorders of the liver. This article presents the concepts and approaches underpinning one of the most exciting and ambitious modeling projects in the field of systems biology and systems medicine. This major multidisciplinary research program is aimed at developing a whole-organ model of the human liver, representing its central physiological functions under normal and pathological conditions The model will be composed of a larger battery of interconnected submodels representing liver anatomy and physiology, integrating processes across hierarchical levels in space, time, and structural organization. In this article, we outline the general architecture of the liver model and present first step taken to reach this ambitious goal.
Collapse
|
13
|
Qu Z, Garfinkel A, Weiss JN, Nivala M. Multi-scale modeling in biology: how to bridge the gaps between scales? PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2011; 107:21-31. [PMID: 21704063 DOI: 10.1016/j.pbiomolbio.2011.06.004] [Citation(s) in RCA: 95] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2011] [Accepted: 06/11/2011] [Indexed: 11/25/2022]
Abstract
Human physiological functions are regulated across many orders of magnitude in space and time. Integrating the information and dynamics from one scale to another is critical for the understanding of human physiology and the treatment of diseases. Multi-scale modeling, as a computational approach, has been widely adopted by researchers in computational and systems biology. A key unsolved issue is how to represent appropriately the dynamical behaviors of a high-dimensional model of a lower scale by a low-dimensional model of a higher scale, so that it can be used to investigate complex dynamical behaviors at even higher scales of integration. In the article, we first review the widely-used different modeling methodologies and their applications at different scales. We then discuss the gaps between different modeling methodologies and between scales, and discuss potential methods for bridging the gaps between scales.
Collapse
Affiliation(s)
- Zhilin Qu
- Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, CA 90095, USA.
| | | | | | | |
Collapse
|
14
|
Mocenni C, Facchini A, Vicino A. Comparison of recurrence quantification methods for the analysis of temporal and spatial chaos. ACTA ACUST UNITED AC 2011. [DOI: 10.1016/j.mcm.2010.04.008] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
15
|
Investigating the robustness of the classical enzyme kinetic equations in small intracellular compartments. BMC SYSTEMS BIOLOGY 2009; 3:101. [PMID: 19814817 PMCID: PMC2778647 DOI: 10.1186/1752-0509-3-101] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/02/2009] [Accepted: 10/08/2009] [Indexed: 11/26/2022]
Abstract
Background Classical descriptions of enzyme kinetics ignore the physical nature of the intracellular environment. Main implicit assumptions behind such approaches are that reactions occur in compartment volumes which are large enough so that molecular discreteness can be ignored and that molecular transport occurs via diffusion. Though these conditions are frequently met in laboratory conditions, they are not characteristic of the intracellular environment, which is compartmentalized at the micron and submicron scales and in which active means of transport play a significant role. Results Starting from a master equation description of enzyme reaction kinetics and assuming metabolic steady-state conditions, we derive novel mesoscopic rate equations which take into account (i) the intrinsic molecular noise due to the low copy number of molecules in intracellular compartments (ii) the physical nature of the substrate transport process, i.e. diffusion or vesicle-mediated transport. These equations replace the conventional macroscopic and deterministic equations in the context of intracellular kinetics. The latter are recovered in the limit of infinite compartment volumes. We find that deviations from the predictions of classical kinetics are pronounced (hundreds of percent in the estimate for the reaction velocity) for enzyme reactions occurring in compartments which are smaller than approximately 200 nm, for the case of substrate transport to the compartment being mediated principally by vesicle or granule transport and in the presence of competitive enzyme inhibitors. Conclusion The derived mesoscopic rate equations describe subcellular enzyme reaction kinetics, taking into account, for the first time, the simultaneous influence of both intrinsic noise and the mode of transport. They clearly show the range of applicability of the conventional deterministic equation models, namely intracellular conditions compatible with diffusive transport and simple enzyme mechanisms in several hundred nanometre-sized compartments. An active transport mechanism coupled with large intrinsic noise in enzyme concentrations is shown to lead to huge deviations from the predictions of deterministic models. This has implications for the common approach of modeling large intracellular reaction networks using ordinary differential equations and also for the calculation of the effective dosage of competitive inhibitor drugs.
Collapse
|