Derivation of nonlinear time-dependent macroscopic conductivity for an electropermeabilization model via homogenization

Resumen Técnico: Derivación de la Conductividad Macroscópica No Lineal Dependiente del Tiempo para un Modelo de Electropermeabilización mediante Homogeneización

Planteamiento del Problema
El artículo aborda el desafío de modelar la electroporación reversible (electropermeabilización) a escala de tejido. Mientras que los mecanismos a nivel celular de creación y cierre de poros están relativamente bien comprendidos, los modelos dinámicos para campos eléctricos a escala de tejido son escasos. Los modelos macroscópicos existentes suelen basarse en suposiciones ad hoc, como prescribir una conductividad de tipo sigmoideo dependiente de la intensidad del campo eléctrico, sin derivar este comportamiento a partir de primeros principios. El problema específico es derivar un modelo macroscópico riguroso para tejido biológico, representado como un medio periódico de células, partiendo de un modelo fenomenológico a nivel celular (Kavian et al. [29]) donde la conductividad de la membrana está acoplada con una ecuación diferencial ordinaria (EDO) para el grado de porosidad. Una dificultad clave es que la conductividad microscópica es constante, pero los datos experimentales sugieren que la conductividad efectiva del tejido es no lineal, dependiente del tiempo y exhibe efectos de memoria.

Metodología
Los autores emplean la teoría de homogeneización periódica para conectar las escalas micro y macro. La metodología procede a través de los siguientes pasos:

1. Análisis Dimensional y Escalamiento: Los autores realizan un análisis dimensional para identificar un parámetro de escalamiento relevante $\varepsilon$, definido como la relación entre el tamaño típico de la célula y el tamaño de la muestra de tejido. Este análisis revela un escalamiento de $\varepsilon^{-1}$ en la condición de transmisión dinámica a través de la membrana celular, lo cual es crucial para capturar el comportamiento asintótico correcto.
2. Formulación del Problema Microscópico: El modelo microscópico consiste en ecuaciones electrostáticas en los dominios intra y extracelulares, acopladas con una condición de transmisión dinámica en las membranas celulares oscilantes. Esta condición involucra un salto en el potencial eléctrico y un término de conductividad superficial no lineal acoplado con una EDO para el grado de porosidad $X$.
3. Bien posed de la Problematica Microscópica: Los autores demuestran la existencia y unicidad de las soluciones al problema microscópico. Un obstáculo técnico significativo es que el término de conductividad superficial no lineal no es ni monótono ni Lipschitz continuo en el sentido estándar. Para superar esto, utilizan una no linealidad modificada (truncamiento) y demuestran una condición de Lipschitz de un solo lado (monotonía generalizada) a lo largo de las soluciones, permitiendo la aplicación de argumentos de monotonía.
4. Expansiones Asintóticas Formales: Se derivan expansiones de dos escalas formales para identificar la estructura del problema límite. Esto revela un sistema acoplado donde el término de orden principal depende de un término corrector que involucra efectos de memoria.
5. Homogeneización Rigurosa: Las expansiones formales se justifican rigurosamente mediante la convergencia de dos escalas. Un desafío mayor es el paso al límite en los términos no lineales definidos en las superficies oscilantes (membranas celulares). Los autores establecen la convergencia fuerte de dos escalas para el salto en el potencial y el grado de porosidad combinando resultados de compacidad con los argumentos de monotonía mencionados anteriormente.
6. Simulación Numérica: Se desarrolla un algoritmo numérico para resolver la ecuación integro-diferencial macroscópica resultante. El algoritmo implica resolver problemas celulares auxiliares (estacionarios y dependientes del tiempo) para computar coeficientes efectivos que dependen de la historia del gradiente de potencial.

Contribuciones Clave

• Derivación Rigurosa de la Conductividad No Lineal: La contribución primaria es la derivación matemática que muestra que una conductividad efectiva no lineal y dependiente del tiempo emerge naturalmente de conductividades microscópicas constantes y un modelo de porosidad dinámica. Esto proporciona una justificación teórica para la dependencia de tipo sigmoideo de la conductividad del tejido respecto a la intensidad del campo eléctrico observada en experimentos.
• Efectos de Memoria: El modelo macroscópico derivado es una ecuación íntegro-diferencial no lineal que exhibe efectos de memoria. La densidad de corriente efectiva depende no solo del campo eléctrico instantáneo, sino también de su historia, mediada por la evolución del grado de porosidad.
• Manejo de No Linealidades No Monótonas: El artículo proporciona un marco riguroso para homogeneizar problemas con condiciones de transmisión dinámica que involucran no linealidades que no son globalmente monótonas o Lipschitz, mediante la demostración de la bien posed de la solución vía monotonía generalizada a lo largo de las trayectorias de la solución.
• Validación Cuantitativa: Los autores computan numéricamente la conductividad efectiva y demuestran que captura la evolución no monótona característica observada en experimentos: una caída inicial (debido a la carga capacitiva de la bicapa lipídica) seguida de un aumento agudo (debido a la formación de poros y el transporte iónico).

Resultos

• Modelo Macroscópico: El problema límite se formula como un sistema que involucra un potencial macroscópico $u_{-1}$ y un corrector $u_0$. La densidad de corriente efectiva viene dada por un término de divergencia que contiene una parte instantánea y una parte dependiente del tiempo que involucra una integral sobre la historia del gradiente de potencial.
• Dinámica de la Conductividad: Las simulaciones numéricas en un entorno bidimensional (tejido entre placas paralelas) confirman que la conductividad efectiva $\sigma_{eff}$ evoluciona con el tiempo. Para campos eléctricos altos, $\sigma_{eff}$ disminuye inicialmente (reflejando la fase de carga capacitiva donde los iones son atraídos hacia la membrana) y luego aumenta bruscamente a medida que se forman los poros.
• Dependencia del Campo: Las simulaciones muestran que la conductividad efectiva en estado estacionario depende de forma no lineal de la intensidad del campo eléctrico aplicado, exhibiendo una tendencia sigmoidea. Esto se alinea con las observaciones experimentales y los modelos fenomenológicos previos, pero aquí se deriva de un fundamento microscópico.
• Convergencia: El algoritmo numérico demuestra una convergencia de segundo orden con respecto al tamaño de la malla y una convergencia de primer orden con respecto a los pasos de tiempo, validando el enfoque computacional.

Significancia y Reivindicaciones
El artículo afirma proporcionar una explicación matemática rigurosa para la dinámica observada experimentalmente de la conductividad del tejido durante la electroporación. Al derivar el modelo macroscópico a partir de un modelo fenomenológico a nivel celular, los autores demuestran que el comportamiento complejo, no lineal y dependiente del tiempo de la conductividad del tejido es una consecuencia directa de la física microscópica de la carga de la membrana y la evolución de los poros. El trabajo cierra la brecha entre los mecanismos celulares y las observaciones a escala de tejido sin recurrir al ajuste empírico de la función de conductividad. Los autores enfatizan que su modelo captura la evolución no trivial de la conductividad efectiva, incluyendo la caída característica y el subsiguiente aumento, lo cual concuerda cualitativamente con los datos experimentales. La significancia radica en la justificación rigurosa de la dependencia no lineal y de forma sigmoidea de la conductividad del tejido respecto al campo eléctrico, ofreciendo una base teórica para los modelos utilizados en la estimulación eléctrica de células biológicas.