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

Résumé technique : Simulation de dynamique moléculaire des effets de la température et de la pression initiales sur le comportement thermique et la dynamique de cavitation d'un tissu tumoral

Énoncé du problème
Les stratégies thérapeutiques basées sur la chaleur pour le cancer, telles que l'hyperthermie et la stimulation mécanique induite par ultrasons, reposent largement sur des mécanismes thermiques, chimiques et mécaniques couplés. Bien que les études macroscopiques aient établi le rôle des champs de pression dans l'amélioration de la pénétration de la chaleur et l'induction de la cavitation (la nucléation, la croissance et l'effondrement de micro-bulles), les mécanismes fondamentaux à l'échelle atomique régissant le couplage thermo-mécanique dans les microenvironnements tumoraux restent mal compris. Plus précisément, les effets combinés des conditions thermodynamiques initiales — à savoir la température et la pression — sur la dynamique de cavitation, la conductivité thermique locale et la redistribution de l'énergie sont insuffisamment explorés. Les modèles existants à l'échelle du continuum et les techniques expérimentales manquent de la résolution nécessaire pour capturer les interactions transitoires à l'échelle atomique, les fluctuations de densité locale et l'évolution précise des mécanismes d'effondrement des bulles sous diverses conditions initiales.

Méthodologie
L'étude emploie des simulations de dynamique moléculaire (MD) utilisant le Large-scale Atomic/Molecular Massively Parallel Simulator (LAMPS) pour étudier un modèle de tissu tumoral à gros grains (coarse-grained).

• Construction du système : Le domaine de calcul (15 × 15 × 15 nm³) est une matrice tumorale de type polymère semi-structurée immergée dans un environnement d'eau TIP3P. Des atomes d'argon dispersés servent de noyaux de cavitation (clusters de gaz) pour initier la formation de bulles.
• Interactions : Les forces interatomiques sont modélisées à l'aide d'un cadre hybride combinant des potentiels de Lennard–Jones (LJ) pour les interactions de van der Waals et des potentiels coulombiques pour les forces électrostatiques. Des interactions liées (étirement de liaison, flexion d'angle, dièdre et termes impropres) sont incluses pour la matrice tumorale.
• Protocole de simulation :
  1. Équilibrage : Le système subit une étape d'équilibrage de 10 ns dans l'ensemble canonique (NVT) en utilisant un thermostat de Nosé–Hoover pour atteindre un état stable à 300 K.
  2. Production : Une phase de production de 10 ns suit. Le système est chauffé à 320 K sur 5 ps à l'aide d'un thermostat de Langevin. Une force sinusoïdale externe dépendante du temps (simulant une excitation ultrasonore) est appliquée à la phase fluide pour les 5 dernières ps.
  3. Variables : L'étude fait varier systématiquement la température initiale (280–310 K) et la pression initiale (0–2 bar).
• Analyse : La conductivité thermique est calculée en utilisant le formalisme de Green–Kubo basé sur les fonctions d'autocorrélation du flux de chaleur. L'évolution structurelle est quantifiée via le rayon de giration de la bulle d'argon, les fonctions de distribution radiale (RDF) et les profils de densité et de contrainte spatiale.

Résultats clés
Les simulations révèlent des réponses distinctes du modèle de tissu tumoral aux variations de la température et de la pression initiales :

• Effet de la température initiale (280–310 K) :

  • Température de saturation : La température d'équilibre finale (saturation) est restée constante à environ 320 K, quelle que soit la température initiale, indiquant une réponse thermique stable indépendante de l'historique thermique initial.
  • Expansion structurelle : L'augmentation de la température initiale a favorisé l'expansion structurelle. Le rayon de giration de la bulle d'argons a augmenté de manière monotone de 42,79 Å (à 280 K) à 42,89 Å (à 310 K), reflétant une mobilité atomique accrue et un confinement effectif réduit.
  • Conductivité thermique : La conductivité thermique a présenté une tendance non monotone, variant entre 0,29 et 0,39 W/m·K. Elle est passée de 0,30 à 0,36 W/m·K (280–290 K), a chuté à 0,29 W/m·K à 300 K, et a atteint un pic de 0,39 W/m·K à 310 K. Ce comportement suggère une compétition entre la mobilité atomique accrue (facilitant le transfert d'énergie) et les contraintes structurelles (limitant l'efficacité à des températures intermédiaires).
  • Contrainte et densité : Les températures plus élevées ont induit une localisation importante des contraintes et une redistribution atomique à travers tout le domaine, bien que l'intégrité structurelle globale ait été maintenue.

• Effet de la pression initiale (0–2 bar) :

  • Température de saturation : Similairement aux variations de température, la température d'équilibre est restée stable à ~320 K pour tous les niveaux de pression.
  • Réponse structurelle : Le rayon de giration a montré une dépendance non monotone vis-à-vis de la pression. Il a culminé à 1 bar (42,94 Å), diminuant légèrement à 0,5 bar (42,87 Å) et 2 bar (42,84 Å). Cela identifie 1 bar comme la condition optimale pour l'expansion de la bulle.
  • Conductivité thermique : La conductivité thermique a également affiché une tendance non monotone, augmentant de 0,29 W/m·K (0 bar) jusqu'à un maximum de 0,38 W/m·K (1 bar), puis diminuant à 0,33 W/m·K (2 bar). Les auteurs attribuent le pic à 1 bar à un équilibre optimal où la réduction de l'espacement interatomique renforce la force d'interaction sans l'excès de compaction structurelle qui restreint la mobilité atomique à des pressions plus élevées (2 bar).

Signification et revendications
L'article affirme combler le fossé entre les observations au niveau du continuum et les mécanismes à l'échelle nanométrique dans la modélisation des tissus tumoraux. En résolvant explicitement les interactions interatomiques, l'étude fournit des informations à l'échelle atomique sur la manière dont les états thermodynamiques initiaux influencent les réponses thermiques induites par la cavitation.

• Nouveauté : Le travail lie systématiquement les états thermodynamiques initiaux (température et pression) avec la dynamique de cavitation et le transfert de chaleur dans un environnement biologiquement pertinent, une combinaison qui n'avait pas été explorée auparavant à l'échelle moléculaire.
• Implications : L'identification de comportements non monotones de la conductivité thermique et de l'expansion structurelle suggère que les stratégies thérapeutiques reposant sur des traitements assistés par pression ou par ultrasons nécessitent une optimisation précise des paramètres. Plus précisément, l'existence d'une pression « optimale » (1 bar dans ce modèle) qui maximise à la fois l'expansion de la bulle et la conductivité thermique implique que l'efficacité thérapeutique peut être améliorée en ajustant ces conditions initiales pour équilibrer l'espacement interatomique et la mobilité atomique.
• Limites : L'étude est confinée à un modèle spécifique à gros grains et à des échelles de temps de simulation (nanosecondes). Elle ne propose pas de nouveaux protocoles expérimentaux ou d'applications cliniques, mais offre des informations fondamentales pour éclairer l'optimisation des thérapies anticancéreuses existantes basées sur la chaleur.