식물 시계의 새로운 통합 접근법

초록

녹조류 Ostreococcus tauri의 시계 유전망을 Bio‑PEPA로 모델링하고, 연속 미분방정식, 확률적 시뮬레이션, 모델 검증을 동시에 적용했다. 분자 수가 적은 시스템에서 연속식 접근법의 한계를 확인하고, 단일 세포가 일정광에서도 자체 진동을 유지할 수 있음을 제시한다. 또한 새벽·황혼 시점을 모델 검증에 포함시켜 주요 단백질·mRNA 농도의 확률분포 변화를 분석하고, 피크 발현 시 변동성이 최소임을 발견했다. 마지막으로 mRNA 분해 속도 변이가 단백질 위상에 미치는 영향을 진화 시스템생물학적 관점에서 탐구해 시계의 돌연변이 내성을 평가한다.

상세 분석

본 연구는 Ostreococcus tauri라는 녹조류의 원시적인 식물 시계 모델을 Bio‑PEPA라는 확률 과정 대수 체계에 구현함으로써, 동일 모델을 세 가지 전혀 다른 분석 파이프라인에 매핑한다는 독창적인 방법론을 제시한다. 첫 번째 파이프라인은 전통적인 연속 미분방정식(ODE) 접근법으로, 대량의 분자 집합을 가정하고 평균 동역학을 기술한다. 그러나 Ostreococcus 시계는 세포당 수백 개 수준의 mRNA와 단백질만을 포함하고 있어, ODE가 가정하는 연속성 근사가 실제 생물학적 변동성을 충분히 포착하지 못한다는 점을 실험적으로 확인한다. 두 번째 파이프라인은 Gillespie 알고리즘 기반의 확률적 시뮬레이션으로, 개별 분자 수준의 무작위성을 직접 재현한다. 이를 통해 일정광(constant light) 조건에서도 단일 세포가 지속적인 진동을 보일 가능성을 제시했으며, 이는 기존에 대량 집단 평균만을 관찰해 온 연구와는 다른 새로운 가설이다. 세 번째 파이프라인은 PRISM 모델 체커를 이용한 모델 검증이다. 여기서는 ‘새벽(dawn)’과 ‘황혼(dusk)’이라는 외부 신호를 논리적 이벤트로 형식화하고, 시간에 따른 주요 종(예: CCA1, TOC1)의 확률분포를 정량적으로 추정한다. 결과는 피크 발현 시점에 변동계수가 최소가 되므로, 실험적 위상 마커로서 피크 시간을 선택하는 것이 가장 효율적임을 뒷받침한다. 마지막으로, 진화 시스템생물학적 분석을 수행해 mRNA 분해 속도(k_deg) 변이가 핵심 단백질 위상(phase)에 미치는 영향을 파라미터 스위프와 민감도 분석으로 조사했다. 시계는 k_deg 변화에 대해 일정 범위 내에서 위상이 크게 변하지 않아, 돌연변이 발생 시에도 기능을 유지할 수 있는 내재적 견고성을 가지고 있음을 보여준다. 전체적으로, 연속‑확률‑형식 검증이라는 삼중 접근법이 각각의 한계를 보완하고, 시계 메커니즘에 대한 다층적 이해를 가능하게 한다는 점이 가장 큰 학술적 기여라 할 수 있다.