초고밀도 초나선 DNA에서의 버블 통계와 위치 특성

초록

본 논문은 초나선 응력이 가해진 이중 가닥 DNA의 열역학적 변성을 Zimm‑Bragg 모델과 Benham 모델을 이용해 전산적으로 분석한다. 전이 행렬과 자체 일관적 선형화 기법을 통해 전체 게놈 규모(10⁶–10⁹ bp)까지 버블의 개방 확률, 크기 분포, 위치 특성을 효율적으로 계산한다. 결과는 무제한 DNA에 비해 초나선 DNA에서 평균 버블 크기가 10²–10³ bp로 크게 증가하고, 실제 유전체에서는 전사 시작점 상류에 큰 버블이 집중되는 경향을 보인다.

상세 분석

이 연구는 두 가지 핵심 모델을 결합한다. 첫 번째는 전통적인 Zimm‑Bragg(ZB) 모델로, 이는 Poland‑Scheraga 모델의 근사 형태이며, 개별 염기쌍(step)마다 닫힌 상태와 개방 상태에 대한 자유에너지를 θ_i 변수(0: 개방, 1: 닫힘)로 표현한다. 전이 행렬 T_i와 T’i를 정의하고, 전체 체계의 분배함수 Z를 트레이스 형태로 전개함으로써 위치별 개방 확률 ⟨θ_i⟩와 버블 크기‑위치 결합 확률 P{i,l}을 O(N) 시간에 정확히 계산한다. 특히, 루프 엔트로피 항 Δg_loop≈20 k_BT를 포함해 작은 버블에 대한 억제 효과를 정량화한다.

두 번째는 Benham 모델로, 초나선 밀도 σ(=α/N, α는 링크 차이)를 외부 제약으로 도입한다. 여기서는 개방된 염기쌍 수 n_o가 초나선 에너지 H_TW(n_o)와 결합해 비선형 효과를 만든다. H_TW는 전단 강성 C≈3.09 k_BT·rad⁻¹·bp⁻¹와 전반적인 초나선 탄성 K≈2220 k_BT/N을 이용해 유도되며, n_o에 대한 로그 항과 1/(βC) 항을 포함한다. Benham 모델은 원래 O(N²) 복잡도를 갖지만, 저자들은 n_o의 변동을 평균값 \bar{n}o 주변에서 1차 전개하고, 효과적인 외부 장 h=∂H_TW/∂n_o|{\bar{n}_o}를 도입한 자체 일관적 선형화(self‑consistent linearization)를 제안한다. 이때 H_eff≈H_ZB+H_TW(\bar{n}_o)+(N-\bar{n}_o)h - h∑θ_i 형태가 되며, h는 자기 일관 방정식 h = (∂H_TW/∂n_o)(h n_o(h))을 통해 수치적으로 해결한다.

동질 서열에 대해서는 작은 σ와 저온 한계에서 해석적 근사식(식 20, 22)을 도출하고, 이론적 예측이 수치 해와 거의 일치함을 확인한다. 이질 서열에 대해서는 뉴턴‑랩슨과 이분법을 결합해 h를 10⁻⁴ 정확도로 10–20회 반복해 구한다. 이렇게 얻은 h를 ZB 전이 행렬에 삽입하면, 전체 유전체에 대해 O(N) 시간에 개별 염기쌍의 개방 확률과 버블 분포를 얻을 수 있다.

실험적으로는 E. coli 전체 게놈(≈4.6 Mbp)을 대상으로 계산했으며, 일반적인 무작위 서열(동일 GC 비율)과 비교했다. 초나선(σ≈−0.06) 조건에서 평균 개방 확률 1−Θ가 무제한 DNA에 비해 10배 이상 증가하고, 평균 버블 길이 L이 10²–10³ bp 수준으로 크게 늘어났다. 특히, 큰 버블(>500 bp)의 발생 확률이 무작위 서열 대비 5~10배 높았다. 이러한 큰 버블은 전사 시작점(TSS) 상류 0–200 bp 구간에 집중되는 경향을 보였으며, 이는 초나선에 의해 전사 개시 부위가 물리적으로 더 쉽게 열리게 된다는 기존 생물학적 가설을 정량적으로 뒷받침한다.

또한, 온도 의존성을 조사한 결과, 초나선이 없는 경우(σ=0)에는 37 °C에서 거의 개방이 없지만, σ=−0.06에서는 동일 온도에서도 1−Θ≈0.02 수준의 지속적인 개방이 존재한다. 이는 초나선이 DNA의 열적 안정성을 크게 낮추어, 생리적 온도에서도 국소적인 변성을 유도한다는 의미다.

전반적으로, 본 논문은 초나선 DNA의 열역학적 변성을 효율적인 전이 행렬 기반 방법과 자체 일관 선형화 기법으로 정량화함으로써, 대규모 게놈 수준에서 버블 통계와 위치 특성을 정확히 예측한다는 점에서 이론 물리·생물학 분야에 중요한 도구를 제공한다.