Stochastically perturbed model of cell electropermeabilization

Titolo: Modello stocasticamente perturbato dell'elettropermeabilizzazione cellulare

Autori: T. Gebäck, O. Misiats, I. Motschan-Armen, I. Pettersson.
Affiliazioni: Chalmers University of Technology & Gothenburg University (Svezia), Virginia Commonwealth University (USA).

---

1. Il Problema e il Contesto

L'elettroporazione (o elettropermeabilizzazione) è un fenomeno transitorio in cui la permeabilità della membrana cellulare aumenta temporaneamente in risposta a brevi impulsi elettrici ad alto voltaggio. Questa tecnica è ampiamente utilizzata in biomedicina, biotecnologia e nella terapia oncologica (ablazione non termica).

I modelli deterministici esistenti (es. [11]) descrivono il fenomeno accoppiando le equazioni elettrostatiche per il potenziale elettrico nei domini intra- ed extracellulare con un'equazione differenziale ordinaria (ODE) non lineare che governa il grado di porosità della membrana. Tuttavia, tali modelli non tengono conto delle fluttuazioni casuali intrinseche, come le variazioni di temperatura o l'incertezza nel campo elettrico applicato.

Obiettivo del lavoro: Introdurre e analizzare matematicamente una versione stocastica del modello deterministico proposto in [11], aggiungendo rumore al sistema per modellare questi effetti casuali, e dimostrare l'esistenza e l'unicità della soluzione per il sistema risultante.

2. Metodologia e Formulazione Matematica

Modello Deterministico di Base

Il modello originale descrive l'evoluzione del salto di potenziale transmembrana $v = [u]$ e del grado di porosità $w$ su una membrana $\Gamma$. Il sistema è composto da:

1. Un'equazione di evoluzione parabolica per $v$, accoppiata con un operatore non locale (mappa di Dirichlet-to-Neumann, $L$) che deriva dalle equazioni di Laplace nel volume cellulare.
2. Un'ODE per $w$, che descrive la cinetica di apertura e chiusura dei pori (con tempi caratteristici $\tau_{ep}$ ed $\tau_{res}$).

Introduzione del Rumore

Gli autori introducono una perturbazione stocastica moltiplicativa e degenerata. Il sistema stocastico (SPDE-ODE accoppiato) è formulato come:
$$\begin{aligned} c_m dv &= -(Lv + S_m(v, w) + \sigma_i \nabla p \cdot n) dt + b(v) dW_t \\ dw &= f(v, w) dt \end{aligned}$$
dove:

• $W_t$ è un processo di Wiener Q-valore in uno spazio di Hilbert.
• $b(v)$ è un operatore di rumore che agisce solo sulla componente SPDE (il potenziale $v$), rendendo il rumore degenerato rispetto al sistema completo.
• Il rumore si propaga alla variabile ODE ($w$) solo indirettamente attraverso l'accoppiamento non lineare nel termine di deriva.

Sfide Tecniche

La principale difficoltà risiede nella natura della non linearità di accoppiamento ($v \cdot w$), che non è né Lipschitziana né monotona in senso classico. Questo impedisce l'applicazione diretta delle teorie standard per le equazioni stocastiche differenziali (SDE/SPDE) basate su operatori monotoni.

Strumenti Matematici Adottati

Per superare queste difficoltà, gli autori utilizzano l'approccio variazionale generalizzato sviluppato da Liu, Röckner e Prévôt:

1. Condizioni Generalizzate: Si dimostrano proprietà di monotonia generalizzata locale e coercività generalizzata per l'operatore non lineare $A(U)$, valide lungo le soluzioni del sistema (e non per tutte le funzioni dello spazio funzionale).
2. Metodo di Galerkin: Si costruiscono approssimazioni finite dimensionali del sistema stocastico.
3. Stime Uniformi: Si derivano stime a priori uniformi rispetto alla dimensione dell'approssimazione, utilizzando il calcolo di Itô e le disuguaglianze di Gronwall adattate al caso stocastico (con fattori esponenziali per gestire la crescita non lineare).
4. Passaggio al Limite: Si dimostra la convergenza delle approssimazioni di Galerkin verso una soluzione variazionale unica.

3. Risultati Principali

1. Esistenza e Unicità (Teorema 3.2):
   È stato dimostrato che il sistema accoppiato SPDE-ODE ammette una soluzione variazionale unica. La soluzione appartiene agli spazi $L^2(\Omega; L^2(0, T; V)) \cap L^2(\Omega; L^\infty(0, T; H))$, garantendo regolarità sufficiente per l'analisi fisica.

   • L'unicità è ottenuta combinando la formula di Itô con una regola del prodotto che introduce un fattore di smorzamento esponenziale, permettendo di gestire la costante di Lipschitz che dipende dalla soluzione stessa.

2. Soluzione "Mild":
   Sotto ipotesi di regolarità aggiuntive sul dominio e sui dati, la soluzione variazionale è dimostrata essere anche una soluzione "mild" (analiticamente debole), collegando il problema alla teoria dei semigruppi di contrazione.

3. Simulazioni Numeriche:
   Gli autori hanno implementato il modello utilizzando il metodo degli elementi finiti (FEM) con il pacchetto FEniCSx, considerando due tipi di rumore:

   • Rumore Additivo: Fluttuazioni globali del potenziale o della tensione applicata.
   • Rumore Moltiplicativo Lineare: Il rumore scala con il valore del potenziale ($b(v) \propto v$).

   Osservazioni dalle simulazioni:

   • Il potenziale di membrana cresce rapidamente, raggiunge un picco e poi decade a causa dell'aumento della porosità ($w$).
   • Sono stati calcolati i medi temporali di $\|v\|^2$ e $\|w\|^2$.
   • La deviazione standard campionaria dei medi temporali decade come $T^{-1/2}$ all'aumentare del tempo finale. Questo comportamento è un forte indicatore numerico dell'esistenza di una misura invariante (ergodicità del sistema), suggerendo che il sistema raggiunge uno stato stazionario statistico.

4. Contributi Chiave

• Generalizzazione del Modello: Estensione di un modello deterministico di elettroporazione al caso stocastico, fondamentale per modellare la variabilità biologica e sperimentale.
• Teoria Matematica Avanzata: Dimostrazione di ben-postezza per un sistema accoppiato SPDE-ODE con rumore degenerato (agente su una sola componente) e non linearità non monotone. Questo riempie un vuoto nella letteratura sulle equazioni stocastiche con condizioni di monotonia generalizzate.
• Analisi Numerica: Fornitura di simulazioni che validano il modello e offrono evidenze numeriche per l'ergodicità, un passo preliminare necessario per lo studio analitico della misura invariante.

5. Significato e Implicazioni

Questo lavoro fornisce un fondamento teorico rigoroso per l'uso di modelli stocastici nell'elettroporazione. La capacità di gestire rumore degenerato e non linearità complesse apre la strada a simulazioni più realistiche dei processi biologici, dove l'incertezza è onnipresente.
La dimostrazione dell'esistenza di una misura invariante (suggerita dai risultati numerici) è cruciale per comprendere il comportamento a lungo termine delle cellule sottoposte a stimolazione elettrica, con potenziali applicazioni nell'ottimizzazione delle terapie tumorali e nella somministrazione di farmaci. Il metodo sviluppato può essere esteso ad altri sistemi biologici accoppiati governati da dinamiche non lineari e stocastiche.