1
|
Hernández JA, Chifflet S, Justet C, Torriglia A. A mathematical model of wound healing in bovine corneal endothelium. J Theor Biol 2023; 559:111374. [PMID: 36460056 DOI: 10.1016/j.jtbi.2022.111374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2022] [Revised: 11/24/2022] [Accepted: 11/25/2022] [Indexed: 11/30/2022]
Abstract
We developed a mathematical model to describe healing processes in bovine corneal endothelial (BCE) cells in culture, triggered by mechanical wounds with parallel edges. Previous findings from our laboratory show that, in these cases, BCE monolayers exhibit an approximately constant healing velocity. Also, that caspase-dependent apoptosis occurs, with the fraction of apoptotic cells increasing with the distance traveled by the healing edge. In addition, in this study we report the novel findings that, for wound scratch assays performed preserving the basal extracellular matrix: i) the healing cells increase their en face surface area in a characteristic fashion, and ii) the average length of the segments of the cell columns actively participating in the healing process increases linearly with time. These latter observations preclude the utilization of standard traveling wave formalisms to model wound healing in BCE cells. Instead, we developed and studied a simple phenomenological model based on a plausible formula for the spreading dynamics of the individual healing cells, that incorporates original evidence about the process in BCE cells. The model can be simulated to: i) obtain an approximately constant healing velocity; ii) reproduce the profile of the healing cell areas, and iii) obtain approximately linear time dependences of the mean cell area and average length of the front active segments per column. In view of its accuracy to account for the experimental observations, the model can also be acceptably employed to quantify the appearance of apoptotic cells during BCE wound healing. The strategy utilized here could offer a novel formal framework to represent modifications undergone by some epithelial cell lines during wound healing.
Collapse
Affiliation(s)
- Julio A Hernández
- Sección Biofísica y Biología de Sistemas, Facultad de Ciencias, Universidad de la República, Iguá s/n esq. Mataojo, 11400 Montevideo, Uruguay.
| | - Silvia Chifflet
- Departamento de Bioquímica, Facultad de Medicina, Universidad de la República, Gral. Flores 2125, 11800 Montevideo, Uruguay
| | - Cristian Justet
- Departamento de Bioquímica, Facultad de Medicina, Universidad de la República, Gral. Flores 2125, 11800 Montevideo, Uruguay
| | - Alicia Torriglia
- Centre de Recherche des Cordeliers, Sorbonne Université, Inserm, Université de Paris, F-75006 Paris, France
| |
Collapse
|
2
|
Murphy RJ, Buenzli PR, Tambyah TA, Thompson EW, Hugo HJ, Baker RE, Simpson MJ. The role of mechanical interactions in EMT. Phys Biol 2021; 18. [PMID: 33789261 DOI: 10.1088/1478-3975/abf425] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Accepted: 03/31/2021] [Indexed: 02/07/2023]
Abstract
The detachment of cells from the boundary of an epithelial tissue and the subsequent invasion of these cells into surrounding tissues is important for cancer development and wound healing, and is strongly associated with the epithelial-mesenchymal transition (EMT). Chemical signals, such as TGF-β, produced by surrounding tissue can be uptaken by cells and induce EMT. In this work, we present a novel cell-based discrete mathematical model of mechanical cellular relaxation, cell proliferation, and cell detachment driven by chemically-dependent EMT in an epithelial tissue. A continuum description of the model is then derived in the form of a novel nonlinear free boundary problem. Using the discrete and continuum models we explore how the coupling of chemical transport and mechanical interactions influences EMT, and postulate how this could be used to help control EMT in pathological situations.
Collapse
Affiliation(s)
- Ryan J Murphy
- Queensland University of Technology, Mathematical Sciences, Brisbane, Australia
| | - Pascal R Buenzli
- Queensland University of Technology, Mathematical Sciences, Brisbane, Australia
| | - Tamara A Tambyah
- Queensland University of Technology, Mathematical Sciences, Brisbane, Australia
| | - Erik W Thompson
- Queensland University of Technology, Institute of Health and Biomedical Innovation, Brisbane, Australia.,Queensland University of Technology, School of Biomedical Sciences, Faculty of Health, Brisbane, Australia.,Translational Research Institute, Brisbane, Australia
| | - Honor J Hugo
- Queensland University of Technology, Institute of Health and Biomedical Innovation, Brisbane, Australia.,Queensland University of Technology, School of Biomedical Sciences, Faculty of Health, Brisbane, Australia.,Translational Research Institute, Brisbane, Australia
| | - Ruth E Baker
- University of Oxford, Mathematical Institute, Oxford, United Kingdom
| | - Matthew J Simpson
- Queensland University of Technology, Mathematical Sciences, Brisbane, Australia
| |
Collapse
|
3
|
Tambyah TA, Murphy RJ, Buenzli PR, Simpson MJ. A free boundary mechanobiological model of epithelial tissues. Proc Math Phys Eng Sci 2020; 476:20200528. [PMID: 33362419 PMCID: PMC7735320 DOI: 10.1098/rspa.2020.0528] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2020] [Accepted: 10/21/2020] [Indexed: 12/12/2022] Open
Abstract
In this study, we couple intracellular signalling and cell-based mechanical properties to develop a novel free boundary mechanobiological model of epithelial tissue dynamics. Mechanobiological coupling is introduced at the cell level in a discrete modelling framework, and new reaction-diffusion equations are derived to describe tissue-level outcomes. The free boundary evolves as a result of the underlying biological mechanisms included in the discrete model. To demonstrate the accuracy of the continuum model, we compare numerical solutions of the discrete and continuum models for two different signalling pathways. First, we study the Rac-Rho pathway where cell- and tissue-level mechanics are directly related to intracellular signalling. Second, we study an activator-inhibitor system which gives rise to spatial and temporal patterning related to Turing patterns. In all cases, the continuum model and free boundary condition accurately reflect the cell-level processes included in the discrete model.
Collapse
Affiliation(s)
| | | | | | - Matthew J. Simpson
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| |
Collapse
|
4
|
Murphy RJ, Buenzli PR, Baker RE, Simpson MJ. Mechanical Cell Competition in Heterogeneous Epithelial Tissues. Bull Math Biol 2020; 82:130. [PMID: 32979100 DOI: 10.1007/s11538-020-00807-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2020] [Accepted: 09/09/2020] [Indexed: 12/11/2022]
Abstract
Mechanical cell competition is important during tissue development, cancer invasion, and tissue ageing. Heterogeneity plays a key role in practical applications since cancer cells can have different cell stiffness and different proliferation rates than normal cells. To study this phenomenon, we propose a one-dimensional mechanical model of heterogeneous epithelial tissue dynamics that includes cell-length-dependent proliferation and death mechanisms. Proliferation and death are incorporated into the discrete model stochastically and arise as source/sink terms in the corresponding continuum model that we derive. Using the new discrete model and continuum description, we explore several applications including the evolution of homogeneous tissues experiencing proliferation and death, and competition in a heterogeneous setting with a cancerous tissue competing for space with an adjacent normal tissue. This framework allows us to postulate new mechanisms that explain the ability of cancer cells to outcompete healthy cells through mechanical differences rather than an intrinsic proliferative advantage. We advise when the continuum model is beneficial and demonstrate why naively adding source/sink terms to a continuum model without considering the underlying discrete model may lead to incorrect results.
Collapse
Affiliation(s)
- Ryan J Murphy
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia.
| | - Pascal R Buenzli
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| | - Ruth E Baker
- Mathematical Institute, University of Oxford, Oxford, UK
| | - Matthew J Simpson
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| |
Collapse
|
5
|
Bubba F, Lorenzi T, Macfarlane FR. From a discrete model of chemotaxis with volume-filling to a generalized Patlak-Keller-Segel model. Proc Math Phys Eng Sci 2020; 476:20190871. [PMID: 32523414 DOI: 10.1098/rspa.2019.0871] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2020] [Accepted: 04/02/2020] [Indexed: 12/26/2022] Open
Abstract
We present a discrete model of chemotaxis whereby cells responding to a chemoattractant are seen as individual agents whose movement is described through a set of rules that result in a biased random walk. In order to take into account possible alterations in cellular motility observed at high cell densities (i.e. volume-filling), we let the probabilities of cell movement be modulated by a decaying function of the cell density. We formally show that a general form of the celebrated Patlak-Keller-Segel (PKS) model of chemotaxis can be formally derived as the appropriate continuum limit of this discrete model. The family of steady-state solutions of such a generalized PKS model are characterized and the conditions for the emergence of spatial patterns are studied via linear stability analysis. Moreover, we carry out a systematic quantitative comparison between numerical simulations of the discrete model and numerical solutions of the corresponding PKS model, both in one and in two spatial dimensions. The results obtained indicate that there is excellent quantitative agreement between the spatial patterns produced by the two models. Finally, we numerically show that the outcomes of the two models faithfully replicate those of the classical PKS model in a suitable asymptotic regime.
Collapse
Affiliation(s)
- Federica Bubba
- Sorbonne Universités, Universités Paris-Diderot, Laboratoire Jacques-Louis Lions, 75005 Paris, France
| | - Tommaso Lorenzi
- School of Mathematics and Statistics, University of St Andrews, St Andrews KY16 9SS, UK.,Department of Mathematical Sciences 'G. L. Lagrange', Dipartimento di Eccellenza 2018-2022, Politecnico di Torino, 10129 Torino, Italy
| | - Fiona R Macfarlane
- School of Mathematics and Statistics, University of St Andrews, St Andrews KY16 9SS, UK
| |
Collapse
|
6
|
Surendran A, Plank MJ, Simpson MJ. Spatial structure arising from chase-escape interactions with crowding. Sci Rep 2019; 9:14988. [PMID: 31628421 PMCID: PMC6800429 DOI: 10.1038/s41598-019-51565-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Accepted: 10/03/2019] [Indexed: 12/17/2022] Open
Abstract
Movement of individuals, mediated by localised interactions, plays a key role in numerous processes including cell biology and ecology. In this work, we investigate an individual-based model accounting for various intraspecies and interspecies interactions in a community consisting of two distinct species. In this framework we consider one species to be chasers and the other species to be escapees, and we focus on chase-escape dynamics where the chasers are biased to move towards the escapees, and the escapees are biased to move away from the chasers. This framework allows us to explore how individual-level directional interactions scale up to influence spatial structure at the macroscale. To focus exclusively on the role of motility and directional bias in determining spatial structure, we consider conservative communities where the number of individuals in each species remains constant. To provide additional information about the individual-based model, we also present a mathematically tractable deterministic approximation based on describing the evolution of the spatial moments. We explore how different features of interactions including interaction strength, spatial extent of interaction, and relative density of species influence the formation of the macroscale spatial patterns.
Collapse
Affiliation(s)
- Anudeep Surendran
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| | - Michael J Plank
- School of Mathematics and Statistics, University of Canterbury, Christchurch, New Zealand.,Te Pūnaha Matatini, A New Zealand Centre of Research Excellence, Auckland, New Zealand
| | - Matthew J Simpson
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia.
| |
Collapse
|
7
|
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
|
8
|
Browning AP, Woodhouse FG, Simpson MJ. Reversible signal transmission in an active mechanical metamaterial. Proc Math Phys Eng Sci 2019; 475:20190146. [PMID: 31423095 PMCID: PMC6694314 DOI: 10.1098/rspa.2019.0146] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2019] [Accepted: 06/11/2019] [Indexed: 01/01/2023] Open
Abstract
Mechanical metamaterials are designed to enable unique functionalities, but are typically limited by an initial energy state and require an independent energy input to function repeatedly. Our study introduces a theoretical active mechanical metamaterial that incorporates a biological reaction mechanism to overcome this key limitation of passive metamaterials. Our material allows for reversible mechanical signal transmission, where energy is reintroduced by the biologically motivated reaction mechanism. By analysing a coarse-grained continuous analogue of the discrete model, we find that signals can be propagated through the material by a travelling wave. Analysis of the continuum model provides the region of the parameter space that allows signal transmission, and reveals similarities with the well-known FitzHugh-Nagumo system. We also find explicit formulae that approximate the effect of the time scale of the reaction mechanism on the signal transmission speed, which is essential for controlling the material.
Collapse
Affiliation(s)
- Alexander P. Browning
- Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Brisbane, Australia
| | | | - Matthew J. Simpson
- Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| |
Collapse
|
9
|
Murphy RJ, Buenzli PR, Baker RE, Simpson MJ. A one-dimensional individual-based mechanical model of cell movement in heterogeneous tissues and its coarse-grained approximation. Proc Math Phys Eng Sci 2019; 475:20180838. [PMID: 31423086 PMCID: PMC6694308 DOI: 10.1098/rspa.2018.0838] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2018] [Accepted: 06/14/2019] [Indexed: 12/21/2022] Open
Abstract
Mechanical heterogeneity in biological tissues, in particular stiffness, can be used to distinguish between healthy and diseased states. However, it is often difficult to explore relationships between cellular-level properties and tissue-level outcomes when biological experiments are performed at a single scale only. To overcome this difficulty, we develop a multi-scale mathematical model which provides a clear framework to explore these connections across biological scales. Starting with an individual-based mechanical model of cell movement, we subsequently derive a novel coarse-grained system of partial differential equations governing the evolution of the cell density due to heterogeneous cellular properties. We demonstrate that solutions of the individual-based model converge to numerical solutions of the coarse-grained model, for both slowly-varying-in-space and rapidly-varying-in-space cellular properties. We discuss applications of the model, such as determining relative cellular-level properties and an interpretation of data from a breast cancer detection experiment.
Collapse
Affiliation(s)
- R. J. Murphy
- Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| | - P. R. Buenzli
- Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| | - R. E. Baker
- Mathematical Institute, University of Oxford, Oxford, UK
| | - M. J. Simpson
- Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| |
Collapse
|