무작위에서 방향성으로 대장균 화학주성의 간소화 모델

초록

본 논문은 대장균(E. coli)의 화학주성(chemotaxis)을 설명하기 위해 기존 5단계 메틸화 수용체 모델을 3단계로 축소한 간소화 모델을 제시한다. 축소된 네트워크에 대해 튜블링(tumbling) 속도의 응답 함수를 해석적으로 도출하고, 약한 공간적 유인물 농도 구배 하에서의 평균 이동 속도(드리프트 속도)를 계산한다. 결과는 실험적 관찰과 일치하며, 모델이 적응(adaptation)과 감도(sensitivity)를 동시에 구현함을 보여준다.

상세 분석

대장균은 화학주성이라는 메커니즘을 통해 영양소가 풍부한 방향으로 이동하거나 독성 물질을 회피한다. 이 과정은 ‘런(run)’이라 불리는 직선 구간과 ‘텀블(tumble)’이라 불리는 방향 전환 단계가 번갈아 나타나는 ‘지그재그’ 움직임으로 구현된다. 런의 지속 시간은 주변 화학물질 농도의 공간 구배에 따라 조절되며, 이는 수용체-신호전달 네트워크에 의해 매개된다. 기존 모델은 수용체의 메틸화 상태를 5단계(04)로 가정하고 복잡한 비선형 미분 방정식으로 기술했지만, 실험적으로는 메틸화 수준이 34 단계만을 주로 이용한다는 점이 알려졌다.

본 연구는 이러한 실험적 사실을 반영해 메틸화 상태를 3단계(0,1,2)로 단순화한다. 각 상태는 활성화 확률 A_i와 비활성화 확률 B_i를 갖으며, 메틸 전이율 k_m와 탈메틸 전이율 k_d가 각각 인접 상태 사이에서 작용한다. 저농도 구배 가정 하에 시스템을 선형화하여 작은 신호에 대한 응답을 분석한다. 핵심은 ‘튜블링 속도 α(t)’가 수용체 활성도 변동에 비례한다는 가정이다. 이를 바탕으로 라플라스 변환을 적용해 α(s)=R(s)·c(s) 형태의 응답 함수 R(s)를 도출한다. 여기서 c(s)는 외부 화학물질 농도 구배의 라플라스 변환이며, R(s)는 네트워크 파라미터(k_m, k_d, 결합 상수 등)에 전적으로 의존한다.

응답 함수는 시간 영역에서 지수적 감쇠와 진동 성분을 포함하는 복합 형태이며, 특히 ‘근완전 적응(near‑perfect adaptation)’을 보장하기 위해 R(0)=0이 성립한다. 이는 장기적으로 일정한 농도에 대해 튜블링 속도가 변하지 않음을 의미한다. 논문은 이 응답 함수를 이용해 약한 구배 G=∇c가 존재할 때 평균 드리프트 속도 v_d를 계산한다. 선형 응답 이론에 따라 v_d = χ·G 로 표현되며, χ는 ‘감도 계수(sensitivity coefficient)’로서 R(s)의 첫 번째 순간(moment)과 런-텀블 주기의 평균값 τ_run, τ_tumble에 의해 결정된다. 구체적으로 χ = (v_0·τ_run^2)/(1+ (k_m+k_d)τ_run)·R′(0) 형태를 취한다. 여기서 v_0는 런 중 속도, R′(0)은 응답 함수의 기울기이다.

수치 시뮬레이션 결과는 이론적 χ값과 실험적으로 측정된 드리프트 속도가 10 % 이내 차이로 일치함을 보여준다. 또한, 메틸화 단계가 3개인 경우에도 ‘적응 시간(≈10 s)’과 ‘감도(≈0.02 µm·s⁻¹·µM⁻¹)’가 기존 5단계 모델과 거의 동일함을 확인한다. 이는 모델 축소가 핵심 동역학을 보존하면서 계산 복잡성을 크게 낮춘다는 중요한 의미를 가진다.

이 연구의 한계는 (1) 약한 구배 가정에만 적용 가능하다는 점, (2) 세포 간 변이와 외부 잡음(noise)을 무시했다는 점, (3) 3단계 메틸화가 실제 생리학적 상황에서 언제든지 충분한지에 대한 실험적 검증이 부족하다는 점이다. 향후 연구에서는 비선형 구배, 다중 화학물질 혼합, 그리고 stochastic simulation을 도입해 모델의 일반성을 검증할 필요가 있다.