대사망 네트워크 충격도 분포의 브랜칭 프로세스 근사 분석

초록

본 논문은 대사망의 구조적 강인성을 평가하기 위한 지표인 충격도(impact degree)의 확률분포를, 브랜칭 프로세스 이론을 이용해 예측하는 방법을 제안한다. 반응‑물질 이분 그래프를 반응 전용 네트워크로 투영하고, 각 반응의 ‘자식 수’를 정의한 뒤 포아송·멱법칙·실제 데이터 기반 세 가지 오프스프링 분포 모델을 적용한다. KEGG에서 추출한 4종(대장균, 바실러스, 효모, 인간) 대사망에 대해 실험적 충격도와 이론적 예측을 비교했으며, 전반적으로 좋은 일치도를 보였지만 사이클 존재와 트리 가정 위반 등 한계도 논의한다.

상세 분석

본 연구는 대사망의 구조적 강인성을 정량화하기 위해 ‘충격도’를 핵심 지표로 채택한다. 충격도는 하나의 반응을 제거했을 때 연쇄적으로 비활성화되는 반응의 총수를 의미하며, 이는 시스템 전반에 퍼지는 ‘연쇄 실패(cascading failure)’와 동일시할 수 있다. 저자들은 먼저 대사망을 이분 그래프(물질 노드 Vc와 반응 노드 Vr)로 정의하고, 각 반응이 활성(1) 혹은 비활성(0) 상태를 갖는 부울 모델을 구축한다. 반응 제거 후 물질·반응 상태를 논리 규칙에 따라 반복 업데이트하면, 최종적으로 비활성화된 반응들의 수가 충격도가 된다.

핵심 기법은 이 과정을 ‘브랜칭 프로세스(branching process)’로 근사하는 것이다. 대사망을 반응‑전용 네트워크로 투영하여, 반응 A의 산물이 반응 B의 기질이 될 경우 A→B 방향의 간선을 만든다. 여기서 ‘자식 수(d_i)’는 해당 반응이 가진 외향 차수(k_out)와 입향 차수(k_in)의 관계에 따라 정의된다. 구체적으로, 입향 차수가 1인 경우에만 외향 차수를 자식 수로 채택하고, 그 외에는 0으로 설정한다. 이는 트리 구조를 가정했을 때, 유일한 전구체 경로를 통해서만 영향이 전파된다는 물리적 직관에 기반한다.

다음으로 저자들은 세 가지 오프스프링 분포 모델을 제시한다.

1. 포아송 모델: 평균 μ를 갖는 포아송 분포를 가정하고, 브랜칭 프로세스의 전체 자식 수(즉, 충격도) r은 Borel 분포를 따른다. μ=1인 임계 상황에서는 P(r)∝r^{-3/2} 형태의 멱법칙이 도출된다.
2. 멱법칙 모델: 실제 대사망이 멱법칙 차수 분포를 보이는 점을 이용해, 오프스프링 d가 P(d)∝d^{-(γ+1)} 형태를 갖는 경우를 분석한다. Saichev 등(2007)의 결과를 인용해, 1<γ<2 구간에서는 P(r)∝r^{-(1+1/γ)}가 얻어지며, γ>2에서는 포아송 경우와 유사한 형태가 된다.
3. 실제 데이터 기반 모델: 경험적으로 측정된 오프스프링 분포를 그대로 사용한다. 확률생성함수 f(s)를 정의하고, 재귀식 F(s)=f(sF(s))를 통해 전체 충격도 분포 F(s)를 구한다. 라그랑주 전개와 미분을 이용해 P(r) = (1/r!) d^{r-1}/ds^{r-1}