대규모 유전자조절망의 진동 조건과 내재 변동성

초록

본 연구는 피드백을 갖는 유전자조절망(GRN)이 내재 변동성 하에서 언제 진동을 유지하거나 소멸하는지를 분석한다. 평균과 두 번째 중심 모멘트를 이용한 2차 모멘트 접근법으로 확률 미분 방정식을 평균‑분산 형태의 ODE로 변환하고, 안정성 및 Hopf 분기 조건을 도출한다. 3‑노드 억제자와 5‑노드 네트워크에 적용한 결과, 시스템 규모가 작아 변동성이 클수록 진동이 억제되며, 일정 규모 이상에서야 지속적인 사이클이 가능함을 확인하였다.

상세 분석

이 논문은 유전자조절망(GRN)의 동역학을 확률적 관점에서 접근한다는 점에서 기존의 결정론적 분석과 차별화된다. 저자들은 ‘두 번째 모멘트 접근법(second‑moment approach)’을 채택해, 각 유전자 발현량 (X_i)의 평균 (\mu_i = \langle X_i\rangle)와 분산·공분산을 나타내는 중앙 2차 모멘트 (C_{ij}= \langle (X_i-\mu_i)(X_j-\mu_j)\rangle)를 독립적인 변수로 만든다. 이때 원래의 마르코프 점프 과정은 Itô‑형식의 확률 미분 방정식으로 기술되며, 평균과 2차 모멘트에 대한 연쇄 법칙을 적용하면 닫힌 형태의 일반 상미분 방정식(ODE) 시스템을 얻는다.

핵심은 이 ODE 시스템의 고정점(stationary point)과 그 주변의 선형화 행렬을 분석해 고유값의 실부와 허부를 조사함으로써, 시스템이 안정 평형(stable equilibrium) 혹은 지속적인 진동(oscillation) 상태에 놓이는지를 판단한다는 점이다. 특히, 고유값이 복소수 쌍으로 실부가 0을 통과하는 Hopf 분기 조건을 명시적으로 도출함으로써, ‘시스템 규모(노드 수)’와 ‘내재 변동성(분산 크기)’ 사이의 정량적 관계를 제시한다.

구체적인 사례로는 전통적인 억제자(repressilator) 3‑노드 모델과, 추가적인 피드백 루프를 포함한 5‑노드 모델을 선택하였다. 두 모델 모두 동일한 형태의 비선형 전사 억제 함수를 사용했으며, 반응 속도 상수와 분해 상수는 동일하게 설정했다. 3‑노드 경우, 노드 수가 작아 평균 분자 수가 수십 개 수준에 머물러 변동성(표준편차/평균 비율)이 30 % 이상에 달한다. 이때 고유값의 실부가 모두 음수이며, Hopf 임계점에 도달하지 못해 진동이 완전히 소멸한다. 반면 5‑노드에서는 평균 분자 수가 증가해 변동성 비율이 15 % 이하로 감소하고, 특정 파라미터 구간에서 고유값의 실부가 양수에서 음수로 바뀌는 Hopf 분기가 관찰된다. 이때 시스템은 작은 진폭의 제한 주기(limit cycle)를 형성하며, 수치 시뮬레이션과 2차 모멘트 해석이 일치한다.

또한 저자들은 ‘임계 최소 규모(critical minimum size)’라는 개념을 도입한다. 이는 변동성에 의해 진동이 억제되는 구간과, 변동성이 충분히 감소해 진동이 유지되는 구간 사이의 경계값이다. 이 임계값은 네트워크 토폴로지(연결 강도, 피드백 루프 수)와 반응 속도 상수의 조합에 따라 달라지며, 분석식으로 근사화할 수 있다. 예를 들어, 억제자 간 결합 강도가 (k)일 때, 임계 노드 수 (N_c)는 (N_c \approx \frac{C}{k^2}) 형태로 표현될 수 있음을 제시한다(여기서 (C)는 시스템 고유의 상수).

결과적으로, 이 연구는 ‘큰 규모일수록 변동성이 평균에 비해 상대적으로 작아지므로, 진동 유지에 유리하다’는 직관을 수학적으로 뒷받침한다. 동시에, 변동성이 큰 소규모 네트워크에서도 진동을 강제로 유지하려면 외부 강제 입력이나 억제 강도 조절 등 추가적인 메커니즘이 필요함을 시사한다. 이러한 통찰은 합성 생물학에서 인공 회로 설계, 그리고 세포 수준에서 자연적인 리듬(예: 세포 주기, circadian rhythm)의 견고성을 이해하는 데 실질적인 가이드라인을 제공한다.