일반화된 운하 규칙을 가진 불린 네트워크의 안정성

초록

본 논문은 유전 조절 모델링에 자주 사용되는 불린 네트워크에서, 입력 중 하나가 특정 ‘운하 상태’를 가질 때 그 입력에 의해 출력이 결정되는 운하 규칙을 일반화한 경우의 동적 안정성을 분석한다. 복잡하고 비랜덤한 네트워크 토폴로지를 고려한 기존의 무질서‑질서 전이 이론을 확장하여, 운하 입력의 비율과 토폴로지의 고유값이 시스템의 전이 임계점을 어떻게 이동시키는지를 수식적으로 도출하고, 수치 시뮬레이션으로 검증한다.

상세 분석

불린 네트워크는 각 노드가 0·1 두 값 중 하나를 취하고, 동기식 업데이트 규칙에 따라 이전 시점의 인접 노드 상태에 의존한다는 점에서 유전 회로의 이산적 모델링에 적합하다. 전통적인 K‑auffman 모델에서는 각 노드가 K개의 무작위 입력을 받아 무작위 Boolean 함수를 적용한다. 그러나 실제 유전 조절에서는 특정 전사인자가 ‘운하 상태’를 가질 때 해당 유전자의 발현이 강제로 결정되는 현상이 관찰되며, 이를 ‘운하 규칙(canalizing rule)’이라 부른다. 운하 규칙은 입력 중 하나가 특정 값을 가질 경우, 나머지 입력들의 상태와 무관하게 출력이 고정되는 특성을 가진다.

본 연구는 이러한 운하 규칙을 일반화한다. 즉, 각 노드 i에 대해 운하 입력의 집합 C_i와 그에 대응하는 운하 상태 s_{i}^{c} (0 또는 1)를 정의하고, 해당 입력이 운하 상태에 있으면 출력이 미리 정해진 값 f_i^{c} 로 고정된다. 운하 입력이 여러 개 존재하거나, 운하 상태가 다중값을 가질 경우도 허용한다. 이러한 일반화는 실제 생물학적 회로에서 복수의 전사인자가 동시에 결정적인 역할을 하는 상황을 포괄한다.

동역학적 안정성을 평가하기 위해 저자들은 평균장(average sensitivity) 개념을 확장하였다. 평균장은 무작위 초기 조건에서 한 단계 업데이트 시 한 노드의 출력이 입력 변동에 의해 바뀔 확률을 의미한다. 운하 규칙이 존재하면 평균장은 운하 입력이 운하 상태에 있을 확률 p_c와 비운하 입력이 변할 때의 민감도 ε에 의해 감소한다. 수식적으로는

λ_eff = λ_0 (1 - p_c) + ε p_c

여기서 λ_0는 전통적인 무작위 Boolean 함수의 평균장(보통 K·2^{-(K-1)})이며, λ_eff는 운하 규칙을 포함한 실제 네트워크의 효과적 평균장이다.

네트워크 토폴로지가 복잡하고 비랜덤일 경우, 평균장만으로는 전이 임계점을 정확히 예측하기 어렵다. 저자들은 인접 행렬 A의 가장 큰 특잇값(또는 스펙트럼 반경) ρ(A)를 도입하여, 전체 시스템의 선형화된 동역학이 x(t+1)=Λ·A·x(t) 형태임을 보였다. 여기서 Λ는 대각 행렬로 각 노드의 평균장을 담고 있다. 시스템이 안정하려면 스펙트럼 반경 ρ(Λ·A) < 1이어야 하며, 이는

ρ(A)·λ_eff < 1

이라는 간단한 조건으로 귀결된다. 따라서 운하 규칙이 증가하면 λ_eff가 감소하고, 동일한 토폴로지에서도 더 큰 ρ(A) 값을 허용하게 된다. 이는 운하 규칙이 네트워크의 ‘무질서‑질서’ 전이를 오른쪽(더 큰 연결도)으로 이동시킨다는 의미이다.

수치 실험에서는 스케일프리, 작은 세계, 그리고 실제 유전 네트워크 토폴로지를 사용해 ρ(A)와 λ_eff의 곱이 1을 초과하는 경우와 미만인 경우를 비교하였다. 운하 입력 비율 p_c를 0에서 0.5까지 변화시켰을 때, 전이 임계점이 명확히 오른쪽으로 이동함을 확인하였다. 특히, 운하 입력이 다중인 경우에도 동일한 선형 조건이 적용됨을 보여주어, 일반화된 운하 규칙이 복합적인 유전 회로에서도 동일한 수학적 프레임워크로 기술될 수 있음을 입증하였다.

이러한 결과는 두 가지 중요한 시사점을 가진다. 첫째, 실제 생물학적 네트워크가 보여주는 높은 안정성은 운하 규칙의 넓은 분포와 토폴로지의 비랜덤성(예: 모듈성, 계층성) 사이의 상호작용으로 설명될 수 있다. 둘째, 네트워크 설계나 합성 생물학에서 운하 입력을 인위적으로 도입하면, 시스템을 보다 견고하게 만들면서도 복잡한 동작을 유지할 수 있는 설계 원칙을 제공한다.