질량작용 반응망에서 스위치 행동을 배제하는 야코비안 기준

초록

본 논문은 질량작용 속도를 갖는 임의의 화학 반응망에 대해, 야코비안 행렬의 행렬식이 영이 되지 않는 경우 네트워크가 주입성(injective) 하며, 따라서 어떤 속도 상수 선택에서도 양의 정상상태가 다중으로 존재할 수 없음을 보이는 새로운 판정법을 제시한다. 또한 완전 개방(open) 네트워크와의 관계, 퇴화 정상상태의 존재 여부, 그리고 P‑행렬과의 연계도 논의한다.

상세 분석

논문은 먼저 화학 반응망을 종(species), 복합(complex), 반응(reaction)이라는 세 집합으로 정의하고, 질량작용 속도법칙에 따라 종 형성 속도 함수 fₖ(c)=∑₍y→y’∈R₎k_{y→y’}c^{y}(y’−y)를 도입한다. 이때 스토이키오메트리 공간 Γ=span{y’−y | y→y’∈R}의 차원 s와 그 직교 보완 Γ^⊥(차원 d=n−s)를 이용해 “총량 보존식” ω·c=const 를 얻는다. 정상상태는 fₖ(c)=0이지만, 스토이키오메트리 클래스 안에서만 비교하면 ω·c가 동일한 두 점 a,b에 대해 fₖ(a)=fₖ(b)=0이면 다중 정상상태가 존재한다는 의미가 된다.

주입성은 “fₖ가 양의 직교 보존식과 결합된 확장 함수 ϕₖ(c)=(ω₁·c,…,ω_d·c,fₖ₁(c),…,fₖ_n(c))에 대해 내부에서 일대일(injective)이다” 로 정의한다. 핵심 정리는 ϕₖ의 야코비안 J(ϕₖ) 행렬식이 모든 양의 c와 모든 양의 속도 상수 k에 대해 영이 아니면 네트워크는 주입성이고, 따라서 어떤 스토이키오메트리 클래스에서도 양의 정상상태가 하나만 존재한다는 것이다. 행렬식은 c와 k에 대한 다항식이며, k에 대해서는 1차, c에 대해서는 다항식 차수가 제한된다. 논문은 이 행렬식의 각 항을 네트워크 구조(반응 스토이키오메트리와 보존식)만으로 완전히 계산할 수 있음을 보여준다.

또한 완전 개방 네트워크(모든 종에 대해 S_i→0 반응이 존재)와 일반 네트워크 사이의 관계를 분석한다. 완전 개방 네트워크의 야코비안 행렬식은 기존 연구(Cracun‑Feinberg)와 동일하게 정의되지만, 일반 네트워크에서는 일부 종에 대한 보존식이 존재하므로 원래 fₖ를 그대로 사용할 수 없고, 보존식으로 대체한 수정 함수 ϕₖ를 사용한다. 저자는 완전 개방 네트워크가 주입성일 경우, 원래 네트워크는 (1) 주입성이면서 모든 정상상태가 비퇴화(non‑degenerate)하거나, (2) 모든 정상상태가 퇴화(degenerate)한다는 두 가지 경우만 남는다는 정리를 증명한다. 퇴화 경우는 보존식과 반응 스토이키오메트리의 특정 선형 종속 구조에 의해 완전히 특징지어진다.

마지막으로 P‑행렬 이론을 연결한다. 기존 연구에서는 −J(fₖ)가 P‑행렬이면 주입성을 보장한다는 결과가 있었는데, 저자는 J(ϕₖ)의 일부 행을 부호 변환하면 P‑행렬이 된다는 것을 보여준다. 따라서 P‑행렬 판정도 본 기준의 충분조건이 된다. 전체적으로 이 논문은 네트워크 구조만으로 다중 양성 정상상태와 퇴화 정상상태를 배제할 수 있는 실용적인 대수적 도구를 제공한다는 점에서 이론 화학·생물학 및 시스템생물학 분야에 큰 의미를 가진다.