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

Resumo Técnico: Simulação de Dinâmica Molecular dos Efeitos da Temperatura e Pressão Iniciais no Comportamento Térmico de Tecido Tumoral e Dinâmica de Cavitação

Definição do Problema
As estratégias terapêuticas baseadas em calor para o câncer, como a hipertermia e a estimulação mecânica induzida por ultrassom, dependem fortemente de mecanismos térmicos, químicos e mecânicos acoplados. Embora estudos macroscópicos tenham estabelecido o papel dos campos de pressão no aumento da penetração de calor e na indução de cavitação (a nucleação, crescimento e colapso de microbolhas), os mecanismos fundamentais em escala atômica que regem o acoplamento térmico-mecânico em microambientes tumorais permanecem mal compreendidos. Especificamente, os efeitos combinados das condições termodinâmicas iniciais — ou seja, temperatura e pressão — sobre a dinâmica de cavitação, a condutividade térmica local e a redistribuição de energia são insuficientemente explorados. Modelos de escala contínua existentes e técnicas experimentais carecem de resolução para capturar interações transientes em escala atômica, flutuações de densidade local e a evolução precisa dos mecanismos de colapso de bolhas sob variações de estados iniciais.

Metodologia
O estudo utiliza simulações de Dinâmica Molecular (MD) usando o Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) para investigar um modelo de tecido tumoral de granulometria grosseira (coarse-grained).

• Construção do Sistema: O domínio computacional (15 × 15 × 15 nm³) consiste em uma matriz tumoral semiestruturada do tipo polimérico imersa em um ambiente de água TIP3P. Átomos de argônio dispersos servem como núcleos de cavitação (aglomerados gasosos) para iniciar a formação de bolhas.
• Interações: As forças interatômicas são modeladas usando uma estrutura híbrida combinando potenciais de Lennard–Jones (LJ) para interações de van der Waals e potenciais Coulombianos para forças eletrostáticas. Interações ligadas (estiramento de ligação, flexão de ângulo, diedro e termos ímpares) estão incluídas para a matriz tumoral.
• Protocolo de Simulação:
  1. Equilibração: O sistema passa por uma etapa de equilibração de 10 ns no ensemble canônico (NVT) usando um termostato de Nosé–Hoover para atingir um estado estável a 300 K.
  2. Produção: Uma corrida de produção de 10 ns se segue. O sistema é aquecido para 320 K ao longo de 5 ps usando um termostato de Langevin. Uma força senoidal externa dependente do tempo (simulando excitação ultrassônica) é aplicada à fase fluida para os 5 ps finais.
  3. Variáveis: O estudo varia sistematicamente a temperatura inicial (280–310 K) e a pressão inicial (0–2 bar).
• Análise: A condutividade térmica é calculada usando o formalismo de Green–Kubo baseado em funções de autocorrelação de fluxo de calor. A evolução estrutural é quantificada via raio de giração da bolha de argônio, funções de distribuição radial (RDF) e perfis de densidade espacial/tensão.

Resultos Principais
As simulações revelam respostas distintas do modelo de tecido tumoral às variações de temperatura e pressão iniciais:

• Efeito da Temperatura Inicial (280–310 K):

  • Temperatura de Saturação: A temperatura de equilíbrio final (saturação) permaneceu constante em aproximadamente 320 K, independentemente da temperatura inicial, indicando uma resposta térmica estável e independente do histórico térmico inicial.
  • Expansão Estrutural: O aumento da temperatura inicial promoveu a expansão estrutural. O raio de giração da bolha de argônio aumentou monotonicamente de 42,79 Å (a 280 K) para 42,89 Å (a 310 K), refletindo maior mobilidade atômica e redução do confinamento efetivo.
  • Condutividade Térmica: A condutividade térmica exibiu uma tendência não monotônica, variando entre 0,29 e 0,39 W/m·K. Ela aumentou de 0,30 para 0,36 W/m·K (280–290 K), caiu para 0,29 W/m·K a 300 K e atingiu o pico de 0,39 W/m·K a 310 K. Este comportamento sugere uma competição entre o aumento da mobilidade atômica (facilitando a transferência de energia) e restrições estruturais (limitando a eficiência em temperaturas intermediárias).
  • Tensão e Densidade: Temperaturas mais altas induziram localização significativa de tensão e redistribuição atômica em todo o domínio, embora a integridade estrutural geral tenha sido mantida.

• Efeito da Pressão Inicial (0–2 bar):

  • Temperatura de Saturação: Semelhante às variações de temperatura, a temperatura de equilíbrio permaneceu estável em ~320 K em todos os níveis de pressão.
  • Resposta Estrutural: O raio de giração mostrou uma dependência não monotônica em relação à pressão. Ele atingiu o pico em 1 bar (42,94 Å), diminuindo ligeiramente a 0,5 bar (42,87 Å) e 2 bar (42,84 Å). Isso identifica 1 bar como a condição ideal para a expansão da bolha.
  • Condutividade Térmica: A condutividade térmica também exibiu uma tendência não monotônica, subindo de 0,29 W/m·K (0 bar) para um máximo de 0,38 W/m·K (1 bar), e declinando para 0,33 W/m·K (2 bar). Os autores atribuem o pico em 1 bar a um equilíbrio ideal onde o espaçamento interatômico reduzido aumenta a força de interação sem a compactação estrutural excessiva que restringe a mobilidade atômica em pressões mais altas (2 bar).

Significância e Alegações
O artigo afirma preencher a lacuna entre observações de nível contínuo e mecanismos em nanoescala no modelagem de tecido tumoral. Ao resolver explicitamente as interações interatômicas, o estudo fornece insights em nível atômico sobre como os estados termodinâmicos iniciais influenciam as respostas térmicas induzidas por cavitação.

• Novidade: O trabalho vincula sistematicamente os estados termodinâmicos iniciais (temperatura e pressão) com a dinâmica de cavitação e a transferência de calor em um ambiente biologicamente relevante, uma combinação não explorada anteriormente em escala molecular.
• Implicações: A identificação de comportamentos não monotônicos na condutividade térmica e na expansão estrutural sugere que as estratégias terapêuticas que dependem de tratamentos assistidos por pressão ou ultrassom requerem otimização precisa de parâmetros. Especificamente, a existência de uma pressão "ideal" (1 bar neste modelo) que maximiza tanto a expansão da bolha quanto a condutividade térmica implica que a eficácia terapêutica pode ser aumentada ao ajustar essas condições iniciais para equilibrar o espaçamento interatômico e a mobilidade atômica.
• Limitações: O estudo está confinado a um modelo específico de granulometria grosseira (coarse-grained) e escalas de tempo de simulação (nanossegundos). Não propõe novos protocolos experimentais ou aplicações clínicas, mas oferece insights fundamentais para informar a otimização de terapias atuais baseadas em calor.