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

기술 요약: 종양 조직의 열적 거동 및 공동 현상 역학에 미치는 초기 온도와 압력 효과에 대한 분자 동역학 시뮬레이션

문제 제기
고온 요법(hyperthermia) 및 초음파 유도 기계적 자극과 같은 열 기반 암 치료 전략은 열적, 화학적, 기계적 메커니즘의 결합된 역할에 크게 의존한다. 거시적 연구를 통해 압력장이 열 침투를 강화하고 공동 현상(cavitation, 미세 기포의 핵 생성, 성장 및 붕괴)을 유도하는 역할이 확립되었으나, 종양 미세환경 내에서 열-기계적 결합을 지배하는 근본적인 원자 단위의 메커σ은 여전히 제대로 이해되지 않고 있다. 특히, 초기 열역학적 조건, 즉 온도와 압력이 공동 현상 역학, 국소 열전도도 및 에너지 재분배에 미치는 결합된 영향은 충분히 탐구되지 않았다. 기존의 연속체 규모 모델과 실험 기법은 다양한 초기 상태 하에서의 과도한 원자 규모 상호작용, 국소 밀도 변동 및 기포 붕괴 메커니즘의 정밀한 진화를 포착할 수 있는 해상도가 부족하다.

방법론
본 연구는 거친 입자(coarse-grained) 종양 조직 모델을 조사하기 위해 LAMMPS(Large-scale Atomic/Molecular Massively Parallel Simulator)를 이용한 분자 동역학(MD) 시뮬레이션을 채택하였다.

• 시스템 구축: 계산 도메인(15 × 15 × 15 nm³)은 TIP3P 물 환경에 담긴 반구조적 고분자 형태의 종양 매트릭스로 구성된다. 공동 현상의 핵(가스 클러스터) 역할을 하기 위해 분산된 아르곤 원자가 사용되었다.
• 상호작용: 원자 간 힘은 반데르발스 상호작용을 위한 레나드-존(Lennard–Jones, LJ) 퍼텐셜과 정전기적 힘을 위한 쿨롱(Coulombic) 퍼텐셜을 결합한 하이브리드 프레임워크를 사용하여 모델링되었다. 종양 매트릭스를 위한 결합 상호작용(결합 신축, 각도 굽힘, 이면각 및 부적절한 항)이 포함되었다.
• 시뮬레이션 프로토콜:
  1. 평형화(Equilibration): 시스템은 300 K에서 안정적인 상태에 도달하기 위해 노제-후버(Nosé–Hoover) 열조절기를 사용하는 정준(NVT) 앙상블에서 10 ns 동안 평형화 단계를 거친다.
  2. 생성(Production): 10 ns의 생성 실행이 뒤따른다. 시스템은 5 ps 동안 란제빈(Langevin) 열조절기를 통해 320 K로 가열된다. 마지막 5 ps 동안에는 유체상에 외부 시간 의존적 사인파 형태의 힘(초음파 자극 시뮬레이션)이 가해진다.
  3. 변수: 본 연구는 초기 온도(280–310 K)와 초기 압력(0–2 bar)을 체계적으로 변화시킨다.
• 분석: 열전도도는 열속 자기상관 함수에 기반한 그린-쿠보(Green–Kubo) 형식을 사용하여 계산된다. 구조적 진화는 아르곤 기포의 회전 반경(radius of gyration), 방사 분포 함수(RDF), 그리고 공간 밀도/응력 프로파일을 통해 정량화된다.

주요 결과
시뮬레이션은 종양 조직 모델이 초기 온도 및 압력 변화에 대해 뚜렷한 반응을 보임을 밝혀냈다:

• 초기 온도(280–310 K)의 영향:

  • 포화 온도: 최종 평형(포화) 온도는 초기 온도와 관계없이 약 320 K로 일정하게 유지되었으며, 이는 초기 열 이력과 무관한 안정적인 열적 응답을 나타낸다.
  • 구조적 팽창: 초기 온도가 높아질수록 구조적 팽창이 촉진되었다. 아르곤 기포의 회전 반경은 42.79 Å(280 K일 때)에서 42.89 Å(310 K일 때)까지 단조 증가하였으며, 이는 향상된 원자 이동성과 감소된 유효 구속을 반영한다.
  • 열전도도: 열전도도는 0.29에서 0.39 W/m·K 사이에서 변하며 비단조적인 경향을 보였다. 온도가 280 K에서 290 K로 상승할 때 0.30에서 0.36 W/m·K로 증가하였고, 300 K에서 0.29 W/m·K로 떨어졌으며, 310 K에서 0.39 W/m·K로 정점에 도달했다. 이러한 거동은 향상된 원자 이동성(에너지 전달 용이)과 구조적 제약(중간 온도에서의 효율 제한) 사이의 경쟁을 시사한다.
  • 응력 및 밀도: 높은 온도는 전체 도메인에 걸쳐 상당한 응력 국부화와 원자 재분포를 유도하였으나, 전반적인 구조적 무결성은 유지되었다.

• 초기 압력(0–2 bar)의 영향:

  • 포화 온도: 온도 변화와 마찬가지로, 평형 온도는 모든 압력 수준에서 ~320 K로 안정적이었다.
  • 구조적 반응: 회전 반경은 압력에 대해 비단조적인 의존성을 보였다. 회전 반경은 1 bar(42.94 Å)에서 정점을 찍었으며, 0.5 bar(42.87 Å)와 2 bar(42.84 Å)에서는 약간 감소하였다. 이는 1 bar가 기포 팽창을 위한 최적의 조건임을 식별한다.
  • 열전도도: 열전도도 역시 비단조적인 경향을 보이며, 0 bar(0.29 W/m·K)에서 1 bar(0.38 W/m·K)로 상승한 후 2 bar(0.33 W/m·K)로 하락하였다. 저자들은 1 bar에서의 정점을 원자 간 간격이 감소하여 상호작용 강도를 높이는 동시에, 높은 압력(2 bar)에서 원자 이동성을 제한하는 과도한 구조적 압착을 일으키지 않는 최적의 균형이 이루어진 결과로 설명한다.

의의 및 주장
본 논문은 종양 조직 모델링에 있어 연속체 수준의 관찰과 나노 규모의 메커니즘 사이의 간극을 메운다고 주장한다. 원자 간 상호작용을 명시적으로 해결함으로써, 본 연구는 초기 열역학적 상태가 공동 현상 유도 열 응답에 어떻게 영향을 미치는지에 대한 원자 수준의 통찰을 제공한다.

• 독창성: 본 연구는 초기 열역학적 상태(온도 및 압력)를 공동 현상 역학 및 열 전달과 체계적으로 연결하였으며, 이는 분자 규모에서 이전에 탐구되지 않았던 조합이다.
• 시사점: 열전도도와 구조적 팽창의 비단조적 거동을 확인한 것은 압력 보조 또는 초음파 기반 치료에 의존하는 치료 전략이 정밀한 파라미터 최적화를 필요로 함을 시사한다. 특히, 기포 팽창과 열전도도를 모두 극대화하는 "최적" 압력(본 모델에서 1 bar)의 존재는, 원자 간 간격과 원자 이동성의 균형을 맞추기 위해 초기 조건을 조절함으로써 치료 효능을 높일 수 있음을 의미한다.
• 한계: 본 연구는 특정 거친 입자 모델과 시뮬레이션 시간 규모(나노초)에 국한되어 있다. 새로운 실험 프로토콜이나 임상 적용을 제안하는 것이 아니라, 기존의 열 기반 암 치료 최적화를 돕기 위한 근본적인 통찰을 제공하는 데 목적이 있다.