이방성 성장에서 나타나는 새로운 부문 형태

초록

이 논문은 콜로니 성장의 이방성 효과를 KPZ와 FKPP 방정식으로 결합한 모델로 분석한다. 파라미터 β와 λ의 비가 1을 벗어나면 새로운 ‘충격파’ 형태의 부문이 나타나며, 이는 이방성 성장의 특징적인 신호이다.

상세 분석

본 연구는 미생물 군집이나 종양과 같은 2차원 확장 현상을 ‘thin‑edge’ 근사법으로 단순화하여, 전면 위치 h(x,t)와 유전형 비율 f(x,t) 두 개의 1차원 변수만으로 기술한다. 전면의 움직임은 수정된 Kardar‑Parisi‑Zhang(KPZ) 방정식(∂ₜh = v + (λ/2)(∂ₓh)² + D_h∂ₓ²h)으로, 여기서 λ는 전면 기울기에 따른 비선형 성장 속도를, v는 평면 전면의 기본 팽창 속도를 나타낸다. 동시에 두 균주 사이의 경쟁은 수정된 Fisher‑Kolmogorov‑Petrovsky‑Piskunov(FKPP) 방정식(∂ₜf = s₀f(1−f) + D_f∂ₓ²f + β∂ₓh ∂ₓf)으로 기술되며, β는 전면 기울기가 유전형 전파에 미치는 영향을 매개한다.

이방성은 전면 속도 v(θ)가 전면의 진행 방향 θ에 의존함으로써 도입된다. θ≈0(가장 느린 성장 방향) 근처에서 v(θ)≈v₀ + (v″(0)/2)θ² 로 전개하면, KPZ 방정식의 비선형 계수 λ가 v₀ → v₀+v″(0) 로 변하고, FKPP 방정식의 β는 v₀(1+ε′(0)) 로 바뀐다. 여기서 ε(θ)≈ε′(0)θ는 이방성에 의해 부문 경계가 전면 법선과 이루는 각도이다. 따라서 isotropic 경우에는 λ=β=v₀가 되지만, 이방성에서는 λ와 β가 서로 독립적인 값이 될 수 있다. 두 파라미터의 비 β/λ가 1이 아닌 경우가 바로 ‘비자명 이방성’이며, 이는 단순한 좌표 재스케일링으로는 제거되지 않는다.

모델의 차원 없는 형태는 D_h/D_f, λ/β, αβ/(s₀D_f) 세 개의 독립 파라미터로 축소된다. 여기서 s₀는 변이의 선택적 이점, α는 두 균주의 팽창 속도 차이, β는 전면 기울기에 의한 전파 편향을 의미한다. 분석 결과, β/λ가 임계값을 초과하면 기존의 부문(전면에 수직으로 뻗는 원형 혹은 둥근 형태)에서 탈피하여, ‘충격파’라 부르는 앞쪽으로 돌출된 kinks가 형성된다. 이 충격파는 피격된 균주가 전면을 앞서며 얇은 파동 형태로 퍼지는 현상으로, 전형적인 isotropic 모델에서는 관찰되지 않는다.

수치 시뮬레이션은 두 가지 수준에서 수행되었다. 첫 번째는 위에서 제시한 연속 방정식들을 격자화한 간소화 모델이며, 두 번째는 반응‑확산 방정식으로 구현한 보다 현실적인 세포 성장 모델이다. 두 모델 모두 β/λ가 1보다 크게 될 때, 전면 앞쪽에 뚜렷한 kinks가 나타나는 것을 확인하였다. 특히, 강한 이방성(β/λ≫1)에서는 충격파가 거의 직선 형태로 전파되어, 피격 균주가 전면을 ‘앞서’는 듯한 시각적 효과를 만든다. 이러한 현상은 실험적으로는 컬라겐 섬유나 화학 구배 등으로 유도된 방향성 성장 환경에서 관찰될 가능성이 있다.

이 연구는 기존의 isotropic 확장 모델이 놓치고 있던 ‘전면 기울기에 의한 경쟁 편향’이라는 메커니즘을 명확히 제시한다. β와 λ를 독립적인 파라미터로 두어 이방성을 정량화함으로써, 실험적 시스템에서 관찰되는 비대칭 부문 형태를 이론적으로 예측하고 해석할 수 있다. 또한, 충격파 형태의 부문은 이방성 성장의 ‘시그니처’로 활용될 수 있어, 미생물 군집이나 종양 조직에서 환경적 이방성을 추정하는 새로운 도구가 될 전망이다.