크로마틴 내 단백질 DNA 결합을 위한 통계역학 격자 모델

초록

통계역학 격자 모델을 이용해 히스톤‑DNA 상호작용과 크로마틴 내 전사인자·다른 염색질 단백질의 경쟁 결합을 정량적으로 기술한다. 이론적 틀을 제시하고, 유전체 전반의 단백질 점유 데이터를 해석하는 기반을 제공한다.

상세 분석

본 논문은 전통적으로 시험관 내 정제된 DNA·단백질 시스템에서 사용되어 온 통계역학 격자 모델을, 최근 고통량 시퀀싱 기반 전사인자·히스톤 점유 데이터가 가능해진 상황에 맞추어 확장한다. 핵심 아이디어는 DNA를 1차원 격자로 간주하고, 각 격자점에 히스톤 옥타머(또는 변형 히스톤)와 전사인자, 기타 크로마틴 결합 단백질이 차지할 수 있는 ‘상태’를 정의하는 것이다. 모델은 결합 자유에너지(ΔG), 결합 상수(Kd), 그리고 격자점 간 상호배제(steric hindrance)를 파라미터화한다. 히스톤이 DNA에 결합하면 약 147 bp 길이의 ‘핵심’ 구간이 차단되어, 그 구간 내 다른 단백질이 결합할 수 없게 되는 ‘배제 효과’를 정확히 반영한다.

논문은 먼저 순수히 히스톤‑DNA 상호작용을 기술하는 ‘단일 격자 모델’을 제시한다. 여기서는 히스톤이 결합·해리하는 평형을 마스터 방정식 형태로 기술하고, 전체 유전체에 걸친 히스톤 밀도 프로파일을 파라미터(예: 평균 결합 자유에너지, 히스톤 간 협동/반협동)로부터 역산한다. 이어서 전사인자·다른 크로마틴 단백질이 히스톤과 경쟁적으로 결합하는 ‘다중 격자 모델’을 구축한다. 이때 전사인자는 히스톤이 차지하지 않은 빈 공간에만 결합 가능하며, 특정 DNA 서열에 대한 친화도(특이적 결합 상수)를 별도로 지정한다. 모델은 전사인자와 히스톤 사이의 ‘협동’(예: 히스톤 변형이 전사인자 친화도를 증가) 혹은 ‘반협동’(예: 히스톤이 전사인자 결합을 억제) 효과를 자유 에너지 항으로 추가할 수 있다.

수학적으로는 전이 행렬(transfer matrix) 방법을 활용해 전체 격자 시스템의 파티션 함수를 효율적으로 계산한다. 파티션 함수는 모든 가능한 결합 배치를 가중합한 것으로, 이를 통해 각 격자점에서 특정 단백질이 차지할 확률(점유율)과 전체 유전체에 걸친 평균 결합 수를 구한다. 또한, Monte Carlo 시뮬레이션을 병행해 전이 행렬 해석의 정확성을 검증한다.

실제 데이터 적용 예시로는 인간 혹은 효모 유전체에서 ChIP‑seq 기반 히스톤 변형(예: H3K27ac)과 전사인자(예: NF‑κB) 점유 데이터를 동시에 입력하고, 모델 파라미터를 최적화해 관측된 점유 패턴을 재현한다. 결과는 히스톤이 높은 밀도로 존재하는 영역에서는 전사인자 점유가 크게 억제되지만, 히스톤 변형이 특정 부위에서 발생하면 전사인자 친화도가 상승해 ‘바이파스’ 현상이 나타남을 보여준다.

이론적 기여는 다음과 같다. 첫째, 히스톤‑DNA 결합을 격자 모델에 정량적으로 포함시켜, 전사인자와의 경쟁·협동을 동시에 다룰 수 있는 통합 프레임워크를 제공한다. 둘째, 전이 행렬을 이용한 해석적 해법으로 대규모 유전체 데이터에도 적용 가능한 계산 효율성을 확보한다. 셋째, 모델 파라미터(결합 자유에너지, 협동 계수 등)를 실험적 점유 데이터와 역연산함으로써, 크로마틴 구조와 전사 조절 메커니즘을 정량적으로 연결한다. 이러한 접근은 향후 유전체 수준에서 전사인자·히스톤 상호작용을 예측하고, 변형 효소 억제제나 크로마틴 재구성 약물의 효과를 시뮬레이션하는 데 활용될 수 있다.