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

Sintesi Tecnica: Derivazione della Conducibilità Macroscopica Non Lineare e Tempo-Dipendente per un Modello di Elettropermeabilizzazione tramite Omogeneizzazione

Enunciato del Problema
Il documento affronta la sfida di modellare l'elettroporazione reversibile (elettropermeabilizzazione) alla scala tissutale. Mentre i meccanismi cellulari di creazione e chiusura dei pori sono relativamente ben compresi, i modelli dinamici per i campi elettrici su scala tissutale sono scarsi. I modelli macroscopici esistenti spesso si affidano ad assunzioni ad hoc, come la prescrizione di una conducibilità di tipo sigmoideo dipendente dall'intensità del campo elettrico, senza derivare tale comportamento da primi principi. Il problema specifico è quello di derivare un rigoroso modello macroscopico per il tessuto biologico, rappresentato come un mezzo periodico di cellule, partendo da un modello fenomenologico a livello cellulare (Kavian et al. [29]) in cui la conducibilità di membrana è accoppiata con un'equazione differenziale ordinaria (ODE) per il grado di porosità. Una difficoltà chiave risiede nel fatto che la conducibilità microscopica è costante, eppure i dati sperimentali suggeriscono che la conducibilità efficace del tessuto sia non lineare, tempo-dipendente e presenti effetti di memoria.

Metodologia
Gli autori utilizzano la teoria dell'omogeneizzazione periodica per colmare il divario tra le scale micro e macro. La metodologia procede attraverso i seguenti passaggi:

1. Analisi Dimensionale e Scaling: Gli autori eseguono un'analisi dimensionale per identificare un parametro di scaling pertinente $\varepsilon$, definito come il rapporto tra la dimensione tipica della cellula e la dimensione del campione di tessuto. Questa analisi rivela uno scaling $\varepsilon^{-1}$ nella condizione di trasmissione dinamica attraverso la membrana cellulare, che è cruciale per catturare il corretto comportamento asintotico.
2. Formulazione del Problema Microscopico: Il modello microscopico consiste in equazioni elettrostatiche nei domini intra ed extracellulare, accoppiate con una condizione di trasmissione dinamica sulle membrane cellulari oscillanti. Questa condizione prevede un salto nel potenziale elettrico e un termine di conducibilità superficiale non lineare accoppiato con un'ODE per il grado di porosità $X$.
3. Ben-posedness del Problema Microscopico: Gli autori dimostrano l'esistenza e l'unicità delle soluzioni al problema microscopico. Un ostacolo tecnico significativo è che il termine di conducibilità superficiale non lineare non è né monotono né Lipschitz continuo nel senso standard. Per superare questo, utilizzano una non linearità modificata (troncamento) e dimostrano una condizione di Lipschitz a lato singolo (monotonia generalizzata) lungo le soluzioni, consentendo l'applicazione di argomenti di monotonia.
4. Espansioni Asintotiche Formali: Vengono derivate espansioni a due scale formali per identificare la struttura del problema limite. Ciò rivela un sistema accoppiato in cui il termine di ordine principale dipende da un termine correttore che coinvolge effetti di memoria.
5. Omogeneizzazione Rigorosa: Le espansioni formali sono giustificate rigorosamente utilizzando la convergenza a due scale. Una sfida principale è il passaggio al limite nei termini non lineari definiti sulle superfici oscillanti (membrane cellulari). Gli autori stabiliscono la forte convergenza a due scale per il salto di potenziale e il grado di porosità combinando risultati di compattezza con gli argomenti di monotonia menzionati sopra.
6. Simulazione Numerica: Viene sviluppato un algoritmo numerico per risolvere la risultante equazione integro-differenziale macroscopica. L'algoritmo prevede la risoluzione di problemi cellulari ausiliari (stazionari e tempo-dipendenti) per calcolare i coefficienti efficaci che dipendono dalla storia del gradiente di potenziale.

Contributi Chiave

• Derivazione Rigorosa della Conducibilità Non Lineare: Il contributo primario è la derivazione matematica che mostra come una conducibilità efficace non lineare e tempo-dipendente emerga naturalmente da conducibilità microscopiche costanti e da un modello di porosità dinamica. Ciò fornisce una giustificazione teorica per la dipendenza sigmoidea della conducibilità tissutale dall'intensità del campo elettrico osservata negli esperimenti.
• Effetti di Memoria: Il modello macroscopico derivato è un'equazione integro-differenziale non lineare che esibisce effetti di memoria. La densità di corrente efficace dipende non solo dal campo elettrico istantaneo, ma anche dalla sua storia, mediata dall'evoluzione del grado di porosità.
• Gestione di Non Linearità Non Monotone: Il lavoro fornisce un quadro rigoroso per omogeneizzare problemi con condizioni di trasmissione dinamiche che coinvolgono non linearità che non sono globalmente monotone o Lipschitz continue, dimostrando la ben-posedness tramite monotonia generalizzata lungo le traiettorie delle soluzioni.
• Validazione Quantitativa: Gli autori calcolano numericamente la conducibilità efficace e dimostrano che essa cattura l'evoluzione non monotona caratteristica osservata negli esperimenti: un calo iniziale (dovuto alla carica capacitiva del doppio strato lipidico) seguito da un aumento rapido (dovuto alla formazione di pori e al trasporto ionico).

Risultati

• Modello Macroscopico: Il problema limite è formulato come un sistema che coinvolge un potenziale macroscopico $u_{-1}$ e un correttore $u_0$. La densità di corrente efficace è data da un termine di divergenza contenente una parte istantanea e una parte tempo-dipendente che coinvolge un integrale sulla storia del gradiente di potenziale.
• Dinamica della Conducibilità: Le simulazioni numeriche in un ambiente bidimensionale (tessuto tra piastre parallele) confermano che la conducibilità efficace $\sigma_{eff}$ evolve nel tempo. Per campi elettrici elevati, $\sigma_{eff}$ inizialmente diminuisce (riflettendo la fase di carica capacitiva in cui gli ioni vengono attratti verso la membrana) e poi aumenta bruscamente con la formazione dei pori.
• Dipendenza dal Campo: Le simulazioni mostrano che la conducibilità efficace a regime dipende in modo non lineare dall'intensità del campo elettrico applicato, esibendo un trend sigmoideo. Ciò si allinea con le osservazioni sperimentali e con i precedenti modelli fenomenologici, ma è qui derivato da una base microscopica.
• Convergenza: L'algoritmo numerico dimostra una convergenza del secondo ordine rispetto alla dimensione della mesh e una convergenza del primo ordine rispetto ai passi temporali, validando l'approccio computazionale.

Significatività e Rivendicazioni
Il lavoro sostiene di fornire una spiegazione matematica rigorosa per la dinamica della conducibilità tissutale durante l'elettroporazione osservata sperimentalmente. Derivando il modello macroscopico da un modello fenomenologico a livello cellulare, gli autori dimostrano che il comportamento complesso, non lineare e tempo-dipendente della conducibilità tissutale è una diretta conseguenza della fisica microscopica della carica della membrana e dell'evoluzione dei pori. Il lavoro colma il divario tra i meccanismi cellulari e le osservazioni a scala tissutale senza ricorrere all'adattamento empirico della funzione di conducibilità. Gli autori sottolineano che il loro modello cattura l'evoluzione non banale della conducibilità efficace, incluso il calo caratteristico e il successivo aumento, il che concorda qualitativamente con i dati sperimentali. La significatività risiede nella rigorosa giustificazione della dipendenza non lineare e sigmoidea della conducibilità tissutale dall'intensità del campo elettrico, offrendo una base teorica per i modelli utilizzati nella stimolazione elettrica di cellule biologiche.