고차원 크로마틴 구조가 제한하는 뉴클레오솜 통계역학

초록

본 연구는 뉴클레오솜이 30 nm 크로마틴 섬유와 그 상위 구조에 의해 어떻게 배열되고 에너지적으로 조절되는지를 설명하기 위해, 유한 크기의 입자와 근접 상호작용을 포함한 열역학 모델을 구축하였다. 고처리량 뉴클레오솜 위치 데이터로부터 1‑입자 밀도와 2‑입자 상호작용을 동시에 추정할 수 있음을 보였으며, 단순히 스테레오적 배제만을 고려한 최소 모델이 전사 영역의 실제 점유 패턴을 재현함을 확인하였다.

상세 분석

이 논문은 뉴클레오솜 배열을 물리적 입자 시스템으로 모델링함으로써, 기존에 주로 서열 기반으로 해석되던 핵심 현상을 새로운 관점에서 재조명한다. 저자들은 뉴클레오솜을 길이 147 bp인 ‘유한 입자’로 보고, 이 입자들 사이에 효과적인 근접 이웃 상호작용(두 입자 사이의 거리 의존적 포텐셜)을 도입하였다. 이때 ‘임의의 DNA 결합 에너지’는 각 위치마다 다른 서열에 의해 결정되는 1‑body 포텐셜로, 전통적인 ‘서열 선호도’와 동일시될 수 있다. 핵심 수식은 통계역학적 퍼스티스(Grand Canonical) 집합을 이용해 입자 수와 위치에 대한 확률분포를 기술하고, 이를 통해 관측된 1‑particle density ρ(i)와 2‑particle correlation g(i,j)를 역으로 추정한다.

특히 저자들은 ρ(i)만을 이용해 1‑body 포텐셜을 추정하는 기존 방법이 2‑body 상호작용을 무시함으로써 실제 에너지 지형을 왜곡시킬 수 있음을 실험적으로 입증하였다. 이를 보완하기 위해, ρ(i)와 g(i,j) 데이터를 동시에 최적화하는 역문제(inverse problem)를 풀어, 두 종류의 포텐셜을 정확히 분리한다. 수치 시뮬레이션에서는 ‘두 몸 상호작용을 무시한 모델’이 서열에 의한 위치 예측에서는 어느 정도 성능을 보이지만, 전사 억제 구역이나 고밀도 영역에서는 큰 오차를 나타낸다. 반면, 전체 모델을 적용하면 실제 in vitro 및 in vivo 데이터와 높은 상관관계를 보이며, 특히 전사체 전반에 걸친 뉴클레오솜 점유율의 ‘파동형’ 패턴을 재현한다.

또 다른 중요한 결과는 ‘최소 모델’이다. 여기서는 1‑body 포텐셜을 거의 무시하고, 오직 스테레오적 배제와 거리 의존적 2‑body 포텐셜만을 사용한다. 이 모델은 전사 개시점 주변에서 뉴클레오솜이 규칙적으로 배열되는 현상을 충분히 설명한다. 즉, 뉴클레오솜 위치는 반드시 서열에 의해 결정되는 것이 아니라, 고차원 크로마틴 구조가 유도하는 물리적 제약—특히 인접 뉴클레오솜 간의 회전 및 굽힘 자유도 제한—에 의해 크게 좌우된다는 점을 강조한다.

이러한 접근법은 기존의 ‘서열 중심’ 모델과 ‘물리적 제약 중심’ 모델 사이의 간극을 메우며, 크로마틴 구조와 뉴클레오솜 배열 사이의 상호작용을 정량적으로 분석할 수 있는 새로운 프레임워크를 제공한다. 향후에는 히스톤 변형, ATP 의존성 리모델링 복합체, 그리고 전사인자 결합 등 동적 요소를 포함시켜, 실제 세포 내 복합적인 조절 메커니즘을 더 정밀히 모델링할 수 있을 것으로 기대된다.