Articulo de referencia

Modelo de Potts celular

En biología computacional , un modelo celular de Potts (CPM, también conocido como modelo de Glazier-Graner-Hogeweg) es un modelo computacional de células y tejidos. Se utiliza ...

En biología computacional , un modelo celular de Potts (CPM, también conocido como modelo de Glazier-Graner-Hogeweg) es un modelo computacional de células y tejidos. Se utiliza para simular el comportamiento celular individual y colectivo, la morfogénesis tisular y el desarrollo del cáncer . El CPM describe las células como objetos deformables con un cierto volumen, que pueden adherirse entre sí y al medio en el que viven. El formalismo puede extenderse para incluir comportamientos celulares como la migración , el crecimiento y la división celular , y la señalización celular . El primer CPM fue propuesto para la simulación de la clasificación celular por François Graner y James A. Glazier como una modificación de un modelo de Potts de Q grande . [ 1 ] El CPM fue popularizado posteriormente por Paulien Hogeweg para el estudio de la morfogénesis. [ 2 ] Aunque el modelo se desarrolló para describir células biológicas , también puede utilizarse para modelar partes individuales de una célula biológica, o incluso regiones de fluido.

Descripción del modelo

Ejemplo esquemático de una estructura reticular utilizada en el modelo celular de Potts.

El CPM consiste en una red euclidiana rectangular , donde cada celda es un subconjunto de sitios de la red que comparten el mismo ID de celda (análogo al espín en los modelos de Potts en física). Los sitios de la red que no están ocupados por celdas son el medio. La dinámica del modelo está regida por una función de energía: el hamiltoniano que describe la energía de una configuración particular de celdas en la red. En un CPM básico, esta energía resulta de la adhesión entre celdas y la resistencia de las celdas a los cambios de volumen. El algoritmo para actualizar el CPM minimiza esta energía .

Para evolucionar el modelo se realizan actualizaciones al estilo de Metrópolis , es decir:

  1. elija un sitio de red aleatorio i
  2. elige un sitio de red vecino aleatorio j para copiar su ID en i .
  3. calcular la diferencia de energía (ΔH{\displaystyle \Delta H}) entre la configuración original y la nueva configuración propuesta.
  4. aceptar o rechazar este evento de copia en función del cambio en la energíaΔH{\displaystyle \Delta H}, de la siguiente manera:
    Si la nueva energía es menor, siempre acepte la copia;
    Si la nueva energía es mayor, acepta la copia con probabilidadmiΔH/T{\displaystyle e^{-\Delta H/T}}(La temperatura de Boltzmann T determina la probabilidad de fluctuaciones energéticamente desfavorables).

El hamiltoniano

El modelo original propuesto por Graner y Glazier contiene células de dos tipos, con diferentes energías de adhesión para células del mismo tipo y células de distinto tipo. Cada tipo de célula también tiene una energía de contacto diferente con el medio, y se supone que el volumen celular se mantiene cercano a un valor objetivo. El hamiltoniano se formula como:

H=i,j vecinosJ(τ(σi),τ(σj))(1δ(σi,σj))+λσi(v(σi)V(σi))2,{\displaystyle {\begin{aligned}H=\sum _{i,j{\text{ vecinos}}}J\left(\tau (\sigma _{i}),\tau (\sigma _{j})\right)\left(1-\delta (\sigma _{i},\sigma _{j})\right)+\lambda \sum _{\sigma _{i}}\left(v(\sigma _{i})-V(\sigma _{i})\right)^{2},\\\end{aligned}}}

donde i , j son sitios de la red, σ i es la celda en el sitio i, τ(σ) es el tipo de celda de la celda σ, J es el coeficiente que determina la adhesión entre dos celdas de tipos τ(σ),τ(σ'), δ es la delta de Kronecker , v(σ) es el volumen de la celda σ, V(σ) es el volumen objetivo y λ es un multiplicador de Lagrange que determina la fuerza de la restricción de volumen.

Las células con un valor J más bajo para su contacto de membrana se adhieren con mayor fuerza entre sí. Por lo tanto, se pueden simular diferentes patrones de clasificación celular variando los valores J.

Extensiones

Con el tiempo, el CPM ha evolucionado desde un modelo específico de clasificación celular hasta un marco general con muchas extensiones, algunas de las cuales están parcial o totalmente fuera de la red. [ 3 ] Se pueden incorporar diversos comportamientos celulares, como la quimiotaxis , la elongación y la haptotaxis , extendiendo el hamiltoniano, H, o el cambio en la energía.ΔH{\displaystyle \Delta H}Se pueden utilizar subredes auxiliares para incluir información espacial adicional, como las concentraciones de sustancias químicas.

Quimiotaxis

En CPM, se puede hacer que las células se muevan en la dirección de una mayor concentración de quimiocinas , aumentando la probabilidad de copiar la ID del sitio j en el sitio i cuando la concentración de quimiocinas es mayor en j . Esto se hace modificando el cambio en la energía.ΔH{\displaystyle \Delta H}con un término que es proporcional a la diferencia de concentración en i y j : [ 2 ]

ΔH=ΔHμ(doidoj){\displaystyle {\begin{aligned}\Delta H'=\Delta H-\mu (C_{i}-C_{j})\\\end{aligned}}}

Dóndeμ{\displaystyle \mu }es la fuerza del movimiento quimiotáctico, ydoi{\displaystyle C_{i}}ydoj{\displaystyle C_{j}}son las concentraciones de la quimiocina en los sitios i y j, respectivamente. El gradiente de quimiocina se implementa típicamente en una red separada de las mismas dimensiones que la red celular.

Modelado multiescala e híbrido mediante CPM

El algoritmo GGH (o CPM) central, que define la evolución de las estructuras a nivel celular, puede integrarse fácilmente con la dinámica de señalización intracelular, la dinámica de reacción-difusión y el modelo basado en reglas para tener en cuenta los procesos que ocurren en escalas de tiempo más bajas (o más altas). [ 4 ] El software de código abierto Bionetsolver puede utilizarse para integrar la dinámica intracelular con el algoritmo CPM. [ 5 ]

Aplicaciones del modelo celular de Potts

El modelo celular de Potts (CPM), también conocido como modelo de Glazier-Graner-Hogeweg (GGH), se utiliza ampliamente en biología computacional para simular sistemas multicelulares. Permite modelar la forma, el movimiento, la adhesión, el crecimiento y las interacciones celulares con el entorno, lo que lo hace idóneo para el estudio de la morfogénesis tisular, el cáncer, la clasificación celular y la formación de patrones. A continuación, se presenta una tabla que resume las principales aplicaciones biológicas del CPM, junto con referencias representativas.

Marcos de trabajo que implementan el modelo celular de Potts

Los siguientes marcos de trabajo permiten a los usuarios aplicar el CPM a la construcción de modelos multicelulares.

Referencias

  1. Graner, François; Glazier, James (1992). "Simulación de la clasificación de células biológicas mediante un modelo de Potts extendido bidimensional" . 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). "Modelado de la morfogénesis: de células individuales a babosas reptantes" . 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 MH; Popławski, Nikodem J.; Swat, Maciej; Glazier, James A. (2007). "El modelo de Glazier-Graner-Hogeweg: extensiones, direcciones futuras y oportunidades para estudios posteriores" . Modelos basados ​​en células individuales en biología y medicina . Matemáticas y biociencias en interacción. pp. 151–167 . doi : 10.1007/978-3-7643-8123-3_7 . ISBN  978-3-7643-8101-1.
  4. Szabó, A; Merks, RM (2013). " Modelado celular de Potts del crecimiento tumoral, la invasión tumoral y la evolución tumoral" . 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). "Integración de la dinámica intracelular mediante CompuCell3D y Bionetsolver: aplicaciones al modelado multiescala del crecimiento e invasión de células cancerosas" . PLOS ONE . 7 (3) e33726. Bibcode : 2012PLoSO...733726A . doi : 10.1371/journal.pone.0033726 . PMC 3312894. PMID 22461894 .  
  • Chen, Nan; Glazier, James A.; Izaguirre, Jesus A.; Alber, Mark S. (2007). "Una implementación paralela del modelo celular de Potts para la simulación de la morfogénesis celular" . Computer Physics Communications . 176 ( 11–12 ): 670–681 . Bibcode : 2007CoPhC.176..670C . doi : 10.1016/j.cpc.2007.03.007 . PMC 2139985. PMID 18084624 .  
  • James Glazier (sitio web profesional)
  • CompuCell3D , un entorno de simulación CPM: Sourceforge
  • SimTK
  • Terreno urbanístico de Notre Dame
  • Modelo de Vida Artificial de morfogénesis multicelular con gradientes generados de forma autónoma para información posicional utilizando el modelo de Potts celular.
  • Autómatas celulares estocásticos