DNA 풀림과 재결합에서의 히스테리시스와 비평형 작업 정리

초록

본 연구는 고정 힘 앙상블에서 양쪽으로 당겨지는 이중 가닥 DNA(dsDNA)의 풀림·재결합 과정에서 발생하는 히스테리시스를 Monte Carlo 시뮬레이션으로 조사한다. 힘을 일정 속도로 증가시켜 위상 경계 위까지 올린 뒤 다시 감소시키는 한 사이클 동안 수행된 작업의 확률분포를 분석하고, 평균 작업이 히스테리시스 면적과 10% 이내의 차이를 보임을 확인한다. 또한 비평형 측정 데이터를 워크 정리를 이용해 반복 적용함으로써 평형 힘‑거리 등온선을 복원하고, 자발적 재결합 과정에서도 동일한 평형·비평형 등온선을 재현한다는 점을 제시한다.

상세 분석

이 논문은 DNA 이중 나선의 역학적 변형을 통계 물리학적 관점에서 접근한다. 고정 힘(정압) 앙상블을 선택함으로써 외부 힘이 직접적으로 시스템에 가해지는 상황을 모델링하고, 힘을 시간에 따라 선형적으로 변화시키는 프로토콜을 설계하였다. 초기 힘 g₀는 풀림이 일어나지 않는 안정 영역에 위치하고, 최대 힘 gₘ은 이론적 위상 경계(풀림‑재결합 전이)를 초과하도록 설정하였다. 힘을 증가시키는 구간과 감소시키는 구간을 연속적으로 수행함으로써 하나의 완전한 사이클을 만든다. 이때 Monte Carlo 시뮬레이션은 각 단계에서 DNA의 미시 상태를 샘플링하고, 작업 W = ∫ g dx(구간) 를 계산한다.

핵심 결과는 작업 분포 P(W) 가 비대칭적인 형태를 보이며, 평균 작업 ⟨W⟩ 가 히스테리시스 루프 면적 A와 거의 일치한다는 점이다. 특히, 당김 속도가 느릴수록 ⟨W⟩‑A 차이는 10% 이하로 수렴하고, 매우 빠른 속도에서는 비평형 효과가 커져 차이가 확대된다. 이는 비평형 열역학에서 알려진 Jarzynski 등식 ⟨e^{−βW}⟩ = e^{−βΔF} 와 연관된다. 저자들은 반복적인 비평형 측정에서 얻은 W 값을 이용해 ⟨e^{−βW}⟩ 를 계산하고, 이를 역으로 자유에너지 차이 ΔF(g) 를 추정한다. 이렇게 얻어진 ΔF(g) 는 직접 평형 Monte Carlo 시뮬레이션으로 구한 힘‑거리 등온선과 거의 일치한다.

또한, 재결합 구간에서 “자발적” 재결합 현상을 관찰하였다. 힘을 감소시킬 때, 시스템은 메타안정 상태를 거쳐 최종적으로 원래의 결합 상태로 돌아가며, 이 과정에서도 비평형 작업 정리를 적용하면 동일한 평형 등온선을 복원할 수 있음을 보였다. 이는 DNA 복제·전사 과정에서 발생하는 실제 물리적 상황을 모델링하는 데 중요한 의미를 가진다.

기술적인 측면에서, 시뮬레이션은 2차원 격자 모델을 기반으로 하며, 각 뉴클레오티드 쌍은 결합 에너지 ε와 외부 힘 g에 의해 결정되는 전이 확률을 가진다. Metropolis 알고리즘을 사용해 온도 T를 일정하게 유지하면서 시스템을 진화시켰다. 당김 속도는 시간당 힘 변화 Δg/Δt 로 정의되었으며, 여러 속도(10⁻⁴ ~ 10⁻¹ ε/k_B) 에 대해 실험을 수행했다. 결과는 속도 의존적인 히스테리시스 루프 형태와 작업 분포의 폭 변화를 명확히 보여준다.

이 연구는 DNA 역학 실험에서 흔히 사용되는 광학 트랩·마이크로피펫 등에서 측정된 비평형 데이터에 워크 정리를 적용해 평형 특성을 추정할 수 있는 실용적인 방법론을 제시한다. 특히, 실험적으로 접근하기 어려운 평형 상태(예: 높은 힘에서의 풀림 상태)를 비평형 경로를 통해 간접적으로 파악할 수 있다는 점이 과학적·기술적 가치를 높인다.