Cellular Potts model

Last updated

In computational biology, a Cellular Potts model (CPM, also known as the Glazier-Graner-Hogeweg model) is a computational model of cells and tissues. It is used to simulate individual and collective cell behavior, tissue morphogenesis and cancer development. CPM describes cells as deformable objects with a certain volume, that can adhere to each other and to the medium in which they live. The formalism can be extended to include cell behaviours such as cell migration, growth and division, and cell signalling. The first CPM was proposed for the simulation of cell sorting by François Graner and James A. Glazier as a modification of a large-Q Potts model. [1] CPM was then popularized by Paulien Hogeweg for studying morphogenesis. [2] Although the model was developed to describe biological cells, it can also be used to model individual parts of a biological cell, or even regions of fluid.

Contents

Model description

Cartoon example of lattice used in Cellular Potts Model. Cpm layout.pdf
Cartoon example of lattice used in Cellular Potts Model.

The CPM consists of a rectangular Euclidean lattice, where each cell is a subset of lattice sites sharing the same cell ID (analogous to spin in Potts models in physics). Lattice sites that are not occupied by cells are the medium. The dynamics of the model are governed by an energy function: the Hamiltonian which describes the energy of a particular configuration of cells in the lattice. In a basic CPM, this energy results from adhesion between cells and resistance of cells to volume changes. The algorithm for updating CPM minimizes this energy.

In order to evolve the model Metropolis-style updates are performed, that is:

  1. choose a random lattice site i
  2. choose a random neighboring lattice site j to copy its ID into i.
  3. calculate the difference in energy () between the original and the proposed new configuration.
  4. accept or reject this copy event based on the change in energy , as follows:
    if the new energy is lower, always accept the copy;
    if the new energy is higher, accept the copy with probability (the Boltzmann temperature T determines the likelihood of energetically unfavorable fluctuations).

The Hamiltonian

The original model proposed by Graner and Glazier contains cells of two types, with different adhesion energies for cells of the same type and cells of a different type. Each cell type also has a different contact energy with the medium, and the cell volume is assumed to remain close to a target value. The Hamiltonian is formulated as:

where i, j are lattice sites, σi is the cell at site i, τ(σ) is the cell type of cell σ, J is the coefficient determining the adhesion between two cells of types τ(σ),τ(σ'), δ is the Kronecker delta, v(σ) is the volume of cell σ, V(σ) is the target volume, and λ is a Lagrange multiplier determining the strength of the volume constraint.

Cells with a lower J value for their membrane contact will stick together more strongly. Therefore, different patterns of cell sorting can be simulated by varying the J values.

Extensions

Over time, the CPM has evolved from a specific model of cell sorting to a general framework with many extensions, some of which are partially or entirely off-lattice. [3] Various cell behaviours, such as chemotaxis, elongation and haptotaxis can be incorporated by extending either the Hamiltonian, H, or the change in energy . Auxiliary sub-lattices may be used to include additional spatial information, such as the concentrations of chemicals.

Chemotaxis

In CPM, cells can be made to move in the direction of higher chemokine concentration, by increasing the probability of copying the ID of site j into site i when the chemokine concentration is higher at j. This is done by modifying the change in energy with a term that is proportional to the difference in concentration at i and j: [2]

Where is the strength of chemotactic movement, and and are the concentration of the chemokine at site i and j, respectively. The chemokine gradient is typically implemented on a separate lattice of the same dimensions as the cell lattice.

Multiscale and hybrid modeling using CPM

Core GGH (or CPM) algorithm which defines the evolution of the cellular level structures can easily be integrated with intracellular signaling dynamics, reaction diffusion dynamics and rule based model to account for the processes which happen at lower (or higher) time scale. [4] Open source software Bionetsolver can be used to integrate intracellular dynamics with CPM algorithm. [5]

Related Research Articles

In a chemical reaction, chemical equilibrium is the state in which both the reactants and products are present in concentrations which have no further tendency to change with time, so that there is no observable change in the properties of the system. This state results when the forward reaction proceeds at the same rate as the reverse reaction. The reaction rates of the forward and backward reactions are generally not zero, but they are equal. Thus, there are no net changes in the concentrations of the reactants and products. Such a state is known as dynamic equilibrium.

In game theory, the Nash equilibrium, named after the mathematician John Nash, is the most common way to define the solution of a non-cooperative game involving two or more players. In a Nash equilibrium, each player is assumed to know the equilibrium strategies of the other players, and no one has anything to gain by changing only one's own strategy. The principle of Nash equilibrium dates back to the time of Cournot, who in 1838 applied it to competing firms choosing outputs.

<span class="mw-page-title-main">Lattice model (physics)</span>

In mathematical physics, a lattice model is a mathematical model of a physical system that is defined on a lattice, as opposed to a continuum, such as the continuum of space or spacetime. Lattice models originally occurred in the context of condensed matter physics, where the atoms of a crystal automatically form a lattice. Currently, lattice models are quite popular in theoretical physics, for many reasons. Some models are exactly solvable, and thus offer insight into physics beyond what can be learned from perturbation theory. Lattice models are also ideal for study by the methods of computational physics, as the discretization of any continuum model automatically turns it into a lattice model. The exact solution to many of these models includes the presence of solitons. Techniques for solving these include the inverse scattering transform and the method of Lax pairs, the Yang–Baxter equation and quantum groups. The solution of these models has given insights into the nature of phase transitions, magnetization and scaling behaviour, as well as insights into the nature of quantum field theory. Physical lattice models frequently occur as an approximation to a continuum theory, either to give an ultraviolet cutoff to the theory to prevent divergences or to perform numerical computations. An example of a continuum theory that is widely studied by lattice models is the QCD lattice model, a discretization of quantum chromodynamics. However, digital physics considers nature fundamentally discrete at the Planck scale, which imposes upper limit to the density of information, aka Holographic principle. More generally, lattice gauge theory and lattice field theory are areas of study. Lattice models are also used to simulate the structure and dynamics of polymers.

The Ising model, named after the physicists Ernst Ising and Wilhelm Lenz, is a mathematical model of ferromagnetism in statistical mechanics. The model consists of discrete variables that represent magnetic dipole moments of atomic "spins" that can be in one of two states. The spins are arranged in a graph, usually a lattice, allowing each spin to interact with its neighbors. Neighboring spins that agree have a lower energy than those that disagree; the system tends to the lowest energy but heat disturbs this tendency, thus creating the possibility of different structural phases. The model allows the identification of phase transitions as a simplified model of reality. The two-dimensional square-lattice Ising model is one of the simplest statistical models to show a phase transition.

The equilibrium constant of a chemical reaction is the value of its reaction quotient at chemical equilibrium, a state approached by a dynamic chemical system after sufficient time has elapsed at which its composition has no measurable tendency towards further change. For a given set of reaction conditions, the equilibrium constant is independent of the initial analytical concentrations of the reactant and product species in the mixture. Thus, given the initial composition of a system, known equilibrium constant values can be used to determine the composition of the system at equilibrium. However, reaction parameters like temperature, solvent, and ionic strength may all influence the value of the equilibrium constant.

In statistical mechanics, the Potts model, a generalization of the Ising model, is a model of interacting spins on a crystalline lattice. By studying the Potts model, one may gain insight into the behaviour of ferromagnets and certain other phenomena of solid-state physics. The strength of the Potts model is not so much that it models these physical systems well; it is rather that the one-dimensional case is exactly solvable, and that it has a rich mathematical formulation that has been studied extensively.

In statistical mechanics, the n-vector model or O(n) model is a simple system of interacting spins on a crystalline lattice. It was developed by H. Eugene Stanley as a generalization of the Ising model, XY model and Heisenberg model. In the n-vector model, n-component unit-length classical spins are placed on the vertices of a d-dimensional lattice. The Hamiltonian of the n-vector model is given by:

<span class="mw-page-title-main">Hubbard model</span> Approximate model used to describe the transition between conducting and insulating systems

The Hubbard model is an approximate model used to describe the transition between conducting and insulating systems. It is particularly useful in solid-state physics. The model is named for John Hubbard.

In statistical mechanics, the two-dimensional square lattice Ising model is a simple lattice model of interacting magnetic spins. The model is notable for having nontrivial interactions, yet having an analytical solution. The model was solved by Lars Onsager for the special case that the external magnetic field H = 0. An analytical solution for the general case for has yet to be found.

<span class="mw-page-title-main">Lattice Boltzmann methods</span> Class of computational fluid dynamics methods

The lattice Boltzmann methods (LBM), originated from the lattice gas automata (LGA) method (Hardy-Pomeau-Pazzis and Frisch-Hasslacher-Pomeau models), is a class of computational fluid dynamics (CFD) methods for fluid simulation. Instead of solving the Navier–Stokes equations directly, a fluid density on a lattice is simulated with streaming and collision (relaxation) processes. The method is versatile as the model fluid can straightforwardly be made to mimic common fluid behaviour like vapour/liquid coexistence, and so fluid systems such as liquid droplets can be simulated. Also, fluids in complex environments such as porous media can be straightforwardly simulated, whereas with complex boundaries other CFD methods can be hard to work with.

<span class="mw-page-title-main">CompuCell3D</span>

CompuCell3D (CC3D) is a three-dimensional C++ and Python software problem solving environment for simulations of biocomplexity problems, integrating multiple mathematical [morphogenesis] models. These include the cellular Potts model (CPM) or Glazier-Graner-Hogeweg model (GGH) which can model cell clustering, growth, division, death, adhesion, and volume and surface area constraints; as well as partial differential equation solvers for modeling reaction–diffusion of external chemical fields and cell type automata for differentiation. By integrating these models CompuCell3D enables modeling of cellular reactions to external chemical fields such as secretion or resorption, and responses such as chemotaxis and haptotaxis.

Paulien Hogeweg is a Dutch theoretical biologist and complex systems researcher studying biological systems as dynamic information processing systems at many interconnected levels. In 1970, together with Ben Hesper, she defined the term bioinformatics as "the study of informatic processes in biotic systems".

In statistical mechanics, the corner transfer matrix describes the effect of adding a quadrant to a lattice. Introduced by Rodney Baxter in 1968 as an extension of the Kramers-Wannier row-to-row transfer matrix, it provides a powerful method of studying lattice models. Calculations with corner transfer matrices led Baxter to the exact solution of the hard hexagon model in 1980.

The Swendsen–Wang algorithm is the first non-local or cluster algorithm for Monte Carlo simulation for large systems near criticality. It has been introduced by Robert Swendsen and Jian-Sheng Wang in 1987 at Carnegie Mellon.

Dynamical mean-field theory (DMFT) is a method to determine the electronic structure of strongly correlated materials. In such materials, the approximation of independent electrons, which is used in density functional theory and usual band structure calculations, breaks down. Dynamical mean-field theory, a non-perturbative treatment of local interactions between electrons, bridges the gap between the nearly free electron gas limit and the atomic limit of condensed-matter physics.

In statistical mechanics, the Griffiths inequality, sometimes also called Griffiths–Kelly–Sherman inequality or GKS inequality, named after Robert B. Griffiths, is a correlation inequality for ferromagnetic spin systems. Informally, it says that in ferromagnetic spin systems, if the 'a-priori distribution' of the spin is invariant under spin flipping, the correlation of any monomial of the spins is non-negative; and the two point correlation of two monomial of the spins is non-negative.

<span class="mw-page-title-main">HPP model</span>

The Hardy–Pomeau–Pazzis (HPP) model is a fundamental lattice gas automaton for the simulation of gases and liquids. It was a precursor to the lattice Boltzmann methods. From lattice gas automata, it is possible to derive the macroscopic Navier-Stokes equations. Interest in lattice gas automaton methods levelled off in the early 1990s, due to rising interest in the lattice Boltzmann methods.

In statistical mechanics, probability theory, graph theory, etc. the random cluster model is a random graph that generalizes and unifies the Ising model, Potts model, and percolation model. It is used to study random combinatorial structures, electrical networks, etc. It is also referred to as the RC model or sometimes the FK representation after its founders Cees Fortuin and Piet Kasteleyn.

The three-state Potts CFT, also known as the parafermion CFT, is a conformal field theory in two dimensions. It is a minimal model with central charge . It is considered to be the simplest minimal model with a non-diagonal partition function in Virasoro characters, as well as the simplest non-trivial CFT with the W-algebra as a symmetry.

James Alexander Glazier is a biophysicist and bioengineer, author, and educator best known for his contributions to the field of multiscale modeling, pattern formation, and morphogenesis in biological systems. Glazier has published numerous articles in leading scientific journals, and his work has been widely recognized within the scientific community. He has also been influential in promoting the use of computational modeling and simulation in the study of complex biological phenomena.

References

  1. Graner, François; Glazier, James (1992). "Simulation of biological cell sorting using a two-dimensional extended Potts model". Phys. Rev. Lett. 69 (13): 2013–7. Bibcode:1992PhRvL..69.2013G. doi:10.1103/PhysRevLett.69.2013. PMID   10046374.
  2. 1 2 Savill, Nicholas J.; Hogeweg, Paulien (1997). "Modelling Morphogenesis: From Single Cells to Crawling Slugs". J. Theor. Biol. 184 (3): 229–235. Bibcode:1997JThBi.184..229S. doi:10.1006/jtbi.1996.0237. hdl:1874/1405. PMID   31940735. S2CID   1611700.
  3. Balter, Ariel; Merks, Roeland M.H.; Popławski, Nikodem J.; Swat, Maciej; Glazier, James A. (2007). "The Glazier-Graner-Hogeweg model: extensions, future directions, and opportunities for further study". Single-Cell-Based Models in Biology and Medicine. Mathematics and Biosciences in Interaction. pp. 151–167. doi:10.1007/978-3-7643-8123-3_7. ISBN   978-3-7643-8101-1.
  4. Szabó, A; Merks, RM (2013). "Cellular potts modeling of tumor growth, tumor invasion, and tumor evolution". Frontiers in Oncology. 3: 87. doi: 10.3389/fonc.2013.00087 . PMC   3627127 . PMID   23596570.
  5. Andasari, Vivi; Roper, Ryan T; Swat, Maciej H; Chaplain, MA (2012). "Integrating intracellular dynamics using CompuCell3D and Bionetsolver: applications to multiscale modelling of cancer cell growth and invasion". PLOS ONE. 7 (3): e33726. Bibcode:2012PLoSO...733726A. doi: 10.1371/journal.pone.0033726 . PMC   3312894 . PMID   22461894.