Articolo metodologico

Un approccio di modellazione computazionale per studiare l'influenza dell'ipertermia sul microambiente tumorale

DOI:

10.3791/65870

⸱

1 dicembre 2023

 ,  , 
1Department of Electrical and Computer Engineering, Kansas State University, 2Department of Interventional Radiology, The University of Texas MD Anderson Cancer Center

In questo articolo

Sommario

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

L'articolo descrive un protocollo per simulare i profili di temperatura transitori e la variazione spazio-temporale accoppiata della pressione del fluido interstiziale a seguito del riscaldamento erogato da un sistema di ipertermia a radiofrequenza dipolare. Il protocollo può essere utilizzato per valutare la risposta dei parametri biofisici che caratterizzano il microambiente tumorale alle tecniche di ipertermia interventistica.

Abstract

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Le proprietà biofisiche del microambiente tumorale differiscono sostanzialmente dai tessuti normali. Una costellazione di caratteristiche, tra cui la diminuzione della vascolarizzazione, la mancanza di drenaggio linfatico e l'elevata pressione interstiziale, riduce la penetrazione delle terapie nei tumori. L'ipertermia locale all'interno del tumore può alterare le proprietà microambientali, come la pressione del fluido interstiziale, portando potenzialmente a miglioramenti nella penetrazione del farmaco. In questo contesto, i modelli computazionali multifisici possono fornire informazioni sull'interazione tra i parametri biofisici all'interno del microambiente tumorale e possono guidare la progettazione e l'interpretazione di esperimenti che testano i bioeffetti dell'ipertermia locale.

Questo articolo descrive un flusso di lavoro passo-passo per un modello computazionale che accoppia equazioni differenziali alle derivate parziali che descrivono la distribuzione della corrente elettrica, il trasferimento di biocalore e la dinamica dei fluidi. L'obiettivo principale è quello di studiare gli effetti dell'ipertermia erogata da un dispositivo bipolare a radiofrequenza sulla pressione del fluido interstiziale all'interno del tumore. Viene presentato il sistema di espressioni matematiche che collegano la distribuzione della corrente elettrica, il trasferimento di biocalore e la pressione del fluido interstiziale, sottolineando i cambiamenti nella distribuzione della pressione del fluido interstiziale che potrebbero essere indotti dall'intervento termico.

Introduzione

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

L'elevata pressione del liquido interstiziale (IFP) è un segno distintivo dei tumori solidi1. La fuoriuscita di liquido nell'interstizio dai vasi sanguigni iperpermeabili è sbilanciata dall'uscita di liquido a causa della compressione delle vene intratumorali e dell'assenza di linfatici 1,2,3. Di concerto con altri parametri biofisici anormali all'interno del microambiente tumorale (TME), tra cui lo stress solido e la rigidità, un elevato IFP mina l'efficacia della somministrazione sistemica e locale di farmaci 4,5,6. La pressione del liquido interstiziale nei tumori solidi varia da 5 mmHg (glioblastoma e melanoma) a 30 mmHg (carcinoma a cellule renali) rispetto a 1-3 mmHg nel tessuto normale2. L'IFP elevato è responsabile dell'aumento del flusso di fluido verso il margine del tumore ed espone le cellule stromali, le cellule infiltrate e altri componenti extracellulari allo stress di taglio 1,4. Le alterazioni meccano-biologiche sostengono un TME immunosoppressivo, ad esempio, aumentando la germinazione endoteliale, che supporta l'angiogenesi, la migrazione e l'invasione delle cellule tumorali, l'espressione del fattore di crescita trasformante β (TGF-β) e l'irrigidimento stromale 7,8,9.

Diversi studi hanno esplorato le terapie basate sull'energia con l'intento di ridurre l'IFP, tra cui gli ultrasuoni a bassa intensità, gli ultrasuoni focalizzati ad alta intensità, i campi elettrici pulsati e le terapie termali 5,10,11. È stato dimostrato che il riscaldamento a temperature comprese tra 40 e 43 °C, indicato come ipertermia lieve, aumenta la perfusione sanguigna tumorale e può quindi contribuire all'espansione delle vene compresse e alla riduzione della pressione vascolare facilitando l'intravaso e il drenaggio del liquido interstiziale11,12. Alcuni studi recenti hanno dimostrato il potenziale dell'ipertermia per ridurre l'IFP e, di conseguenza, per facilitare la distribuzione di farmaci o mezzi di contrasto all'interno di un tumore13,14. Questi studi mostrano anche un aumento dell'infiltrazione delle cellule T dopo ipertermia rispetto ai gruppi di controllo senza trattamento13.

I risultati promettenti degli esperimenti in vivo su piccoli animali motivano ulteriori studi che impiegano approcci computazionali per far progredire la comprensione di come i parametri fisici all'interno del TME siano influenzati dagli interventi fisici 4,15,16,17. I risultati dei modelli computazionali possono integrare gli studi sperimentali in vivo per scoprire la relazione causa-effetto alla base del riscaldamento locale (o di altre fonti di energia esterne) e dell'IFP. Ciò può essere particolarmente istruttivo date le sfide con la misurazione delle variazioni spaziali nell'IFP con trasduttori di pressione basati su catetere e ago, che in genere forniscono misurazioni puntuali 9,16,18,19. Nel contesto della somministrazione di farmaci, la comprensione dei principali meccanismi biofisici è essenziale per definire il protocollo di riscaldamento appropriato e la finestra temporale per l'iniezione del farmaco per aumentare la probabilità di un'efficace distribuzione del farmaco. Le informazioni quantitative in termini di cambiamenti nelle caratteristiche biofisiche del TME, incluso ma non limitato all'IFP, potrebbero anche fornire informazioni sull'interpretazione della risposta immunologica (ad esempio, l'infiltrazione delle cellule T) agli stimoli esterni.

Presentiamo un protocollo per la modellazione computazionale delle modifiche mediate termicamente ai profili IFP tumorali. In particolare, il protocollo descrive in dettaglio come modellare un apparato personalizzato per piccoli animali per la somministrazione di terapia termica controllata con corrente a radiofrequenza, simulare profili di temperatura transitori dopo il riscaldamento e accoppiare simulazioni fluidodinamiche per calcolare la variazione spazio-temporale dell'IFP tumorale in risposta alla terapia termica. Questo modello rispecchia le caratteristiche essenziali della configurazione sperimentale che abbiamo utilizzato in un modello di tumore sottocutaneo (McArdle RH7777, ATCC) in un precedente studio sperimentale20.

La Figura 1 mostra il modello computazionale che abbiamo implementato per calcolare le variazioni indotte termicamente dell'IFP in un tumore circondato da tessuto normale. Un paio di aghi ipodermici inseriti nel tumore sono modellati per fornire riscaldamento con corrente a radiofrequenza a 500 kHz. Nel dominio tumorale si assume un materiale poroso, composto da due fasi: la fase solida rappresenta la matrice extracellulare solida e la fase fluida rappresenta il liquido interstiziale. Nel caso di una variazione di pressione o di una deformazione della matrice derivante da uno stimolo esterno, ad esempio un aumento della temperatura, i componenti solidi e fluidi si riorganizzano. Ciò provoca il movimento del liquido interstiziale attraverso la matrice solida extracellulare 16,17,21.

Dalla teoria della poroelasticità, il tensore di sollecitazione S (Pa) (equazione [1]) è la combinazione del termine elastico che descrive la variazione di volume della componente solida rispetto alle condizioni iniziali, e di un termine poroso che descrive la sollecitazione indotta dalla pressione idrostatica della componente fluida.

figure-introduction-1(1)

Dove, λ, μ (Pa) sono i parametri di Lamé, E è il tensore della deformazione, e è il tensore della deformazione volumetrica, Pi (Pa) è la pressione del fluido interstiziale (I è la matrice identità). Si assumono condizioni di stato stazionario per la componente solida sotto sollecitazione poroelastica, il che significa che le componenti del tensore di sollecitazione sono ortogonali, figure-introduction-2.

La Figura 2 mostra il sistema di equazioni matematiche implementato nel modello poroelastico descritto e l'interazione tra i componenti del modello multifisico presentato. Il flusso di lavoro delle simulazioni computazionali comprende:

Equazioni dei problemi elettrici. La soluzione delle equazioni del problema elettrico fornisce la fonte di calore RF mediata nel tempo Q (riscaldamento Joule). A tal fine, un'approssimazione quasi statica delle equazioni di Maxwell viene utilizzata per calcolare la distribuzione del campo elettrico medio nel tempo E (V/m) (Figura 2, blocco 1).

Equazioni dei problemi termici. La soluzione dell'equazione del biocalore di Pennes (Figura 2, blocco 2) fornisce la variazione spaziale e temporale della temperatura T (°C) come risultato della fonte di calore (Q) legata all'energia elettromagnetica assorbita, del riscaldamento passivo legato alla conduzione termica dei tessuti (figure-introduction-3), e dell'effetto dissipatore di calore della perfusione sanguigna tissutale (cWb(T) (T - Tb)). Il termine dissipatore di calore approssima lo scambio di calore tra il sangue che scorre nel microvascolare e il tessuto adiacente dove viene assorbita l'energia elettromagnetica. L'equazione del trasferimento di calore include anche il termine di avvezione (figure-introduction-4), che descrive la variazione della temperatura causata dal movimento del fluido interstiziale attraverso la matrice extracellulare del modello poroelastico. Tuttavia, questo termine ha un impatto trascurabile sul profilo di temperatura rispetto agli altri meccanismi responsabili del cambiamento di temperatura.

Equazioni di problemi fluidodinamiche. L'equazione di conservazione della massa (Figura 2, blocco 3) combinata con la legge di Darcy (Figura 2, blocco 4) fornisce come output la variazione spaziale e temporale della pressione del fluido interstiziale Pirisultante dall'equilibrio tra la sorgente (figure-introduction-5) e il pozzo (figure-introduction-6 ) del fluido. Il termine di pressione transitoria sul lato sinistro dell'equazione di conservazione della massa, figure-introduction-7 , descrive il riarrangiamento dei componenti fluidi e solidi nel materiale poroelastico. Ciò è causato dalla variazione della pressione del fluido interstiziale, Pi, guidata dalla variazione della pressione vascolare Pvin funzione della temperatura.

La differenza tra la pressione vascolare (Pv) e la pressione del fluido interstiziale (Pi) è la fonte del fluido che scorre attraverso la matrice extracellulare. Il termine sink è legato alla differenza di pressione tra i vasi linfatici (PL) e lo spazio interstiziale (Pi). Nei tessuti normali, la pressione nel sistema vascolare linfatico (~ -6-0 mmHg) è fino a due volte inferiore alla pressione del fluido interstiziale13. Questa differenza di pressione garantisce l'efficacia dei vasi linfatici nel drenare l'eccesso di liquido che fuoriesce dalla parete dei vasi sanguigni nell'interstizio. Per il modello tumorale qui presentato, abbiamo trascurato il contributo del sistema linfatico 4,16,22.

Le espressioni matematiche dalle equazioni da (2) a (5) sono utilizzate per descrivere la dipendenza dalla temperatura della conducibilità elettrica e termica dei tessuti e della perfusione sanguigna dei tessuti23,24. Due diversi modelli matematici sono utilizzati per descrivere la dipendenza dalla temperatura della perfusione sanguigna nel dominio del tessuto normale e in quello del tessuto tumorale, rispettivamente24,25. I modelli mostrano che la perfusione sanguigna aumenta con la temperatura fino a nove volte rispetto al basale nel tessuto normale e solo circa due volte il valore basale nel dominio tumorale. Per entrambi i modelli, l'aumento della perfusione sanguigna è limitato alle temperature all'interno dell'intervallo di ipertermia lieve (inferiore a 45 °C). Vale la pena ricordare che le espressioni matematiche, le equazioni (4) e (5), non descrivono completamente i meccanismi alla base delle variazioni dipendenti dalla temperatura nella perfusione sanguigna nei due diversi tipi di tessuto. Tuttavia, aiutano a rappresentare la limitata perfusione che tipicamente caratterizza il microambiente tumorale rispetto ai tessuti normali.

figure-introduction-8(2)

figure-introduction-9(3)

figure-introduction-10 (4)

figure-introduction-11 (5)

figure-introduction-12(6)

figure-introduction-13(7)

In questo studio, abbiamo utilizzato le equazioni (6) e (7) per modellare la pressione vascolare in funzione della perfusione sanguigna sia per i modelli di tessuto normale che per quelli di tessuto tumorale26. Dalle equazioni (4) e (5), la velocità del flusso sanguigno può essere espressa come il rapporto tra la perfusione sanguigna e la densità del sangue. La relazione tra flusso sanguigno e pressione vascolare è ben stabilita in letteratura3: la velocità del flusso sanguigno e la resistenza geometrica (o conducibilità, Lp) del sistema vascolare determinano la differenza di pressione all'interno del vaso sanguigno. La pressione vascolare può essere espressa in funzione della temperatura (Equazioni (6) e (7)), facendo leva su questa relazione e sul modello temperatura-dipendente della perfusione sanguigna (Equazioni (4) e (5)).

L'implementazione del flusso di lavoro computazionale (Figura 2) e le proprietà dipendenti dalla temperatura dei modelli tissutali sono descritte in dettaglio nella sezione seguente. Tutte le proprietà dei materiali e le loro descrizioni e valori di base (cioè a temperatura corporea) sono elencati nella Tabella 1. Consulta la Tabella dei materiali per i dettagli su COMSOL Multiphysics installato su un computer utilizzato per implementare questo protocollo computazionale. Il problema elettrico è stato modellato utilizzando il modulo AC/DC; il trasferimento di biocalore è stato modellato utilizzando la fisica del trasferimento di calore; e il problema della fluidodinamica è stato modellato utilizzando l'interfaccia Matematica.

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Protocollo

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

1. Costruire il modello di un sistema bipolare a radiofrequenza

  1. Passaggi preliminari per impostare l'interfaccia
    1. Avviate COMSOL Multiphysics e cliccate su Model Wizard.
    2. Selezionare 3D come dimensione spaziale.
    3. Seleziona il modulo di fisica AC/DC | Campi e correnti elettriche | Correnti elettriche.
    4. Seleziona il modulo di trasferimento del calore | Trasferimento di calore nei solidi.
    5. Seleziona il modulo Matematica | Interfacce PDE | Forma a coefficiente PDE.
    6. Seleziona Studio | Dipendente dal tempo. Fare clic su Fine.
    7. Una volta visualizzata l'area di lavoro Comsol:
      1. Seleziona Multifisica | Riscaldamento elettromagnetico. Con questo passaggio, la densità di perdita di potenza elettromagnetica viene automaticamente accoppiata come fonte di calore per l'equazione di trasferimento del biocalore.
        NOTA: Se la multifisica non viene visualizzata automaticamente, specificare manualmente la fonte di calore elettromagnetica (mostrata in COMSOL come densità di perdita volumetrica). Per maggiori dettagli su come aggiungere la fonte di calore, vedere la sezione "Fisica", passaggio 2 "Configurazione per il problema termico".
      2. Seleziona Studio dalla barra multifunzione superiore | Fasi di studio | Transitorio di frequenza.
  2. Definite le geometrie. Dalla barra multifunzione superiore, selezionare Geometria, quindi:
    1. Definite due coni con le quote elencate nella Tabella 2.
    2. Posizionare i coni alla distanza indicata nella Tabella 2 (distanziati tra loro). Questi due coni modelleranno i due aghi ipodermici utilizzati per costruire il sistema RF bipolare.
    3. Duplicare i due coni precedenti per modellare l'isolamento degli aghi; modificare la dimensione del cono in base alle dimensioni riportate in Tabella 2.
    4. Selezionate un cilindro (altezza, hm e diametro dm) per modellare la massa del muscolo posizionato a z = - 9 mm (x = 0, y = 0). I valori di ciascuna dimensione sono elencati nella Tabella 2.
    5. Selezionare un cilindro (altezza, hs e diametro ds) per modellare lo strato sottile di pelle posizionato a z = 4 mm (x = 0, y = 0). I valori di ciascuna dimensione sono elencati nella Tabella 2.
    6. Selezionare una sfera (diametro, dt) per modellare il tumore sottocutaneo posizionato a z = -0,5 mm (x = 0, y = 0). La dimensione del tumore è elencata nella Tabella 2.
    7. Per facilitare la selezione delle geometrie nelle fasi successive del protocollo, si consiglia quanto segue:
      1. Dalla barra multifunzione Geometria , selezionare Operazioni virtuali | Formare domini compositi.
      2. Selezionare tutti i domini relativi alla parte elettricamente conduttiva degli aghi per creare una geometria composita.
      3. Ripetere la stessa procedura per creare domini compositi per le geometrie di isolamento dell'ago.
  3. Definire le proprietà dei modelli biologici tissutali.
    NOTA: I passaggi seguenti descrivono la procedura per implementare le espressioni matematiche descritte dalle equazioni (2)-(7).
    1. Dal nodo Componente, fare clic con il pulsante destro del mouse per selezionare Definizioni.
    2. In Funzioni, selezionare Analitica.
      1. Specificare il nome della funzione (ad esempio, k_muscle o sigma_muscle) e digitare l'espressione matematica coerente con l'Eq. 2).
      2. Specificare la temperatura (T) come argomento.
      3. Specificare le Unità della funzione: S/m nel caso di conducibilità elettrica.
      4. Ripetere i passaggi precedenti da 1 a 3 per implementare l'Eq. 3, modificando l'unità di conseguenza (ad esempio, W/(m·K) per la conducibilità termica).
      5. Specificare l'unità per l'argomento: K (kelvin) per la temperatura. In Parametri grafico, specificare l'intervallo di valori dell'argomento della funzione (ad esempio, temperatura). Per seguire questo protocollo, utilizzare un intervallo di 33-100 °C (306,15-373,15 K).
      6. Ripetere i passaggi precedenti da 1 a 5 per aggiungere le funzioni dipendenti dalla temperatura delle conducibilità elettriche (Eq. 2) e termiche (Eq. 3) per ciascun modello di tessuto (ad esempio, muscolo, pelle e tumore) utilizzando i valori nominali elencati nella Tabella 1 (il tessuto normale si riferisce sia al muscolo che alla pelle).
    3. In Funzioni, selezionare A tratti per implementare le equazioni (4)-(7):
      1. Specificare il nome della funzione.
      2. Specificare la temperatura (T) come argomento della funzione.
      3. Digitare l'espressione matematica per ogni intervallo di temperatura coerente con le equazioni (4)-(7).
      4. Ripetere i passaggi precedenti da 1 a 3 per aggiungere le funzioni dipendenti dalla temperatura della perfusione sanguigna e della pressione vascolare per ciascun modello di tessuto utilizzando i valori nominali elencati nella Tabella 1 (il tessuto normale si riferisce sia al muscolo che alla pelle).
  4. Assegnare le proprietà del materiale ai componenti geometrici.
    1. Dal nodo del componente, selezionare Materiali.
    2. Seleziona i materiali bianchi per includere tessuto normale, tessuto tumorale, sangue, PTFE e acciaio inossidabile.
    3. Abilitare la selezione manuale e selezionare l'entità geometrica corrispondente al materiale specificato.
      1. Il tessuto normale è associato a geometrie che modellano il muscolo e la pelle.
      2. Il tumore e i tessuti sanguigni sono associati alla geometria del tumore.
      3. Il materiale PTFE è associato alle geometrie che modellano l'isolante dell'ago.
      4. Il materiale in acciaio inossidabile è associato alle geometrie dei coni che modellano il terreno e gli aghi attivi.
    4. Per le conducibilità elettriche e termiche dipendenti dalla temperatura 23, digitare il nome scelto della funzione e l'argomento correlato (ad esempio, T) che appare nel nodo Definizioni.
    5. Per le proprietà del materiale che non dipendono dalla temperatura, fare riferimento ai valori di base27 elencati nella Tabella 1.
      NOTA: Ci affidiamo alla teoria poroelastica per calcolare la pressione 16,17,26. I passaggi seguenti mostrano come le proprietà di un materiale poroso possono essere assegnate a un dominio specifico.
    6. Da Materiali, selezionare Altri materiali | Materiale poroso.
    7. Fare clic con il pulsante destro del mouse su Materiale poroso per selezionare i componenti fluidi e solidi . Selezionare Nodo Fluido e in Proprietà fluido selezionare Sangue (definito nei passaggi precedenti). Selezionare Nodo solido e in Proprietà solide selezionare Tumore (definito nei passaggi precedenti). Nel nodo Solido, specificare la frazione di volume definita come θS (Tabella 1).
    8. Abilitare la selezione manuale e selezionare l'entità geometrica corrispondente al materiale specificato. Per seguire questo protocollo, si supponga che solo la regione tumorale sia un dominio poroelastico.
  5. Mesh
    1. In Nodo Mesh selezionare Dimensione e selezionare una mesh più fine predefinita.
    2. Aggiungi la funzione Tetraedrico libero nel nodo Mesh . Questo passaggio consente di ottenere una mesh raffinata nelle aree critiche.
      NOTA: Per questo modello, abbiamo identificato i bordi del tumore e l'estremità distale dei modelli di ago ipodermico come aree critiche.
    3. Seleziona le geometrie di interesse e personalizza la dimensione massima (0,25 mm) e minima dell'elemento in modo che il componente più piccolo (ad esempio, la punta dell'ago) sia discretizzato da almeno quattro elementi di mesh (la mesh completa è composta da 1.487.828 elementi).

2. Fisica

  1. Configurazione per il problema elettrico
    NOTA: I passaggi seguenti forniscono informazioni su come impostare i parametri per calcolare la distribuzione del campo elettrico (Figura 2, blocco 1) che fornirà la fonte di calore a radiofrequenza (Q).
    1. Fare clic con il pulsante destro del mouse sul nodo Correnti elettriche.
    2. Per le condizioni al contorno elettriche mostrate nella Figura 3A, selezionare Terminale e Terra come contorni.
      1. Per Terminale, selezionare manualmente l'estremità prossimale (in alto) di uno dei due aghi. L'ago identificato fornirà l'alimentazione in ingresso.
      2. In Terminale, selezionare Potenza e specificare il valore in base al protocollo energetico desiderato. Per seguire questo protocollo, selezionare 0,5 W per l'ipertermia lieve sulla base di esperimenti preliminari ex vivo 20.
      3. Selezionare Rettifica e selezionare manualmente la superficie prossimale del secondo ago. Questo ago fungerà da elettrodo di ritorno per il percorso di ritorno della corrente elettrica.
      4. Applicare l'isolamento elettrico alla superficie esterna rimanente del modello.
  2. Configurazione per il problema termico
    NOTA: I passaggi seguenti mostrano come includere le funzioni di perfusione sanguigna dipendenti dalla temperatura (equazioni 4 e 5) nell'equazione di trasferimento del biocalore per modellare il dissipatore di calore causato dal flusso sanguigno.
    1. Selezionare il nodo Trasferimento di calore nei solidi e specificare 33 °C come valore iniziale della temperatura.
    2. Per modellare l'effetto del dissipatore di calore dovuto al flusso sanguigno, fare clic con il pulsante destro del mouse su Trasferimento di calore nei solidi, aggiungere il dominio Fonte di calore e selezionare la geometria in cui deve essere considerato l'effetto del dissipatore di calore (ad esempio, tumore e tessuto normale ). Seleziona Sorgente Generale | Definito dall'utente in cui è possibile digitare l'espressione per il dissipatore di calore.
    3. Per le condizioni al contorno termiche mostrate nella Figura 3B, fare clic con il pulsante destro del mouse su Trasferimento di calore, aggiungere Flusso di calore come condizione al contorno e specificare le superfici esterne a cui viene applicato il flusso di calore. Selezionare Flusso di calore convettivo come tipo di flusso. Per il coefficiente di trasmittanza termica, utilizzare h = 15 W/(m2 · K) modellare il meccanismo dello scambio termico naturale tra la pelle e l'aria28. Specificare la temperatura esterna. Utilizzare T = 20 °C per modellare la temperatura ambiente nell'ambiente di laboratorio.
  3. Preparazione per il problema fluidodinamico
    NOTA: I passaggi seguenti descrivono come implementare l'equazione di conservazione della massa illustrata in Figura 2 (Blocco 3) e come può essere collegato alla variazione della temperatura.
    1. Selezionate il nodo PDE forma coefficiente e specificate Pressione (Pressure ) come variabile dipendente. In questa fase, l'unità Pascal (Pa) viene assegnata automaticamente.
      NOTA: Una volta calcolata la simulazione, i risultati possono essere visualizzati e/o esportati utilizzando l'unità scelta. Presentiamo i risultati utilizzando l'unità mmHg per coerenza con la letteratura (vedere la sezione dei risultati rappresentativi).
    2. Specificare l'unità di conduttanza del fluido 1/s come quantità del termine di origine.
    3. Definisci il nome per identificare la variabile (Pi greco, pressione del fluido interstiziale in questo studio).
    4. Fare clic con il pulsante destro del mouse sul nodo PDE Coefficient Form e selezionare il dominio Coefficient Form . Specificare l'entità geometrica a cui si riferisce l'equazione (tumore). Ripetere gli stessi passaggi e selezionare il tessuto rimanente (tessuto normale) a cui verrà applicata una PDE diversa.
    5. Per il modello tumorale, specificare i seguenti coefficienti e termini per ottenere l'equazione di conservazione della massa (Figura 2 blocco 3): coefficiente di diffusione Kidel tumore (Tabella 1); coefficiente figure-protocol-1 di smorzamento ); Termine figure-protocol-2 di origine . Per il modello tumorale, trascurare il contributo del sistema linfatico. Impostare tutti gli altri coefficienti su zero.
    6. Per il modello di tessuto normale, specificare i seguenti coefficienti e termini per ottenere l'equazione di conservazione della massa (Figura 2 blocco 3): coefficiente di diffusione Ki del tessuto normale (Tabella 1); coefficiente figure-protocol-3 di smorzamento ; Termine figure-protocol-4 di origine . Per considerare il tessuto normale come un tessuto funzionante normale, si consideri il contributo del sistema linfatico. Impostare tutti gli altri coefficienti su zero.
    7. Per creare il collegamento con la simulazione elettromagnetico-termica, esprimere la pressione vascolare Pv in funzione della temperatura (per mezzo della variabile di perfusione sanguigna, vedi Equazioni 6 e 7).
    8. Fare clic con il pulsante destro del mouse su Coefficient Form PDE e selezionare Valori iniziali. Selezionare il dominio geometrico (tumore) e ripetere lo stesso passaggio per il modello di tessuto normale (tessuto normale). Specificare Pi0per il tumore e il tessuto normale in base ai valori riportati nella Tabella 1.
    9. Per le condizioni al contorno relative allo studio fluidodinamico, mostrate nella Figura 3C, fare clic con il pulsante destro del mouse su Coefficient Form PDE e selezionare Condizioni al contorno di Dirichlet. Selezionare la superficie esterna del dominio tissutale normale e assegnare il valore Pi0corrispondente al tessuto normale (Tabella 1).

3. Esegui le simulazioni e visualizza i risultati

NOTA: Come ultimo passaggio prima del calcolo, specificare il tempo (simulando la durata della procedura) e la frequenza operativa:

  1. Selezionare Frequenza-Transitorio dal nodo Studio.
    1. Specificare l'unità di tempo.
    2. Da Tempi di uscita, selezionare l'intervallo (sul lato destro) e specificare 0 s come inizio, 5 s come passo e 900 s come stop.
    3. Impostare la frequenza su 500e3 Hz.
  2. Selezionare Calcola per eseguire le simulazioni.
  3. Per la visualizzazione dei risultati, selezionare Set di dati nel nodo Risultati.
    1. Fare clic con il pulsante destro del mouse per selezionare un piano di taglio per definire il piano da utilizzare per visualizzare le distribuzioni 2D (ad esempio, zx-plane a y = 0).
    2. Fare clic con il pulsante destro del mouse per selezionare un punto di taglio nel volume 3D per visualizzare la variazione di un parametro nel tempo.
  4. Da Risultati sulla barra multifunzione superiore,
    1. Selezionare il gruppo di grafici 2D per visualizzare la distribuzione bidimensionale di una variabile (ad esempio, la temperatura) su uno dei piani identificati nei passaggi precedenti.
    2. Selezionare il gruppo di grafici 1D per visualizzare i risultati 1D (ad esempio, la pressione nel tempo) nel punto o in più punti identificati nei passaggi precedenti.
      NOTA: Il tempo necessario per eseguire le simulazioni con le impostazioni descritte in questo protocollo è di circa 2,5 ore.

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Risultati

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

La distribuzione omogenea dell'elevata pressione del fluido interstiziale all'interno del tumore e un calo ai valori normali (0-3 mmHg) alla periferia sono i tratti distintivi della TME. La Figura 4 e la Figura 5 mostrano le condizioni iniziali (t = 0 min) di temperatura (A), pressione del fluido interstiziale (B) e velocità del fluido (C). Prima di iniziare il riscaldamento, quando la temperatura iniziale è di 33 °C, il valore della pressione del fluido interst...

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Discussione

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Presentiamo un protocollo di modellazione computazionale per accoppiare simulazioni elettrico-termiche transitorie con simulazioni fluidodinamiche per studiare l'impatto dell'ipertermia RF sui profili di pressione dei fluidi termici e interstiziali nei tumori. L'aspetto chiave è la costruzione di un flusso di lavoro numerico in grado di catturare la relazione esistente tra temperatura e pressione vascolare, che a sua volta guida le variazioni della pressione del fluido interstiziale.

Abbiamo u...

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Dichiarazioni

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Gli autori non hanno conflitti di interesse da rivelare.

Ringraziamenti

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Lo studio è stato sostenuto da sovvenzioni della National Science Foundation (n. 2039014) e del National Cancer Institute (R37CA269622).

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Materiali

Elenco dei materiali utilizzati in questo articolo
NomeAziendaNumero di catalogoCommenti
COMSOL Multiphysics (v. 6.0)COMSOL AB, Stoccolma, SveziaSoftware utilizzato per implementare il flusso di lavoro computazionale descritto nel protocollo
Dell 1.8.0, Intel(R) Core(TM) di 11a generazione i7-11850H @ 2.50GHz, 2496 Mhz, 8 Core(s), 16 Processore(i) logico(i), 32 GBDell Inc.Laptop utilizzato per eseguire simulazioni computazionali
RAM

Riferimenti

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Nia, H. T., Munn, L. L., Jain, R. K. Physical traits of cancer. Science. 370 (6516), 546-556 (2020).
  2. Heldin, C. -H., Rubin, K., Pietras, K., Östman, A. High interstitial fluid pressure - an obstacle in cancer therapy. Nature Reviews Cancer. 4 (10), 806-813 (2004).
  3. Jain, R. K. Determinants of tumor blood flow: a review. Cancer Research. 48, 2641-2658 (1988).
  4. Stylianopoulos, T., Munn, L. L., Jain, R. K. Reengineering the physical microenvironment of tumors to improve drug delivery and efficacy: from mathematical modeling to bench to bedside. Trends in Cancer. 4 (4), 292-319 (2018).
  5. Sheth, R. A., Hesketh, R., Kong, D. S., Wicky, S., Oklu, R. Barriers to drug delivery in interventional oncology. Journal of Vascular and Interventional Radiology. 24 (8), 1201-1207 (2013).
  6. Chauhan, V. P., Stylianopoulos, T., Boucher, Y., Jain, R. K. Delivery of molecular and nanoscale medicine to tumors: transport barriers and strategies. Annual Review of Chemical and Biomolecular Engineering. 2 (1), 281-298 (2011).
  7. Li, R., et al. Interstitial flow promotes macrophage polarization toward an M2 phenotype. Molecular Biology of the Cell. 29 (16), 1927-1940 (2018).
  8. Stine, C. A., Munson, J. M. Autologous gradient formation under differential interstitial fluid flow environments. Biophysica. 2 (1), 16-33 (2022).
  9. Provenzano, P. P., et al. Enzymatic targeting of the stroma ablates physical barriers to treatment of pancreatic ductal adenocarcinoma. Cancer Cell. 21 (3), 418-429 (2012).
  10. Pal, K., Sheth, R. A. Engineering the tumor immune microenvironment through minimally invasive interventions. Cancers. 15 (1), 196(2022).
  11. Dunne, M., Regenold, M., Allen, C. Hyperthermia can alter tumor physiology and improve chemo- and radio-therapy efficacy. Advanced Drug Delivery Reviews. 163-164, 98-124 (2020).
  12. Vaupel, P., et al. From localized mild hyperthermia to improved tumor oxygenation: physiological mechanisms critically involved in oncologic thermo-radio-immunotherapy. Cancers. 15 (5), 1394(2023).
  13. Stapleton, S., et al. Radiation and heat improve the delivery and efficacy of nanotherapeutics by modulating intratumoral fluid dynamics. ACS Nano. 12 (8), 7583-7600 (2018).
  14. Li, Q., Zhou, Y., Zhang, F., McGregor, H., Yang, X. Radiofrequency hyperthermia enhances locally delivered oncolytic immuno-virotherapy for pancreatic adenocarcinoma. CardioVascular and Interventional Radiology. 45 (12), 1812-1821 (2022).
  15. Mpekris, F., et al. Combining microenvironment normalization strategies to improve cancer immunotherapy. Proceedings of the National Academy of Sciences. 117 (7), 3728-3737 (2020).
  16. Netti, P. A., Baxter, L. T., Boucher, Y., Skalak, R., Jam, R. K. Time-dependent behavior of interstitial fluid pressure in solid tumors: implications for drug delivery. Cancer Research. 15 (55), 5451-5458 (1995).
  17. Andreozzi, A., Iasiello, M., Netti, P. A. Effects of pulsating heat source on interstitial fluid transport in tumour tissues. Journal of The Royal Society Interface. 17 (170), 612-626 (2020).
  18. Leunig, M., Goetz, A. E., Messmer, K. Interstitial fluid pressure in solid tumors following hyperthermia: possible correlation with therapeutic response. Cancer Research. 52, 487-490 (1992).
  19. Muñoz, N. M., et al. Immune modulation by molecularly targeted photothermal ablation in a mouse model of advanced hepatocellular carcinoma and cirrhosis. Scientific Reports. 12 (1), 14449(2022).
  20. Bottiglieri, A., et al. RF-hyperthermia to modulate tumor interstitial fluid pressure: an in vivo pilot study. 38th Annual Society for Thermal Medicine Meeting. , (2023).
  21. Baxter, L. T., Jain, R. K. Transport of fluid and macromolecules in tumors. I. Role of interstitial pressure and convection. Microvascular Research. 37 (1), 77-104 (1989).
  22. Stapleton, S., et al. A mathematical model of the enhanced permeability and retention effect for liposome transport in solid tumors. PLoS ONE. 8 (12), 1-10 (2013).
  23. Rossmann, C., Haemmerich, D. Review of temperature dependence of thermal properties, dielectric properties, and perfusion of biological tissues at hyperthermic and ablation temperatures. Critical Reviews in Biomedical Engineering. 42 (6), 467-492 (2014).
  24. Song, C. W., Lokshina, A., Rhee, J. G., Patten, M., Levitt, S. H. Implication of blood flow in hyperthermic treatment of tumors. IEEE Transactions on Biomedical Engineering. 31 (1), 9-16 (1984).
  25. Tompkins, D. T., et al. Temperature-dependent versus constant-rate blood perfusion modelling in ferromagnetic thermoseed hyperthermia: results with a model of the human prostate. International Journal of Hyperthermia. 10 (4), 517-536 (1994).
  26. Andreozzi, A., Iasiello, M., Netti, P. A. A thermoporoelastic model for fluid transport in tumour tissues. Journal of The Royal Society Interface. 16 (154), 0030-0046 (2019).
  27. Hasgall, P. A., et al. IT'IS Database for thermal and electromagnetic parameters of biological tissues. , (2022).
  28. Cavagnaro, M., et al. Influence of the target tissue size on the shape of ex vivo microwave ablation zones. International Journal of Hyperthermia. 31 (1), 48-57 (2015).
  29. Munson, J., Shieh, A. Interstitial fluid flow in cancer: implications for disease progression and treatment. Cancer Management and Research. 19 (6), 317-328 (2014).
  30. Muñoz, N. M., et al. Influence of injection technique, drug formulation and tumor microenvironment on intratumoral immunotherapy delivery and efficacy. Journal for ImmunoTherapy of Cancer. 9 (2), 0018-0027 (2021).
  31. Swartz, M. A., Lund, A. W. Lymphatic and interstitial flow in the tumour microenvironment: linking mechanobiology with immunity. Nature Reviews Cancer. 12 (3), 210-219 (2012).
  32. Mehta, A., Oklu, R., Sheth, R. A. Thermal ablative therapies and immune checkpoint modulation: can locoregional approaches effect a systemic response. Gastroenterology Research and Practice. 2016, 1-11 (2016).
  33. Song, C. W., Park, H., Griffin, R. J. Improvement of tumor oxygenation by mild hyperthermia. Radiation Research. 155 (4), 515-528 (2001).
  34. Dewhirst, M. W., Oleson, J. R., Kirkpatrick, J., Secomb, T. W. Accurate three-dimensional thermal dosimetry and assessment of physiologic response are essential for optimizing thermoradiotherapy. Cancers. 14 (7), 1701(2022).

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Ristampe e permessi

Richiedi il permesso di riutilizzare il testo o le figure di questo articolo JoVE

Richiedi permesso

Tag

Trattamento ipertermicopressione del fluido interstizialetrasferimento di calore biologicoipertermia a radiofrequenzadinamica dei fluidiperfusione sanguignaequazioni differenziali alle derivate parzialidispositivo a radiofrequenza bipolare

Articoli correlati