양성 회로와 이산 동역학 시스템의 고정점 상한

초록

이 논문은 정수 구간들의 데카르트 곱 위에서 정의된 이산 동역학 시스템에 대해, 상호작용 그래프의 양성 회로 구조만을 이용해 고정점(또는 어트랙터)의 최대 개수를 제한하는 일반적인 상한을 제시한다. 기존의 토마스 추측(양성 회로가 없으면 고정점이 하나 이하)과 부울 네트워크에 대한 이전 결과들을 모두 포함·확장한다.

상세 분석

본 연구는 n개의 유한 정수 구간 X₁,…,Xₙ의 데카르트 곱 X=∏₁ⁿ Xᵢ 위에 정의된 임의의 지도 F:X→X에 대해, 그 동역학적 특성인 고정점 수를 그래프 이론적 관점에서 분석한다. 핵심 도구는 ‘이산 야코비안 행렬’과 이를 기반으로 정의된 지역 상호작용 그래프 G_F(x,v)이다. 여기서 (x,v)∈X′는 현재 상태 x와 한 단계 변화를 나타내는 방향 벡터 v∈{−1,1}ⁿ이며, x+v∈X이어야 한다. 행렬 원소 f_{ij}(x,v)=f_i(x+v_j e_j)−f_i(x)·v_j은 변수 x_j가 증가(또는 감소)할 때 f_i가 어떻게 변하는지를 이산적으로 측정한다. 이 값의 부호에 따라 G_F(x,v)에 양성·음성 간선이 부여된다. 논문은 기존 정의에 ‘임계값( threshold )’ 조건을 추가해 G_F(x,v)′을 정의함으로써, 그래프가 실제 비동기 전이 그래프 Γ(F)와 완전히 일치하도록 보강한다. 이렇게 함으로써 모든 지역 그래프가 전역 그래프 G(F)=⋃_{(x,v)}E(G_F(x,v))의 부분 그래프가 된다.

양성 회로는 회로 내 음성 간선의 개수가 짝수인 경로로 정의되며, 회로에 포함된 모든 정점 집합 I가 ‘양성 피드백 정점 집합’이면, 각 양성 회로가 최소 하나의 정점을 I와 교차한다는 의미다. 주요 정리(정리 3, 정리 5)는 다음과 같다. 임의의 I⊆{1,…,n}이 모든 지역 그래프 G_F(x,v)의 양성 회로에 대해 피드백 정점 집합이라면, 고정점(또는 어트랙터)의 개수는 ∏{i∈I}|X_i| (정리 3) 혹은 ∑{i∈I}|T_i(G_F)|+1 (정리 5) 이하로 제한된다. 여기서 T_i(G_F) 는 i번째 변수와 연관된 임계값 t=x_i+v_i/2 중, 해당 회로에 실제로 기여하는 값들의 집합이다. 특히 X가 {0,1}ⁿ인 부울 경우 |X_i|=2이므로 기존의 2^{|I|} 상한을 재현한다.

증명은 I의 크기에 대한 귀납법으로 전개된다. 기본 단계에서는 I=∅이면 모든 지역 그래프가 양성 회로를 갖지 않으므로 토마스 추측에 따라 고정점이 하나 이하임을 보인다. 귀납 단계에서는 I에 포함된 최소 인덱스(예: 1)를 선택하고, X₁을 T₁(G_F) 로 정의된 임계값에 따라 구간 P 로 분할한다. 각 구간마다 F를 제한·조정한 새로운 지도 ˜F를 구성하고, ˜F의 지역 그래프가 원래 그래프의 부분 그래프임을 보임으로써 귀납 가정을 적용한다. 이 과정에서 T_i(˜F)⊆T_i(G_F)임을 이용해 상한이 유지됨을 확인한다. 최종적으로 모든 가능한 구간에 대해 어트랙터 수가 위의 식을 초과할 수 없음을 증명한다.

이 결과는 두 가지 중요한 의미를 가진다. 첫째, 양성 회로의 존재 여부만으로도 고정점 수에 대한 강력한 제한을 제공하므로, 실제 생물학적 네트워크에서 상호작용 부호만 알면 다중 안정성 가능성을 정량적으로 평가할 수 있다. 둘째, 부울 네트워크를 포함한 모든 이산 모델에 적용 가능하므로, 기존의 특수 경우(예: 토마스 추측, Aracena‑Goles‑Demongeot 결과)를 일관되게 일반화한다. 또한 비동기 전이 그래프 Γ(F)를 사용함으로써, 연속적인 동시 업데이트 모델이 아닌 실제 유전자 발현의 비동기적 특성을 더 잘 반영한다는 점에서도 의의가 크다.