DNA 상태 변동이 자기조절 유전자 단일세포 동역학에 미치는 영향

초록

동적 평균장 이론을 도입해 유전자 발현의 단일세포 확률적 동역학을 분석하였다. DNA 상태 변동의 비평형·비정상(adabatic) 특성을 명시적으로 고려함으로써 단일세포 수준의 두 시점 상관함수와 응답함수를 유도하였다. 이 방법을 자기조절 유전자 모델에 적용한 결과, 완화 시간이 비정상적으로 증가하거나 상관이 진동적으로 감쇠하는 등 다양한 동역학 현상이 예측되었다. 단백질 수의 상관과 응답의 비율로 정의되는 유효 “온도”는 DNA 상태 전이가 빈번할 때는 작게, 전이가 드물어질수록 크게 나타나며, 후자의 경우 잡음이 크게 강화됨을 시사한다.

상세 요약

이 논문은 유전자 발현 과정에서 DNA가 전사활성 상태와 억제 상태 사이를 전이하는 ‘DNA 상태 변동’이 단일세포 수준의 동역학에 미치는 영향을 정량적으로 분석한다. 기존의 모델들은 보통 DNA 전이를 매우 빠른(adiabatic) 과정으로 가정하거나, 평균적인 전이율만을 사용해 확률적 변동을 간과하는 경우가 많았다. 그러나 실제 세포 내에서는 전사인자 결합·해리, 크로마틴 구조 변화 등으로 인해 DNA 상태 전이가 비교적 느리게 일어나며, 이는 단백질 생산·분해 과정과 시간척도가 겹쳐 비평형 현상을 초래한다.

저자들은 이러한 비정상(non‑adiabatic) 상황을 포착하기 위해 ‘동적 평균장(dynamical mean‑field) 이론’을 구축한다. 핵심 아이디어는 DNA 상태를 이진 변수(s=0,1)로 두고, 단백질 수 n(t)를 연속적인 확률 변수로 취급한 뒤, 마코프 과정으로 기술되는 전체 마스터 방정식을 두 변수의 결합 확률분포 P(s,n,t)에 대해 전개한다. 여기서 평균장은 DNA 상태 전이율 k_on, k_off와 단백질 합성·분해율 γ, β 등을 포함한다. 평균장 근사를 적용하면 고차 상관을 무시하고 1‑점 및 2‑점 함수에 대한 폐쇄식 연립 미분방정식을 얻을 수 있다.

이 식들을 풀어 두 시점 상관함수 C(t,τ)=⟨δn(t)δn(t+τ)⟩와 응답함수 R(t,τ)=δ⟨n(t+τ)⟩/δh(t) (h는 외부 자극) 를 구한다. 중요한 점은 C와 R가 전이율에 따라 전혀 다른 형태를 보인다는 것이다. 전이가 빠를수록(adiabatic limit) DNA 상태는 거의 평균화되어 유전자가 효과적으로 ‘연속적’인 생산원으로 동작한다. 이 경우 C는 단순한 지수 감쇠를 보이며, R과의 비율인 유효 온도 T_eff=C/R는 작은 상수값에 수렴한다. 반면 전이가 느릴 때는 DNA가 장시간 특정 상태에 머무르면서 ‘스위치’ 역할을 하게 되고, 단백질 수는 급격히 상승·감소하는 펄스 형태를 만든다. 이때 C는 진동적 감쇠(oscillatory decay)를 보이며, 완화 시간 τ_relax이 비정상적으로 길어지는 현상이 나타난다. 또한 T_eff는 시간에 따라 크게 변동하고, 평균적으로는 크게 증가한다. 이는 DNA 전이가 드물어질수록 시스템이 ‘비평형 잡음’에 더 민감해짐을 의미한다.

자기조절 유전자 모델에 적용한 결과는 두 가지 흥미로운 현상을 강조한다. 첫째, 전이율이 특정 구간을 지날 때 완화 시간이 급격히 증가하는 ‘critical slowing down’ 현상이 관찰된다. 이는 전이율이 전이율과 단백질 분해율 사이의 공명 조건에 근접할 때 발생한다. 둘째, 전이율이 매우 낮은 경우 상관함수가 주기적인 진동을 보이며 점차 감쇠하는데, 이는 DNA가 장시간 억제 상태에 머물다 갑자기 활성화되면서 단백질 수가 ‘펄스’처럼 급증하고 다시 감소하는 과정이 반복되기 때문이다.

이러한 결과는 실험적으로 단일세포 형광 추적이나 마이크로플루이딕스 시스템을 이용해 DNA 전이율을 조절(예: 전사인자 농도 변화, 크로마틴 변형제 투여)함으로써 검증할 수 있다. 특히 ‘유효 온도’ 개념은 비평형 시스템에서 잡음 강도를 정량화하는 새로운 지표로 활용 가능하며, 세포 내 신호전달 회로가 어떻게 외부 교란에 대해 민감도와 견고성을 조절하는지를 이해하는 데 기여한다.