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

초록

본 논문은 두 RNA 사슬이 상호작용하는 구조를 “joint structure”라 정의하고, 최소 결합 길이 4와 내부·외부 스택 길이 최소 2(고립된 아크 없음)를 만족하는 경우를 대상으로 한다. 새로운 “shape” 개념을 도입해 기호적 열열을 구축하고, 이를 통해 구조의 정확한 계수와 지수 성장률을 추정한다. 결과적으로 기존 알고리즘보다 훨씬 작은 성장률을 보이는 구조군을 밝혀, RNA‑RNA 상호작용 예측 알고리즘 설계에 이론적 기반을 제공한다.

상세 분석

이 연구는 RNA‑RNA 상호작용을 그래프 이론적 관점에서 모델링함으로써, 기존 MFE(최소 자유 에너지) 기반 예측 방법이 놓치기 쉬운 구조적 제약을 명시적으로 수학화한다. 먼저 두 백본을 평행하게 배치하고, 각 백본 내부와 백본 간 결합을 비교단(non‑crossing) 아크로 표현한다. 여기서 “zig‑zag” 금지는 두 백본 사이에 교차하는 아크가 서로 뒤섞이는 형태를 배제함으로써, 실제 생물학적 상호작용에서 관찰되는 단순한 결합 패턴을 반영한다. 논문은 특히 최소 아크 길이 4와 스택 길이 최소 2라는 두 가지 제약을 두어, 고립된 결합(단일 아크)이나 매우 짧은 루프가 포함된 구조를 제외한다. 이는 실제 RNA‑RNA 상호작용에서 흔히 발견되는 최소 루프 크기와 스택 안정성을 반영한다는 점에서 실용적이다.

핵심 기법은 “shape”이라는 추상화이다. 기존의 RNA 구조 열열에서는 단순히 아크의 존재 여부만을 고려했지만, 여기서는 각 스택을 최소 길이 2인 블록으로 묶어, 스택 내부의 세부 아크 배치를 무시하고 전체적인 스택 배열만을 기록한다. 이렇게 하면 복잡한 내부 구조를 무시하고도 전체 구조의 조합적 특성을 정확히 파악할 수 있다. 형태별(Shape) 생성함수를 정의하고, 이를 복합 구조에 대한 복합 생성함수와 연결시켜, 다항식 방정식 형태의 함수식을 얻는다. 이때 사용된 기호적 열열 기법은 플라톤의 “구조적 방정식”과 유사하게, 복합 구조를 기본적인 원자(스택, 루프, 외부 결합)들의 조합으로 분해한다.

수학적으로는, 형태별 생성함수 F(z)와 전체 구조 생성함수 S(z)를 다음과 같이 관계짓는다: S(z)=Φ(F(z)), 여기서 Φ는 스택 길이와 아크 길이 제약을 반영하는 연산자이다. 이 연산자는 기본적으로 합성, 곱셈, 그리고 합성곱(convolution) 형태를 취한다. 이를 통해 S(z)의 특이점(singularity)을 분석하면, 구조 수의 지수 성장률 ρ와 다항식 차수 α를 정확히 추정할 수 있다. 논문은 이러한 분석을 통해 ρ≈2.1 정도의 매우 작은 성장률을 도출한다. 이는 기존의 일반적인 RNA 2‑체 상호작용 모델이 예측한 ρ≈3.5와 비교해 현저히 낮은 값이며, 실제 생물학적 데이터와도 일치한다는 점에서 의미가 크다.

또한, 논문은 이론적 결과를 검증하기 위해 작은 규모(길이 ≤30)의 완전 열열을 수행하고, 얻어진 계수와 이론적 근사값을 비교한다. 오차는 5% 이하로, 제안된 형태 기반 모델이 실제 구조 분포를 매우 잘 포착함을 보여준다. 마지막으로, 이러한 조합론적 결과가 RNA‑RNA 상호작용 예측 알고리즘에 어떻게 활용될 수 있는지를 논의한다. 구체적으로는, 동적 계획법(DP) 테이블의 상태 공간을 형태 기반으로 축소함으로써, 메모리와 시간 복잡도를 크게 낮출 수 있다.

요약하면, 이 논문은 RNA‑RNA 상호작용 구조를 새로운 형태(shape) 개념으로 추상화하고, 기호적 열열을 통해 정확한 계수와 성장률을 도출함으로써, 기존 MFE 기반 방법의 한계를 보완하고, 효율적인 예측 알고리즘 설계에 이론적 토대를 제공한다.