[논문 리뷰] Analysis of Reaction Network Systems Using Tropical Geometry

Satya Swarup Samal, Dima Grigoriev외 2명|arXiv (Cornell University)|2015. 08. 24.

Polynomial and algebraic computation참고 문헌 18인용 수 4

한 줄 요약

이 논문은 다항식 또는 유리 함수 속도 법칙을 가진 반응 네트워크 시스템을 분석하기 위해 토포로지컬 기하학 기반의 방법을 제안한다. 이 방법은 차등 방정식에서 반대 부호의 지배적 단항식이 균형을 이루는 '토로피컬 평형'을 식별함으로써 작동한다. 이 접근법은 이러한 평형을 분지로 묶으며, 최소 분지는 가장 단순한 준안정 상태를 나타내고, 이들의 연결성을 계산하여 가능한 역학적 전이를 추론한다. 이를 통해 모델 단순화와 열화상도 및 그래프 구조를 통한 느린 변수 식별이 가능하다.

ABSTRACT

We discuss a novel analysis method for reaction network systems with polynomial or rational rate functions. This method is based on computing tropical equilibrations defined by the equality of at least two dominant monomials of opposite signs in the differential equations of each dynamic variable. In algebraic geometry, the tropical equilibration problem is tantamount to finding tropical prevarieties, that are finite intersections of tropical hypersurfaces. Tropical equilibrations with the same set of dominant monomials define a branch or equivalence class. Minimal branches are particularly interesting as they describe the simplest states of the reaction network. We provide a method to compute the number of minimal branches and to find representative tropical equilibrations for each branch.

연구 동기 및 목표

• 다항식 또는 유리 함수 속도 함수를 가진 복잡한 반응 네트워크 시스템을 분석하기 위한 새로운 방법을 개발하기 위해.
• 반대 부호의 지배적 단항식이 균형을 이루는 토포로지컬 평형을 통해 반응 네트워크에서 준안정 상태를 식별하기 위해.
• 토로피컬 평형을 공통된 지배적 단항식 집합을 공유하는 바탕으로 분지로 묶고, 최소 분지는 가장 단순한 역학적 상태를 나타내도록 하기 위해.
• 최소 분지 간의 연결성을 계산하여 준안정 상태 간의 가능한 전이를 추정하기 위해.
• 검증 및 모델 단순화 통찰을 위해 Biomodels 데이터베이스의 실제 생화학 모델에 이 방법을 적용하기 위해.

제안 방법

• 다항식 또는 유리 함수 속도 함수를 사용한 질량 작용 속도 법칙을 이용해 반응 네트워크를 수립하기 위해.
• 작은 매개변수 ε을 통한 스케일링을 도입하여 종 농도의 주요 크기를 할당함으로써 점근적 분석을 가능하게 하기 위해.
• 각 종의 미분 방정식에서 반대 부호의 지배적 단항식이 균형을 이루는 경우를 토포로지컬 평형으로 정의하기 위해.
• 해결 집합을 나타내기 위해 토포로지컬 초곡면의 유한한 교차를 구성함으로써 토포로지컬 전예약을 구축하기 위해.
• 공통된 지배적 단항식 집합을 공유하는 평형을 기반으로 평형을 분지로 묶으며, 최소 분지는 기초 성분으로 간주된다.
• 분지 간 포함 관계와 연결 관계를 각각 방향성 있는 그래프와 무방향 그래프로 표현하기 위해.

실험 결과

연구 질문

• RQ1 다항식 역학을 가진 반응 네트워크에 대해 토포로지컬 평형을 체계적으로 계산하는 방법은 무엇인가?
• RQ2 토로피컬 평형의 해 공간의 구조는 어떻게 되며, 이를 의미 있는 분지로 분해할 수 있는가?
• RQ3 최소 분지는 생화학적 시스템에서 불변 다양체와 느린 변수와 어떻게 관련이 있는가?
• RQ4 최소 분지의 연결성은 무엇이며, 이는 준안정 상태 간의 가능한 역학적 전이에 어떤 정보를 제공하는가?
• RQ5 이 방법은 Biomodels 데이터베이스의 실제 생화학 모델에 효과적으로 적용될 수 있는가?

주요 결과

• Biomodels 데이터베이스의 모델들에서 ε = 1/5일 때 최소 분지 수는 1에서 17까지 변동하였고, 타당한 해 시스템 총 수는 최대 9,847에 이르렀다.
• 최소 분지 수와 타당한 전체 시스템 수의 비율은 항상 낮았으며, 이는 최소 분지가 해의 희박하지만 구조적으로 중요한 부분집합을 대표하고 있음을 시사한다.
• 스케일링된 순서(μi − ai)의 열화상도는 종을 빠른 동역학과 느린 동역학으로 자연스럽게 군집화함을 보여주었으며, 이는 모델 단순화에서 느린 변수로 사용될 수 있음을 지지한다.
• 모든 테스트된 모델(BIOMD0000000035, 0000000072)에서 최소 분지 그래프는 단일 연결 성분을 가졌으며, 준안정 상태 간 전이에 대해 단일 지배적 경로가 존재함을 시사한다.
• 분지 간 포함 관계를 나타내는 방향성 그래프는 풍부한 계층적 구조를 드러내었으며, 다양한 평형 해의 복잡한 관계를 시사한다.
• 이 방법은 토포로지컬 평형을 통해 불변 다원체 후보를 성공적으로 식별하였으며, 최소 분지는 최소 불변 다원체에 해당할 가능성이 있으며, 이는 향후 연구를 위한 가설로 남겨진다.