p53 부정 피드백 루프의 부울 네트워크 동기화와 확률적 순환

초록

이 논문은 p53‑Mdm2 등 핵심 부정 피드백 회로를 부울 네트워크로 모델링하고, 결정론적·확률적 두 버전을 비교한다. 외부 손상 신호가 주어지면 시스템은 저농도 정상 상태에서 벗어나 제한된 상태 순환(리미트 사이클)으로 수렴하며, 이 순환은 잡음 수준이 크게 변해도 안정성을 유지한다. 마코프 연쇄의 순환 이론을 적용해 확률적 모델의 주요 순환을 정량화하고, 전력 스펙트럼의 피크 주기가 해당 순환의 상태 수와 일치함을 확인한다. 주요 순환과 그 외 순환 사이의 순환 강도 차이가 ‘확률적 동기화’를 야기한다는 점을 강조한다.

상세 분석

본 연구는 p53 경로의 핵심 부정 피드백 루프, 즉 p53‑Mdm2, p53‑Wip1, ATM‑Chk2 등 세 개의 작은 회로를 선택하고, 각각을 0·1 이진 상태로 표현한 부울 네트워크를 구축하였다. 결정론적 모델에서는 모든 노드를 동기식으로 업데이트하며, 각 노드의 논리 함수를 실제 생화학적 상호작용(전사·번역·분해 등)의 정성적 특성을 반영하도록 설계하였다. 예컨대 p53은 DNA 손상 신호가 존재하면 1로 전환되고, Mdm2은 p53이 활성화될 때만 1이 되는 식으로 정의된다. 이러한 규칙은 Boolean 함수의 진리표를 통해 명시적으로 제시되며, 전체 네트워크는 2ⁿ(여기서 n은 노드 수)개의 가능한 상태 공간을 갖는다.

확률적 모델은 결정론적 업데이트에 잡음을 도입한 형태로, 각 노드가 자신의 논리 함수에 따라 전이할 확률을 β(역온도) 파라미터를 통해 조절한다. β가 무한대에 가까울수록 결정론적 전이와 동일해지고, β가 0에 가까워지면 완전 무작위 전이가 된다. 이때 전체 시스템은 유한 마코프 연쇄로 기술되며, 전이 행렬 P는 각 상태 i에서 j로 전이될 확률 P_{ij}로 구성된다.

연구진은 마코프 연쇄의 순환 이론(circulation theory)을 활용해 상태 전이 그래프를 순환(사이클)들의 집합으로 분해하였다. 순환은 각 사이클 C에 대해 순환 강도 γ(C) = Π_{(i→j)∈C} P_{ij} - Π_{(j→i)∈C} P_{ji} 로 정의되며, 이는 해당 사이클이 순방향으로 흐르는 정도를 정량화한다. 수치 실험 결과, β가 중간값(예: 5~10)에서도 γ가 가장 큰 하나의 주요 사이클이 존재했으며, 이 사이클은 결정론적 모델에서 관찰되는 리미트 사이클과 정확히 동일한 상태 순서를 가진다. 즉, 확률적 모델의 ‘지배 순환’이 결정론적 리미트 사이클의 자연적인 확장임을 확인하였다.

또한, 시간 시계열 데이터를 푸리에 변환하여 전력 스펙트럼을 구했을 때, 주요 피크 주기가 주요 순환의 상태 수와 일치함을 발견했다. 이는 순환 강도가 압도적으로 큰 사이클이 시스템의 동기화된 진동을 주도한다는 물리적 직관과 일치한다. 잡음 수준을 크게 증가시켜도(β→0) 주요 순환의 γ가 완전히 사라지지는 않으며, 대신 주변 순환들의 γ가 상대적으로 증가해 동기화 정도가 감소한다는 점을 통해 시스템의 강인성을 정량적으로 평가하였다.

생물학적 해석 측면에서, 부정 피드백이 p53을 저농도 정상 상태로 유지하면서도 손상 신호가 들어오면 빠르게 고농도 진동 상태(리미트 사이클)로 전환한다는 기존 실험 결과와 모델이 일치한다. 특히, 잡음(세포 내 변동성)에도 불구하고 이러한 전이와 진동이 유지된다는 점은 세포가 외부 스트레스를 견디면서도 정확한 결정 메커니즘을 보유하고 있음을 시사한다.

요약하면, 본 논문은 부울 네트워크와 마코프 순환 이론을 결합해 p53 부정 피드백 루프의 결정론적·확률적 동역학을 통합적으로 분석했으며, 주요 순환과 전력 스펙트럼 사이의 정량적 연결고리를 제시함으로써 ‘확률적 동기화’ 현상의 메커니즘을 명확히 밝혔다.