[논문 리뷰] Analysis of the Conradi-Kahle Algorithm for Detecting Binomiality on Biological Models

Alexandru Iosif, Hamid Rahkooy|arXiv (Cornell University)|2019. 12. 14.

Computational Drug Discovery Methods참고 문헌 12인용 수 4

한 줄 요약

이 논문은 생물학적 모델에서 유도된 다항식 아이디얼의 이진형성(binomiality)을 탐지하기 위한 Conradi-Kahle 알고리즘을 분석한다. 이 알고리즘은 Macaulay2와 Maple에 구현되어 있으며, BioModels 레포지토리의 20개 생물모델에 대해 그뢰브너 기저(Gröbner bases) 계산과 성능을 비교한다. 주요 발견은 Conradi-Kahle 알고리즘이 대부분의 경우 1초 이내로 매우 빠르게 작동하지만, 이질적인 아이디얼에서의 완전성 부족으로 인해 오직 3개의 모델에서만 정확하게 이진형성을 식별한다는 것이다. 반면, 더 느리지만 600초 이내에 종료되는 경우에 한해 그뢰브너 기저는 비이진형성을 탐지할 수 있다.

ABSTRACT

We analyze the Conradi-Kahle Algorithm for detecting binomiality. We present experiments using two implementations of the algorithm in Macaulay2 and Maple on biological models and assess the performance of the algorithm on these models. We compare the two implementations with each other and with Gröbner bases computations up to their performance on these biological models.

연구 동기 및 목표

• 생물학적 모델의 다항식 아이디얼에서 이진형성을 탐지하기 위한 Conradi-Kahle 알고리즘의 효율성과 효과성을 분석하는 것.
• Macaulay2와 Maple에 구현된 Conradi-Kahle 알고리즘의 성능을 실제 생물학적 모델에서 그뢰브너 기저 계산과 비교하는 것.
• 이질적인 아이디얼에 적용되었을 때 Conradi-Kahle 알고리즘의 한계를 조사하는 것. 생물학적 반응 네트워크에서 흔한 이질적 아이디얼을 고려한다.
• 그뢰브너 기저가 합리적인 시간 내에 종료되지 않을 경우에도 알고리즘이 생화학적 네트워크의 평형 상태 아이디얼에서 이진형성을 신뢰성 있게 결정할 수 있는지 평가하는 것.

제안 방법

• 20개의 생물모델은 BioModels 레포지토리에서 가져왔으며, 속도 계수 매개변수에 대한 유리함수 체 위의 이질적 다항식 아이디얼이다.
• 알고리즘은 계수 행렬의 기저 행렬 기반의 축소된 행기저형을 반복적으로 적용하고, 3개 이상의 비영원소를 가진 행이 있는지 검사함으로써 이진형성을 테스트한다.
• 기준으로서 그뢰브너 기저 계산을 사용하며, Eisenbud와 Sturmfels의 보조정리 1.2를 적용하여 그뢰브너 기저 감소를 통한 이진형성 테스트를 수행한다.
• Macaulay2와 Maple에 모두 알고리즘을 구현하고 실행 시간을 기록하며, 두 시스템 간의 결과를 비교한다.
• 원래의 이질적 아이디얼의 동차화된 형태에 알고리즘을 적용하며, 동차 아이디얼에 대해서만 완전성이 보장된다는 점을 인지한다.
• 성능 평가 시에는 CPU 시간과 이진형성 결정 성공률을 측정하며, 그뢰브너 기저의 경우 600초 타임아웃을 설정한다.

실험 결과

연구 질문

• RQ1 생물학적 모델의 이질적 다항식 아이디얼에서 이진형성을 탐지할 때 Conradi-Kahle 알고리즘의 속도와 정확도는 어떻게 되는가?
• RQ2 실제 생물학적 모델에서 Conradi-Kahle 알고리즘이 그뢰브너 기저보다 얼마나 뛰어나거나 열등한가?
• RQ3 효율성은 높지만 대부분의 생물모델에서 이진형성을 결정하지 못하는 이유는 무엇인가?
• RQ4 Conradi-Kahle 알고리즘은 비이진형 아이디얼을 탐지할 수 있는가, 아니면 오직 이진형성을 확인하는 데만 가능한가?
• RQ5 Macaulay2와 Maple에 구현된 Conradi-Kahle 알고리즘의 계산 성능은 생물학적 모델에서 어떻게 비교되는가?

주요 결과

• 모든 생물모델에 대해 Conradi-Kahle 알고리즘이 1초 이내로 실행되며, Macaulay2에서는 0.01~0.7초, Maple에서는 0.02~6초의 실행 시간을 기록한다.
• 알고리즘은 오직 3개의 생물모델(9, 262, 264)에서만 이진형성을 정확히 식별하며, Bin (C-K) 열에 'Yes'를 반환하고 나머지 모든 모델에 대해서는 'No'를 반환한다.
• 그뢰브너 기저 계산은 20개의 생물모델 중 10개에서만 600초 이내에 종료되며, 가장 오래 걸린 것은 Maple에서 biomodel 153에 대해 500초였다.
• biomodel 46, 187, 243, 335의 경우 그뢰브너 기저는 비이진형성을 정확히 탐지했지만, Conradi-Kahle 알고리즘은 이질적 아이디얼에서의 불완전성으로 인해 이를 탐지하지 못했다.
• 양쪽 방법이 성공한 경우, 실행 시간은 유사한 편이었다: biomodel 264의 경우 C-K (M2)는 0.7초, GB (M2)는 2초였고, C-K (Maple)는 0.03초, GB (Maple)는 0.04초였다.
• Macaulay2와 Maple의 구현 간에 Conradi-Kahle 알고리즘 또는 그뢰브너 기저 계산에서 성능 차이가 없었다.