소프트 타원 접촉 포텐셜을 이용한 DNA 거시 모델링

초록

본 논문은 기존 지질 이중층 모델에 사용된 소프트 타원 접촉 포텐셜(ECP)을 DNA의 염기쌍 상호작용에 적용한 새로운 거시(코스 그레인) 모델을 제시한다. 베딩·다이헤드럴·배제 부피·전기·용매 효과를 포함한 여섯 가지 포텐셜을 결합해 20염기쌍 두 가닥을 시뮬레이션하고, 온도 구배에 따른 하이브리다이제이션 전이와 실험적 용융곡선을 재현하였다.

상세 분석

이 연구는 DNA의 복잡한 3차원 구조와 열역학적 거동을 효율적으로 포착하기 위해, 입자 간 방향 의존성을 내재한 타원형 접촉 포텐셜(ECP)을 핵심 상호작용으로 채택한 점이 가장 큰 특징이다. 기존의 구형 기반 모델들은 방향성을 별도의 각도 의존 포텐셜으로 보정해야 했지만, ECP는 입자의 회전 행렬(α̂, β̂)과 접근 방향(r̂)을 직접 입력함으로써 염기 평면의 평면성 및 스택킹을 자연스럽게 구현한다. 특히, ϵ와 σ를 각각 방향 의존 함수로 정의하고, ν=2, μ=1 파라미터를 적용해 Gay‑Berne 형태를 변형함으로써 염기 간 결합 깊이와 최소 접근 거리를 정확히 조절한다.

베딩(V_bend)과 다이헤드럴(V_dihedral) 포텐셜은 각각 2차 조화식과 코사인 형태로 정의되어, 당‑염기‑당 골격의 굽힘과 비틀림 자유도를 제한한다. 여기서는 하드 제약(RATTLE) 알고리즘을 이용해 결합 길이를 고정함으로써 고주파 진동을 배제하고, 저주파 구조 변화를 포착하는 데 집중한다. 배제 부피는 WCA 형태로 구현했으며, 염기‑염기 상호작용을 제외한 모든 입자 쌍에 적용한다. 전기적 상호작용은 Debye‑Hückel 식을 사용해 용액 내 이온 효과를 근사했으며, 100 mM Na⁺ 조건에서 약 1 nm의 Debye 길이를 적용하였다. 용매 효과는 3SPN 모델에서 차용한 Morse형 V_solv을 통해 당‑당 간 간접적인 끌어당김을 제공, 이는 실험적 DNA 직경(≈20 Å)과 일치하도록 파라미터화되었다.

시뮬레이션 설정은 Fortran 90 기반으로, 20염기씩 두 가닥을 격자에 초기 배치하고, 온도 T* = 0.9에서 가열 후 ΔT = 0.02씩 냉각하며 10⁶ 스텝씩 평형을 잡았다. Langevin 열조절기를 사용해 translational·rotational 자유도를 동시에 잡으며, 시간 스텝은 0.005 (무차원)으로 설정했다. 평면형 타원 입자를 사용함에도 불구하고, 평면‑수직(θ = 90°)과 평면‑평면(θ = 0°) 배향에 따라 서로 다른 포텐셜 곡선을 보이며, 이는 염기 스택킹과 수소 결합을 구분하는 데 기여한다.

결과적으로 온도 T* ≈ 0.48에서 평균 염기쌍 결합률(Φ)이 급격히 감소하고, 특정 열용량(C_v)이 피크를 보이며 1차 상전이를 나타낸다. 시뮬레이션에서 얻은 용융 곡선은 실험 데이터와 sigmoid 함수로 잘 맞추어졌으며, 이는 모델이 실제 DNA의 열적 변성을 재현함을 의미한다. 다만, 파라미터 최적화 과정이 상세히 기술되지 않았고, 염기 서열 다양성·이온 농도 변화에 대한 검증이 부족한 점은 향후 보완이 필요하다. 또한, ECP의 계산 비용이 구형 Lennard‑Jones 대비 상승할 가능성이 있으므로, 대규모 시스템(수천 염기) 적용 시 효율성 평가가 요구된다.