Travelling Waves in a Mathematical Model for Oncolytic Virotherapy

Título: Ondas de Travessia em um Modelo Matemático para Viroterapia Oncolítica

Autores: Negar Mohammadnejad e Thomas Hillen
Instituição: Universidade de Alberta, Canadá

1. Problema e Contexto

A viroterapia oncolítica (OVT) é uma estratégia promissora de tratamento de câncer que utiliza vírus geneticamente modificados para infectar e destruir seletivamente células tumorais. Um desafio central nesta área é entender se o vírus, quando administrado como monoterapia, consegue se espalhar através do tumor e erradicá-lo antes que o próprio tumor invada os tecidos circundantes.

O problema matemático abordado neste trabalho é a análise de um sistema de equações de reação-difusão não cooperativo de três componentes, que descreve a dinâmica espacial entre:

1. Células cancerígenas não infectadas ($C$);
2. Células cancerígenas infectadas ($I$);
3. Vírus livres (virions, $V$).

O modelo em questão, baseado em trabalhos anteriores de Baabdulla et al., é não cooperativo (não monotônico), o que impede a aplicação direta de teorias clássicas de ondas de travessia desenvolvidas para sistemas cooperativos ou monotônicos. O objetivo principal é estabelecer rigorosamente a existência de soluções de ondas de travessia (travelling waves) que representem a invasão do vírus no tecido tumoral e determinar a velocidade mínima de propagação.

2. Metodologia

Os autores utilizam uma abordagem analítica rigorosa baseada na teoria de equações diferenciais parciais (EDP) e análise funcional:

• Não Dimensionalização: O sistema original foi transformado em um sistema adimensional para simplificar a análise, mantendo a estrutura biológica essencial.
• Transformação de Variáveis: Para estudar a invasão do vírus no estado de "câncer não tratado", os autores realizaram uma mudança de variáveis ($B = 1 - C$), deslocando o equilíbrio correspondente ao tumor puro para a origem $(0,0,0)$.
• Formulação de Ondas de Travessia: O problema foi convertido em um sistema de Equações Diferenciais Ordinárias (EDOs) não lineares de segunda ordem, assumindo uma solução da forma $(B(\xi), I(\xi), V(\xi))$ onde $\xi = x + ct$.
• Construção de Soluções Superiores e Inferiores (Upper and Lower Solutions): Devido à natureza não cooperativa do sistema, os autores construíram manualmente soluções superiores e inferiores que satisfazem as condições de fronteira. Isso envolveu uma análise detalhada das relações de dispersão (eigenvalores) na borda da onda (frente de invasão).
• Teorema do Ponto Fixo de Schauder: Uma vez estabelecidas as soluções superior e inferior, os autores definiram um operador integral $T$ em um espaço de Banach apropriado. Eles demonstraram que este operador é contínuo, compacto e mapeia um conjunto convexo, fechado e limitado em si mesmo. Pelo Teorema do Ponto Fixo de Schauder, a existência de um ponto fixo (que corresponde à solução da onda de travessia) é garantida.
• Análise de Velocidade Mínima: A velocidade mínima $\bar{c}$ foi determinada analisando a curva de dispersão $S(\rho)$ derivada da linearização do sistema na frente da onda.

3. Contribuições Chave

• Existência Rigorosa: O trabalho fornece a primeira prova matemática rigorosa da existência de ondas de travessia positivas para este modelo específico de três componentes não cooperativo.
• Determinação da Velocidade Mínima: Os autores identificam um valor crítico de velocidade mínima $\bar{c}$ tal que ondas de travessia positivas existem para todas as velocidades $c \geq \bar{c}$.
• Identificação de Regimes "Neblinosos" (Foggy Regimes): O método de construção das soluções superior e inferior revela que, em certas regiões paramétricas (dependendo de parâmetros como a taxa de infecção, taxa de lise e difusão), a existência de ondas não pode ser concluída com a técnica atual. O artigo destaca essas zonas de incerteza como novas questões matemáticas a serem resolvidas.
• Análise de Estabilidade Linear: Uma análise detalhada dos autovalores da matriz Jacobiana na frente da onda foi realizada para garantir que as soluções construídas sejam biologicamente relevantes (não oscilatórias na frente da onda).

4. Resultados Principais

• Teorema Principal (Teorema 3.6): Sob condições específicas nos parâmetros (principalmente relacionadas à relação entre a taxa de liberação viral $\theta$ e a taxa de morte celular $\gamma$), o sistema admite soluções de ondas de travessia que conectam o estado de tumor não infectado $(1, 0, 0)$ a um estado de coexistência ou oscilação, permanecendo limitado no infinito.
• Velocidade Crítica: A velocidade mínima $\bar{c}$ é definida como o mínimo da curva de velocidade $S(\rho)$ sobre um intervalo admissível de taxas de decaimento $\rho$.
• Simulações Numéricas: As simulações realizadas em MATLAB (usando o solver pdepe) confirmam a formação de frentes de invasão. A comparação entre as velocidades numéricas observadas e a velocidade teórica mínima $\bar{c}$ mostra uma boa concordância na maioria dos casos, embora discrepâncias ocorram nas regiões paramétricas onde a prova analítica é menos restritiva (os regimes "neblinosos").
• Dinâmica Oscilatória: O modelo captura transições para dinâmicas oscilatórias (órbitas periódicas) atrás da frente de onda quando o parâmetro de bifurcação $\theta$ ultrapassa um certo limiar (bifurcação de Hopf), o que é consistente com observações experimentais e modelos anteriores.

5. Significado e Impacto

• Validação Teórica: O trabalho valida matematicamente a capacidade do vírus de se propagar como uma onda de invasão em tumores, fornecendo uma base teórica sólida para estratégias de OVT.
• Otimização de Parâmetros: Ao identificar os parâmetros que influenciam a velocidade mínima da onda, o estudo sugere quais características virais (como infectividade ou tamanho do broto viral) devem ser otimizadas para garantir que o vírus se espalhe mais rápido do que o tumor cresce.
• Desafios Futuros: A identificação de regimes onde a existência de ondas permanece ambígua ("foggy regimes") aponta para a necessidade de desenvolver novas técnicas analíticas para sistemas não cooperativos mais complexos.
• Relevância Clínica: Os resultados reforçam a ideia de que a OVT como monoterapia pode não ser suficiente para erradicação total em todos os cenários, destacando a importância de entender a dinâmica espacial e a necessidade de terapias combinadas (ex: com imunoterapia) para superar as limitações de propagação viral.

Em resumo, o artigo é uma contribuição significativa para a matemática aplicada na oncologia, unindo análise funcional avançada com modelagem biológica para elucidar os mecanismos de propagação de vírus oncolíticos.