[논문 리뷰] Reachability in Biochemical Dynamical Systems by Quantitative Discrete Approximation

Luboš Brim, J. Fabriková외 2명|arXiv (Cornell University)|2011. 07. 29.

Gene Regulatory Network Analysis참고 문헌 20인용 수 8

한 줄 요약

이 논문은 일반 미분 방정식(OEDs)로 모델링된 생화학적 동역학 시스템을 위한 매개변수화된 정량적 이산 근사(QDAA) 방법을 제안한다. 이 방법은 유한한 상태 수로 제한된 근사치를 통해 정확한 도달 가능성 분석을 가능하게 하며, 직사각형 추상화에서 발생하는 잘못된 전이를 줄이기 위해 국소 벡터장의 행동을 통합함으로써 도달 가능한 종 농도에 대한 더 날카운 봉인을 달성한다. 기존 방법보다 더 좁은 범위를 제공한다.

ABSTRACT

In this paper a novel computational technique for finite discrete approximation of continuous dynamical systems suitable for a significant class of biochemical dynamical systems is introduced. The method is parameterized in order to affect the imposed level of approximation provided that with increasing parameter value the approximation converges to the original continuous system. By employing this approximation technique, we present algorithms solving the reachability problem for biochemical dynamical systems. The presented method and algorithms are evaluated on several exemplary biological models and on a real case study. This is a full version of the paper published in the proceedings of CompMod 2011.

연구 동기 및 목표

• 직사각형 추상화의 한계를 해결하기 위해, 직사각형 경계에서 벡터장이 과도하게 근사화되어 발생하는 과도한 잘못된 행동을 줄이기 위한 것이다.
• 매개변수를 증가시킬수록 원래 연속 시스템에 수렴하는 매개변수화된 이산 근사 기법을 개발하여 확장 가능하고 정확한 도달 가능성 계산을 가능하게 하기 위한 것이다.
• 초기 조건과 매개변수의 불확실성 하에서 상태 공간의 도달 가능성에 대한 전역 분석을 가능하게 함으로써 수치 시뮬레이션의 계산 비용을 줄이는 효율적인 대안을 제공하기 위한 것이다.
• 모델 체킹에서 잘못된 반례의 정도를 정량화하고 줄이기 위해 근사 품질 측도를 제공하기 위한 것이다.
• 미래의 생물학적 시스템 모델 체킹 도구에 이 방법을 통합하여 복잡하고 고차원적인 동역학 시스템의 분석 능력을 향상시키기 위한 것이다.

제안 방법

• 연속 상태 공간을 직사각형으로 분할하여 유한한 이산 근사를 구성하지만, 이는 직사각형 내외 경계에서의 벡터장 행동에 대한 정량적 측도를 통합함으로써 향상된다.
• 근사 수준을 제어하는 매개변수화된 정밀화 과정을 도입하여 매개변수를 증가시킬수록 원래 ODE 시스템에 수렴하도록 보장한다.
• 도달 가능한 집합의 정밀도를 높이기 위해 후행 시뮬레이션을 사용하며, 분할 밀도를 증가시키지 않아도 종 농도의 간격 범위를 개선한다.
• 직사각형 전이 시스템 추상화를 적용하지만, 전이 규칙을 국소 벡터장 역학에 맞게 수정하여 잘못된 전이를 줄인다.
• 알고리즘 2를 사용하여 도달 가능한 집합의 과도한 근사를 계산하며, 이는 정방향 및 역방향 시뮬레이션을 통합하여 상태 변수의 범위를 좁힌다.
• 프로토타입 소프트웨어 구현을 통해 소규모 모델과 *E. coli*에서의 암모니아 수송의 실제 사례 연구에 대해 방법을 평가하였으며, 결과를 수치 시뮬레이션과 직사각형 추상화와 비교하였다.

실험 결과

연구 질문

• RQ1 생화학적 ODE 시스템의 도달 가능성 분석에서 잘못된 전이를 줄일 수 있는 매개변수화된 이산 근사 방법을 개발할 수 있는가?
• RQ2 제안된 QDAA 방법은 기존의 표준 직사각형 추상화에 비해 도달 가능한 집합 추정의 정확도를 어떻게 향상시키는가?
• RQ3 분할 밀도를 증가시키지 않아도 후행 시뮬레이션은 도달 가능한 간격을 얼마나 잘 개선하고 과도한 근사를 줄일 수 있는가?
• RQ4 전통적인 추상화 기법보다 더 좁은 범위의 종 농도(예: NH₄⁺in 및 NH₃in)를 제공하면서도 계산적으로 실현 가능한가?
• RQ5 QDAA의 근사 품질은 선택된 매개변수 값과 어떻게 관련되어 있으며, 진정한 연속적 행동으로 수렴하는가?

주요 결과

• 암모니아 수송 모델에서 QDAA는 NH₄⁺in의 상한을 5.4×10⁻⁴로 추정하였으며, 이는 직사각형 추상화의 추정치 5.4×10⁻⁴과 유사하지만, 과도한 근사화가 줄어들어 수치 시뮬레이션 결과보다 한 단계 더 가까운 정확도를 보였다.
• QDAA 방법은 NH₃in에 대해 [10⁻⁸, 1.1×10⁻⁶]의 도달 가능한 간격을 도출하였으며, 이는 초깃값의 범위와 일치하며 시뮬레이션 데이터와도 일관성을 보였다.
• NH₄⁺in의 경우 QDAA 결과인 [10⁻⁶, 10⁻⁵]는 수치 시뮬레이션 결과보다 기존의 직사각형 추상화에 비해 한 단계 더 가까운 정밀도를 보이며 정밀도 향상을 입증하였다.
• 이 방법은 특히 직사각형 경계에서 국소 벡터장 역학을 통합함으로써 도달 가능성 분석에서 잘못된 행동을 성공적으로 줄였으며, 더 미세한 분할을 요구하지 않았다.
• 알고리즘 2에서의 후행 시뮬레이션은 NH₃in 및 NH₄⁺in과 같이 초깃값이 제한된 종의 경우에 특히 범위를 좁히는 데 뚜렷한 기여를 하였다.
• 프로토타입 구현은 QDAA가 실제 생물학적 모델에 효과적으로 적용될 수 있으며, 형식적 검증과 시스템 생물학 분석에 적합한 신뢰할 수 있는 과도한 근사를 제공할 수 있음을 보여주었다.