Artículo de método

Un enfoque de modelado computacional para investigar la influencia de la hipertermia en el microambiente tumoral

DOI:

10.3791/65870

⸱

1 de diciembre de 2023

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

En este artículo

Resumen

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

El artículo describe un protocolo para simular los perfiles de temperatura transitoria y la variación espacio-temporal acoplada de la presión del fluido intersticial después del calentamiento suministrado por un sistema de hipertermia de radiofrecuencia dipolar. El protocolo se puede utilizar para evaluar la respuesta de los parámetros biofísicos que caracterizan el microambiente tumoral a las técnicas de hipertermia intervencionista.

Resumen

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

Las propiedades biofísicas del microambiente tumoral difieren sustancialmente de las de los tejidos normales. Una constelación de características, que incluyen disminución de la vascularización, falta de drenaje linfático y presión intersticial elevada, disminuye la penetración de los tratamientos en los tumores. La hipertermia local dentro del tumor puede alterar las propiedades microambientales, como la presión del líquido intersticial, lo que puede conducir a mejoras en la penetración del fármaco. En este contexto, los modelos computacionales multifísicos pueden proporcionar información sobre la interacción entre los parámetros biofísicos dentro del microambiente tumoral y pueden guiar el diseño y la interpretación de experimentos que prueben los efectos biológicos de la hipertermia local.

Este artículo describe un flujo de trabajo paso a paso para un modelo computacional que acopla ecuaciones diferenciales parciales que describen la distribución de corriente eléctrica, la transferencia de biocalor y la dinámica de fluidos. El objetivo principal es estudiar los efectos de la hipertermia administrada por un dispositivo de radiofrecuencia bipolar sobre la presión del líquido intersticial dentro del tumor. Se presenta el sistema de expresiones matemáticas que relacionan la distribución de corriente eléctrica, la transferencia de biocalor y la presión del fluido intersticial, enfatizando los cambios en la distribución de la presión del fluido intersticial que podrían ser inducidos por la intervención térmica.

Introducción

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

La presión elevada del líquido intersticial (IFP) es una característica distintiva de los tumores sólidos1. La fuga de líquido hacia el intersticio desde los vasos sanguíneos hiperpermeables se ve desequilibrada por la salida de líquido debido a la compresión de las venas intratumorales y a la ausencia de linfáticos 1,2,3. Junto con otros parámetros biofísicos que son anormales dentro del microambiente tumoral (TME), incluyendo el estrés sólido y la rigidez, la IFP elevada socava la eficacia de la administración de fármacos tanto sistémica como local 4,5,6. La presión del líquido intersticial en los tumores sólidos oscila entre 5 mmHg (glioblastoma y melanoma) y 30 mmHg (carcinoma de células renales) en comparación con 1-3 mmHg en el tejido normal2. La IFP alta es responsable de aumentar el flujo de líquido hacia el margen del tumor y expone las células estromalas, las células infiltradas y otros componentes extracelulares al estrés de cizallamiento 1,4. Las alteraciones mecanobiológicas sostienen un TME inmunosupresor, por ejemplo, al aumentar la brotación endotelial, que apoya la angiogénesis, la migración e invasión de células cancerosas, la expresión transformante del factor de crecimiento-β (TGF-β) y el endurecimiento del estroma 7,8,9.

Varios estudios han explorado las terapias basadas en la energía con la intención de disminuir la PFI, incluyendo el ultrasonido de baja intensidad, el ultrasonido focalizado de alta intensidad, los campos eléctricos pulsados y las terapias térmicas 5,10,11. Se ha demostrado que el calentamiento a temperaturas en el rango de 40-43 °C, denominado hipertermia leve, aumenta la perfusión sanguínea tumoral y, por lo tanto, puede contribuir a la expansión de las venas comprimidas y a la reducción de la presión vascular al facilitar la intravasación y el drenaje del líquido intersticial11,12. Algunos estudios recientes han demostrado el potencial de la hipertermia para reducir la PIIF y, en consecuencia, para facilitar la distribución de fármacos o agentes de contraste dentro de un tumor13,14. Estos estudios también muestran un aumento de la infiltración de células T después de la hipertermia en comparación con los grupos control sin tratamiento13.

Los prometedores resultados de los experimentos in vivo con animales pequeños motivan la realización de nuevos estudios que emplean enfoques computacionales para avanzar en la comprensión de cómo los parámetros físicos dentro del TME se ven afectados por las intervenciones físicas 4,15,16,17. Los resultados de los modelos computacionales pueden complementar los estudios experimentales in vivo para descubrir la relación causa-efecto que subyace al calentamiento local (u otras fuentes de energía externas) y el IFP. Esto puede ser particularmente instructivo dados los desafíos con la medición de las variaciones espaciales en la IFP con transductores de presión basados en catéteres y agujas, que generalmente proporcionan mediciones puntuales 9,16,18,19. En el contexto de la administración de fármacos, la comprensión de los mecanismos biofísicos clave es esencial para definir el protocolo de calentamiento adecuado, así como la ventana de tiempo para la inyección de fármacos, a fin de aumentar la probabilidad de una distribución eficaz del fármaco. La información cuantitativa en términos de cambios en las características biofísicas del TME, incluido, entre otros, el IFP, también podría proporcionar información sobre la interpretación de la respuesta inmunológica (por ejemplo, infiltración de células T) a estímulos externos.

Presentamos un protocolo para el modelado computacional de cambios mediados térmicamente en los perfiles IFP tumorales. Específicamente, el protocolo detalla cómo modelar un aparato personalizado para animales pequeños para administrar terapia térmica controlada con corriente de radiofrecuencia, simular perfiles de temperatura transitoria después del calentamiento y combinar simulaciones de dinámica de fluidos para calcular la variación espacio-temporal de la IFP tumoral en respuesta a la terapia térmica. Este modelo refleja las características esenciales de la configuración experimental que utilizamos en un modelo de tumor subcutáneo (McArdle RH7777, ATCC) en un estudio experimental previo20.

En la Figura 1 se muestra el modelo computacional que implementamos para calcular los cambios inducidos térmicamente en el IFP en un tumor rodeado de tejido normal. Se modela un par de agujas hipodérmicas insertadas en el tumor para administrar calor con corriente de radiofrecuencia a 500 kHz. Se supone un material poroso en el dominio tumoral, compuesto por dos fases: la fase sólida representa la matriz extracelular sólida, y la fase fluida representa el líquido intersticial. En el caso de un cambio de presión o una deformación de la matriz resultante de un estímulo externo, por ejemplo, un aumento de la temperatura, los componentes sólidos y fluidos se reorganizan. Esto provoca el movimiento del líquido intersticial a través de la matriz sólida extracelular 16,17,21.

A partir de la teoría de la poroelasticidad, el tensor de tensión S (Pa) (ecuación [1]) es la combinación del término elástico que describe el cambio en el volumen del componente sólido en relación con las condiciones iniciales, y un término poroso que describe la tensión inducida por la presión hidrostática del componente fluido.

figure-introduction-1(1)

Donde, λ, μ (Pa) son los parámetros de Lamé, E es el tensor de deformación, e es el tensor de deformación volumétrica, Pi (Pa) es la presión del fluido intersticial (I es la matriz identidad). Se asumen condiciones de estado estacionario para el componente sólido bajo tensión poroelástica, lo que significa que los componentes del tensor de tensión son ortogonales, figure-introduction-2.

La Figura 2 muestra el sistema de ecuaciones matemáticas implementado en el modelo poroelástico descrito y la interacción entre los componentes del modelo multifísico presentado. El flujo de trabajo de las simulaciones computacionales incluye:

Ecuaciones de problemas eléctricos. La solución de las ecuaciones del problema eléctrico proporciona la fuente de calor de RF promediada en el tiempo Q (calentamiento Joule). Con este fin, se utiliza una aproximación cuasiestática a las ecuaciones de Maxwell para calcular la distribución del campo eléctrico promediado en el tiempo E (V/m) (Figura 2, bloque 1).

Ecuaciones de problemas térmicos. La solución de la ecuación de Pennes sobre el biocalor (Figura 2, bloque 2) proporciona la variación espacial y temporal de la temperatura T (°C) como resultado de la fuente de calor (Q) ligada a la energía electromagnética absorbida, el calentamiento pasivo ligado a la conducción térmica de los tejidos (figure-introduction-3), y el efecto disipador de calor de la perfusión sanguínea tisular (cWb(T) (T - Tb)). El término del disipador de calor se aproxima al intercambio de calor entre la sangre que fluye en la microvasculatura y el tejido adyacente donde se absorbe la energía electromagnética. La ecuación de transferencia de calor también incluye el término de advección (figure-introduction-4), que describe el cambio en la temperatura causado por el movimiento del fluido intersticial a través de la matriz extracelular del modelo poroelástico. Sin embargo, este término tiene un impacto insignificante en el perfil de temperatura en comparación con los otros mecanismos responsables del cambio de temperatura.

Ecuaciones de problemas de dinámica de fluidos. La ecuación de conservación de masa (Figura 2, bloque 3) combinada con la ley de Darcy (Figura 2, bloque 4) da como resultado la variación espacial y temporal de la presión intersticial del fluido Piresultante del equilibrio entre la fuente (figure-introduction-5) y el sumidero (figure-introduction-6 ) del fluido. El término de presión transitoria en el lado izquierdo de la ecuación de conservación de masa, figure-introduction-7 , describe el reordenamiento de los componentes fluido y sólido en el material poroelástico. Esto se produce por la variación de la presión del líquido intersticial, Pi, impulsada por la variación de la presión vascular Pven función de la temperatura.

La diferencia entre la presión vascular (Pv) y la presión del líquido intersticial (Pi) es la fuente del líquido que fluye a través de la matriz extracelular. El término sumidero está relacionado con la diferencia de presión entre los vasos linfáticos (PL) y el espacio intersticial (Pi). En el tejido normal, la presión en la vasculatura linfática (~ -6-0 mmHg) es hasta dos veces menor que la presión del líquido intersticial13. Esta diferencia de presión asegura la eficacia de los vasos linfáticos para drenar el exceso de líquido que se extravasa de la pared de los vasos sanguíneos al intersticio. Para el modelo tumoral que aquí se presenta, no se tuvo en cuenta la contribución del sistema linfático 4,16,22.

Las expresiones matemáticas de las ecuaciones (2) a (5) se utilizan para describir la dependencia de la temperatura de la conductividad eléctrica y térmica de los tejidos y la perfusión sanguínea de los tejidos23,24. Se utilizan dos modelos matemáticos diferentes para describir la dependencia de la temperatura de la perfusión sanguínea en los dominios de tejido normal y tumoral, respectivamente24,25. Los modelos muestran que la perfusión sanguínea aumenta con la temperatura hasta nueve veces en comparación con la línea de base en el tejido normal y solo aproximadamente dos veces el valor de referencia en el dominio tumoral. Para ambos modelos, el aumento de la perfusión sanguínea se limita a las temperaturas dentro del rango de hipertermia leve (por debajo de 45 °C). Vale la pena mencionar que las expresiones matemáticas, ecuaciones (4) y (5), no describen completamente los mecanismos subyacentes a los cambios dependientes de la temperatura en la perfusión sanguínea en los dos tipos diferentes de tejido. Sin embargo, ayudan a representar la perfusión limitada que suele caracterizar el microambiente tumoral en comparación con los tejidos normales.

figure-introduction-8(2)

figure-introduction-9(3)

figure-introduction-10 (4)

figure-introduction-11 (5)

figure-introduction-12(6)

figure-introduction-13(7)

En este estudio, utilizamos las ecuaciones (6) y (7) para modelar la presión vascular en función de la perfusión sanguínea tanto para modelos normales como de tejido tumoral26. A partir de las ecuaciones (4) y (5), la tasa de flujo sanguíneo se puede expresar como la relación entre la perfusión sanguínea y la densidad sanguínea. La relación entre el flujo sanguíneo y la presión vascular está bien establecida en la literatura3: la tasa de flujo sanguíneo y la resistencia geométrica (o conductividad, Lp) de la vasculatura determinan la diferencia de presión dentro del vaso sanguíneo. La presión vascular puede expresarse en función de la temperatura (Ecuaciones (6) y (7)), aprovechando esta relación y el modelo de perfusión sanguínea dependiente de la temperatura (Ecuaciones (4) y (5)).

La implementación del flujo de trabajo computacional (Figura 2) y las propiedades dependientes de la temperatura de los modelos de tejidos se describen en detalle en la siguiente sección. Todas las propiedades del material y sus descripciones y valores de referencia (es decir, a temperatura corporal) se enumeran en la Tabla 1. Consulte la Tabla de Materiales para obtener detalles sobre COMSOL Multiphysics instalado en una computadora utilizada para implementar este protocolo computacional. El problema eléctrico se modeló utilizando el módulo AC/DC; La transferencia de biocalor se modeló utilizando la física de transferencia de calor; y el problema de dinámica de fluidos se modeló utilizando la interfaz de Matemáticas.

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Protocolo

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

1. Construir el modelo de un sistema de radiofrecuencia bipolar

  1. Pasos preliminares para configurar la interfaz
    1. Inicie COMSOL Multiphysics y haga clic en Model Wizard.
    2. Seleccione 3D como dimensión espacial.
    3. Seleccione el módulo AC/DC Physics | Campos y corrientes eléctricas | Corrientes eléctricas.
    4. Seleccione el módulo de transferencia de calor | Transferencia de calor en sólidos.
    5. Seleccione el módulo de Matemáticas | Interfaces PDE | Forma de coeficiente PDE.
    6. Seleccionar Estudio | Depende del tiempo. Haga clic en Listo.
    7. Una vez que aparezca el espacio de trabajo de Comsol:
      1. Seleccione Multifísica | Calentamiento electromagnético. Con este paso, la densidad de pérdida de potencia electromagnética se acopla automáticamente como fuente de calor para la ecuación de transferencia de biocalor.
        NOTA: Si Multiphysics no aparece automáticamente, especifique manualmente la fuente de calor electromagnética (mostrada en COMSOL como Densidad de pérdida volumétrica). Para obtener más detalles sobre cómo agregar la fuente de calor, consulte la sección 'Física', paso 2 'Configuración para el problema térmico'.
      2. Seleccione Estudiar en la cinta superior | Pasos del estudio | Frecuencia transitoria.
  2. Defina las geometrías. En la cinta superior, seleccione Geometría y, a continuación:
    1. Defina dos conos con las dimensiones enumeradas en la Tabla 2.
    2. Coloque los conos a la distancia indicada en la Tabla 2 (espaciados dentre sí ). Estos dos conos modelarán las dos agujas hipodérmicas utilizadas para construir el sistema de RF bipolar.
    3. Duplique los dos conos anteriores para modelar el aislamiento de las agujas; modificar el tamaño del cono de acuerdo con las dimensiones reportadas en la Tabla 2.
    4. Seleccione un cilindro (altura, hm y diámetro dm) para modelar la masa de músculo colocada en z = - 9 mm (x = 0, y = 0). Los valores de cada dimensión se enumeran en la Tabla 2.
    5. Seleccione un cilindro (altura, hs y diámetro ds) para modelar la capa delgada de piel colocada a z = 4 mm (x = 0, y = 0). Los valores de cada dimensión se enumeran en la Tabla 2.
    6. Seleccione una esfera (diámetro, dt) para modelar el tumor subcutáneo colocado a z = -0,5 mm (x = 0, y = 0). El tamaño del tumor se indica en el Cuadro 2.
    7. Para facilitar la selección de las geometrías en los siguientes pasos del protocolo, recomendamos lo siguiente:
      1. En la cinta Geometría , seleccione Operaciones virtuales | Formar dominios compuestos.
      2. Seleccione todos los dominios relacionados con la parte conductora de electricidad de las agujas para crear una geometría compuesta.
      3. Repita el mismo procedimiento para crear dominios compuestos para las geometrías de aislamiento de la aguja.
  3. Definir las propiedades de los modelos biológicos de tejidos.
    NOTA: Los siguientes pasos describen el procedimiento para implementar las expresiones matemáticas descritas por las ecuaciones (2)-(7).
    1. En el nodo Componente, haga clic con el botón derecho para seleccionar Definiciones.
    2. En Funciones, seleccione Analítico.
      1. Especifique el nombre de la función (por ejemplo, k_muscle o sigma_muscle) y escriba la expresión matemática coherente con la Ec. 2).
      2. Especifique la temperatura (T) como argumento.
      3. Especifique las unidades de la función: S/m en el caso de conductividad eléctrica.
      4. Repita los pasos anteriores del 1 al 3 para implementar la Ec. 3, modificando la unidad en consecuencia (es decir, W/(m·K) para la conductividad térmica).
      5. Especifique la unidad para el argumento: K (kelvin) para la temperatura. En Parámetros de trazado, especifique el rango de valores del argumento de función (es decir, temperatura). Para seguir este protocolo, utilice un rango de 33-100 °C (306,15-373,15 K).
      6. Repita los pasos anteriores del 1 al 5 para sumar las funciones dependientes de la temperatura de las conductividades eléctricas ( Ec. 2) y térmicas (Ec. 3) para cada modelo de tejido (es decir, músculo, piel y tumor) utilizando los valores nominales enumerados en la Tabla 1 (el tejido normal se refiere tanto al músculo como a la piel).
    3. En Funciones, seleccione A partes para implementar las ecuaciones (4)-(7):
      1. Especifique el nombre de la función.
      2. Especifique la temperatura (T) como argumento de la función.
      3. Escriba la expresión matemática para cada intervalo de temperatura consistente con las ecuaciones (4)-(7).
      4. Repita los pasos anteriores del 1 al 3 para sumar las funciones dependientes de la temperatura de la perfusión sanguínea y la presión vascular para cada modelo de tejido utilizando los valores nominales enumerados en la Tabla 1 (el tejido normal se refiere tanto al músculo como a la piel).
  4. Asigne propiedades de material a los componentes de geometría.
    1. En el nodo componente, seleccione Materiales.
    2. Seleccione materiales en blanco que incluyan tejido normal, tejido tumoral, sangre, PTFE y acero inoxidable.
    3. Habilite la selección manual y seleccione la entidad geométrica correspondiente al material especificado.
      1. El tejido normal se asocia con geometrías que modelan el músculo y la piel.
      2. El tumor y los tejidos sanguíneos se asocian con la geometría del tumor.
      3. El material de PTFE se asocia con las geometrías que modelan el aislante de la aguja.
      4. El material de acero inoxidable se asocia con las geometrías de los conos que modelan el suelo y las agujas activas.
    4. Para las conductividades eléctricas y térmicas dependientes de la temperatura 23, escriba el nombre elegido de la función y el argumento relacionado (es decir, T) que aparece en el nodo Definiciones.
    5. Para conocer las propiedades del material que no dependen de la temperatura, consulte los valores de referencia27 enumerados en la Tabla 1.
      NOTA: Nos basamos en la teoría poroelástica para calcular la presión 16,17,26. Los siguientes pasos muestran cómo las propiedades de un material poroso se pueden asignar a un dominio específico.
    6. En Materiales, seleccione Más materiales | Material poroso.
    7. Haga clic con el botón derecho en Material poroso para seleccionar componentes fluidos y sólidos . Seleccione Nodo Fluido y, en Propiedades del fluido, seleccione Sangre (definida en los pasos anteriores). Seleccione Nodo sólido y, en Propiedades sólidas, seleccione Tumor (definido en los pasos anteriores). En el nodo Sólido, especifique la fracción de volumen definida como θS(Tabla 1).
    8. Habilite la selección manual y seleccione la entidad geométrica correspondiente al material especificado. Para seguir este protocolo, suponga que solo la región tumoral es un dominio poroelástico.
  5. Mallado
    1. En el nodo Malla, seleccione Tamaño y seleccione una malla más fina predefinida.
    2. Agregue la función Tetraédrico libre en el nodo Malla. Este paso permite una malla refinada en las áreas críticas.
      NOTA: Para este modelo, identificamos los bordes del tumor y el extremo distal de los modelos de aguja hipodérmica como áreas críticas.
    3. Seleccione las geometrías de interés y personalice el tamaño máximo (0,25 mm) y mínimo del elemento de manera que el componente más pequeño (por ejemplo, la punta de la aguja) sea discretizado por al menos cuatro elementos de malla (la malla completa consta de 1.487.828 elementos).

2. Física

  1. Configuración para el problema eléctrico
    NOTA: Los siguientes pasos proporcionan información sobre cómo configurar los parámetros para calcular la distribución del campo eléctrico (Figura 2, bloque 1) que proporcionará la fuente de calor de radiofrecuencia (Q).
    1. Haga clic con el botón derecho en el nodo Corrientes eléctricas.
    2. Para las condiciones de contorno eléctrico que se muestran en la Figura 3A, seleccione Terminal y Tierra como límites.
      1. En Terminal, seleccione manualmente el extremo proximal (en la parte superior) de una de las dos agujas. La aguja identificada proporcionará la potencia de entrada.
      2. En Terminal, seleccione Potencia y especifique el valor según el protocolo de energía deseado. Para seguir este protocolo, seleccione 0,5 W para una hipertermia leve basada en experimentos preliminares ex vivo 20.
      3. Seleccione Tierra y seleccione manualmente la superficie proximal de la segunda aguja. Esta aguja actuará como un electrodo de retorno para la trayectoria de retorno de la corriente eléctrica.
      4. Aplique aislamiento eléctrico a la superficie externa restante del modelo.
  2. Configuración para el problema térmico
    NOTA: Los siguientes pasos muestran cómo incluir las funciones de perfusión sanguínea dependientes de la temperatura (Ecuaciones 4 y 5) en la ecuación de transferencia de biocalor para modelar el disipador de calor causado por el flujo sanguíneo.
    1. Seleccione el nodo Transferencia de calor en sólidos y especifique 33 °C como valor inicial de la temperatura.
    2. Para modelar el efecto del disipador de calor debido al flujo sanguíneo, haga clic con el botón derecho en Transferencia de calor en sólidos, agregue el dominio de la fuente de calor y seleccione la geometría donde se debe considerar el efecto del disipador de calor (es decir, tumor y tejido normal ). Seleccione la fuente general | Definido por el usuario donde se puede escribir la expresión para el disipador de calor.
    3. Para las condiciones de contorno térmico que se muestran en la Figura 3B, haga clic con el botón derecho en Transferencia de calor, agregue Flujo de calor como condición de contorno y especifique las superficies externas a las que se aplica el flujo de calor. Seleccione Flujo de calor convectivo como tipo de flujo. Para el coeficiente de transferencia de calor, utilice h = 15 W/(m2 · K) modelar el mecanismo de intercambio natural de calor entre la piel y el aire28. Especifique la temperatura externa. Utilice T = 20 °C para modelar la temperatura ambiente en el entorno de laboratorio.
  3. Configuración para el problema de la dinámica de fluidos
    NOTA: Los siguientes pasos describen cómo implementar la ecuación de conservación de masa ilustrada en Figura 2 (Bloque 3) y cómo se puede relacionar con la variación de la temperatura.
    1. Seleccione el nodo PDE de forma de coeficiente y especifique Presión como variable dependiente. En esta etapa, la unidad Pascal (Pa) se asigna automáticamente.
      NOTA: Una vez calculada la simulación, los resultados pueden visualizarse y/o exportarse utilizando la unidad elegida. Presentamos los resultados utilizando la unidad mmHg para mantener la consistencia con la literatura (ver la sección de resultados representativos).
    2. Especifique la unidad de conductancia del fluido 1/s como cantidad del término de origen.
    3. Defina el nombre para identificar la variable (Pi, presión del fluido intersticial en este estudio).
    4. Haga clic con el botón derecho en el nodo PDE de Formulario de coeficiente y seleccione el dominio Formulario de coeficiente . Especifique la entidad geométrica a la que se refiere la ecuación (tumor). Repita los mismos pasos y seleccione el tejido restante (tejido normal) al que se aplicará un PDE diferente.
    5. Para el modelo tumoral, especificar los siguientes coeficientes y términos para obtener la ecuación de conservación de masa (Figura 2, bloque 3): coeficiente de difusión Kidel tumor (Tabla 1); Coeficiente figure-protocol-1 de amortiguación ); Término figure-protocol-2 de origen . Para el modelo tumoral, hay que descuidar la contribución del sistema linfático. Establezca todos los demás coeficientes iguales a cero.
    6. Para el modelo de tejido normal, especificar los siguientes coeficientes y términos para obtener la ecuación de conservación de masa (Figura 2 bloque 3): coeficiente de difusión Ki del tejido normal (Tabla 1); Coeficiente figure-protocol-3 de amortiguación ; Término figure-protocol-4 de origen . Para considerar el tejido normal como un tejido que funciona normalmente, considere la contribución del sistema linfático. Establezca todos los demás coeficientes iguales a cero.
    7. Para crear el vínculo con la simulación electromagnético-térmica, exprese la presión vascular Pv en función de la temperatura (mediante la variable de perfusión sanguínea, ver Ecuaciones 6 y 7).
    8. Haga clic con el botón derecho en Formulario de coeficiente PDE y seleccione Valores iniciales. Seleccione el dominio geométrico (tumor) y repita el mismo paso para el modelo de tejido normal (tejido normal). Especifique Pi0para el tumor y el tejido normal de acuerdo con los valores enumerados en la Tabla 1.
    9. Para las condiciones de contorno relacionadas con el estudio de fluidodinámica, que se muestran en la Figura 3C, haga clic con el botón derecho en Forma de coeficiente PDE y seleccione Condiciones de contorno de Dirichlet. Seleccione la superficie externa del dominio de tejido normal y asigne el valor de Pi0correspondiente al tejido normal (Tabla 1).

3. Ejecute las simulaciones y muestre los resultados

NOTA: Como último paso antes de calcular, especifique el tiempo (simulando la duración del procedimiento) y la frecuencia de operación:

  1. Seleccione Frecuencia-Transitorio en el nodo Estudio.
    1. Especifique la (s) unidad(es) de tiempo.
    2. En Tiempos de salida, seleccione el rango (en el lado derecho) y especifique 0 s como inicio, 5 s como paso y 900 s como parada.
    3. Establezca la frecuencia en 500e3 Hz.
  2. Seleccione Proceso para ejecutar las simulaciones.
  3. Para la visualización de los resultados, seleccione Conjuntos de datos en el nodo Resultados.
    1. Haga clic con el botón derecho para seleccionar un plano de corte y definir el plano que se utilizará para visualizar distribuciones 2D (por ejemplo, plano zx en y = 0).
    2. Haga clic con el botón derecho para seleccionar un punto de corte en el volumen 3D y mostrar la variación de un parámetro a lo largo del tiempo.
  4. En Resultados en la cinta superior,
    1. Seleccione el grupo de trazado 2D para visualizar la distribución bidimensional de una variable (por ejemplo, temperatura) en uno de los planos identificados en los pasos anteriores.
    2. Seleccione el grupo de trazado 1D para visualizar los resultados 1D (por ejemplo, la presión a lo largo del tiempo) en el punto o varios puntos identificados en los pasos anteriores.
      NOTA: El tiempo para ejecutar las simulaciones con la configuración descrita en este protocolo es de aproximadamente 2,5 h.

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Resultados

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

La distribución homogénea de una alta presión de líquido intersticial dentro del tumor y una caída a los valores normales (0-3 mmHg) en la periferia son características distintivas del TME. Las figuras 4 y 5 muestran las condiciones iniciales (t = 0 min) de temperatura (A), presión intersticial del fluido (B) y velocidad del fluido (C). Antes de iniciar el calentamiento, cuando la temperatura inicial es de 33 °C, el valor de la presión del líquido intersticial d...

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Discusión

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

Presentamos un protocolo de modelado computacional para acoplar simulaciones electrotérmicas transitorias con simulaciones fluidodinámicas para estudiar el impacto de la hipertermia de RF en los perfiles de presión de fluidos térmicos e intersticiales en tumores. El aspecto clave está en la construcción de un flujo de trabajo numérico capaz de capturar la relación existente entre la temperatura y la presión vascular, que a su vez impulsa los cambios en la presión del líquido intersticial.

Util...

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Divulgaciones

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

Los autores no tienen conflictos de intereses que revelar.

Agradecimientos

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

El estudio fue financiado por subvenciones de la Fundación Nacional de Ciencias (n.º 2039014) y el Instituto Nacional del Cáncer (R37CA269622).

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Materiales

Lista de materiales utilizados en este artículo
NombreEmpresaNúmero de catálogoComentarios
COMSOL Multiphysics (v. 6.0)COMSOL AB, Estocolmo, SueciaSoftware utilizado para implementar el flujo de trabajo computacional descrito en el protocolo
Dell 1.8.0, Intel(R) Core(TM) i7-11850H de 11.ª generación @ 2.50GHz, 2496 Mhz, 8 núcleo(s), 16 procesador(es) lógico(s), 32 GB de RAMDell Inc. Portátil utilizado para ejecutar simulaciones computacionales

Referencias

$$\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).

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Reimpresiones y permisos

Solicitar permiso para reutilizar el texto o las figuras de este artículo de JoVE

Solicitar permiso

Etiquetas

Tratamiento de hipertermiapresi n del fluido intersticialtransferencia de biocalorhipertermia por radiofrecuenciadin mica de fluidosperfusi n sangu neaecuaciones diferenciales parcialesdispositivo de radiofrecuencia bipolar

Artículos relacionados