[논문 리뷰] Generalized Mass-Action Systems and Positive Solutions of Polynomial Equations with Real and Symbolic Exponents

Stefan C. Müller, Georg Regensburger|arXiv (Cornell University)|2014. 06. 25.

Protein Structure and Dynamics참고 문헌 48인용 수 24

한 줄 요약

이 논문은 실수 및 기호적 반응 속도 차수를 가진 일반화된 마스터 반응계에 화학 반응망 이론을 확장하여, 스토이키오메트릭 및 반응 속도 차수 부분공간에 대한 부호 벡터 조건을 사용해 양의 평형 상태의 존재성과 유일성을 증명한다. 이는 버치의 정리의 일반화이며, 오리엔티드 매트로이드 이론과 이항방정식을 통해 다평형성에 대한 알고리즘 기반 기준을 제공한다.

ABSTRACT

Dynamical systems arising from chemical reaction networks with mass action kinetics are the subject of chemical reaction network theory (CRNT). In particular, this theory provides statements about uniqueness, existence, and stability of positive steady states for all rate constants and initial conditions. In terms of the corresponding polynomial equations, the results guarantee uniqueness and existence of positive solutions for all positive parameters. We address a recent extension of CRNT, called generalized mass-action systems, where reaction rates are allowed to be power-laws in the concentrations. In particular, the (real) kinetic orders can differ from the (integer) stoichiometric coefficients. As with mass-action kinetics, complex balancing equilibria are determined by the graph Laplacian of the underlying network and can be characterized by binomial equations and parametrized by monomials. In algebraic terms, we focus on a constructive characterization of positive solutions of polynomial equations with real and symbolic exponents. Uniqueness and existence for all rate constants and initial conditions additionally depend on sign vectors of the stoichiometric and kinetic-order subspaces. This leads to a generalization of Birch's theorem, which is robust with respect to certain perturbations in the exponents. In this context, we discuss the occurrence of multiple complex balancing equilibria. We illustrate our results by a running example and provide a MAPLE worksheet with implementations of all algorithmic methods.

연구 동기 및 목표

• 기존 질량작용 속도 법칙을 넘어서 실수 및 기호적 반응 속도 차수를 가진 시스템으로 화학 반응망 이론(CRNT)을 일반화하는 것.
• 모든 속도 상수와 초기 조건에 대해 일반화된 질량작용계에서 양의 평형 상태의 존재성과 유일성을 보장하는 조건을 설정하는 것.
• 스토이키오메트릭 부분공간과 반응 속도 차수 부분공간 간의 부호 벡터 교차를 통해 다평형성을 특성화하는 것.
• 컴퓨터 대수를 사용하여 이러한 조건을 검증하는 알고리즘적 방법을 제공하고, 연속적인 예시를 통해 설명하는 것.
• 기호적 지수를 가진 일반화된 다항식 시스템으로 버치의 정리를 확장하여, 지수의 변화에 대해 안정성을 보장하는 것.

제안 방법

• 스토이키오메트릭 계수와는 다를 바 있는 실수값의 반응 속도 차수를 가진 반응 복합체를 사용하여 일반화된 질량작용계를 수학적으로 정의한다.
• 반응망의 그래프 라플라시안에서 유도된 이항방정식을 통해 복합 균형 평형 상태를 모델링한다.
• 스토이키오메트릭 부분공간 $S$와 반응 속도 차수 부분공간 $ ilde{S}$의 부호 벡터를 사용하여 양의 해의 존재성과 유일성을 판단한다.
• 오리엔티드 매트로이드 이론을 적용하여 스토이키오메트릭 행렬과 반응 속도 차수 행렬의 치로토프를 비교하고, 부호 벡터의 호환성을 확보한다.
• 행렬의 소수와 정부호 계산을 사용하여 다평형성 조건 $\sigma(S) \cap \sigma(\tilde{S}^\perp) \neq \{0\}$을 검증하는 알고리즘을 구축한다.
• 모든 방법을 재현 가능하고 실용적인 응용을 위해 MAPLE 워크시트에 구현한다.

실험 결과

연구 질문

• RQ1 반응 속도 차수와 속도 상수에 어떤 조건이 성립할 경우 일반화된 질량작용계가 모든 스토이키오메트릭 호환성 클래스에서 유일한 양의 평형 상태를 갖는가?
• RQ2 실수 및 기호적 지수를 가진 다항방정식의 양의 해가 알고리즘적으로 어떻게 존재하고 유일한지를 결정할 수 있는가?
• RQ3 일반화된 질량작용계가 다수의 복합 균형 평형 상태를 가질 수 있는 경우는 언제이며, 이러한 현상을 가능하게 하는 구조적 조건은 무엇인가?
• RQ4 기호적 지수를 가진 시스템에 대해 일반화된 버치의 정리는 어떻게 적용되며, 어떤 변화가 해의 유일성을 유지하는가?
• RQ5 부호 벡터와 오리엔티드 매트로이드 구조는 일반화된 화학 반응망에서 다평형성을 특성화하는 데 어떤 역할을 하는가?

주요 결과

• 연속적인 예시에서 $a, b, c > 0$ 이고 $a < c$ 라면, 스토이키오메트릭 부분공간과 반응 속도 차수 부분공간의 치로토프가 일치하므로 모든 양의 스토이키오메트릭 호환성 클래스에 대해 유일한 복합 균형 평형 상태가 존재한다.
• 만약 $a > c$ 이고 $a, b, c > 0$ 이라면, $\sigma(S) \cap \sigma(\tilde{S}^\perp) \neq \{0\}$ 이면 다수의 복합 균형 평형 상태가 존재할 가능성 있다.
• $\sigma(S) = \sigma(\tilde{S})$ 이고 $(+,\dots,+) \in \sigma(S^\perp)$ 라면, 모든 속도 상수에 대해 양의 평형 상태의 존재성과 유일성이 보장된다.
• 지수의 소규모 변화에 대해 일반화된 버치의 정리는 부호 벡터의 부분공간 호환성이 유지되는 한 안정적으로 성립한다.
• MAPLE 워크시트에 구현된 알고리즘 프레임워크는 일반화된 다항식 시스템에 대해 유일성과 다평형성 조건을 체계적으로 검증할 수 있게 한다.
• 일반화된 질량작용계에서 다평형성은 스토이키오메트릭 부분공간과 반응 속도 차수 부분공간의 부호 벡터가 비자명하게 교차할 경우 발생하며, 이는 질량작용계의 결과를 일반화한다.