DNA 용매 장벽을 포함한 메조스코픽 모델의 열·기계 특성 연구

초록

본 연구는 용매와의 상호작용을 고려한 Peyrard‑Bishop‑Dauxois (PBD) 모델을 이용해 DNA의 열적·기계적 거동을 분석한다. 균일한 A‑T 서열의 용융 전이를 통해 모델 파라미터를 최적화하고, 버블 형성 및 안정성을 조사하였다. 이후 P5 프로모터 서열에 대해 분자 동역학 시뮬레이션을 수행하고, 궤적에 대한 주성분 분석(PCA)을 적용해 주요 동적 모드를 추출하였다. PCA 결과는 기존 생물학적 실험과 일치함을 확인했으며, 용매 장벽이 DNA 거동에 미치는 영향을 정량적으로 규명하였다.

상세 요약

본 논문은 기존 PBD 모델에 용매와의 상호작용을 나타내는 ‘solvation barrier’를 추가함으로써 DNA 이중나선의 열적 변성(멜팅)과 기계적 변형을 보다 현실적으로 모사한다. 먼저, 균일한 A‑T 서열을 대상으로 온도 구배 시뮬레이션을 수행해 멜팅 온도(T_m)를 측정하고, 실험값과 일치하도록 결합 상수(D), 비선형 스프링 상수(k), 그리고 용매 장벽 높이(V_solv)를 튜닝한다. 파라미터 최적화 과정에서 용매 장벽은 온도 상승에 따라 결합이 급격히 약해지는 비선형 함수를 취하며, 이는 실제 물 환경에서 수소 결합이 물 분자에 의해 차단되는 현상을 반영한다.

파라미터가 확정된 후, 시뮬레이션을 통해 열적 거품(bubble)의 형성 확률 분포와 평균 수명을 계산한다. 결과는 온도가 T_m 근처에서 거품이 국소적으로 발생하고, 용매 장벽이 존재할 경우 거품의 크기가 작아지고 수명이 짧아지는 경향을 보인다. 이는 용매가 DNA 내부에 물리적 ‘완충’ 역할을 하여 결합 파괴를 억제한다는 가설을 뒷받침한다.

다음 단계에서는 실제 전사 조절에 중요한 P5 프로모터 서열을 대상으로 장시간(수백 나노초) 분자 동역학 궤적을 수집한다. 수집된 3차원 좌표 데이터에 대해 주성분 분석(PCA)을 적용하면, 첫 번째와 두 번째 주성분이 각각 전체 서열의 ‘열적 팽창/수축 모드’와 ‘국소적 거품 개방 모드’에 해당함을 확인한다. 특히, 전사 시작점 근처(‑10, ‑35 박스)에서 두드러진 변동이 관찰되며, 이는 기존 실험에서 보고된 전사 인자 결합 부위와 높은 상관성을 가진다.

마지막으로, 용매 장벽을 제거한 순수 PBD 모델과 비교했을 때, 용매를 포함한 모델이 거품 발생 위치와 빈도에서 더 정확한 생물학적 패턴을 재현한다는 점을 강조한다. 이는 DNA의 열·기계적 거동을 예측할 때 용매 효과를 무시하면 중요한 정보를 놓칠 수 있음을 시사한다.