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

본 논문은 유전자의 발현 수준을 정수 구간으로 모델링하고, 그 동역학을 비동기적 업데이트 규칙에 따라 정의된 지도 F:X→X 로 기술한다. 서론에서는 유전자 네트워크 분석에서 상호작용 그래프가 핵심적인 역할을 함을 강조하고, 기존 연구(Thomas의 추측, Richard‑Comet 정리, Aracena‑Goles‑Demongeot 정리)의 한계를 지적한다. 특히, 기존 결과들은 부울(0/1) 네트워크에 국한되거나, 전역 그래프에 양성·음성 간선이 동시에 존재할 때 적용이 어려운 점이 있었다. 2절에서는 비동기 상태 전이 그래프 Γ(F)를 정의한다. 각 정점은 시스템의 전체 상태 x∈X이며, 한 정점에서 다른 정점으로의 간선은 어떤 한 유전자의 농도가 목표값 f_i(x) 쪽으로 한 단계(±1) 이동할 수 있음을 의미한다. 고정점은 Γ(F)에서 외부로 나가는 간선이 없는 정점이며, 어트랙터는 더 큰 의미에서 탈출이 불가능한 최소 트랩 도메인으로 정의된다. 이는 생물학적 의미에서 안정된 세포 상태와 다중 안정성을 포괄한다. 3절에서는 이산 야코비안 행렬 F′(x,v)와 이를 기반으로 한 지역 상호작용 그래프 G_F(x,v) 를 도입한다. f_{ij}(x,v)는 변수 x_j가 한 단계 변할 때 f_i가 어떻게 변하는지를 측정하며, 그 부호에 따라 양성·음성 간선이 부여된다. 기존 정의에 ‘임계값’ 조건을 추가해 G_F(x,v)′을 정의함으로써, 그래프가 실제 전이 가능성(Γ(F)와의 일치)을 반영하도록 보강한다. 전역 상호작용 그래프 G(F)는 모든 (x,v) 에 대한 지역 그래프의 간선 합집합이며, 양성·음성 간선이 동시에 존재할 수 있다. 4절은 본 논문의 핵심 정리들을 제시한다. 정리 4 (Richard‑Comet)에서는 모든 지역 그래프가 양성 회로를 갖지 않을 경우 Γ(F) 가 유일한 어트랙터를 가진다는 기존 결과를 재언급한다. 정리 5 (주요 결과)에서는 임의의 피드백 정점 집합 I (모든 지역 그래프의 양성 회로가 I와 교차) 에 대해 어트랙터 수가 ∑_{i∈I} |T_i(G_F)| + 1 이하임을 보인다. 여기서 T_i(G_F)는 i번째 변수와 연관된 임계값 t=x_i+v_i/2 중, 실제 양성 회로에 기여하는 값들의 집합이다. 정리 3은 정리 5의 특수 경우로, |T_i(G_F)| 를 |X_i| 로 대체해 고정점 수가 ∏_{i∈I}|X_i| 이하임을 제시한다. 부울 경우 X_i={0,1} 이므로 ∏|X_i|=2^{|I|} 가 되며, 이는 Aracena‑Goles‑Demongeot 결과와 일치한다. 증명은 I의 원소 개수에 대한 귀납법으로 전개된다. 기본 단계(I=∅)에서는 모든 지역 그래프가 양성 회로를 갖지 않으므로 토마스 추측에 따라 고정점이 하나 이하임을 보인다. 귀납 단계에서는 I에 포함된 최소 인덱스(예: 1)를 선택하고, X₁을 T₁(G_F) 로 정의된 임계값에 따라 구간 P 로 분할한다. 각 구간 Y∈P에 대해 F를 수정한 ˜F 를 정의하고, ˜F의 지역 그래프가 원래 그래프의 부분 그래프임을 보인다. 이때 T_i(˜F)⊆T_i(G_F) 임을 이용해 귀납 가정에 의해 어트랙터 수가 ∑_{i∈I\{1\}}|T_i|+1 이하임을 얻고, 구간의 개수 |P|=|T₁|+1 를 곱해 최종 상한을 도출한다. 증명 과정에서 비동기 전이와 임계값 조건이 핵심적인 역할을 하며, 전역 그래프 G(F)와 지역 그래프 G_F(x,v) 사이의 포함 관계를 정교히 이용한다. 마지막으로 5절에서는 결과의 생물학적 의미와 Thomas의 논리적 방법과의 연관성을 논의한다. 양성 회로는 양성 피드백 루프를 의미하며, 이는 다중 안정성(세포 분화 등)의 근본 메커니즘으로 해석될 수 있다. 본 정리는 회로 구조만 알면 고정점 수를 정량적으로 제한할 수 있음을 보여, 실험적으로 상호작용 부호만 확보된 경우에도 네트워크의 동적 가능성을 평가할 수 있는 강력한 도구를 제공한다. 또한, 비동기 업데이트 모델을 사용함으로써 실제 유전자 발현의 비동시성을 반영하고, 기존 동시 업데이트 모델보다 현실적인 분석이 가능함을 강조한다.