종양 혈관의 꼬불꼬불 형성: 형태‑탄성 모델

초록

본 논문은 종양 모세관을 성장하는 초탄성 원통으로 모델링하고, 주변 조직을 선형 탄성 스프링 기초로 가정한다. 성장 텐서는 축방향으로만 비균일하게 적용되며, 이를 통해 축성장 파라미터 γ가 일정 임계값을 초과하면 원통이 토루스(굽힘) 형태로 전이한다. 선형 안정성 분석은 증분 변형을 유한 변형 위에 겹쳐 풀고, Stroh 형식과 임피던스 매트릭스 방법으로 임계 모드와 임계 γ를 계산한다. 비선형 후버클링은 혼합 유한요소법으로 시뮬레이션하였다. 결과는 주변 조직의 탄성 강성 μ_k와 기하학적 비율 α_R이 토루스 발생에 중요한 역할을 함을 보여준다.

상세 분석

이 연구는 종양 혈관의 비정상적인 고곡률(토루스)을 형태‑탄성 이론으로 설명한다. 기본 가정은 모세관 벽이 불압축성 네오후크(neo‑Hookean) 초탄성 물질이며, 성장은 변형 그래디언트 F를 성장 텐서 G와 탄성 변형 F_e의 곱으로 분해하는 곱분해법(multiplicative decomposition)으로 기술된다. 여기서 G=diag(1,1,γ) 로 설정해 축방향 성장만을 고려한다. 기본 축대칭 해는 반지름 r(R)와 압축성 라그랑지 승수 p(r)으로 명시되며, 내부·외부 경계조건은 무전단(내부)와 선형 스프링(외부)으로 정의된다.

선형 안정성 분석에서는 기존의 유한 변형 위에 작은 증분 변위 δu를 도입하고, 네오후크 물질의 즉시 탄성 모듈러스 텐서 A₀를 이용해 증분 피올라‑키르히호프 응력 δP₀를 전개한다. 연속성 방정식과 운동량 균형을 결합한 PDE 시스템을 Stroh 형식으로 변환하면 6×6 행렬 N이 등장하고, 이를 Riccati 형태의 임피던스 매트릭스 Z(r)에 대한 미분 방정식으로 정리한다. 초기 조건 Z(r_i)=0과 외부 경계조건 (T+μ_kR₀I)U=0을 적용해 det(Z+μ_kR₀I)=0이 만족될 때 비자명 해가 존재함을 확인한다. 이 조건을 수치적으로 풀어 γ와 파동수 k, 원주 방향 모드 m에 대한 임계 곡선을 도출한다.

결과는 μ_k=0(즉, 자유‑버클링)일 때는 고전적인 Euler‑버클링과 동일하게 γ≈1에서 m=1, k→0인 장축 모드가 가장 먼저 불안정해짐을 보여준다. 그러나 μ_k>0인 경우, 즉 주변 조직이 스프링 기초로 작용하면 장축 파동수 k가 작아질수록 γ_cr이 급격히 증가한다. 이는 주변 조직의 저항이 축방향 압축을 억제해 토루스 발생을 지연시키는 메커니즘을 의미한다. 또한 기하학적 비율 α_R=R_i/R_o가 작을수록(즉, 얇은 벽) 불안정성 임계가 낮아져 토루스가 쉽게 발생한다.

비선형 후버클링 단계에서는 혼합 유한요소법을 사용해 성장 γ를 점진적으로 증가시키며, 완전 비선형 평형을 계산한다. 시뮬레이션 결과는 초기 선형 분석에서 예측된 모드와 일치하면서도, 큰 γ에서 복합적인 파동수와 다중 모드가 결합된 복잡한 토루스 형태가 나타난다. 이는 실제 종양 혈관이 보이는 비정상적인 곡률과 관류 저항 증가를 설명하는 데 충분히 현실적인 메커니즘을 제공한다.

전반적으로 이 논문은 성장‑구속 상호작용, 물질 비선형성, 그리고 주변 조직의 탄성 기초가 결합된 형태‑탄성 모델을 통해 종양 모세관의 토루스 형성을 정량적으로 예측하고, 치료 약물 전달 효율을 저해하는 물리적 원인을 이론적으로 규명한다.