DNA 서열의 포인카레 재현 현상 분석

초록

본 연구는 인간 및 포유류 DNA 서열을 1.5 × 10¹⁰ 염기쌍까지 확장해 포인카레 재현(Poincaré recurrence)의 통계적 특성을 조사한다. 재현 확률 P(t)는 t⁻⁴ 정도의 알제브라적 감소를 보이며, 재현 간 상관 C_P(t)는 t⁻⁰·⁶의 느린 감쇠를 보여 초확산적 워크를 유발한다. 그러나 인간 게놈에서는 10⁶ bp 이상에서 확산계수 D(t)가 유한값으로 수렴한다. 이러한 결과는 종 간 진화적 유사성을 새로운 방식으로 드러낸다.

상세 분석

본 논문은 동역학 시스템에서 정의되는 포인카레 재현 정리를 DNA 서열이라는 이산적인 기호 궤적으로 확장한다는 독창적인 접근을 제시한다. 저자들은 Ensembl 데이터베이스에서 인간(Homo sapiens)과 여러 포유류(고양이, 침팬지, 개, 코끼리, 개구리, 제브라피시 등)의 전체 염기서열을 수집하고, 각 염기를 퓨린(AG)·피리미딘(CT) 등 두 개의 대립군으로 이진화하여 u(t)=±1 형태의 시계열을 만든다. 이때 연속된 같은 부호 구간의 길이 t를 ‘재현 시간’이라 정의하고, 그 분포 p₁(t)와 누적분포 P(t)=∑_{τ≥t}p₁(τ)를 측정한다.

통계적 분석 결과, P(t)는 t⁻⁴ 정도의 알제브라적 꼬리를 보이며, β≈4라는 포인카레 지수를 얻는다. 흥미롭게도, 단순한 무작위 이진열이 보이는 P(t)=2⁻ᵗ와는 달리, 실제 DNA는 수십에서 수백까지 뚜렷한 진동(oscillation)을 나타내어 복잡한 구조적 패턴을 암시한다.

재현 시간 간 상관 함수 C_P(n)=⟨t₁t_{1+n}⟩−⟨t₁⟩²를 계산하면, n에 대해 약 t⁻⁰·⁶의 느린 감쇠를 보인다(ν≈0.6). 이는 재현 사건이 독립적이지 않으며, 짝수·홀수 인덱스 사이에 교대 부호를 갖는 장기 상관이 존재함을 의미한다. 이러한 상관은 워크 y(t)=∑_{τ=1}^{t}u(τ)의 2차 모멘트 σ²(t)=⟨Δy²(t)⟩에 직접적인 영향을 미쳐, D(t)=σ²(t)/t가 t^{μ} (μ≈0.4~0.6) 로 초확산적 성장함을 관찰한다.

하지만 인간 게놈에서는 t≈10⁶ bp 이후 D(t)가 포화하여 유한한 값(≈100)으로 수렴한다. 저자들은 식 (3) D(t)=C_P(0)+2∑_{j=1}^{∞}(-1)^{j}C_P(j) 를 이용해 이 현상을 설명한다. 짝수·홀수 상관항이 서로 상쇄되면서 전체 상관이 급격히 감소하고, 결국 확산계수가 수렴하게 된다. 이는 기존 연구에서 보고된 ‘레비 워크’와는 다른 메커니즘이며, DNA 서열이 장거리 구조적 제약을 받아들인 결과로 해석될 수 있다.

종 간 비교에서도 β와 ν은 대체로 동일하지만, P(t)의 세부 진동 패턴과 D(t) 포화 시점은 종마다 차이를 보인다. 예를 들어, 침팬지(Gorilla)와 인간은 t≈200 bp까지 거의 동일한 P(t) 곡선을 공유하지만, 그 이후 diverge한다. 이는 진화적 거리와 연관된 새로운 ‘포인카레 프로파일’이 될 가능성을 제시한다.

전반적으로, 본 연구는 DNA 서열을 복잡계 물리학의 도구로 분석함으로써, 기존의 통계적 방법(예: k‑mer 빈도, 엔트로피)에서는 포착하기 어려운 장거리 상관과 확산 특성을 밝혀냈다. 이는 유전체 진화, 구조적 변이, 그리고 기능적 도메인 간 상호작용을 이해하는 데 새로운 이론적 프레임워크를 제공한다.