패시비티 기반 상호연결 시스템의 안정성 판정과 생화학 반응망 적용

초록

본 논문은 생화학 반응망에서 영감을 얻은 비선형 상호연결 시스템 군에 대한 안정성 검사를 제시한다. 주요 결과 중 하나는 부분 시스템들의 패시비티 특성, 네트워크의 연결 구조, 그리고 연결 항의 부호 정보를 모두 포함하는 ‘소산성 행렬(dissipativity matrix)’의 대각 안정성으로부터 전체 네트워크의 전역 점근적 안정성을 판정한다는 것이다. 이 안정성 검사는 저자들의 이전 연구에서 제시된 순환 네트워크에 대한 ‘시컨트 기준(secant criterion)’을 포괄하며, 그래프 형태로 표현되는 일반적인 연결 구조로 확장한다. 두 번째 주요 결과는 상태 곱(term)들을 포함하도록 허용함으로써, 순환 시스템을 포함한 보다 넓은 모델 클래스에 적용 가능하도록 한다. 새로운 안정성 기준은 MAPK(미트겐 활성화 단백질 키네이스) 캐스케이드 모델과 대사망을 모티프로 한 분기형 연결 구조에 적용 사례를 제시한다. 마지막으로, 동일한 시스템으로 구성된 구획(compartment) 시스템에 확산 항이 존재할 때 안정성이 얼마나 강인한지를 다루는 추가 결과도 포함한다.

상세 요약

이 논문은 현대 시스템생물학에서 핵심적인 문제인 복잡한 생화학 네트워크의 동적 안정성을 수학적으로 검증하는 새로운 프레임워크를 제시한다. 기존의 ‘시컨트 기준’은 순환형 구조에 한정되어 있었으며, 각 반응 단계가 단일 입력‑단일 출력 형태로 모델링될 때만 적용 가능했다. 그러나 실제 세포 내 대사 경로나 신호 전달 경로는 다중 입력‑다중 출력, 분기·합류 구조, 그리고 상태 변수 간의 곱셈적 상호작용을 포함하는 복합적인 형태를 띤다. 저자들은 이러한 현실을 반영하기 위해 두 가지 혁신적인 요소를 도입한다. 첫째, 각 서브시스템을 ‘패시비티’라는 에너지 흐름 관점에서 분석하고, 그 패시비티 정도를 행렬 형태의 ‘소산성 행렬’에 정량화한다. 이 행렬은 (i) 서브시스템의 입력‑출력 패시비티 지표, (ii) 연결 그래프의 인접성, (iii) 연결 신호의 부호(양성·음성) 정보를 동시에 담고 있다. 행렬이 대각적으로 안정(diagonal stability)하면, 즉 적절한 양의 대각 행렬 D가 존재해 D·M + Mᵀ·D가 음정(definite negative)임을 보이면 전체 네트워크는 전역적으로 점근 안정성을 갖는다. 이는 Lyapunov‑함수 기반의 전통적 방법보다 계산적으로 효율적이며, 그래프 이론과 행렬 분석을 결합함으로써 구조적 직관을 제공한다. 둘째, 상태 곱(term) 즉, x_i·x_j 형태의 비선형 결합을 허용하도록 확장함으로써, 효소 촉매 반응식(Michaelis‑Menten)이나 복합 형성 반응 등 실제 생화학 모델에 흔히 등장하는 비선형성을 그대로 보존한다. 이 확장은 기존의 선형화 기반 안정성 검증이 놓치는 중요한 동적 특성을 포착한다. 논문은 두 가지 실제 사례를 통해 이론의 적용 가능성을 검증한다. 첫 번째는 MAPK 캐스케이드 모델로, 여러 단계의 인산화·탈인산화 반응이 순차적으로 진행되는 전형적인 신호 전달 경로이다. 여기서 저자들은 소산성 행렬을 구성하고, 대각 안정성을 확인함으로써, 파라미터 변화에도 불구하고 전체 캐스케이드가 안정적으로 동작함을 보였다. 두 번째는 대사망을 모티프로 한 분기형 네트워크로, 하나의 전구체가 여러 하위 경로로 분기되는 구조를 갖는다. 이 경우에도 동일한 행렬 기반 검증이 가능함을 보여, 그래프 구조가 복잡해도 일관된 안정성 판단이 가능함을 입증한다. 마지막으로, 동일한 서브시스템을 여러 구획에 복제하고 각 구획 사이에 확산(디퓨전) 항을 추가한 경우, 소산성 행렬에 확산 효과를 반영한 보정 항을 포함시켜도 대각 안정성이 유지되는 조건을 제시한다. 이는 조직 수준에서 세포 간 물질 교환이 존재할 때도 시스템이 안정적일 수 있음을 의미한다. 전체적으로 이 논문은 패시비티와 그래프 이론을 결합한 새로운 수학적 도구를 제공함으로써, 복잡한 생화학 네트워크의 설계·분석·제어에 있어 중요한 이정표가 된다.