DNA 삼중위치 모델로 본 장·단 온도와 힘에 의한 변성 연구

초록

본 논문은 Knotts·Rathore·Schwartz·de Pablo가 제안한 3‑site coarse‑grain DNA 모델을 약간의 베이스 페어링 파라미터 조정으로 장쇄 DNA의 열·기계 변성을 정량적으로 재현하도록 개선하였다. 수정된 모델은 동질 서열의 용융 온도를 실험값과 통계 모델값에 맞추고, 이질 서열에서는 단계적 변성을, 상온에서의 인장 파괴 힘은 약 10 pN 수준을 제공한다. 또한 긴 동질 서열의 열 변성이 1차 상전이임을 확인하고, 파괴 힘의 임계 지수 σ≈0.70을 얻었다. 모델은 DNA의 이중 나선 구조와 메이저·마이너 그루브를 현실적으로 재현하면서, 베이스 페어를 끊는 에너지와 힘이 기대 범위에 있음을 보여준다.

상세 분석

본 연구는 기존 KRSdP 모델이 갖는 3 site‑per‑nucleotide 구조를 유지하면서, 베이스 페어링 에너지 항을 실제 DNA의 열·기계 거동에 맞게 재조정한 점이 핵심이다. 원 논문에서는 모든 상보적 염기쌍에 대해 bp V 항을 적용했으나, 이는 하나의 염기가 두 개 이상의 상보 염기와 결합하는 비물리적 상황을 초래한다. 저자들은 이를 방지하기 위해 bp V를 동일한 i 쌍에만 적용하도록 제한했고, 그 결과 두 가닥 간 결합이 다소 약해져 용융 온도가 낮아지는 문제가 드러났다. 이를 보완하기 위해 AT와 GC의 결합 에너지를 각각 3.90 kcal·mol⁻¹, 4.37 kcal·mol⁻¹ 로 상승시켰으며, 이렇게 조정된 파라미터는 480 bp 길이의 동질 서열(…AAAAA… 및 …GGGGG…)에 대해 각각 335 K와 374 K의 용융 온도를 재현한다. 이는 50 mM Na⁺ 조건에서 통계적 폴리머 모델이 예측한 값(336.3 K, 373.1 K)과 거의 일치한다.

열 변성 실험에서는 절반이 열린 초기 상태를 설정하고 온도에 따라 개방·재결합 속도를 측정함으로써 임계 온도를 정의하였다. 온도와 속도 사이의 선형 관계는 전이의 1차 성격을 시사한다. 또한, Langevin 동역학 시뮬레이션에서 사용된 마찰 계수 γ=5 ns⁻¹가 실제 DNA의 확산 속도와 차이가 있음을 인식하고, 3D 확산 계수를 통해 실험적 γ≈5×10¹¹ s⁻¹(≈500 ns⁻¹) 정도가 필요함을 제시하였다.

기계 변성에서는 일정 온도(298 K)에서 외부 힘을 가해 DNA를 인장했으며, 변성 전이점에서의 임계 힘이 약 10 pN 수준으로 측정되었다. 이는 실험적 단일분자 힘분광법 결과와 일치한다. 또한, 임계 힘과 온도 사이의 관계를 분석한 결과, F_c∝(T_c−T)^σ 형태의 스케일링을 보이며 σ≈0.70이라는 비정수 지수를 얻었다. 이는 기존 1차 전이 모델에서 기대되는 σ=0.5와 차이가 있어, 모델이 포함하는 비선형 탄성 및 스택 상호작용이 전이 거동에 영향을 미침을 시사한다.

전반적으로, 모델은 DNA의 나선 구조와 주요·보조 그루브를 정확히 재현하면서, 열·기계 변성에 필요한 에너지와 힘을 실험적 범위 내에 두었다. 이는 DNA‑단백질 상호작용, 특히 특정 결합 부위에서의 국소 개방 현상을 고해상도 시뮬레이션으로 탐구할 수 있는 강력한 도구가 된다.