RNA 구조와 의사결절의 조합론

초록

본 논문은 의사결절(pseudoknot)을 포함한 RNA 구조의 조합적 특성을 연구한다. k-비교단(k‑noncrossing) RNA 의사결절 구조를 최대 상호 교차 호의 집합으로 분류하고, 원형 RNA에 대한 열거식도 제시한다. 특히 3‑비교단 구조와 전통적인 2‑비교단(이차 구조)에서 각각 4항 및 2항 재귀식을 도출하였다. 또한 2‑아크(길이 2인 호)를 제외한 제한된 k‑비교단 구조의 일반적 열거식을 제공한다.

상세 분석

논문은 RNA 2차 구조를 그래프 이론의 호(arc) 모델로 표현하고, 의사결절을 허용하는 경우 호들의 교차 패턴이 복잡해지는 점에 주목한다. 기존 연구는 주로 비교단(k‑noncrossing) 조건 하에서 호들의 최대 교차 수를 제한했지만, 저자들은 “최대 상호 교차 호 집합(maximal mutually intersecting arcs)”이라는 새로운 분류 체계를 도입한다. 이는 동일한 교차 그룹에 속하는 호들을 하나의 블록으로 묶어, 구조 전체를 블록들의 순열과 내부 배치로 해석할 수 있게 한다. 이러한 접근은 생성함수(generating function)를 구성할 때 블록별 독립성을 활용해 복합적인 교차 패턴을 조합론적으로 분해하는 데 유리하다.

핵심 수학적 도구는 매개변수화된 다항식과 비선형 재귀식이다. 3‑비교단 구조에 대해서는 4항 재귀식이 도출되는데, 이는 이전에 알려진 3‑항 혹은 2‑항 재귀식보다 더 높은 차수의 상호작용을 포착한다. 구체적으로, n‑길이 RNA에 대한 구조 수 a_n은 a_{n‑1}, a_{n‑2}, a_{n‑3}, a_{n‑4}의 선형 결합으로 표현되며, 각 항의 계수는 Catalan 수와 그 변형인 Motzkin 수, 그리고 k‑비교단 제한에 따른 조정 인자를 포함한다. 이 재귀식은 동적 프로그래밍 구현 시 O(n) 시간 복잡도로 모든 구조를 열거할 수 있게 해준다.

또한, 원형 RNA(폐쇄형 서열)의 경우 호의 시작점과 끝점이 순환적으로 연결되므로, 선형 경우와는 다른 경계 조건이 필요하다. 저자는 원형 구조를 선형 구조에 “wrap‑around” 제약을 추가하는 방식으로 모델링하고, 이에 대한 생성함수를 기존 선형 생성함수에 적절한 대칭 항을 더함으로써 얻는다. 결과적으로 원형 3‑비교단 구조의 열거식은 선형 경우와 동일한 형태의 재귀식을 가지지만, 초기 조건이 달라져 전체 카운트가 약간 증가한다.

제한된 k‑비교단 구조, 즉 2‑아크를 배제한 경우는 생물학적 의미가 있다. 2‑아크는 최소 루프 길이가 3인 경우에 해당하며, 실제 RNA에서 흔히 관찰되지 않는다. 이를 배제함으로써 생성함수에 (1‑x^2)와 같은 차단 인자를 곱해 주면, 기존 k‑비교단 생성함수에서 불필요한 항을 제거할 수 있다. 저자는 이 과정을 일반 k에 대해 귀납적으로 증명하고, 최종적으로는 모든 k에 대해 통합된 폐쇄형 식을 제시한다. 이 식은 k가 증가함에 따라 급격히 복잡해지는 구조 공간을 효율적으로 추정할 수 있게 해준다.

전체적으로 논문은 조합론, 대수적 생성함수, 그리고 재귀적 구조 카운팅을 결합해 RNA 의사결절 구조의 정확한 열거식을 제공한다. 이는 기존의 확률적 모델이나 시뮬레이션 기반 접근법에 비해 이론적 근거가 명확하고, 알고리즘 구현 시 메모리와 시간 효율성을 크게 향상시킬 수 있다. 또한, 제시된 4항 및 2항 재귀식은 향후 RNA 설계와 구조 예측 알고리즘에 직접 적용 가능하며, 생물학적 데이터와의 정량적 비교에도 활용될 여지가 크다.