연속시간 확률 이진 모델링으로 보는 세포 동역학

초록

본 논문은 세포 내 분자 상태를 이진값(활성/비활성)으로 단순화하고, 각 상호작용을 연속시간 마코프 과정으로 기술함으로써, 전통적인 화학반응식 대신 유한 개의 선형 미분방정식으로 전체 집단의 동역학을 예측하는 새로운 정량 모델링 프레임워크를 제시한다.

상세 분석

이 연구는 기존의 화학반응식 모델링이 요구하는 정확한 반응계수, 화학양론, 절대 농도 등의 상세 정보가 부족한 상황에서, 세포 집단의 이질성과 내재된 확률성을 직접적으로 반영할 수 있는 대안을 제시한다. 핵심 아이디어는 각 분자를 ‘활성(T)’ 또는 ‘비활성(F)’이라는 두 가지 상태만으로 표현하고, 이들 사이의 전이율을 연속시간 확률 과정(마코프 프로세스)으로 모델링한다는 점이다. 이렇게 하면 전체 시스템은 2^N개의 가능한 상태(여기서 N은 노드 수)로 정의되지만, 전이율 행렬이 선형이므로 마스터 방정식 형태의 유한 차원 선형 미분방정식으로 정확히 기술할 수 있다.

논문은 먼저 두 노드(A→B)와 상호 억제(A↔B) 시스템을 통해 전이율 행렬을 구성하고, 초기 조건에 따라 해석적 해를 도출한다. 여기서 얻은 결과는 전통적인 ODE 기반 모델이 예측하는 단일 안정상태와 달리, 집단 수준에서는 두 안정상태가 동시에 존재하는 이질적 분포를 자연스럽게 설명한다. 또한, 부정 피드백 회로를 포함한 2노드 진동 모델을 분석하면서, 개별 세포는 진동하지만 전이율이 지수분포를 따르기 때문에 집단 평균에서는 급격히 감쇠되는 ‘데코히런스’ 현상이 발생함을 보여준다. 이는 마코프 가정(지수 대기시간)이 실제 세포 내 타이밍 변동성을 어떻게 평균화하는지를 직관적으로 설명한다.

다음으로 Boolean 연산(AND, OR 등)을 확장하여 다중 노드 네트워크에도 적용 가능함을 제시한다. 특히, TNF‑NF‑κB 신호전달 경로를 6개의 노드와 13개의 전이율 파라미터만으로 재구성한 사례는, 기존 33종 분자와 110개의 반응계수를 필요로 했던 전통 모델에 비해 훨씬 간결하면서도 WT, IκB KO, A20 KO와 같은 유전자 결손 실험 결과를 정량적으로 재현한다는 점에서 실용성을 강조한다.

이 프레임워크의 장점은 (1) 선형 방정식으로 해석이 용이해 수치적 불안정성이나 발산 문제를 회피한다, (2) 파라미터가 ‘시간 스케일’과 ‘활성/비활성 시그널 강도’ 두 종류만으로 제한돼 실험 데이터와의 매핑이 직관적이다, (3) 세포 간 이질성을 초기 조건에 그대로 반영함으로써 집단 평균이 아닌 확률 분포 자체를 예측한다는 점이다. 반면, 마코프 가정에 의한 지수 대기시간은 실제 생물학적 과정에서 관측되는 비지수적 대기시간(예: 멀티스텝 촉진)과 차이가 있을 수 있으며, 이진화된 상태표현이 복합적인 농도 의존성을 충분히 포착하지 못한다는 한계도 존재한다. 향후에는 비마코프 전이, 다중값(다중 상태) 확장, 그리고 세포 간 상호작용을 전이율 행렬에 확률적으로 결합하는 방법을 모색함으로써 현재 모델의 적용 범위를 넓힐 수 있을 것이다.