양성 피드백이 세포막에서 활성 신호분자 클러스터와 전파 속도에 미치는 영향

초록

본 연구는 양성 피드백과 느린 확산이 결합될 때 Ras와 유사한 신호분자가 세포막에서 클러스터를 형성하고, 그 경계가 일정한 속도로 전파되는 메커니즘을 단순 모델과 Monte Carlo 시뮬레이션으로 분석한다. 클러스터 전파 전선은 KPZ 방정식과 일치하는 거칠기를 보이며, 임계 핵 크기와 성장 법칙이 확산계수와 피드백 강도에 의존함을 확인하였다.

상세 분석

이 논문은 세포막에서 일어나는 신호전달을 Ras‑GTP 전환 과정에 빗대어, 세 종류의 입자(Z = Ras‑GDP, X = Ras‑GTP, Y = SOS)와 세 가지 반응(기본 활성화, 양성 피드백 활성화, SOS 탈착)으로 구성된 최소 모델을 제시한다. 핵심 가정은 막면에서의 확산이 반응 속도와 비슷하거나 느리다는 점이며, 이는 실제 막단백질의 확산계수가 0.01–1 µm²/s 수준인 점과 일치한다. 모델은 격자 기반의 kinetic Monte Carlo 알고리즘으로 구현되었으며, 입자 간의 하드코어 배제와 동일한 격자점당 최대 하나의 입자 제한을 둔다. 반응 확률은 실제 물리적 속도 상수와 격자 간격·시간 스케일을 연결하는 변환식(k₁ = 2k₁^MC dl₀^d 등)으로 정의된다.

시뮬레이션 결과는 차원에 따라 전혀 다른 동역학을 보인다. 1차원에서는 Y 입자가 X 양쪽에 갇히면 영구적으로 비활성화돼 시스템이 정지 상태에 빠진다. 반면 2차원에서는 초기의 무작위 Z‑Y 혼합에서 X 입자가 점차 생성되고, 양성 피드백에 의해 일정 크기 이상의 클러스터가 ‘임계 핵’으로 성장한다. 핵의 크기 l_c는 확산계수 D와 피드백 반응 상수 k₂, 그리고 평균 농도 ρ_Y·ρ_Z에 따라 l_c ≈ √(D/(k₂ρ_Yρ_Z)) 로 추정되지만, 실제 시뮬레이션에서는 비선형성 및 통계적 변동 때문에 이 식과는 약간의 편차를 보인다.

클러스터 성장 속도는 전선의 전파 속도 v와 직접 연결된다. 전선은 일정한 평균 속도로 전진하면서 KPZ(1+1) 차원의 스케일링 지수를 갖는 거칠기를 나타낸다. 즉, 전선 높이 w(L,t) ∼ t^β, β≈1/3, 그리고 공간 상관 길이 ξ ∼ t^{1/z}, z≈3/2 로서 KPZ universality class와 일치한다. 이는 양성 피드백이 단순한 확산을 넘어 비선형 성장 메커니즘을 제공함을 의미한다.

또한, 구조인자 S(q,t)를 이용해 임계 핵의 파수 q_c를 추정하고, q_c가 확산계수와 Y·Z 농도에 따라 변하는 것을 확인하였다. 전반적인 X 입자 비율 f_X(t)은 초기 급격한 상승 후 포화하는 S자형 곡선을 보이며, 이는 고정된 핵 수 n이 성장하면서 영역 A(t) ∝ πR(t)²가 선형적으로 확장된다는 가정 하에 df_X/dt = k n dA/dt (1‑f_X) 형태의 방정식으로 근사될 수 있다. 실제 데이터는 R(t) ∝ t 형태의 선형 성장과 어느 정도 일치하지만, 클러스터 간 병합, 비원형 전선, 그리고 Y 입자의 소멸(k₃) 등이 오차를 유발한다.

결론적으로, 양성 피드백과 느린 확산의 결합은 세포막 수준에서 신호분자의 국소적 클러스터링과 빠른 전파를 동시에 가능하게 하며, 이러한 현상은 KPZ와 같은 비평형 통계 물리학적 프레임워크로 기술될 수 있음을 보여준다. 이는 Ras‑SOS 시스템뿐 아니라, 막 근처에서 피드백이 작용하는 다양한 신호전달 네트워크에 보편적으로 적용될 수 있다.