Travelling Waves in a Mathematical Model for Oncolytic Virotherapy

A continuación presento un resumen técnico detallado del artículo "Travelling Waves in a Mathematical Model for Oncolytic Virotherapy" (Olasas viajeras en un modelo matemático para la viroterapia oncolítica), escrito por Negar Mohammadnejad y Thomas Hillen.

1. Planteamiento del Problema

La viroterapia oncolítica (OVT) es una estrategia prometedora que utiliza virus genéticamente modificados para infectar y destruir selectivamente células tumorales. Un desafío central en este campo es determinar si el virus, administrado como monoterapia, puede propagarse a través del tumor y erradicarlo antes de que las células cancerosas invadan el tejido circundante.

El problema matemático abordado en este trabajo es el análisis de la existencia de soluciones de ondas viajeras en un modelo de reacción-difusión no cooperativo de tres componentes. El modelo describe la interacción espacial entre:

• $C$: Células cancerosas no infectadas.
• $I$: Células cancerosas infectadas.
• $V$: Concentración de viriones libres (virus).

El sistema es no cooperativo (no monótono), lo que significa que la matriz Jacobiana de las cinéticas contiene entradas no diagonales negativas. Esto impide la aplicación de la teoría estándar de ondas viajeras para sistemas cooperativos o monótonos, que es común en la literatura. El objetivo es demostrar rigurosamente la existencia de ondas viajeras positivas que conectan el estado de tumor no tratado con un estado de invasión viral, y determinar la velocidad mínima de propagación ($\bar{c}$) necesaria para que estas ondas existan.

2. Metodología

Los autores emplean un enfoque analítico riguroso basado en la teoría de ecuaciones diferenciales parciales (EDP) y análisis funcional:

1. Modelo y No-dimensionalización:

   • Se parte del modelo de Baabdulla et al. [2] y se aplica una no-dimensionalización para simplificar los parámetros.
   • Se introduce un cambio de variable $B = 1 - C$ para trasladar el equilibrio del tumor no tratado (donde $C=1$) al origen $(0,0,0)$, facilitando el análisis de la invasión hacia un estado "libre de enfermedad".

2. Formulación de la Onda Viajera:

   • Se asume una solución de la forma $(B, I, V)(x, t) = (B, I, V)(\xi)$ con $\xi = x + ct$.
   • Esto transforma el sistema de EDPs en un sistema de ecuaciones diferenciales ordinarias (EDOs) de segundo orden no lineales.

3. Análisis de la Relación de Dispersión:

   • Se linealiza el sistema alrededor del equilibrio trivial $(0,0,0)$ en el frente de la onda.
   • Se derivan las relaciones de dispersión para encontrar la velocidad de la onda $S(\rho)$ en función de la tasa de decaimiento espacial $\rho$.
   • Se identifica la velocidad mínima teórica $c_m = \inf_{\rho>0} S(\rho)$.

4. Construcción de Soluciones Superiores e Inferiores (Upper and Lower Solutions):

   • Dado que el sistema no es cooperativo, no se puede usar el teorema de comparación estándar. En su lugar, los autores construyen explícitamente funciones de solución superior ($\bar{B}, \bar{I}, \bar{V}$) y solución inferior ($\underline{B}, \underline{I}, \underline{V}$).
   • Estas funciones se definen mediante combinaciones de exponenciales y constantes, ajustadas para satisfacer desigualdades diferenciales específicas que "encuadran" la solución real.
   • Se imponen condiciones estrictas sobre los parámetros (como $a$, $\gamma$, $\theta$) para asegurar que las componentes de los vectores propios asociados sean positivas, lo cual es crucial para la validez de las soluciones.

5. Teorema del Punto Fijo de Schauder:

   • Se define un operador integral $T$ basado en la variación de parámetros.
   • Se construye un conjunto $\Gamma$ en un espacio de Banach ponderado, acotado por las soluciones superior e inferior.
   • Se demuestra que el operador $T$ es continuo, compacto y deja invariante el conjunto $\Gamma$.
   • La aplicación del Teorema del Punto Fijo de Schauder garantiza la existencia de una solución fija, que corresponde a la onda viajera buscada.

3. Contribuciones Clave

• Extensión de la teoría a sistemas no cooperativos: El trabajo supera las limitaciones de la teoría de sistemas monótonos, aplicando técnicas de soluciones superior/inferior a un modelo de viroterapia de tres componentes con estructura no monótona.
• Determinación de la velocidad mínima crítica ($\bar{c}$): Se establece rigurosamente que las ondas viajeras positivas existen para todas las velocidades $c \geq \bar{c}$. El valor de $\bar{c}$ no es simplemente el mínimo de la curva de dispersión lineal ($c_m$), sino que depende de las intersecciones de múltiples curvas ($S(\rho)$, $H_1(\rho)$, $H_3(\rho)$, $H_4(\rho)$) que surgen de las condiciones de las soluciones superior e inferior.
• Identificación de "Regímenes Nebulosos" (Foggy Regimes): El análisis revela regiones paramétricas donde la existencia de ondas viajeras no puede determinarse concluyentemente con los métodos actuales. Esto surge cuando las curvas de velocidad se intersectan de manera compleja, dejando un intervalo de incertidumbre entre la velocidad mínima teórica y la velocidad numérica observada.

4. Resultados Principales

• Teorema de Existencia (Teorema 3.6): Bajo ciertas condiciones de los parámetros biológicos (específicamente $\theta > a\gamma$ y restricciones en la relación entre las curvas de velocidad), el sistema admite soluciones de ondas viajeras no negativas y no triviales para cualquier velocidad $c \geq \bar{c}$.
• Comportamiento Asintótico: Las soluciones conectan el estado de tumor puro $(1, 0, 0)$ en $\xi \to -\infty$ con un estado acotado (que puede ser un equilibrio de coexistencia o una órbita periódica si hay bifurcación de Hopf) en $\xi \to \infty$.
• Validación Numérica:
  • Se realizaron simulaciones numéricas en MATLAB utilizando el solver pdepe.
  • Los resultados numéricos muestran que las velocidades de onda observadas ($c_N$) coinciden estrechamente con el mínimo de la curva de dispersión lineal ($c_m$) en la mayoría de los casos, aunque a veces difieren ligeramente de la velocidad teórica estricta $\bar{c}$ debido a las restricciones de las soluciones superior/inferior.
  • Se observó la formación de frentes de invasión y dinámicas oscilatorias detrás del frente, dependiendo de la magnitud del parámetro de estallido viral ($\theta$).

5. Significado e Impacto

• Fundamento Teórico: Este trabajo proporciona una base matemática rigurosa para entender la propagación espacial de virus oncolíticos, llenando un vacío en la literatura que se había centrado principalmente en simulaciones numéricas o modelos cooperativos.
• Implicaciones Clínicas: La identificación de la velocidad mínima de propagación es crucial para el diseño de tratamientos. Sugiere que para que la viroterapia sea efectiva como monoterapia, los parámetros virales (como la tasa de infección o el tamaño del estallido) deben ser lo suficientemente altos para superar la velocidad de invasión del tumor.
• Nuevas Direcciones de Investigación: La identificación de las "regiones nebulosas" donde la existencia de ondas es ambigua plantea preguntas matemáticas abiertas. Sugiere que se necesitan nuevas técnicas analíticas o refinamientos en la construcción de soluciones superior/inferior para cerrar estas brechas y predecir con mayor precisión los umbrales de éxito terapéutico en condiciones biológicas complejas.

En resumen, el artículo combina análisis matemático avanzado con simulaciones numéricas para demostrar que, bajo condiciones específicas, los virus oncolíticos pueden propagarse como ondas viajeras estables a través de un tumor, ofreciendo una herramienta teórica para optimizar estrategias de tratamiento contra el cáncer.