성장 전이 종양 구체의 조직 고유 형태 역학

초록

본 연구는 3차원 정점 모델에 그래프 기반 폴리헤드 분할 알고리즘을 도입해 세포 분열을 구현하고, 표면 장력·생존‑괴사 경계 장력·증식 리밍 두께·활성 잡음이라는 네 가지 내재 기계 인자가 종양 구체의 구형 유지와 돌출형( lobulated) 변형을 어떻게 결정하는지 체계적으로 규명한다.

상세 분석

논문은 먼저 기존 3D 정점 모델이 세포 수 증가를 다루지 못하는 한계를 지적하고, 이를 극복하기 위해 Graph Vertex Model(GVM)에 폴리헤드 분할을 위한 그래프 변환 규칙을 추가한다. 분할 평면을 임의로 설정하고, 평면이 교차하는 모든 면·에지를 식별한 뒤, 패턴 매칭을 통해 해당 부분 그래프를 추출하고 노드·관계의 삭제·생성을 동시에 수행함으로써 두 자식 세포를 생성한다. 이 과정은 경계 세포와 내부 세포 모두에 적용 가능하도록 설계돼, 완전한 3D 토폴로지 변화를 자연스럽게 구현한다.

기계적 에너지 함수 W는 표면·내부 면적 A_mn에 장력 Γ_mn을 곱한 항과 부피 변형을 억제하는 압축 탄성항 κ_V(V_m‑V_0)^2 로 구성된다. 표면 장력 Γ와 살아‑괴사 경계(LNI) 장력 Γ_LNI 를 독립적으로 조정함으로써 종양 외부와 내부 인터페이스의 물리적 차이를 모델링한다. 영양소 확산·소모를 고려한 반응‑확산 방정식으로 영양소 농도 c_l을 층 번호 l에 따라 기하급수적으로 감소하도록 설정하고, 임계 농도 c_TH 를 초과하는 세포만이 분열하도록 한다. 이때 증식 리밍 두께 λ는 c_TH 와 직접 연관돼, λ가 클수록 증식 가능한 층이 늘어나며, 이는 실제 종양에서 관찰되는 살아있는 외피와 괴사 중심의 구조를 재현한다.

시뮬레이션은 초기 50세포 구체에서 시작해 2000세포까지 성장시키며, 각 파라미터 조합에 대해 50번의 독립 실행을 수행한다. 형태 지표로는 구형도( reduced volume v = 6√π V / A^{3/2})와 단면 불규칙도( reduced perimeter S = P_cs / √(4π A_cs))를 사용한다. 결과는 Γ_LNI 가 높을 때 LNI가 강한 장력 경계가 되어 표면 파동을 억제하고 구형을 유지한다는 것을 보여준다. 반면 Γ_LNI 가 낮고 λ 가 중간 정도일 때, LNI와 표면 사이의 기계적 결합이 약해져 내부 성장 압력이 표면에 전달되고, 이는 국소 팽창과 돌출을 야기해 v 가 감소하고 S 가 증가하는 lobulated 형태를 만든다. λ 가 매우 크면 전체 구체가 균일히 성장해 LNI 효과가 사라지면서 형태가 다시 평탄해진다.

또한, 활성 잡음( active tension fluctuations )을 τ_relaxation 시간 스케일로 모델링해 조직의 점탄성 특성을 도입하였다. 잡음 강도가 클수록 세포 간 장력 변동이 커져, 작은 불안정이 증폭돼 표면에 불규칙성이 급격히 증가한다. 이는 실제 종양에서 관찰되는 비정상적인 세포-세포 접촉과 ECM와의 동적 상호작용을 물리적으로 설명한다.

종합적으로, 논문은 “표면·LNI 장력, 증식 리밍 두께, 활성 잡음”이라는 세 축이 종양 구체의 형태 전이를 결정한다는 ‘조직 고유 형태 역학 삼위일체’를 제시한다. 이 삼요소는 서로 보완적으로 작용해, 특정 조합에서는 종양이 매끄러운 구형을 유지하고, 다른 조합에서는 조기 전이와 연관된 돌출형 형태로 변한다. 이러한 물리적 메커니즘은 기존의 분자·생물학적 설명을 보완하며, 치료적 표적( 예: 장력 조절, 성장 억제, 기계적 잡음 감소) 제시에도 기여한다.