DNA 언징 역학 모델링과 실험 설계 분석

초록

본 논문은 외부 힘에 의해 DNA가 풀리는 과정을 동역학적으로 모델링한다. 포크의 서열 의존 이동, 광학·자기 트랩의 비드, 그리고 풀린 단일가닥 및 연결 고리의 고분자 거동을 모두 포함한 통합 시뮬레이션을 제시한다. 고전적인 평형 가정 모델과 비교해 비평형 효과와 장비 파라미터가 해석에 미치는 영향을 정량화한다.

상세 분석

이 연구는 DNA 언징 실험에서 흔히 간과되던 비평형 고분자 동역학을 체계적으로 포함한 모델을 구축하였다. 먼저, 포크가 이동하는 자유에너지 지형을 염기쌍 간 결합 에너지 g₀(bᵢ,bᵢ₊₁)와 ssDNA가 신장될 때의 작업 g_ss(f)를 이용해 총 자유에너지 G(n;B)=∑_{i=1}^{n}g₀을 정의하고, 이를 기반으로 포크의 전이율을 Kramers 이론 형태로 기술한다. 고분자 사슬은 Gaussian Hook, 수정된 자유 관절 사슬(FJC), 그리고 확장 가능한 웜‑라이크 체인(WLC) 모델을 각각 ssDNA와 dsDNA에 적용해 힘‑연장 관계 l(f)를 얻는다. 특히, 수정된 FJC 모델에 탄성 파라미터 γ_ss=800 pN·nm을 도입해 실험 데이터와의 적합성을 높였다.

시스템 전체는 두 개의 광학 트랩(또는 하나의 자기 트랩)과 dsDNA 핸들, 그리고 풀린 ssDNA 두 가닥으로 구성된다. 각 요소는 스프링 상수 k₁, k₂, k_m^{ss}, k_m^{ds} 로 표현되며, 평형 상태에서는 힘 균형식 k₁·x₁ = w’{ds}=w’{ss}=k₂·(X−x₄)=\bar f 로부터 전단력을 계산한다. 포크가 열리면서 n이 증가하면 ssDNA 길이가 늘어나고, 이에 따라 전체 강성 k_eff는 1/k_eff = 1/k₁ + 1/k_m^{ds}·N_ds + 2/(k_m^{ss}·(N₀^{ss}+n)) + 1/k₂ 로 감소한다. 이 감소는 포크 이동 속도와 잡음 스펙트럼에 직접적인 영향을 미친다.

시뮬레이션은 고정 힘(fixed‑force)과 고정 변위(fixed‑extension) 두 모드에서 수행되었다. 고정 힘 모드에서는 포크가 자유 에너지 최소값 사이를 홉핑하는 확률이 힘에 민감하게 변하며, 특히 임계 힘 f_c≈15.9 pN 근처에서 두 최소 사이 장벽이 12 k_BT 정도로 낮아져 반복적인 스위칭이 관찰된다. 고정 변위 모드에서는 트랩 강성이 낮을수록 포크에 전달되는 힘 변동이 감소하고, 결과적으로 전이율이 평형값에 근접한다는 점을 확인하였다. 또한, 비드와 핸들의 열 잡음이 전체 시스템의 신호‑대‑잡음 비율에 미치는 영향을 정량화했으며, n이 커질수록 ssDNA의 열 변동이 1/√n 비율로 감소함을 보여준다.

마지막으로, 기존 문헌에서 가정한 “폴리머는 즉시 평형에 도달한다”는 전제를 완화함으로써, 실험 장비의 시간 해상도(10 kHz 이하)와 공간 해상도(0.1 nm)에서 관측 가능한 비평형 효과를 명시적으로 예측한다. 이는 차세대 단일분자 시퀀싱 기술이 요구하는 정확한 힘‑시간 프로파일을 설계하는 데 필수적인 통찰을 제공한다.