RNA 상호작용 구조의 조합론적 분석

초록

본 논문은 두 비암호화 RNA가 형성할 수 있는 결합 구조를 수학적으로 모델링하고, 그 구조들의 개수를 생성함수와 특이점 분석을 통해 정확히 추정한다. 인트라분자 결합은 가짜결절(pseudoknot)이 없고, 인터분자 결합은 교차하지 않으며, ‘지그재그’ 형태를 배제한 ‘공동 구조(joint structure)’를 정의한다. 이를 바탕으로 명시적 재귀식과 비대칭적인 비율을 도출하고, 대수적 근사식을 통해 구조 수의 지수적 성장률을 제시한다.

상세 분석

논문은 먼저 RNA‑RNA 상호작용을 ‘다이어그램’ 형태로 시각화하고, 각 파트너의 내부 결합을 가짜결절이 없는 이중 비교차(non‑crossing) 구조, 즉 이차원 평면에 겹치지 않게 배치된 완전 매칭으로 모델링한다. 인터분자 결합 역시 교차를 허용하지 않으며, 두 파트너 사이에 발생할 수 있는 ‘지그재그’ 구성—즉, 한 파트너의 결합이 다른 파트너의 결합 사이에 끼어 들어 복잡한 꼬리를 형성하는 경우—을 엄격히 배제한다. 이러한 제약은 구조 공간을 크게 제한하지만, 동시에 정확한 조합론적 계산을 가능하게 만든다.

저자들은 먼저 단일 RNA의 2차원 비교차 구조에 대한 기존 결과(카탈란 수와 그 변형)를 확장하여, 두 사슬이 상호작용하는 경우의 생성함수 G(x,y)를 정의한다. 여기서 x와 y는 각각 첫 번째와 두 번째 RNA의 길이를 나타내는 변수이며, G는 두 사슬의 길이쌍 (n,m)에 대한 구조 수 a_{n,m}의 이중 거듭제곱급수이다. 제약 조건을 반영하기 위해 G는 기본적인 비교차 구조의 생성함수 C(z)와 결합 매트릭스 M(x,y)를 곱하는 형태로 전개된다. 특히, M은 인터분자 결합이 교차하지 않도록 하는 ‘Dyck‑path’와 유사한 구조를 갖으며, 이는 Catalan‑type 항을 포함한다.

특이점 분석 단계에서는 G의 가장 가까운 특이점(반경 r)과 그 특이점의 종류(알게브라적, 로그형 등)를 파악한다. 저자들은 C(z)의 특이점이 z=1/4임을 이용해, M(x,y)의 특이점이 x·y=1/4와 연관됨을 보인다. 복합 특이점 이론을 적용해, a_{n,m}는 대수적 형태 a·n^{-3/2}·b^{n}·c^{m}·(1+o(1)) 로 근사될 수 있음을 증명한다. 여기서 b와 c는 각각 파트너 사슬의 성장 상수이며, a는 특이점의 잔여값에 의해 결정된다.

또한, 논문은 a_{n,m}에 대한 명시적 재귀식을 제시한다. 재귀식은 작은 길이의 구조 수를 기반으로, 새로운 결합을 추가하거나 기존 결합을 확장하는 두 가지 경우로 나뉜다. 구체적으로, a_{n,m}=∑{i≤n,j≤m} a{i,j}·a_{n−i,m−j}·C_{i+j} 형태이며, 여기서 C_k는 k‑길이 단일 사슬의 비교차 구조 수이다. 이 재귀식은 동적 프로그래밍 구현을 통해 효율적인 계산이 가능함을 보여준다.

결과적으로, 저자들은 공동 구조의 수가 단일 사슬 구조보다 빠르게 성장하지만, 여전히 다항식 차수와 지수적 성장 상수가 명확히 규정된다는 점을 강조한다. 이러한 정량적 이해는 MFE(최소 자유 에너지) 기반 예측 알고리즘의 설계와 복잡도 분석에 직접적인 영향을 미친다.