크로마틴 구조에 의한 이중 상호작용으로 인한 뉴클레오솜 배열의 통계역학

초록

본 연구는 뉴클레오솜이 1차원 배열을 이루는 과정에서, 크로마틴 섬유 형성에 따른 회전 설정 보존이 1011 bp 간격의 링크어 길이 이산화를 유도한다는 점을 밝혀낸다. 유한 크기의 입자와 효과적인 두-체 상호작용을 포함한 열역학 모델을 구축하고, 대규모 뉴클레오솜 위치 데이터로부터 1체와 2체 에너지를 역추정한다. 결과는 두-체 힘이 전장적인 1011 bp 주기의 핵심 원동력이며, 전사 영역에서 관찰되는 뉴클레오솜 정렬을 설명한다.

상세 요약

이 논문은 뉴클레오솜이 DNA에 결합한 히스톤 옥타머로서, 서로 간의 거리와 회전 설정이 크로마틴 섬유의 고차 구조 형성에 결정적인 역할을 한다는 가설을 정량적으로 검증한다. 저자들은 1차원 격자 모델에 유한 크기의 입자를 배치하고, 입자 간에 거리 의존적인 두-체 상호작용 U(d)와 외부 포텐셜 V(i)를 도입한 통계역학 프레임워크를 제시한다. 핵심은 ‘회전 설정 보존’이라는 제약으로, DNA의 10.5 bp/회전이라는 구조적 특성 때문에 특정 링크어 길이(10~11 bp)의 자유에너지 장벽이 낮아진다. 이를 두-체 상호작용으로 수학적으로 표현하면, 인접 뉴클레오솜 사이의 거리 d가 10·n ± 1 bp( n∈ℕ)일 때 U(d) 가 최소가 되며, 이는 실험적으로 관측되는 주기성을 재현한다.

역추정 과정에서는 대규모 MNase‑seq 혹은 ATAC‑seq 데이터에서 얻은 뉴클레오솜 점유 확률 p(i)와 쌍점유 확률 p₂(i,j)를 이용한다. 베이즈 추정과 최대우도 방법을 결합해, V(i)와 U(d)를 동시에 최적화한다. 특히, 두-체 상호작용을 무시하면 10~11 bp 주기가 급격히 사라지고, 전사 시작점 주변의 정렬도 불규칙해지는 것을 확인했다. 이는 기존에 ‘1체 포텐셜만으로 설명 가능한’ 모델이 실제 크로마틴 환경을 충분히 포착하지 못한다는 강력한 증거다.

또한, in vitro와 in vivo 데이터 모두에서 동일한 두-체 파라미터가 적용 가능함을 보여, 세포 내외부에서 동일한 물리적 메커니즘이 작동함을 시사한다. 전사 억제 구역에서는 강한 외부 포텐셜이 존재하지만, 두-체 상호작용이 여전히 주기적인 배열을 유지시킨다. 이러한 결과는 뉴클레오솜 재배열이 전사 조절에 미치는 영향을 물리학적 관점에서 재해석할 수 있는 기반을 제공한다.