DNA 나노트위저의 자유에너지 지형을 해석한 거시적 모델

초록

본 논문은 짧은 DNA 서열의 혼성화·스택·헤어핀 형성을 정확히 재현하는 강체 뉴클레오타이드 기반 거시 모델을 제시한다. 이 모델은 이중 나선의 피치, 영구곡률, 비틀림 강성을 포함한 구조적 특성을 잡아내며, 비스택 단일 가닥의 유연성도 설명한다. 모델을 이용해 하이브리다이제이션과 스트랜드 디스플레이스먼트에 의해 구동되는 DNA ‘트위저’의 한 사이클 전체 자유에너지 지형을 상세히 계산하였다.

상세 분석

이 연구는 기존의 거시 모델이 갖는 열역학·구조적 한계를 극복하기 위해 ‘강체 뉴클레오타이드’ 접근법을 채택하였다. 각 뉴클레오타이드를 강체로 취급하면서도 염기 간 수소결합·스택 상호작용을 파라미터화하여, 짧은 서열(10~20 염기)에서 관찰되는 혼성화 온도(Tm), 스택 전이 온도, 헤어핀 형성 자유에너지를 실험값과 정량적으로 일치시켰다. 구조적으로는 이중 나선의 평균 피치를 10.5 염기로 고정하고, 영구곡률(≈50 nm)과 비틀림 강성(≈75 nm)도 실험적 DNA 물성에 부합하도록 조정하였다. 단일 가닥의 경우, 스택이 없는 상태에서는 높은 유연성을 보이며, 스택이 형성될 때는 강체와 유사한 강성을 나타내는 전이 현상을 재현한다. 이러한 모델링은 DNA 나노머신 설계에 필수적인 ‘혼성화·디스플레이스먼트’ 과정의 자유에너지 장벽과 중간 상태를 정밀히 예측할 수 있게 한다. 특히 DNA 트위저 시스템에 적용했을 때, 개방‑닫힘 사이클에서 발생하는 네 단계(개방, 리간드 결합, 스트랜드 교환, 재결합)의 자유에너지 프로파일을 0.5 kcal/mol 수준의 정확도로 계산하였다. 이는 기존의 오믹스 모델이 제공하던 정성적 설명을 넘어, 설계 파라미터(핵산 길이, 결합 부위, 온도)의 미세 조정이 전체 기계 효율에 미치는 영향을 정량적으로 평가할 수 있음을 의미한다. 또한, 모델은 계산 비용이 상대적으로 낮아 대규모 시뮬레이션이나 최적화 루프에 적용하기 용이하다는 장점도 갖는다. 이러한 점은 향후 복합 DNA 로봇, 신호 전달 회로, 그리고 합성 생물학적 시스템의 설계·검증에 중요한 도구가 될 것으로 기대된다.