정규와 포화 RNA 2차 구조의 점근적 특성

초록

본 논문은 RNA 2차 구조 중 정규(canonical)와 포화(saturated) 구조의 조합론적 수량을 점근적으로 분석한다. 정규 구조는 비틀림을 최소화해 계산 속도가 빨라지는 특징을, 포화 구조는 더 이상 염기쌍을 추가할 수 없는 최대 매칭을 의미한다. 저자들은 각각의 구조군에 대해 생성함수와 특이점 분석을 수행해, 전체 구조의 수 $1.104366 n^{-3/2}2.618034^n$에 대비해 정규 구조는 $c_1 n^{-3/2}\alpha^n$, 포화 구조는 $c_2 n^{-3/2}\beta^n$ 형태의 성장률을 보임을 보였다. 또한 포화 구조의 평균 염기쌍 수는 $0.337361n$, quasi‑random 포화 구조는 $0.340633n$이며, 모든 경우에 상태 밀도는 정규분포에 수렴한다는 결론을 제시한다.

상세 분석

논문은 먼저 RNA 2차 구조를 비교적 엄격한 정의(비교적 겹치지 않는 쌍, 최소 루프 길이 3 등) 하에 형식 언어로 모델링하고, 이를 바탕으로 정규(canonical)와 포화(saturated) 구조를 각각 정의한다. 정규 구조는 모든 스템(stem)이 최소 길이 2를 갖고, 스템 사이에 최소 하나의 비스테프(bulge) 혹은 루프가 존재하도록 제한한다. 이러한 제약은 실제 바이오인포매틱스 파이프라인, 특히 Vienna RNA Package에서 동적 계획법(DP) 테이블의 상태 수를 크게 감소시켜 계산 시간을 단축한다는 점에서 실용적 의미가 있다. 저자들은 정규 구조의 생성함수를 $C(z)=\sum_{n\ge0}c_n z^n$ 로 두고, 스템과 루프의 조합 규칙을 이용해 비선형 방정식 $C(z)=z+z^2C(z)+z^3C(z)^2$ 형태로 정리한다. 여기서 $z$는 염기 수에 대한 지표 변수이며, 특이점 $z_0$를 찾기 위해 알게된 근사값을 뉴턴법으로 수치해석한다. 특이점 근처에서 $C(z)\sim K(1-z/z_0)^{-1/2}$ 형태가 되므로, 전형적인 플라네르-아시멜로프 정리를 적용하면 $c_n\sim \kappa n^{-3/2}z_0^{-n}$ 를 얻는다. 계산 결과 $z_0^{-1}= \alpha \approx 2.54027$, $\kappa\approx0.927$ 로, 전체 구조 성장률 $2.618034$ 보다 약간 낮은 값을 보인다. 이는 정규 구조가 전체 구조의 약 $78%$ 정도를 차지함을 의미한다.

포화 구조는 “더 이상 염기쌍을 추가할 수 없는” 최대 매칭을 의미한다. 이를 조합적으로 기술하기 위해 저자들은 “포화 스템-루프”라는 기본 블록을 정의하고, 이 블록이 서로 겹치지 않도록 배치하는 문제를 트리 구조의 생성함수 $S(z)$ 로 전환한다. 포화 구조의 생성함수는 $S(z)=z+z^2S(z)+z^3S(z)^2+z^4S(z)^3$ 와 같은 고차 비선형 방정식으로 귀결된다. 특이점 분석을 통해 $z_0^{-1}= \beta \approx 1.69562$, $\kappa’\approx0.323954$ 를 얻으며, 이는 스템-루프만을 고려한 고전적인 $2^{n-2}$ 성장률에 비해 현저히 낮은 지수이다. 즉, 포화 구조는 전체 구조 중 극히 작은 비율을 차지한다.

또한 저자들은 평균 염기쌍 수를 구하기 위해 $C’(z)$, $S’(z)$ 를 이용한 기대값 연산을 수행한다. 포화 구조의 경우 $\mathbb{E}