활성 에피제네틱 전환의 시공간 동기화와 도메인 형성

초록

DNA 메틸화와 탈메틸화 효소가 동시에 작용해 발생하는 활발한 전환 주기가 염색질 구조에 의해 장거리·비상호작용으로 연결될 때, 위상 고정(domain) 형태의 동기화 영역이 형성되고, 유전체 이질성에 의해 그 성장(코어싱)이 제한된다. 저자들은 최소 모델을 제시하고, 분석·시뮬레이션을 통해 동기화 정도와 도메인 크기를 예측했으며, 단일세포 시퀀싱 데이터로 정성적 검증하였다.

상세 분석

본 논문은 DNA 메틸화·탈메틸화 사이클을 ‘생화학적 오실레이터’로 모델링하고, 각 CpG 자리의 위상을 ϕ_i 로 표현한다. 오실레이터는 고유 전이율 ω_i 와 다른 자리들의 위상에 의존하는 상호작용 항 k_i(ϕ) 로 구성된다. 염색질 접힘에 의해 장거리 상호작용이 발생하며, 이는 거리 |i‑j|^‑λ 형태의 감쇠와 지수적 컷오프를 포함한다. 특히, 메틸화된 부위가 DNMT3를 더 많이 모집해 인근의 비메틸화 부위 전이를 촉진하는 비상호작용(non‑reciprocal) 특성을 f₁(ϕ)·f₂(ϕ′) 형태로 수식화하였다.

연속극한으로 전이하면 마스터 방정식은 Langevin 형태의 stochastic Kuramoto 방정식(식 3)으로 변환된다. 여기서 J는 결합 강도, λ는 접촉 확률 지수, m은 ψ₂(메틸화) 단계에 있는 CpG 비율이다. 전역 동기화 정도는 Kuramoto order parameter r 로 정의되며, 위상 밀도 ρ(θ,ω,z,t) 의 연속 방정식(식 5‑6)으로 기술된다.

분석에서는 먼저 유전체 이질성(ω_i 분포)을 무시하고 동질 오실레이터 군을 고려한다. 수치적으로 deterministic 부분을 풀어보면, 초기 무작위 위상에서 제한된 장거리 상호작용에 의해 ‘phase‑locked domain’이 형성되고, 시간이 지남에 따라 도메인이 점차 합쳐져 전체적인 위상 고정에 도달한다. 도메인 성장 속도는 J에 비단조적으로 의존한다(작은 J에서는 느리게, 중간 J≈1에서 가장 빠르게, 큰 J에서는 다시 둔화).

다음으로 강한 무질서(유전체 위치 변이) 상황을 강‑무질서 재규격화 방법으로 분석한다. 여기서는 ω_i 가 지수분포 g(ω)=μ⁻¹e^{‑μω} 로 주어지고, 상호작용에 의해 유도되는 평균 메틸화 비율 m 이 식 (8‑9) 로 자기일관적으로 결정된다. 강한 결합(J→∞)에서는 m→2/3, 약한 결합(J→0)에서는 m→1/2 로 수렴한다. 이론적 예측은 다양한 N, Ω 값을 갖는 stochastic lattice 시뮬레이션과 일치한다(그림 3a).

마지막으로 실제 단일세포 bisulfite‑sequencing 데이터( C, 5mC, 5hmC 구분)와 비교한다. 데이터에서 관찰된 메틸화 변동성 및 지역별 위상 동기화 정도는 모델이 예측한 부분 동기화와 제한된 도메인 크기와 일치한다는 점을 제시한다.

핵심 통찰은 다음과 같다. (1) 염색질 구조가 제공하는 장거리·비상호작용은 전역 동기화를 억제하고, 대신 국소적으로 위상이 고정된 도메인을 만든다. (2) 유전체 이질성(ω_i 분포)은 도메인 형성을 촉진하면서도 전체 동기화 수준을 제한한다. (3) 결합 강도 J가 중간값일 때 도메인 성장 속도가 최적이며, 이는 세포가 필요에 따라 빠르게 지역적 동기화를 조절할 수 있음을 시사한다. (4) 모델은 실험적 데이터와 정성적으로 일치해, DNA 메틸화 전환이 세포 분화 조절에 활용될 수 있는 새로운 물리적 메커니즘을 제시한다.