세포 소기관 내 효소 반응식의 견고성 재검토

초록

이 논문은 작은 세포 소기관에서 효소 반응이 전통적인 대수식(맥락식)으로는 정확히 설명되지 않음을 보이고, 분자 수가 적은 경우와 물질 운반 방식(확산·소포 매개)에 따른 내재적 잡음이 반응 속도에 미치는 영향을 정량화한다. 200 nm 이하의 미세 구획에서는 기존 모델과 수백 퍼센트 차이가 발생할 수 있음을 제시한다.

상세 분석

본 연구는 효소 반응을 마스터 방정식으로 기술한 뒤, 대사 정상 상태(steady‑state) 가정 하에 평균값과 분산을 포함하는 mesoscopic rate equation을 도출한다. 핵심은 두 가지 비전통적 요인을 명시적으로 모델링한 점이다. 첫째, 소기관 부피가 작아 분자 복제수가 수십 개 수준으로 감소하면, 확률적 플럭투에이션이 반응 속도에 비선형적으로 결합한다. 이를 위해 시스템 크기 확장법(system‑size expansion)을 적용해 1/V 차수까지의 교정항을 포함시켰으며, 결과 식은 기존 미카엘리스‑멘텐(Michaelis‑Menten) 식에 잡음‑보정항을 더한 형태가 된다. 둘째, 기질이 단순 확산이 아니라 소포·과립에 의해 운반될 경우, 입자 도착 간격이 포아송 과정이 아닌 복합 점프 프로세스로 모델링된다. 이 경우 입자 도착률이 시간에 따라 변동하고, 입자 집합이 한 번에 다수 도착하는 ‘burst’ 현상이 발생한다. 저자들은 이러한 burst 현상을 반영해 입자 도착 확률 분포를 지수‑감쇠형으로 가정하고, 그 결과 효소‑기질 복합체 형성 확률이 급격히 변동함을 보였다. 특히 경쟁 억제제(I) 존재 시, 억제제와 기질이 동시에 소포에 탑재될 확률이 낮아 억제 효과가 과소평가되며, 이는 전통적인 억제 상수(Ki) 추정에 큰 오차를 초래한다. 수치 시뮬레이션은 부피가 50 nm에서 200 nm 사이일 때 반응 속도가 기존 모델 대비 150 %~300 % 차이 나는 것을 확인했으며, 부피가 무한대로 커질 경우(즉, V→∞) 잡음‑보정항이 사라져 기존 미카엘리스‑멘텐 식으로 수렴한다. 이러한 결과는 세포 내 소기관(예: 시냅스 소포, 리소좀, 골지체 등)에서 효소‑기질 상호작용을 ODE 기반 네트워크 모델로 단순화하는 것이 위험할 수 있음을 경고한다.