새로운 알파벳 의존형 형태 전이 무작위 RNA 정렬

초록

연구진은 무작위 RNA 사슬에서 염기 종류 수 c 가 변함에 따라 2 차 구조의 완전성에 급격한 전이가 일어난다 고 밝혀냈다 c 가 임계값 이하일 때 모든 염기가 짝을 이루어 f 가 1 에 수렴하고 c 가 임계값 초과일 때는 일부 염기가 짝을 못 이루어 f 가 1 미만이 된다 임계값은 2 이상 4 이하로 제한되며 수치 해석은 약 2.7 로 제시한다

상세 분석

본 논문은 무작위 이형 중합체 RNA 를 고전적인 카키드 형태의 2 차 구조 모델에 매핑하여 온도가 0 에 가까워질 때 형성되는 결합된 염기 비율 f 를 알파벳 크기 c 의 함수로 분석한다 저자들은 먼저 c=2 인 경우 모든 서열이 완전 매칭을 가질 수 있음을 구성적 알고리즘을 통해 증명한다 이 알고리즘은 연속된 동일 염기 쌍을 찾아 결합하고 남은 부분을 교대로 매칭하는 방식으로 진행된다 이를 통해 f 가 1 에 수렴함을 보인다 다음으로 c 가 커질수록 완전 매칭을 지원하는 서열의 수가 전체 서열 수에 비해 지수적으로 감소한다는 카탈란 수와 모츠키 경로의 조합을 이용한 조합론적 추정을 제시한다 특히 c<4 일 때는 Dyck 경로에 해당하는 완전 매칭이 전체 경우의 수를 압도하지만 c>4 일 때는 Motzkin 경로가 지배적으로 나타나며 이때는 일정 비율의 공백(갭)이 필연적으로 존재한다 저자들은 이론적 상한을 이용해 f(c) 의 상한식을 도출하고, 수치 시뮬레이션을 통해 n=200 길이의 서열에 대해 f(c) 를 측정한다 결과는 이론적 상한보다 낮으며 c 가 2.7 정도에서 급격히 감소하는 전이를 보인다 또한 베르누이 매칭 모델을 도입해 c 를 실수값으로 연장한 경우에도 전이가 약 2.7 에서 발생함을 확인한다 마지막으로 생물학적 함의를 논의하면서 실제 RNA 가 네 종류의 염기로 구성된 점이 이 임계값 근처에 위치함을 강조한다 이는 RNA 세계 가설과도 연관되어 알파벳 크기가 진화적 최적화에 기여했을 가능성을 시사한다