단백질 상호작용 네트워크 위상 모델링 두 단계 접근법

초록

본 논문은 단백질‑단백질 상호작용 네트워크의 위상을 재현하기 위해 ‘선호적 소멸’과 ‘유사성 기반 연결’ 두 메커니즘을 결합한 두 단계 모델을 제안한다. 하나의 자유 파라미터 α만을 조정함으로써 실제 생물 데이터와 동일한 정점 수, 엣지 수, 차수 분포, 클러스터 크기 분포, 클러스터링 계수 및 평균 최단 경로 길이를 얻는다. 결과적으로 단백질 네트워크는 스케일프리보다는 작은 세계(small‑world) 특성을 보이며, 모델은 다양한 미생물 및 인간 데이터에 잘 맞는다.

상세 분석

이 연구는 기존의 무작위 혹은 선호적 부착 모델이 단백질‑단백질 상호작용 네트워크의 복합적인 위상 특성을 설명하지 못한다는 점에 착안한다. 저자들은 초기 완전 그래프(모든 정점이 서로 연결된 상태)에서 시작해 두 가지 확률적 규칙을 순차적으로 적용한다. 첫 번째 단계는 ‘선호적 소멸(preferential depletion)’이다. 각 반복에서 무작위 정점 i를 선택하고, 그 정점의 인접 엣지 eij 중 하나를 선택해 제거한다. 제거 확률 p_{i,j}=p_j/N_i는 인접 정점 j의 차수 k_j에 의존하며, p_j=(k_j-α)/k_j (k_j>1) 로 정의된다. α>0이면 고차수 정점의 엣지가 보존될 확률이 높아져, 저차수 정점이 먼저 소멸되는 ‘가난한 자는 가난하게 남는다’ 현상이 발생한다. α=0이면 완전 무작위 소멸이 된다. 두 번째 단계는 ‘유사성 기반 연결(similarity)’이다. 정점 i와 j를 무작위로 선택하고, 공통 이웃 수 N_c(i,j)를 이용해 연결 확률 p_{i,j}=