플라나리아 장기의 자가조직 성장·분기·전신 스케일링 모델

초록

플라나리아의 장은 몸 크기에 따라 길이, 분기 수, 간격 등이 전형적인 멱법칙으로 스케일링한다. 저자들은 장 골격을 정량화하고, 형태 불안정과 모폴젠 확산‑반응을 결합한 2차원 인터페이스 모델을 제시한다. 모델은 곡률 억제, 방향장, 조직 성장률을 포함해 인터페이스 속도를 조절하고, 실험에서 관찰된 스케일링 지수를 재현한다. 이를 통해 복잡한 분기 구조가 자가조직적으로 형성되고, 전신 성장에 따라 전신적인 스케일링 법칙을 따르는 메커니즘을 제시한다.

상세 분석

본 논문은 두 가지 축을 중심으로 깊이 있게 전개된다. 첫 번째 축은 플라나리아 장기의 정량적 측정이다. 저자들은 전체 장 골격을 1픽셀 폭의 토폴로지 보존 스켈레톤으로 추출하고, 개별 분기 길이(ℓₙ), 분기 간 거리(d), 분기 각도(ψₙ) 등을 측정하였다. 이를 바탕으로 전체 장 길이(L_gut), 평균 분기 길이(ℓ), 분기 수(N), 장 면적(A_b) 등과 몸의 가로·세로 길이(L_x, L_y)와의 관계를 로그‑선형 회귀로 분석했다. 결과는 L_gut ∝ A_b^{0.75}, ℓ ∝ A_b^{0.19}, N ∝ A_b^{0.59}, d ∝ A_b^{0.24}와 같은 멱법칙 형태의 스케일링을 보였으며, 작은 개체에서는 주요 분기의 기하학적 편향으로 인해 약간의 비선형성이 관찰되었다. 또한 분기 방향은 cos²ψ의 이중 피크(수직·수평)로 나타났으며, 크기에 무관하게 일정한 방향성 분포를 유지한다는 점이 흥미롭다.

두 번째 축은 이론 모델이다. 저자들은 장의 성장 과정을 “in” 영역(장 조직)과 “out” 영역(기타 조직)을 구분하는 움직이는 인터페이스(R)로 추상화한다. 인터페이스는 법선 방향으로 ∂ₜR = (v_n + u_n) n 로 이동하며, 여기서 v_n은 모폴젠 농도와 곡률에 의해 조절되는 속도, u_n은 조직 전체의 팽창 속도이다. 구체적으로 v_n = χ(c,m) – βκ 로 정의하고, χ(c,m) = Γ(c)Θ(m) 로 분해한다. Γ(c) = v₀ – γc 로 선형 억제(또는 활성) 관계를, Θ(m) = 1 – δ(1 – n·m) 로 외부 방향장(m)의 영향을 반영한다. 곡률 항 βκ는 양의 곡률(돌출부)일수록 성장 속도를 감소시켜 인터페이스를 안정화한다. 모폴젠은 반응‑확산‑대류 방정식 ∂ₜc_i + u·∇c_i = D∇²c_i – (k_i + ∇·u)c_i + s_i 로 기술되며, 인터페이스에서 연속성과 연속적인 기울기를 만족한다.

이 모델을 성장하는 도메인(플라나리아 전체 크기 증가) 위에 수치적으로 구현하면, 초기 작은 원형 “장”이 모폴젠 농도 구배에 의해 국소적으로 빠르게 팽창하고, 곡률 억제와 방향장에 의해 불안정이 증폭되어 가지가 분기한다. 파라미터 γ, β, δ, D, k_i 등을 실험에서 얻은 스케일링 지수에 맞추어 조정하면, 시뮬레이션 결과가 L_gut, ℓ, N, d의 멱법칙 지수를 정량적으로 재현한다. 특히, 조직 성장률 u가 시간에 따라 선형적으로 증가할 경우, 인터페이스 불안정이 지속적으로 재점화되어 무한히 복잡한 분기 구조가 생성되는 것이 관찰된다. 이는 실제 플라나리아가 영양 상태에 따라 급격히 크기가 변하면서도 장 구조가 자동으로 재조정되는 현상을 설명한다.

핵심적인 통찰은 다음과 같다. (1) 형태 불안정은 단순히 물리적 현상이 아니라, 모폴젠 농도와 곡률, 방향성 신호가 결합된 복합 피드백 루프에 의해 조절된다. (2) 전신 성장(도메인 팽창)은 인터페이스 동역학에 시간‑스케일을 제공하여, 동일한 물리‑화학 파라미터라도 성장 단계에 따라 다른 스케일링 법칙을 나타낼 수 있다. (3) 모델은 특정 조직(플라나리아 장)뿐 아니라, 폐, 혈관, 유선 등 다른 계층적 분기 기관에도 적용 가능함을 시사한다.