Stochastically perturbed model of cell electropermeabilization

1. Problema e Contexto

O artigo aborda o fenômeno da eletroporação (ou eletropermeabilização reversível), que é o aumento transitório da permeabilidade da membrana celular induzido por pulsos elétricos de alta tensão e curta duração. Embora a eletroporação seja amplamente utilizada em biomedicina, biotecnologia e tratamento de câncer, a modelagem matemática existente é predominantemente determinística.

O problema central deste trabalho é a falta de modelos que incorporem efeitos aleatórios inerentes ao processo biológico e físico, tais como flutuações de temperatura ou incertezas no campo elétrico aplicado. O objetivo é estender o modelo fenomenológico determinístico proposto anteriormente por Kavian et al. [11] para um contexto estocástico, formulando um sistema acoplado de Equações Diferenciais Estocásticas Parciais (SPDEs) e Equações Diferenciais Ordinárias (ODEs).

2. Metodologia

Formulação Matemática

O modelo descreve a evolução do potencial elétrico transmembrana ($v$) e o grau de porosidade da membrana ($w$).

• Domínio: O sistema é definido na interface $\Gamma$ (membrana celular), acoplado às equações eletrostáticas nos domínios intra e extracelular.
• Equações:
  • Uma SPDE não linear para o potencial $v$, contendo um termo de ruído multiplicativo e degenerado (atuando apenas na componente da SPDE).
  • Uma ODE não linear para a porosidade $w$, que evolui dependendo de $v$ e de tempos característicos de eletroporação ($\tau_{ep}$) e selamento ($\tau_{res}$).
• Acoplamento: As equações são acopladas de forma não linear através de um termo de condutividade superficial $S_m(v, w)$ e da função de evolução de $w$, que depende de $v$. A não linearidade principal é do tipo $v \cdot w$.

Abordagem Analítica

A prova de existência e unicidade da solução enfrenta desafios significativos devido às propriedades da não linearidade:

• Desafios: As não linearidades não são Lipschitz contínuas nem monótonas globalmente, o que impede o uso de teoremas clássicos de operadores monotônicos.
• Técnica: Os autores adotam uma abordagem variacional baseada no método de Galerkin.
• Referência Teórica: Seguem a metodologia desenvolvida por Liu e Rößner [19], que generaliza as condições de monotonicidade e coercividade. Isso permite tratar coeficientes localmente monotônicos e ruído multiplicativo.
• Estratégia de Prova:
  1. Construção de aproximações de Galerkin em dimensão finita.
  2. Estabelecimento de estimativas uniformes (limites superiores) para as aproximações, utilizando fórmulas de Itô e desigualdades de Gronwall adaptadas ao contexto estocástico (com fatores de peso exponencial para lidar com a monotonicidade generalizada).
  3. Passagem ao limite para obter uma solução variacional.
  4. Prova de unicidade da solução (unicidade de trajetória) utilizando um argumento de produto e a fórmula de Itô com um fator de peso específico para cancelar o termo de crescimento da monotonicidade.

3. Principais Contribuições

1. Formulação Estocástica: Introdução de um modelo de eletroporação que incorpora ruído estocástico (aditivo e multiplicativo) diretamente na membrana celular, capturando flutuações térmicas e incertezas no campo aplicado.
2. Prova de Bem-Postura (Well-posedness): Estabelecimento rigoroso da existência e unicidade de uma solução variacional para o sistema acoplado SPDE-ODE.
   • O ruído é multiplicativo e degenerado (age apenas na equação do potencial, não diretamente na porosidade). A mistura (mixing) nas variáveis da ODE ocorre indiretamente através do acoplamento não linear no termo de deriva (drift).
3. Superação de Limitações Clássicas: Demonstração de que o sistema satisfaz as condições de monotonicidade generalizada e crescimento polinomial necessárias para a teoria de Liu e Rößner, mesmo que as não linearidades não sejam Lipschitz globais.
4. Solução Mild: Prova de que a solução variacional também é uma solução "mild" (suave) sob condições de regularidade adequadas do domínio e dos dados.

4. Resultados Numéricos

Os autores realizaram simulações numéricas utilizando o método de elementos finitos (FEniCSx) e o método de Euler semi-implícito para resolver o sistema estocástico. Foram testados dois cenários de ruído:

• Ruído Aditivo: Representa flutuações temporais globais ou incertezas na tensão aplicada.
• Ruído Multiplicativo Linear: Onde a intensidade do ruído depende do valor do potencial ($v$).

Observações das Simulações:

• Comportamento Temporal: O potencial elétrico cresce rapidamente até um pico e depois decai conforme a porosidade da membrana aumenta (efeito de selamento).
• Ergodicidade e Medida Invariante: Foram calculadas médias temporais do quadrado do potencial e da porosidade.
  • As médias temporais de trajetórias individuais convergem.
  • O desvio padrão amostral das médias temporais decai com a taxa $T^{-1/2}$ à medida que o tempo final aumenta.
• Indicação de Medida Invariante: O decaimento do desvio padrão sugere numericamente a existência de uma medida invariante (estado estacionário) para o sistema, embora a prova analítica disso seja deixada para trabalhos futuros.

5. Significado e Impacto

Este trabalho é significativo por preencher uma lacuna importante na modelagem de processos bioelétricos. Ao demonstrar a bem-postura matemática de um modelo estocástico acoplado com ruído degenerado e não linearidades complexas, o artigo fornece uma base teórica sólida para:

• Previsão de Incertezas: Permitir que pesquisadores quantifiquem a variabilidade nos resultados da eletroporação devido a fatores aleatórios.
• Otimização de Terapias: Melhorar o planejamento de tratamentos de ablação por eletroporação e entrega de genes, onde a consistência e a previsibilidade são cruciais.
• Avanço Teórico: Oferecer um exemplo de aplicação bem-sucedida da teoria de SPDEs com monotonicidade generalizada a sistemas biológicos acoplados, abrindo caminho para a análise de outros modelos complexos de interação célula-ambiente.

Em resumo, o artigo combina análise matemática rigorosa com simulação numérica para validar um modelo mais realista e robusto da eletroporação celular, considerando a natureza inerentemente estocástica dos fenômenos biológicos.