내재적 변동이 지연 시스템에 미치는 영향과 유전자 발현 모델 적용

초록

본 논문은 시스템 크기의 역수 전개를 이용해 내재적 잡음이 포함된 지연 반응 시스템을 분석한다. 결정론적 모델이 진동을 보이지 않는 파라미터 영역에서도, 확률적 효과가 자발적 진동을 유발함을 보이며, 그 진동의 파워 스펙트럼을 해석적으로 도출하고 시뮬레이션과 일치함을 확인한다. 이론을 단순 1차원 모델에 적용한 뒤, 자동 억제 피드백을 갖는 Hes1 mRNA‑단백질 시스템에 적용해 생물학적 의미를 탐구한다.

상세 분석

이 연구는 지연 미분 방정식에 내재적 잡음이 결합된 경우를 다루기 위해, 전통적인 마스터 방정식에서 시스템 크기 N에 대한 역수 전개( van Kampen 확장)를 수행한다. 먼저, 지연이 포함된 반응 네트워크를 확률적 전이율로 정의하고, 평균장( deterministic limit )을 얻기 위해 N→∞ 한계에서 라그랑지 방정식을 도출한다. 이때 결정론적 해는 고정점에 수렴하거나, 파라미터에 따라 제한된 주기적 궤도를 가질 수 있다. 그러나 1/N 차수의 변동항을 보존하면, 선형화된 푸아송 잡음이 지연된 선형 시스템에 결합되어 복소 고유값을 갖는 유효 동역학을 만든다. 특히, 고정점 주변에서의 플럭투에이션이 복소 고유진동수를 갖게 되면, 시스템은 ‘내재적 스토캐스틱 진동(stochastic oscillations)’을 보이며, 이는 결정론적 모델이 안정적인 고정점에 머무는 경우에도 나타난다.

저자들은 이 현상을 정량화하기 위해 푸아송 잡음의 푸리에 변환을 수행하고, 지연 시간 τ와 반응 속도 파라미터 (예: 전사·번역·분해율) 사이의 관계를 포함한 파워 스펙트럼 S(ω)를 해석적으로 계산한다. 스펙트럼은 분모에 (ω²+… ) 형태의 항과 지연에 의해 발생하는 코사인·사인 항이 결합된 복합 구조를 가지며, 특정 파라미터 조합에서 피크가 형성된다. 이 피크의 중심 주파수와 폭은 시스템 크기 N에 반비례적으로 좁아지며, N이 충분히 클 경우 이론적 예측과 Gillespie‑type 지연 시뮬레이션이 거의 일치한다는 점을 실험적으로 입증한다.

이론적 프레임워크는 1차원 ‘toy model’—예를 들어, X → X+1 (생성)와 X → ∅ (소멸) 반응에 지연 τ가 포함된 형태—에 적용되어, 지연이 없는 경우와 비교했을 때 잡음에 의한 진동이 어떻게 나타나는지를 명확히 보여준다. 이어서, 실제 생물학적 시스템인 Hes1 mRNA‑단백질 네트워크에 모델을 확장한다. Hes1은 전사 억제 피드백을 통해 자신의 mRNA 생산을 지연된 방식으로 억제한다. 이 시스템은 일반적으로 낮은 복제 수(N≈10–100)와 긴 지연(τ≈30 분) 때문에 큰 변동성을 보이며, 기존 결정론적 모델은 안정적인 고정점을 예측한다. 그러나 본 논문의 확장된 stochastic delay 이론을 적용하면, 실험적으로 관찰되는 약 2 시간 주기의 진동이 내재적 잡음에 의해 자연스럽게 발생함을 보여준다.

핵심적인 통찰은 다음과 같다. 첫째, 지연은 시스템의 고유주파수를 정의하지만, 잡음은 그 주파수를 실제 진동으로 전환시킨다. 둘째, 시스템 크기가 작을수록 (즉, 분자 수가 적을수록) 잡음의 상대적 강도가 커져 스펙트럼 피크가 더 뚜렷해진다. 셋째, 파워 스펙트럼의 형태는 지연 시간과 반응 속도에 대한 민감도를 제공하므로, 실험 데이터와 비교해 파라미터 추정에 활용될 수 있다. 마지막으로, 이 접근법은 단순 1차원 모델을 넘어, 다변량 지연 네트워크에도 일반화 가능함을 시사한다.