[논문 리뷰] Dynamical Equivalence and Linear Conjugacy of Chemical Reaction Networks: New Results and Methods

Matthew D. Johnston, David Siegel외 1명|arXiv (Cornell University)|2011. 11. 06.

Computational Drug Discovery Methods참고 문헌 29인용 수 13

한 줄 요약

이 논문은 질량작용 속도 법칙 하에서 조밀한, 희박한, 가역적인, 세부 균형 잡힌, 복합 균형 잡힌 선형 콘주게이트 반응 네트워크를 계산하기 위해 혼합정수선형계획법(MILP) 방법을 도입한다. 조밀한 선형 콘주게이트 네트워크가 복합체 집합에 의해 유일하게 결정됨을 증명하고, 평형 선택이 복합체 또는 세부 균형 잡힌 실현의 존재성에 영향을 주지 않음을 보이며, 이는 속도 상수가 사전에 알려지지 않은 상태에서 반응 메커니즘의 구조 기반 분석을 가능하게 한다.

ABSTRACT

In the first part of this paper, we propose new optimization-based methods for the computation of preferred (dense, sparse, reversible, detailed and complex balanced) linearly conjugate reaction network structures with mass action dynamics. The developed methods are extensions of previously published results on dynamically equivalent reaction networks and are based on mixed-integer linear programming. As related theoretical contributions we show that (i) dense linearly conjugate networks define a unique super-structure for any positive diagonal state transformation if the set of chemical complexes is given, and (ii) the existence of linearly conjugate detailed balanced and complex balanced networks do not depend on the selection of equilibrium points. In the second part of the paper it is shown that determining dynamically equivalent realizations to a network that is structurally fixed but parametrically not can also be written and solved as a mixed-integer linear programming problem. Several examples illustrate the presented computation methods.

연구 동기 및 목표

• 특정 역학적 성질인 가역성, 세부 균형, 복합 균형 등을 갖는 선호되는 선형 콘주게이트 반응 네트워크 구조를 최적화 기반으로 계산하기 위한 방법을 개발한다.
• 기존의 MILP 프레임워크를 확장하여 복합체, 반응, 복합체의 최소 및 최대 수를 갖는 네트워크를 구조적 제약 조건 하에서 처리할 수 있도록 한다.
• 네트워크 구조는 고정되어 있으나 속도 상수가 알려지지 않은 경우의 역학적 동치성 문제를 다루며, 구조적 분석만 가능하게 한다.
• 조밀한 선형 콘주게이트 네트워크의 유일성과 복합체/세부 균형 잡힘 존재성이 다양한 평형점에서 불변하는 이론적 결과를 증명한다.
• 매개변수 없이도 구조적 분석이 가능한 역학적 동치성 및 콘주게이트성의 적용 범위를 넓히기 위해 반응 네트워크 이론을 확장한다.

제안 방법

• 구조적 희박성과 제약 조건을 활용하여 선형 콘주게이트 네트워크 탐색 문제를 혼합정수선형계획법(MILP) 문제로 수립한다.
• 스토이히오메트릭 행렬과 복합체 행렬 표현을 사용하여 네트워크의 구조와 역학을 표현하며, 속도 상수를 최적화 변수로 포함한다.
• 선형 등식 및 부등식 조건을 통해 약한 가역성, 가역성, 세부 균형, 복합 균형을 강제하는 제약 조건을 적용한다.
• 비자명한 선형 콘주게이트를 통한 변환 기반 접근법을 도입하여 네트워크 간 역학을 유지하면서 동적 행동을 보존한다.
• 양의 대각선 상태 변환을 활용하여 역학적 동치성과 콘주게이트성을 연결하고, 평형에 영향을 받지 않는 분석을 가능하게 한다.
• 속도 상수가 고정되지 않은 경우에도 MILP 프레임워크를 확장하여 특정 운동학적 매개변수에 의존하지 않는 구조적 분석이 가능하도록 한다.

실험 결과

연구 질문

• RQ1 최적화를 통해 특정 역학적 성질(예: 가역성, 세부 균형)을 갖는 가장 희박하고 가장 조밀한 선형 콘주게이트 네트워크를 계산할 수 있는가?
• RQ2 가장 조밀한 선형 콘주게이트 네트워크의 구조는 복합체 집합과 양의 대각선 상태 변환에 의해 유일하게 결정되는가?
• RQ3 선형 콘주게이트 복합 균형 잡힌 또는 세부 균형 잡힌 네트워크의 존재성은 평형점의 선택에 의존하는가?
• RQ4 속도 상수가 알려지지 않은 상태에서 네트워크 구조만 알려졌을 때 역학적 동치성이 확인되고 계산될 수 있는가?
• RQ5 질량작용 속도 법칙 하에서 역학을 유지하지만 반응 구조가 다른 대체 실현을 체계적으로 생성할 수 있는 방법은 무엇인가?

주요 결과

• 조밀한 선형 콘주게이트 네트워크는 속도 상수에 관계없이 복합체 집합과 양의 대각선 상태 변환에 의해 유일하게 결정된다.
• 선형 콘주게이트 세부 균형 잡힌 및 복합 균형 잡힌 네트워크의 존재성은 평형점의 선택과 무관하여, 임의의 양의 평형에서의 구축이 가능하다.
• 모든 선형 콘주게이트 네트워크의 가중치 없는 방향성 반응 그래프는 가장 조밀한 네트워크의 부분 그래프이므로, 구조적 계층성을 보장한다.
• 약한 가역성, 가역성, 복합체/세부 균형 잡힘 조건 하에서 희박 및 조밀 실현을 찾는 문제는 MILP로 공식화할 수 있으며, 이는 효율적인 계산을 가능하게 한다.
• 속도 상수가 고정되지 않은 경우에도 복합체 수가 최소 또는 최대인 동적 동치 실현의 계산은 여전히 MILP 문제로 남아 있어, 구조적 분석이 가능하다.
• 동적 동치 실현과 동일한 구조를 갖는 복합 균형 네트워크는 속도 상수 재할당을 통해 구성 가능하며, 이는 네트워크 설계에서 매개변수 유연성을 보여준다.