세포 간 신호분자 교환의 동기화 모델링

초록

본 논문은 두 박테리아 세포 사이의 신호분자 교환을 로지스틱 차분 방정식으로 모델링하고, 환경 파라미터(농도 c, 로지스틱 계수 r, 수용체 친화도 p)의 변동이 동기화와 혼돈에 미치는 영향을 Lyapunov 지수와 Cross‑Sample Entropy로 분석한다.

상세 분석

이 연구는 세포 간 신호전달을 “내·외부 환경 파라미터에 의해 변조되는 비선형 발진 시스템”으로 간주하고, 두 개의 상호 연결된 로지스틱 맵을 도입한다. 각 맵은 내부 신호분자 농도 x 또는 y를 나타내며, 업데이트 규칙은
x_{t+1}=r x_t(1‑x_t)+c p y_t, y_{t+1}=r y_t(1‑y_t)+c p x_t
와 같이 표현된다. 여기서 r은 전통적인 로지스틱 성장률로, 환경이 신호분자 증식에 제공하는 전반적 여건을 의미한다. c는 두 세포 사이의 물리적·화학적 결합 강도(예: 확산 거리, 매질 점도)이며, p는 수용체·운반체의 친화도(단백질 무질서와 내재 노이즈에 의해 변동)이다.

모델링 단계에서 저자들은 파라미터 공간을 체계적으로 탐색한다. 먼저 r=3.95(혼돈 영역 근접)에서 c를 0‒1 범위로 스캔하고, p를 0, 0.1, 0.2, 0.3, 0.4 등 다섯 단계로 변화시켰다. 각 파라미터 조합에 대해 Lyapunov 지수를 계산해 시스템이 발산(양의 지수)하거나 수렴(음의 지수)하는지를 판정하였다. 결과는 c≈0.4를 경계로 두 개의 구역이 형성됨을 보여준다. 낮은 농도(c<0.4)에서는 주로 비동기 상태가 유지되며, 간헐적인 동기화 윈도우가 존재한다. 반면 c>0.4에서는 대부분의 경우 Lyapunov 지수가 음수가 되어 완전 동기화가 확보된다. p가 증가할수록 동기화 구역이 확대되는 경향이 관찰되었으며, 이는 친화도가 높을수록 수용체-신호 결합이 강해져 상호 작용이 강화된다는 생물학적 직관과 일치한다.

동기화 정도를 정량화하기 위해 Cross‑Sample Entropy(Cross‑SampEn)를 도입하였다. Cross‑SampEn은 두 시계열 간 비동기성을 측정하는 지표로, 값이 작을수록 동기화가 강함을 의미한다. 저자들은 c가 0.4를 초과할 때 Cross‑SampEn이 급격히 감소함을 보고, 이는 Lyapunov 분석과 일관된 결과임을 확인했다. 또한 p가 0.6 이상이면 거의 모든 c 구간에서 낮은 Cross‑SampEn을 보였으며, 이는 친화도가 충분히 높을 경우 환경 변동에 대한 내성이 크게 향상된다는 것을 시사한다.

이 모델의 강점은 복잡한 생화학적 네트워크를 최소한의 변수( r, c, p )로 축소하면서도, 비선형 동역학과 확률적 교란을 동시에 고려한다는 점이다. 그러나 몇 가지 한계도 존재한다. 첫째, 실제 세포는 다수의 신호분자를 동시에 방출·수용하므로 2차원 맵으로의 축소가 모든 상호작용을 포착하지 못한다. 둘째, 파라미터 p를 고정된 확률값으로 다루었지만, 실제 수용체 친화도는 시간에 따라 동적으로 변할 수 있다(예: 포스트‑번역 변형). 셋째, 로지스틱 성장 모델은 포화 효과와 자가 억제 메커니즘을 단순화하므로, 고농도에서의 비선형 억제 현상을 과소평가할 가능성이 있다.

향후 연구에서는 (1) 다세포 네트워크로 확장하여 클러스터링 및 파동 전파 현상을 탐색하고, (2) stochastic differential equations를 도입해 연속적인 잡음 항을 명시적으로 모델링하며, (3) 실험적 데이터(예: 형광 라벨링을 통한 신호분자 농도 측정)와의 정량적 검증을 수행함으로써 모델의 생물학적 타당성을 강화할 수 있다.