[논문 리뷰] Tutorial of numerical continuation and bifurcation theory for systems and synthetic biology

M. G. Blyth, Ludovic Renson외 1명|arXiv (Cornell University)|2020. 08. 12.

Gene Regulatory Network Analysis인용 수 4

한 줄 요약

이 가이드는 비선형 역학을 연구하는 데 유용한 계산 도구로 수치적 연속성과 분기 분석을 소개한다. 이론적 배경을 설명하고 뉴런 모델을 활용한 응용 사례를 제시하며, 프로그래밍 경험, 성능 요구사항, 시스템 복잡도에 따라 적합한 소프트웨어 도구를 선택하는 데 도움을 주며, 생물학적 모델링에서 재현 가능성과 사용성에 중점을 둔다.

ABSTRACT

Mathematical modelling allows us to concisely describe fundamental principles in biology. Analysis of models can help to both explain known phenomena, and predict the existence of new, unseen behaviours. Model analysis is often a complex task, such that we have little choice but to approach the problem with computational methods. Numerical continuation is a computational method for analysing the dynamics of nonlinear models by algorithmically detecting bifurcations. Here we aim to promote the use of numerical continuation tools by providing an introduction to nonlinear dynamics and numerical bifurcation analysis. Many numerical continuation packages are available, covering a wide range of system classes; a review of these packages is provided, to help both new and experienced practitioners in choosing the appropriate software tools for their needs.

연구 동기 및 목표

• 비선형 역학에 대한 배경 지식이 없는 시스템 및 합성 생물학 연구자들이 수치적 연속성과 분기 분석을 쉽게 접근할 수 있도록 하는 것.
• 분기 분석이 뉴런 동역학에서의 번개 발화 시작과 같은 중요한 전이를 드러내는 방법을 보여주는 것.
• 사용자 전문성, 성능, 시스템 유형에 따라 도구 선택에 도움이 되는 기존 수치적 연속성 소프트웨어 패키지 간의 실용적 비교를 제공하는 것.
• 해석적으로 해결하기 어려운 비선형 미분방정식을 분석하기 위한 계산 방법의 사용을 장려하는 것.
• 스크립팅 기능과 생물학적 응용에 특화된 도구상자 기능을 강조하여 재현 가능하고 확장 가능한 연구를 지원하는 것.

제안 방법

• 비선형 역학의 기초 개념, 예를 들어 상태 변수, 시간에 따른 진화, 그리고 시스템의 정성적 행동(예: 평형점, 진동, 카오스)을 소개한다.
• 매개변수 변화에 따른 시스템 동역학의 정성적 변화로 분기를 정의하며, 식량망 모델에서의 캐나드 폭발과 같은 예시로 설명한다.
• 전도도 기반 뉴런 모델을 사용하여 수치적 연속성이 분기 경로를 추적하고 번개 발화 행동을 설명하는 방법을 시연한다.
• 15개 이상의 수치적 연속성 소프트웨어 패키지를 검토하며, 시스템 유형(ODE, PDE, 맵, 하이브리드, 비연속 시스템)과 인터페이스 유형(GUI, 스크립팅, 라이브러리 기반)으로 분류한다.
• 사용성, 성능, 분기 탐지 능력, 확장성 측면에서 핵심 도구(XPPAUTO, MatCont, PyDSTool, CoCo)를 비교한다.
• 재현 가능성과 협업을 위해 스크립팅의 중요성을 강조하며, 시스템 생물학 응용에 특화된 도구상자 기능을 갖춘 PyDSTool의 장점을 설명한다.

실험 결과

연구 질문

• RQ1 수치적 연속성은 뉴런 활동을 기술하는 비선형 생물학적 모델에서 분기를 탐지하고 분석하는 데 어떻게 활용될 수 있는가?
• RQ2 주요 수치적 연속성 소프트웨어 패키지 간의 주요 차이점은 무엇인가? 사용성, 성능, 기능 측면에서 비교해보라.
• RQ3 기본적인 프로그래밍 경험만 있는 연구자와 고급 프로그래밍 능력을 가진 연구자에게 가장 적합한 소프트웨어 도구는 무엇인가?
• RQ4 다양한 소프트웨어 패키지는 주기 궤도나 하이브리드 시스템에서의 슬라이딩 분기와 같은 복잡한 분기 구조를 어떻게 다루는가?
• RQ5 수치적 연속성 도구는 시스템 생물학에서 합성 유전자 네트워크의 설계 및 분석을 어떻게 향상시킬 수 있는가?

주요 결과

• 수치적 연속성은 표준 시뮬레이션으로는 포착할 수 없는 비선형 생물학적 모델에서 분기를 탐지할 수 있으며, 예를 들어 뉴런 모델에서 번개 발화로의 전이를 포함한다.
• XPPAUTO와 PyDSTool는 그래픽 사용자 인터페이스와 단계별 평면 분석 도구를 제공하므로, 프로그래밍 경험이 적은 사용자에게 적합하다.
• 고성능 계산과 복잡한 분기 추적을 필요로 하는 경우 MatCont가 가장 적합하다. 다양한 유형의 분기를 탐지하고 추적할 수 있기 때문이다.
• PyDSTool과 CoCo는 확장성과 대규모 워크플로우 통합이 필요한 고급 사용자에게 이상적이다. 특히 PyDSTool는 시스템 생물학 응용에 특화된 도구상자 기능을 제공한다.
• 스크립팅 인터페이스를 갖춘 소프트웨어(예: PyDSTool, CoCo)는 재현 가능하고 협업이 가능한 연구를 지원하며, 분석의 적응 및 자동화를 용이하게 한다.
• 이 가이드는 PDE 기반 모델에 적합한 PDECONT와 pde2path, PDE 및 비연속 시스템에 각각 특화된 WAVETRAIN과 SlideCont를 식별하여, 가용 도구의 다양성을 부각시킨다.