DNA 결함 부위의 염기쌍 호흡에서 자기조직화 임계현상

초록

본 연구는 12염기쌍 길이의 DNA에 중앙 결함을 도입하고, 전구체 전자구조와 확률 미분 방정식(SDE) 모델을 이용해 염기쌍 호흡을 시뮬레이션하였다. 푸리에 변환 결과, 거리 변동이 1/f 형태의 핑크 노이즈를 보이며, 스펙트럼 기울기가 -1에 근접하는 전력법칙을 나타냈다. 이는 염기쌍 호흡이 자기조직화 임계현상(self‑organized criticality, SOC)이며, DNA 전이(멜팅)의 핵생 단계가 2차 상전이와 유사함을 시사한다.

상세 분석

이 논문은 DNA 이중 나선의 미세한 구조 변동, 즉 ‘염기쌍 호흡(breathing)’ 현상이 복잡계 물리학에서 다루는 자기조직화 임계현상(SOC)의 전형적인 예임을 실증적으로 보여준다. 연구자는 12염기쌍으로 구성된 짧은 DNA 서열에 중앙에 결함(예: T·A 대신 G·C 결함)을 삽입하고, 두 가지 독립적인 접근법을 적용하였다. 첫 번째는 AMBER 패키지를 이용한 전통적인 전자구조 기반 전원분자역학(MD) 시뮬레이션으로, 원자 수준에서 힘장 파라미터와 용매 효과를 모두 고려한다. 두 번째는 각 염기쌍 사이의 거리 변수를 연속적인 확률 미분 방정식(SDE) 형태로 모델링한 것으로, 마찰항과 백색 잡음(화이트 노이즈) 항을 포함한다. 두 모델 모두 시간에 따라 기록된 염기쌍 간 거리 데이터를 푸리에 변환(Fourier Transform)하여 스펙트럼을 얻었다. 흥미롭게도, 두 경우 모두 저주파 영역에서 파워 스펙트럼이 $S(f)\propto f^{-\alpha}$ 형태를 보였으며, $\alpha$ 값이 0.9~1.1 사이로 -1에 매우 근접했다. 이는 전통적인 백색 잡음($\alpha=0$)이나 갈색 잡음($\alpha=2$)과는 구별되는 핑크 노이즈($\alpha\approx1$) 특성을 의미한다. 특히, SDE 모델에서 입력 잡음은 순수한 백색 잡음이었음에도 불구하고, 시스템 자체가 비선형 복원력과 마찰을 통해 잡음 스펙트럼을 변형시켜 핑크 노이즈를 생성한다는 점은 비선형 동역학이 임계 상태를 자율적으로 유도한다는 강력한 증거다. 또한, 결함 부위 주변의 염기쌍이 가장 큰 스펙트럼 진폭을 보이며, 이는 결함이 ‘핵생(nucleation)’ 역할을 수행해 전체 이중 나선이 열리기 전 단계에서 국소적인 불안정성을 증폭시킨다는 물리적 해석을 가능하게 한다. 이러한 결과는 DNA 멜팅이 연속적인 2차 상전이와 유사하게, 작은 국소 변동이 임계점에서 전역적인 구조 변화를 촉발한다는 기존 가설을 정량적으로 뒷받침한다.