Polymer physics of nuclear organization and function

왜 읽었나?

Amitai & Holcman, Physics Reports 678, 1–83 (2017)

Chromatin dynamics 쪽으로는 가장 잘 정리된 리뷰 논문이라는 임박사님의 추천천

핵심 질문: 세포 핵 내부의 DNA/크로마틴 구조와 동역학을 폴리머 물리학으로 어떻게 기술할 수 있는가?

PDF: Polymer physics of nuclear organization and function.pdf


0. 논문의 목적과 범위

이 리뷰는 핵 구조와 기능을 연결짓기 위해 세 가지 축을 다룬다:

  1. 폴리머 모델 — DNA/크로마틴의 물리적 기술
  2. 수학적 분석 — 루핑, MFET 등 확률론적 공식 유도
  3. 실험 데이터 해석 — Hi-C, 단일입자추적(SPT) 데이터와의 연결

1. 핵 구조 개요 (Section 2)

다중 스케일 구조

스케일구조
2 nmDNA 이중나선
11 nm뉴클레오솜 (히스톤 8량체에 감긴 DNA)
30 nm크로마틴 파이버
300 nm루프 도메인
700 nm ~ 1400 nm응축 염색체

연구자들이 주목한 DNA 특성들

  • 이상확산(Anomalous diffusion): MSD , (서브확산). 거의 모든 크로마틴 좌위에서 관찰됨.
  • 크로마틴 영역(Chromosomal territories): 각 염색체는 핵 내부에서 뚜렷한 공간적 영역을 점유.
  • 자발적 DNA 루핑: 멀리 떨어진 두 유전자 좌위가 물리적으로 만나 유전자 발현을 조절 (Lac operon 등).
  • 텔로미어 클러스터링: 효모 핵에서 텔로미어가 2~6개의 클러스터를 형성 (Sir4, yKu 단백질 매개).
  • DSB(이중가닥 파손) 수리: 파손된 좌위가 핵 내부를 탐색하며 相동 서열을 찾음.
  • 핵 봉투와 유전자 조절: 핵 공극(nuclear pore)과의 상호작용이 유전자 발현 조절.

주요 실험 기법

기법내용
Live-cell imaging / SPTGFP-LacI-LacO 시스템으로 단일 좌위를 수십 nm, 수백 ms 해상도로 추적
3C / 4C / 5C / Hi-C화학 결찰(ligation)로 공간적으로 가까운 DNA 쌍의 접촉 빈도 측정. Hi-C는 전체 게놈 저해상도 지도
Super-resolution microscopy (STORM 등)크로마틴 응축/이완을 나노미터 해상도로 시각화
Brownian dynamics simulation폴리머 모델을 기반으로 실험 통계 재현 및 예측
Monte-Carlo simulation루프 구성 샘플링, Hi-C 데이터 역설계

2. Poymer models

독립된 노트로 작성: Polymer Models

3. 이상확산 (Anomalous Diffusion) (Section 3.10)

이상확산 공식

체제
서브확산 (subdiffusion)
정상 확산
슈퍼확산
탄도운동 (ballistic)

다양한 폴리머 모델의 이상확산 지수

모델
Rouse
Zimm (유체역학 포함)
Reptation
-polymer
실험값 (크로마틴 좌위, 효모)

연속시간 랜덤워크 (CTRW)

대기 시간 분포 (0 < a < 1)일 때 서브확산. 분수 확산 방정식:

분수 브라운 운동 (FBM)

Hurst 지수 를 가진 메모리 커널:

이상확산 지수:


4. 폴리머 루핑과 MFET (Section 4)

왜 루핑이 중요한가?

두 유전자 좌위가 물리적으로 만나는 것이 유전자 발현, 수리의 선결 조건. 루핑 속도 = MFET의 역수.

4.1 First Encounter Time 정의

두 모노머 , 간 첫 만남 시간:

4.2 Szabo-Schulten-Schulten 접근

end-to-end 거리 의 평균 MFET (소형 극한):

이 작으면 , 이 크면 으로 지배.

4.3 Fokker-Planck 고유값 전개

MFET는 FP 연산자의 첫 번째 고유값의 역수:

3차원 자유 공간 Rouse 폴리머의 :

MFET 점근 공식 (자유 공간, 3D):

MFET 점근 공식 (자유 공간, 2D):

4.4 가둔 도메인에서의 MFET (MFETC)

핵 내부 가둠을 조화 퍼텐셜로 근사. 총 에너지:

조화 퍼텐셜 강도 를 반경 에 맞춰 보정:

MFETC 공식:

중요 특성: 에서도 MFETC는 유한값으로 수렴. 루핑 시간 분포는 Poisson 분포로 잘 근사됨.

4.5 β-폴리머의 MFET

에서: (Rouse의 와 비교)


5. Hi-C 데이터 해석과 폴리머 모델 (Sections 5–6)

Hi-C란?

화학 결찰로 공간적으로 가까운 DNA 쌍을 포착하여 게놈 전체 접촉 빈도 행렬 산출. 크로마틴 조직의 “평균 사진”.

접촉 빈도의 폴리머 예측

Rouse 폴리머의 경우 접촉 빈도 (genomic 거리 의 함수):

CC 데이터 해석 모델들

5.1 Fractal Globule 모델

비평형 상태의 비얽힘(non-entangled) 구조. DNA가 자기 유사적(auto-similar) 규칙으로 구역을 채움. 현재는 크로마틴 구조를 설명하기에 불충분한 것으로 평가됨.

5.2 β-폴리머 모델

접촉 확률:

를 데이터에 피팅하여 크로마틴 내재적 장거리 상호작용의 강도 추정.

5.3 String Binder and Switch (SBS) 모델

  • 크로마틴을 자기 회피 사슬로 모델링
  • 일부 모노머( 비율)가 결합 분자(binders)와 연결 가능
  • 결합자 농도 에서 상전이(붕괴) 발생
  • CTCF 단백질과 같은 chromatin 조절 분자 표현

5.4 Co-polymer 모델

에피게놈 상태를 반영한 이종 폴리머. 에너지:

5.5 루프 돌출 (Loop Extrusion) 모델

  • Condensin/cohesin 효소가 능동적으로 크로마틴 루프를 밀어냄
  • Loop-Extruding Factor(LEF)가 염색체를 따라 이동
  • TAD(위상적 연합 도메인)와 루프 도메인 형성을 설명
  • 저농도: 불완전한 루프, 고농도: 연속적인 루프 배열

5.6 TAD 재구성 모델 (Shukron & Holcman)

Hi-C 행렬의 접촉 빈도 감쇠 지수 를 맞추도록 랜덤 connector와 고정 connector 배치.


6. 텔로미어 조직 (Section 7)

효모 핵의 텔로미어

  • 32개 텔로미어 → 2~6개 클러스터 형성 (평균 4개씩)
  • 핵 봉투 표면에서 확산·집합·해산 반복

만남 속도

텔로미어 만남 속도 (핵 표면 위의 2D Brownian particle로 coarse-graining):

해산 모델

개 텔로미어 클러스터의 해산 속도:

실험값: s → 두 텔로미어 결합의 평균 수명

공존 확률:

( s, s)


7. 단일입자추적 (SPT) 분석 (Section 8)

실험 방법

GFP-LacI/LacO 시스템 → 특정 좌위를 수십 nm 해상도로 추적. 시간 해상도 수백 ms.

핵심 통계량

자기상관 함수 (Cross-correlation):

MSD:

제약 길이 (Constraint Length):

유효 확산계수:

이상확산 지수 변화 원인들

원인효과
Rouse polymer
β-polymer ()
핵 회전/외부 드리프트 증가 (~)
진동력 증가 ()
DSB 후 크로마틴 이완 증가

외력 추출 공식

관찰 모노머 , 상호작용 모노머 사이의 유효 스프링 상수:

SPT로부터 추정:


8. DNA 수리의 물리학 (Section 9)

이중가닥 파손(DSB) 수리 과정

  1. NHEJ (비상동 말단 결합): 두 파손 말단을 직접 재연결
  2. HR (상동 재조합): 상동 서열을 탐색하여 정확히 수리 (rate-limiting step!)

물리적 특성 변화 (DSB 후)

  • 제약 길이 증가 → 더 넓은 핵 공간 탐색
  • 유효 스프링 상수 감소 → 상호작용 완화
  • 이상확산 지수 증가 → 크로마틴 이완 (에 가까워짐)

β-폴리머로 해석한 크로마틴 이완

: 고응축, : 이완. DSB 후 증가 = 크로마틴 이완 = 탐색 가속.

결론: 크로마틴 이완이 탐색 시간을 4~8배 단축.

DSB의 핵 봉투 재배치

수리 불가능한 DSB → 핵 공극(nuclear pore) 또는 Mps3(SUN domain 단백질)로 이동. 전위(translocation) 방지 보호 기전으로 해석.


9. 모델 선택 가이드 (Section 10.1)

상황권장 모델
Rouse
β-polymer
외력 + Rouse (해석 어려움)
Hi-C 해석 (전체 구조)SBS + 랜덤 루프
텔로미어 역학Coarse-grained Brownian particles on surface
DSB 탐색β-polymer + 이완 모델

10. 미해결 문제들 (Section 10.2)

  • 복수 태그 좌위의 상관관계로 역학 추출
  • Hi-C 데이터의 최적 폴리머 모델?
  • SPT와 Hi-C 통계의 통합
  • DSB 수리 단백질의 동역학 모델링
  • 전사인자의 프로모터 탐색 과정 재방문
  • NHEJ 전위 확률 계산
  • 핵 손상 후 뉴클레오솜 재배치 모델
  • 암세포 Hi-C 데이터의 특이적 특성 탐색

관련 노트