코어그레인 DNA 모델의 구조·역학·열역학 특성 분석

초록

본 논문은 Ouldridge et al. (2010)에서 제시된 코어그레인 DNA 모델을 기반으로, 사슬 연결, 배제 부피, 염기 스태킹, 수소 결합 등 효과적인 상호작용을 통해 단일 가닥 스태킹, 이중 나선 형성, 헤어핀 구조 등 다양한 DNA 현상을 재현한다. 실험 데이터와 비교해 단일 가닥 영구 길이, 이중 나선의 굽힘·비틀림·신장 강성, 그리고 말단, 버블, 내부 루프와 같은 복합 구조에 대한 열역학적 거동을 정량화한다. 모델이 실제 DNA의 물리적 특성을 얼마나 정확히 포착하는지 평가하고, 스태킹과 프레이링이 이중 나선 형성 전이에서 차지하는 역할을 탐구한다.

상세 분석

이 연구는 기존 Ouldridge et al. (2010) 모델을 확장·보완하여, DNA의 구조·역학·열역학을 통합적으로 평가한다. 모델은 네 가지 기본 상호작용을 사용한다. 첫째, 체인 연결을 담당하는 펜틀리톤-볼츠만 형태의 결합 포텐셜은 골격의 연속성을 보장한다. 둘째, 배제 부피는 입자 간의 겹침을 방지하며, Lennard‑Jones‑type 잠재력을 통해 적절한 거리에서 강한 반발을 구현한다. 셋째, 염기 스태킹은 인접 염기 사이의 방향성 있는 포텐셜로, 온도 의존적인 스태킹 에너지를 부여해 단일 가닥의 강직성을 재현한다. 넷째, 수소 결합은 상보적 염기쌍 사이에 비방향성·거리 의존적인 포텐셜을 적용해 이중 나선의 형성을 촉진한다.

열역학적 검증에서는 단일 가닥 스태킹 전이 온도(Tm)를 시뮬레이션하고, 실험값과의 차이를 2–3 K 이내로 최소화했다. 이는 스태킹 포텐셜의 온도 스케일링 파라미터가 적절히 조정된 결과이며, 단일 가닥의 영구 길이(Lp)가 약 1.5 nm(실험값 1.4–1.6 nm)으로 재현되었다. 이중 나선의 굽힘 강성은 영구 길이 50 nm 정도로, 실험적 DNA(≈50 nm)와 일치한다. 비틀림 강성은 토크-각도 관계를 통해 추정했으며, 약 2 × 10⁻¹⁹ J·rad⁻¹로, 실제 DNA의 비틀림 모듈러스와 비슷한 수준이다. 신장 강성은 1000 pN·nm⁻¹ 정도로, 실험적 단일 사슬 DNA의 1000–1200 pN·nm⁻¹ 범위와 부합한다.

복합 구조에 대한 테스트에서는 말단(덩글링) 염기, 버블(돌출 염기), 내부 루프를 포함한 시퀀스를 설계하고, 각각의 자유 에너지 변화를 계산했다. 결과는 말단이 존재할 경우 결합 자유 에너지가 약 0.5 kcal·mol⁻¹ 감소하고, 버블이 클수록 전이 온도가 낮아지는 경향을 보였다. 특히, 스태킹이 약화된 부위에서 프레이링 현상이 두드러져, 이중 나선의 전이 폭이 넓어지고 전이 온도가 감소한다는 점을 확인했다. 이러한 현상은 실험적 DNA 열역학 데이터와 정량적으로 일치한다.

전반적으로 모델은 복잡한 DNA 구조와 동역학을 저렴한 계산 비용으로 재현하면서도, 실험과의 차이를 수 kBT 수준으로 억제한다는 장점을 가진다. 다만, 고농도 이온 환경이나 특정 금속 이온 결합을 포함한 경우에는 포텐셜 파라미터의 추가 보정이 필요할 것으로 보인다.