DNA 버블 호흡으로 안정성 파라미터 추정 시뮬레이티드 어닐링 활용

초록

본 논문은 DNA 이중 나선의 국소 개구(버블) 호흡 데이터를 이용해 염기쌍 간 수소 결합 및 인접 염기 스태킹 자유 에너지를 추정하는 방법을 제시한다. 시간 시리즈로부터 얻은 버블 크기 분포를 목표 함수로 삼아 시뮬레이티드 어닐링(SA) 최적화를 수행함으로써, 잡음이 섞인 데이터에서도 기존 실험값과 일치하는 높은 정확도의 파라미터를 복원한다.

상세 분석

이 연구는 DNA의 국소 변성 현상인 ‘버블’이 열적 요동에 의해 개폐되는 동역학을 통계 물리학적으로 모델링한 뒤, 그 모델에 내재된 자유 에너지 파라미터를 역추정하는 새로운 접근법을 제시한다. 먼저, 염기쌍별 수소 결합 자유 에너지 ε_hb(AT), ε_hb(GC)와 인접 염기쌍 간 스태킹 자유 에너지 ε_st를 정의하고, 버블 형성 시 두 개의 스태킹 접촉이 끊어지는 ‘링 팩터’ ξ를 도입한다. 버블의 크기 m와 왼쪽 경계 x_L에 대해, Boltzmann 인자를 이용한 파티션 함수 Z(x_L,m) = ξ′(1+m)^{‑c}∏{x=x_L+1}^{x_L+m}u_hb(x) ∏{x=x_L+1}^{x_L+m+1}u_st(x) 을 도출하고, 여기서 c 는 고분자 루프의 엔트로피 지수(3차원 자가회피 사슬에 대해 c≈1.76)이다. 이 식을 통해 특정 위치와 크기의 버블이 존재할 평형 확률 P_eq(x_L,m) 을 구한다.

동역학적으로는 버블 크기가 한 염기씩 열리고 닫히는 마스터 방정식 ∂P/∂t = W P 를 사용한다. 실제 시간 시리즈는 Gillespie 알고리즘을 통해 시뮬레이션했으며, 실험에서는 형광 라벨을 이용해 특정 염기 주변이 모두 개방될 때만 신호 I(t)=1을 기록하도록 가정한다. 이렇게 얻은 I(t) 로부터 버블 크기와 위치의 히스토그램을 만들고, 이를 P_eq 와 비교한다.

역추정 단계에서는 비용 함수 C = ∑_{i=1}^{M}