시간 지연과 영역 크기가 만든 2차원 패턴 형성의 복합 동역학

초록

본 연구는 유전자 발현 시간 지연과 도메인 크기가 반응‑확산 시스템의 투링 패턴 형성에 미치는 영향을 조사한다. 수정된 Schnakenberg 모델인 Ligand Internalisation(LI) 모델을 2차원으로 확장하여 선형 안정성 분석과 수치 시뮬레이션을 수행하였다. 시간 지연이 증가하면 패턴이 나타나는 시간이 거의 선형적으로 늘어나며, 영역 크기에 대해서는 패턴 형성 시간이 비단조적(non‑monotonic)으로 변한다. 다양한 초기 조건(무작위 섭동 및 불안정 매니폴드)과 두 파라미터(시간 지연, 영역 크기) 차트가 제시된다.

상세 분석

본 논문은 기존 1차원 연구에서 관찰된 시간 지연(τ)과 패턴 형성 시간 사이의 선형 관계를 2차원 영역으로 일반화한다는 점에서 의미가 크다. LI 모델은 Schnakenberg 모델에 활성자(u)의 내부화 과정을 지연항(τ)으로 도입한 변형형이며, 반응항에 ˆu(t‑τ), ˆv(t‑τ) 형태로 과거 상태가 포함된다. 모델식(2)은 두 확산 계수 d_u=0.01, d_v=0.2을 고정하고, 도메인 크기 L_x, L_y를 비정규화된 스케일링 파라미터로 사용한다.

선형 안정성 분석에서는 균일 평형 E*=(a+b, b/(a+b)^2)를 기준으로 작은 섭동 ζ, η를 도입하고, 지연을 포함한 선형화된 시스템(6)을 유도한다. 고유모드 형태 e^{λt}cos(k_xπx/L_x)cos(k_yπy/L_y) 를 가정하면 특성 방정식(9)‑(10)이 얻어지며, 여기서 e^{−λτ} 항이 지연의 핵심 효과를 담당한다. 파라미터 집합 (a, b, d_u, d_v, τ, L_x, L_y) 에 대해 실수부 최대값 α(τ) 를 정의하고, α>0이면 투링 불안정이 발생한다는 전통적인 기준을 유지한다.

분석 결과는 다음과 같다. 첫째, τ가 증가하면 λ의 실수부가 감소하여 α가 점진적으로 작아지며, 결국 α가 0을 넘는 최소 τ값이 존재한다. 이는 패턴이 나타나기까지 걸리는 시간이 τ에 비례한다는 선형 관계를 수치적으로 확인한다. 둘째, L_x와 L_y를 변화시켰을 때 α는 단조적으로 변하지 않는다. 특정 영역 크기에서는 가장 큰 성장률을 보이다가, 크기가 더 커지거나 작아지면 성장률이 감소한다는 비단조적(역 U자형) 곡선을 만든다. 이는 파동수(k_x, k_y)와 도메인 고유주기가 일치하거나 불일치할 때 발생하는 모드 선택 현상과 연관된다.

또한 저자는 초기 조건에 따라 패턴 전이 경로가 달라짐을 보여준다. 무작위 섭동은 대부분의 경우 가장 불안정한 모드가 빠르게 지배하게 하여 전형적인 스트립이나 점 패턴을 생성한다. 반면, 불안정 매니폴드에서 추출한 초기 조건은 특정 모드가 장기적으로 유지되도록 하여, 동일한 파라미터 집합에서도 서로 다른 패턴(예: 격자형, 원형)으로 수렴한다. 이러한 초기 조건 의존성은 실제 생물학적 시스템에서 초기 세포 배치나 외부 자극이 최종 형태에 미치는 영향을 시사한다.

마지막으로, 두 파라미터(τ, L) 차트는 투링 영역을 시각화한다. τ가 작을 때는 넓은 (a, b) 영역에서 불안정이 발생하지만, τ가 커질수록 불안정 영역이 축소되고, 특정 L값에서는 불안정 영역이 다시 확대되는 현상이 관찰된다. 이는 시간 지연과 공간 규모가 상호 보완적으로 작용해 패턴 형성의 임계 조건을 조절한다는 중요한 통찰을 제공한다.

요약하면, 이 연구는 시간 지연과 도메인 크기가 투링 패턴의 발생 시점과 형태를 복합적으로 조절한다는 사실을 정량적·정성적으로 입증했으며, 2차원 LI 모델을 통한 분석 프레임워크는 향후 복잡한 조직 형성 모델에 적용 가능성을 열어준다.