Molecular dynamics simulation of the effect of initial temperature and initial pressure on thermal behavior and cavitation dynamics in tumor tissue

Resumen Técnico: Simulación de Dinámica Molecular de los Efectos de la Temperatura y Presión Iniciales en el Comportamiento Térmico del Tejido Tumoral y la Dinámica de Cavitación

Planteamiento del Problema
Las estrategias terapéuticas basadas en calor para el cáncer, como la hipertermia y la estimulación mecánica inducida por ultrasonido, dependen en gran medida de mecanismos térmicos, químicos y mecánicos acoplados. Si bien los estudios macroscópicos han establecido el papel de los campos de presión en la mejora de la penetración del calor y la inducción de la cavitación (la nucleación, crecimiento y colapso de microburbujas), los mecanismos fundamentales a nivel atómico que gobiernan el acoplamiento térmico-mecánico en los microentornos tumorales siguen estando insuficientemente comprendidos. Específicamente, los efectos combinados de las condiciones termodinámicas iniciales —es decir, la temperatura y la presión— sobre la dinámica de cavitación, la conductividad térmica local y la redistribución de energía no han sido explorados suficientemente. Los modelos existentes a escala de continuo y las técnicas experimentales carecen de la resolución necesaria para capturar las interacciones transitorias a escala atómica, las fluctuaciones locales de densidad y la evolución precisa de los mecanismos de colapso de las burbujas bajo estados iniciales variables.

Metodología
El estudio emplea simulaciones de Dinámica Molecular (MD) utilizando el Simulador de Masa Masivamente Paralelo Atómico/Molecular a Gran Escala (LAMMPS) para investigar un modelo de tejido tumoral de grano grueso.

• Construcción del Sistema: El dominio computacional (15 × 15 × 15 nm³) consiste en una matriz tumoral de tipo polímero semiestructurada sumergida en un entorno de agua TIP3P. Átomos de argón dispersos sirven como núcleos de cavitación (conglomerados de gas) para iniciar la formación de burbujas.
• Interacciones: Las fuerzas interatómicas se modelan mediante un marco híbrido que combina potenciales de Lennard–Jones (LJ) para las interacciones de Van der Waals y potenciales de Coulomb para las fuerzas electrostáticas. Se incluyen interacciones ligadas (estiramiento de enlaces, flexión de ángulos, diedros y términos impropios) para la matriz tumoral.
• Protocolo de Simulación:
  1. Equilibrio: El sistema se somete a una etapa de equilibrado de 10 ns en el ensamble canónico (NVT) utilizando un termostato de Nosé–Hoover para alcanzar un estado estable a 300 K.
  2. Producción: Le sigue una carrera de producción de 10 ns. El sistema se calienta a 320 K durante 5 ps mediante un termostato de Langevin. Se aplica una fuerza sinusoidal externa dependiente del tiempo (que simula la excitación ultrasónica) a la fase fluida durante los últimos 5 ps.
  3. Variables: El estudio varía sistemáticamente la temperatura inicial (280–310 K) y la presión inicial (0–2 bar).
• Análisis: La conductividad térmica se calcula mediante el formalismo de Green–Kubo basado en funciones de autocorrelación de flujo de calor. La evolución estructural se cuantifica mediante el radio de giro de la burbuja de argón, funciones de distribución radial (RDF) y perfiles de densidad/esfuerzo espacial.

Resultos Clave
Las simulaciones revelan respuestas distintas del modelo de tejido tumoral ante las variaciones de la temperatura y la presión iniciales:

• Efecto de la Temperatura Inicial (280–310 K):

  • Temperatura de Saturación: La temperatura de equilibrio final (saturación) se mantuvo constante en aproximadamente 320 K independientemente de la temperatura inicial, lo que indica una respuesta térmica estable e independiente de la historia térmica inicial.
  • Expansión Estructural: El aumento de la temperatura inicial promovió la expansión estructural. El radio de giro de la burbuja de argón aumentó monótonamente de 42.79 Å (a 280 K) a 42.89 Å (a 310 K), lo que refleja una mayor movilidad atómica y una reducción del confinamiento efectivo.
  • Conductividad Térmica: La conductividad térmica mostró una tendencia no monótona, variando entre 0.29 y 0.39 W/m·K. Aumentó de 0.30 a 0.36 W/m·K (280–290 K), cayó a 0.29 W/m·K a 300 K y alcanzó un máximo de 0.39 W/m·K a 310 K. Este comportamiento sugiere una competencia entre la mejora de la movilidad atómica (que facilita la transferencia de energía) y las restricciones estructurales (que limitan la eficiencia a temperaturas intermedias).
  • Esfuerzo y Densidad: Las temperaturas más altas indujeron una localización significativa del esfuerzo y una redistribución atómica en todo el dominio, aunque se mantuvo la integridad estructural general.

• Efecto de la Presión Inicial (0–2 bar):

  • Temperatura de Saturación: De manera similar a las variaciones de temperatura, la temperatura de equilibrio se mantuvo estable en ~320 K en todos los niveles de presión.
  • Respuesta Estructural: El radio de giro mostró una dependencia no monótona de la presión. Alcanzó su máximo a 1 bar (42.94 Å), disminuyendo ligeramente a 0.5 bar (42.87 Å) y 2 bar (42.84 Å). Esto identifica a 1 bar como la condición óptima para la expansión de la burbuja.
  • Conductividad Térmica: La conductividad térmica también mostró una tendencia no monótona, ascendiendo de 0.29 W/m·K (0 bar) a un máximo de 0.38 W/m·K (1 bar), para luego descender a 0.33 W/m·K (2 bar). Los autores atribuyen el pico a 1 bar a un equilibrio óptimo donde la reducción del espaciamiento interatómico mejora la fuerza de interacción sin la excesiva compactación estructural que restringe la movilidad atómica a presiones más altas (2 bar).

Significación y Reivindicaciones
El artículo afirma cerrar la brecha entre las observaciones a nivel de continuo y los mecanismos a nanoescala en el modelado de tejidos tumorales. Al resolver explícitamente las interacciones interatómicas, el estudio proporciona perspectivas a nivel atómico sobre cómo los estados termodinámicos iniciales influyen en las respuestas térmicas inducidas por la cavitación.

• Novedad: El trabajo vincula sistemáticamente los estados termodinámicos iniciales (temperatura y presión) con la dinámica de cavitación y la transferencia de calor en un entorno biológicamente relevante, una combinación no explorada previamente a la escala molecular.
• Implicaciones: La identificación de comportamientos no monótonos en la conductividad térmica y la expansión estructural sugiere que las estrategias terapéuticas que dependen de tratamientos asistidos por presión o ultrasonido requieren una optimización precisa de los parámetros. Específicamente, la existencia de una presión "óptima" (1 bar en este modelo) que maximiza tanto la expansión de la burbuja como la conductividad térmica implica que la eficacia terapéutica puede mejorarse ajustando estas condiciones iniciales para equilibrar el espaciamiento interatómico y la movilidad atómica.
• Limitaciones: El estudio está limitado a un modelo específico de grano grueso y escalas de tiempo de simulación (nanosegundos). No propone nuevos protocolos experimentales ni aplicaciones clínicas, sino que ofrece perspectivas fundamentales para informar la optimización de las terapias actuales basadas en calor.