세포 내 생화학 네트워크의 두 모멘트 접근법 심층 탐구

초록

본 논문은 화학 마스터 방정식(CME) 기반의 stochastic 모델링에서 평균과 공분산을 동시에 추적하는 두 모멘트 근사(2MA)를 확장한다. 비정형 반응식과 상대 농도 개념을 도입해 기존 모델의 한계를 보완하고, 효모 세포주기 모델에 적용해 MPF 재설정 현상이 평균‑분산 결합에 의해 어떻게 주기 시간 클러스터링을 일으키는지 설명한다.

상세 분석

두 모멘트 근사(2MA)는 CME에서 파생된 확률적 동역학을 평균(1차 모멘트)과 공분산(2차 모멘트) 수준에서 동시에 기술한다는 점에서 기존의 화학 라그랑지안 방정식(CLE)이나 선형 잡음 근사(LNA)와 근본적으로 차별된다. CLE는 연속적인 확률 미분 방정식으로 평균과 분산을 독립적으로 다루는 경향이 있으며, LNA는 시스템이 큰 복제수에 가까울 때만 정확성을 보장한다. 반면 2MA는 평균과 분산이 비선형적으로 결합되는 경우에도 적용 가능하도록 설계되었으며, 이는 특히 반응 차수가 높은 경우나 비정형(비원소) 반응이 다수 포함된 생화학 네트워크에서 중요한 역할을 한다.

저자들은 기존 2MA 유도 과정을 확장하여 두 가지 주요 제한을 해소하였다. 첫째, 비원소 반응을 하나의 복합 반응식으로 모델링하면서도 해당 반응의 스토키오메트리와 속도 상수를 정확히 반영하도록 수식화하였다. 이는 실제 실험에서 여러 원소 반응이 하나의 효소 촉매 단계로 근사되는 상황을 정량적으로 설명한다. 둘째, 절대 농도가 아닌 상대 농도(예: 세포 부피 대비 농도)를 사용함으로써 세포 성장이나 분열 과정에서 부피 변화가 모델에 미치는 영향을 자연스럽게 포함시켰다. 이러한 확장은 2MA를 보다 일반적인 생물학적 시스템에 적용할 수 있게 만든다.

논문의 핵심 검증 사례는 효모(시스테스) 세포주기 모델이다. 기존 ODE 기반 모델은 MPF(미토시스 촉진 인자)의 급격한 상승과 감소를 평균값만으로 기술했으며, 실험에서 관찰되는 주기 시간의 다중 클러스터링 현상을 설명하지 못했다. 2MA를 적용한 결과, MPF 평균값 근처에서 분산이 급증하고, 이때 평균‑분산 결합 항이 비선형적으로 작용해 MPF가 여러 차례 “리셋”되는 현상이 포착되었다. 즉, G2/M 전이 구간에서 평균값이 임계치를 넘을 때마다 분산이 크게 늘어나 시스템이 불안정해지고, 이는 곧 MPF 농도의 급격한 감소(리셋)로 이어진다. 이러한 리셋 현상이 반복되면서 세포주기의 길이가 서로 다른 몇 개의 뚜렷한 그룹으로 나뉘는 것이 실험 데이터와 일치한다.

또한 저자들은 2MA의 수치적 안정성을 검증하기 위해 Gillespie 알고리즘 기반의 정확한 stochastic 시뮬레이션과 비교하였다. 평균값과 공분산 모두에서 오차가 5% 이내에 머물렀으며, 특히 변동성이 큰 구간에서 LNA보다 현저히 높은 정확도를 보였다. 이는 2MA가 비선형 결합 효과를 포착함으로써 실제 stochastic 동역학을 보다 정밀하게 근사한다는 강력한 증거가 된다.

결론적으로, 이 연구는 2MA가 비정형 반응과 상대 농도라는 현실적인 제약을 포함하면서도, 복잡한 세포 주기와 같은 다중 안정 상태 시스템을 해석하는 데 유용한 도구임을 입증한다. 향후 다른 세포 신호 전달 경로나 대사 네트워크에도 동일한 접근법을 적용하면, 잡음이 기능적 역할을 하는 메커니즘을 정량적으로 규명할 수 있을 것으로 기대된다.