[논문 리뷰] Reduction of chemical systems by delayed quasi-steady state assumptions

Tomáš Vejchodský, Radek Erban외 1명|arXiv (Cornell University)|2014. 06. 17.

Gene Regulatory Network Analysis참고 문헌 28인용 수 6

한 줄 요약

이 논문은 복잡한 화학계를 단순화하기 위해 빠른 종의 동역학을 지연 미분 방정식으로 대체하는 지연 준정적상가정(D-QSSA)을 소개한다. 이는 기존 준정적상가정(QSSA)보다 정확도를 향상시키기 위한 새로운 방법이다. D-QSSA는 빠른 종이 정적상에 도달하는 데 소요되는 시간 지연을 고려하여 유한차원 근사값을 제공하며, 오차는 지수적으로 감소한다. 이는 Hes1 및 세포주기 시스템과 같은 진동하는 모델에서 기존 QSSA보다 오차가 5배 이상 낮아, 뛰어난 성능을 보인다.

ABSTRACT

Mathematical analysis of mass action models of large complex chemical systems is typically only possible if the models are reduced. The most common reduction technique is based on quasi-steady state assumptions. To increase the accuracy of this technique we propose delayed quasi-steady state assumptions (D-QSSA) which yield systems of delay differential equations. We define the approximation based on D-QSSA, prove the corresponding error estimate, and show how it approximates the invariant manifold. Then we define a class of well mixed chemical systems and formulate assumptions enabling the application of D-QSSA. We also apply the D-QSSA to a model of Hes1 expression and to a cell-cycle model to illustrate the improved accuracy of the D-QSSA with respect to the standard quasi-steady state assumptions.

연구 동기 및 목표

• 진동하는 시스템에서 특히 임시 동역학을 포착하지 못하는 기존 준정적상가정(QSSA)의 부정확성을 해결한다.
• 빠른 종이 즉각적으로 평형에 도달한다는 가정을 하는 QSSA의 한계를 극복하기 위해 시간 지연을 도입한다.
• 식별 가능한 빠른 종과 느린 종을 가진 균일한 혼합 질량작용 법칙 모델에 D-QSSA를 적용할 수 있는 체계적 프레임워크를 개발한다.
• 특정 조건 하에서 D-QSSA 근사의 오차가 지수적으로 0으로 수렴함을 보여주는 엄밀한 오차 추정을 제공한다.
• 실제 생물학적 모델, 특히 Hes1 발현 및 세포주기 조절에서 표준 QSSA가 진동을 유지하지 못하는 데서 나타나는 D-QSSA의 효과를 입증한다.

제안 방법

• D-QSSA를 빠른 ODE를 지연 미분 방정식(DDE)으로 대체하는 감소 기법으로 정의하여, 빠른 종이 준정적상에 도달하는 데 소요되는 시간 지연을 모델링한다.
• 빠른 부분계에서 가장 큰 속도 상수의 역수로부터 유도된 일정한 지연을 사용하여 D-QSSA 근사를 수립한다.
• D-QSSA의 오차 추정을 증명하여, 특정 조건 하에서 순시 오차가 시간이 지남에 따라 지수적으로 감소함을 보여준다.
• 먼저 선형 ODE 시스템에 D-QSSA를 적용한 후, 빠른 종이 직접적으로 서로 영향을 주지 않는다고 가정하여 비선형 질량작용 시스템으로 확장한다.
• 초기 조건의 일정한 연장 기법을 사용하여 무한차원 상태공간을 유지하면서도 유한차원 상태공간을 유지함으로써 감소 목표를 보존한다.
• 두 가지 생물학적 모델에서 방법을 검증하였다: 9변수 Hes1 진동자와 3변수 세포주기 모델로, 각각 2개와 1개의 DDE로 감소시켰다.

실험 결과

연구 질문

• RQ1 준정적상가정에 시간 지연을 도입하면 화학동역학에서 모델 감소의 정확도가 향상될 수 있는가?
• RQ2 D-QSSA 근사의 오차는 시간이 지남에 따라 어떻게 변화하는가? 그리고 분석적으로 경계하거나 추정할 수 있는가?
• RQ3 특히 진동하는 시스템에서 D-QSSA는 어떤 종류의 화학계에서 표준 QSSA를 능가하는가?
• RQ4 D-QSSA는 표준 QSSA에서 잃어버리는 진동 행동을 유지할 수 있는가, 특히 피드백 조절 시스템에서?
• RQ5 다양한 지연이 D-QSSA의 정확도에 미치는 영향은 무엇이며, 실생활에서 일정한 지연 근사를 신뢰성 있게 사용할 수 있는가?

주요 결과

• D-QSSA는 ODE 시스템을 지연 미분 방정식(DDE) 시스템으로 감소시키며, 일정한 역사 초기화를 통해 유한차원 상태공간을 유지한다.
• D-QSSA 근사의 오차는 시간이 지남에 따라 지수적으로 감소하며, 선형 시스템에 대해 엄밀한 상한이 유도되었다.
• Hes1 모델에서 D-QSSA는 표준 QSSA에서 약 30% 수준이었던 주기 오차를 약 1–2%로 감소시켜, 진동 정확도 향상에 있어 뚜렷한 개선을 보였다.
• 세포주기 모델에서는 표준 QSSA가 파괴하는 지속적인 진동을 D-QSSA가 성공적으로 복원하여, 진동이 지연된 부정성 피드백에 의해 유도됨을 시사한다.
• D-QSSA를 통해 유도된 지연은 이전 연구에서 세밀한 분석을 통해 도출된 값과 일치하여, 이 방법의 기존 결과와의 일관성을 검증한다.
• 수치 실험 결과, 시간에 따라 변하는 지연을 일정한 근사로 대체해도 정확한 결과를 얻을 수 있었지만, 최적의 지연 선택은 민감하며 향후 추가 연구가 필요하다.