Computational Oncology of Chemotaxis-Driven Tumour--Immune Spatial Patterning and Stability

Resumen Técnico: Oncología Computacional del Patrón Espacial y la Estabilidad de la Quimiotaxis en la Interacción Tumor–Inmune

Planteamiento del Problema
La progresión de los tumores sólidos se caracteriza por una heterogeneidad espacial, que incluye patrones de infiltración inmunitaria, exclusión inmunitaria y necrosis. Si bien los modelos de reacción-difusión han descrito durante mucho tiempo el crecimiento tumoral y la dinámica inmunitaria, el papel específico de la quimiotaxis mediada por quimiocinas en la generación de patrones espacialmente heterogéneos (como la transición entre tumores "calientes" y "fríos") requiere una investigación matemática y numérica rigurosa. El problema central abordado es la formulación y el análisis de un sistema mínimo de tumor-inmune-quimiocina que acopla el crecimiento logístico tumoral, la muerte mediada por el sistema inmunitario, la producción de quimiocinas y la migración inmunitaria dirigida. El desafío radica en identificar los umbrales de estabilidad que separan la coexistencia homogénea de la formación de patrones impulsada por la quimiotaxis, asegurando al mismo tiempo que las simulaciones numéricas preserven las restricciones biológicas, como la no negatividad de las densidades celulares y los límites de tejido de flujo nulo.

Metodología
Los autores desarrollan un marco de oncología computacional basado en un sistema de reacción-difusión-quimiotaxis definido en un dominio tisular acotado $\Omega \subset \mathbb{R}^d$ con condiciones de contorno de flujo nulo.

1. Formulación Matemática:

   • Variables: El sistema rastrea la densidad de células tumorales ($T$), la densidad de células efectoras inmunitarias ($E$) y la concentración de quimiocinas ($A$).
   • Mecanismos: Las células tumorales crecen logísticamente y son eliminadas por las células inmunitarias mediante contacto ($ET$). Las quimiocinas son producidas por las células tumorales y los contactos $ET$, reclutando células inmunitarias y guiando su migración mediante quimiotaxis ($-\chi \nabla \cdot (E \nabla A)$).
   • No dimensionalización: El sistema se nondimensionaliza para reducir el espacio de parámetros, resultando en un sistema con variables $u$ (tumor), $v$ (inmune) y $w$ (quimiocina). Los parámetros adimensionales clave incluyen la sensibilidad quimiotáctica ($\xi$), el suministro inmunitario ($\sigma_0, \sigma_1$) y las tasas de decaimiento ($\delta, \ell$).

2. Esquema Numérico:

   • Discretización: Se emplea un método de volúmenes finitos conservativo. Los flujos difusivos se aproximan mediante diferencias centradas.
   • Tratamiento de la Quimiotaxis: El flujo quimiotáctico se discretiza utilizando un esquema de upwind (corriente arriba) para garantizar la robustez y la positividad en regímenes con gradientes pronunciados.
   • Integración Temporal: Se utiliza un esquema de Diferenciación de BDF (Backward Differentiation Formula) implícito (de primer o segundo orden) para manejar la rigidez derivada de los términos de difusión y reacción.
   • Diagnósticos: El esquema incorpora el monitoreo estricto de positividad, diagnósticos de balance de masa (verificando las leyes de conservación discretas) y cálculos de residuo post-hoc para distinguir la formación genuina de patrones de los artefactos numéricos.

3. Enfoque Analítico:

   • Bien posed de la solución: Los autores establecen la solvabilidad clásica local y demuestran que las soluciones permanecen no negativas. Derivan estimaciones a priori, incluyendo un límite uniforme en la densidad tumoral y estimaciones de masa global para las variables inmunitarias y de quimiocina.
   • Análisis de Equilibrio: El artículo identifica equilibrios espacialmente homogéneos, específicamente el estado libre de tumor y los estados de coexistencia.
   • Estabilidad Lineal: Se realiza un análisis de estabilidad lineal alrededor del equilibrio de coexistencia. Los autores derivan una relación de dispersión modo-dependiente para determinar cómo la sensibilidad quimiotáctica afecta la estabilidad de las perturbaciones de longitud de onda finita.

Resultados Clave

1. Umbral Libre de Tumor: El equilibrio libre de tumor ($u^*=0$) es linealmente estable si el suministro inmunitario basal supera la tasa de pérdida inmunitaria ($\sigma_0 > \delta$). Crucialmente, los autores demuestran que la quimiotaxis no influye en este umbral específico; la estabilidad del estado libre de tumor está determinada únicamente por el equilibrio entre el reclutamiento inmunitario y el decaimiento.
2. Equilibrios de Coexistencia: La existencia de un estado de coexistencia espacialmente homogéneo se reduce a encontrar la raíz de una ecuación escalar $F(u)=0$. El artículo establece las condiciones bajo las cuales existe un equilibrio de coexistencia positivo (por ejemplo, cuando $\sigma_0 < \delta$ y ciertos parámetros de reclutamiento son positivos).
3. Inestabilidad Impulsada por la Quimiotaxis: El análisis de estabilidad lineal del estado de coexistencia revela que la quimiotaxis actúa como un acoplamiento amplificado por el número de onda. Mientras que el modo espacialmente uniforme ($k=0$) permanece inalterado por la quimiotaxis, los modos de longitud de onda finita pueden volverse inestables si la sensibilidad quimiotáctica $\xi$ excede un valor crítico $\xi_c$. Esto conduce a una inestabilidad de tipo Turing donde emergen patrones espacialmente heterogéneos a partir de un estado homogéneo.
4. Verificación Numérica: El esquema de volúmenes finitos propuesto verifica con éxito los umbrales teóricos. Los experimentos numéricos confirman los modos inestables dominantes predichos, los mapas de sensibilidad, y la preservación de la positividad y el balance de masa. Los diagnósticos de residuo confirman que los patrones observados son consistentes con el sistema de EDP (ecuaciones diferenciales parciales).

Significancia y Reivindicaciones
El artículo afirma contribuir a la oncología computacional de tres maneras específicas:

1. Formulación del Modelo: Formula un modelo de reacción-difusión-quimiotaxis mínimo y analíticamente tratable que acopla explícitamente la muerte tumoral mediada por el sistema inmunitario, el reclutamiento mediado por quimiocinas y la migración dirigida.
2. Umbrales Analíticos: Deriva umbrales interpretables para el control libre de tumor y la inestabilidad de longitud de onda finita impulsada por la quimiotaxis, vinculando los mecanismos de retroalimentación biológica directamente con los criterios de estabilidad y formación de patrones.
3. Fiabilidad Numérica: Implementa un esquema de volúmenes finitos conservativo con características específicas (flujo quimiotáctico de upwind, tratamiento de contorno de flujo nulo y diagnósticos rigurosos) diseñadas para preservar las restricciones biológicas de no negatividad y conservación de masa.

Los autores posicionan este trabajo como un puente entre el análisis matemático y la simulación numérica, proporcionando un marco para investigar cómo el reclutamiento mediado por quimiocinas y la migración quimiotáctica moldean la organización espacial de las interacciones tumor-inmune. El estudio se presenta como una "caricatura" de la transición biológica entre la coexistencia homogénea y la heterogeneidad impulsada por la quimiotaxis (el díctico de tumores "calientes/fríos"), ofreciendo una base mecánica para comprender los patrones de exclusión e infiltración inmunitaria sin pretender modelar resultados clínicos específicos o proponer nuevas terapias experimentales.