다중계층 부울 네트워크로 보는 세포 상태 역학

초록

본 논문은 에피제네틱, 유전자, 전사인자 등 서로 다른 세포 구성요소를 명시적으로 구분하고, 이들 간 상호작용의 유형을 반영한 고차원 부울 네트워크(Higher‑Order Boolean Network, HOBN)를 제안한다. 규제 상호작용과 그 강도를 행렬로 표현하고, 해당 행렬의 고유값 문제를 풀어 시스템의 안정성을 분석한다. 안정·혼돈 두 가지 전통적 구분 외에, 세포 상태가 구성요소들의 독립적 연구만으로 설명 가능한 ‘단순’ 상태와, 복합적인 상호작용을 고려해야 하는 ‘복합’ 상태를 구분한다. 마지막으로, 더 높은 수준과 고차 동역학을 포함하도록 모델을 확장하는 방법을 제시한다.

상세 요약

이 연구는 기존 부울 네트워크 모델이 세포 내 다양한 조절 층을 동일한 노드 집합에 혼합함으로써, 각 층의 상대적 기여도를 정량화하기 어려웠던 점을 지적한다. 이를 해결하기 위해 저자들은 세 가지 기본 구성요소—에피제네틱 마크, 유전자 발현, 전사인자 활동—를 각각 별도의 서브네트워크로 정의하고, 이들 사이의 상호작용을 ‘제어’와 ‘피드백’ 두 종류로 구분한다. 수학적으로는 각 서브네트워크의 상태를 0‑1 벡터로 나타내고, 전체 시스템을 3×3 블록 행렬 M 으로 표현한다. 행렬 원소 M_{ij} 는 i‑층이 j‑층에 미치는 평균적인 논리적 영향력과 연결 강도를 합산한 값이며, 이는 랜덤 부울 함수의 평균 감도와 연결 확률을 곱한 형태로 계산된다. 시스템의 동역학은 x(t+1)=F(M·x(t)) 와 같이 비선형 함수 F 에 의해 진행되지만, 평균적인 동작을 선형화하면 고유값 λ 가 주요 안정성 지표가 된다. |λ|<1이면 모든 작은 교란이 소멸해 안정(ordered) 상태에 머무르고, |λ|>1이면 교란이 증폭돼 혼돈(chaotic) 영역에 진입한다. 흥미로운 점은, λ이 1에 근접한 경우에도 행렬의 구조에 따라 ‘단순’ 상태와 ‘복합’ 상태가 구분된다는 것이다. 즉, λ<1이라도 M 의 비대각원소가 크게 작용하면 서로 다른 층 간의 상호 의존성이 강해져, 개별 층을 독립적으로 분석해도 전체 세포 상태를 예측할 수 없는 복합적 거동이 나타난다. 반대로, 비대각원소가 약하면 각 층을 별도로 연구해도 전체 동역학을 재구성할 수 있다. 이러한 구분은 기존의 ‘안정‑혼돈’ 이분법을 넘어, 세포 운명 결정 메커니즘을 더 세밀히 파악할 수 있는 새로운 틀을 제공한다. 또한, 저자들은 고차 상호작용(예: 삼중 상호작용)이나 다중 레벨(예: 조직 수준)까지 확장 가능한 일반화된 행렬 M 구조를 제시함으로써, 복잡한 생물학적 네트워크를 계층적·다중 차원적으로 모델링할 수 있음을 보여준다.