열역학 구조 합의 원리로 미세RNA 성숙 위치와 구조 예측 및 보존 서브그룹 분류

초록

본 연구는 miRNA의 열역학적 특성과 2차 구조 보존성을 이용해 성숙 miRNA 위치를 예측하고, 종간 보존 서브그룹을 기능 없이 분류한다. 학습 기반 알고리즘으로 44종 2780개의 miRNA 전구체에서 22 bp 성숙 부위를 79.4 % 정확도로 예측했으며, 동물에서는 루프 중심과의 정규화 거리, 식물에서는 ss‑RNA 꼬리와의 정규화 거리가 보존된 절단 신호임을 제시한다.

상세 분석

이 논문은 miRNA 성숙 과정에서 열역학적 안정성(ΔG)과 2차 구조적 특징이 중요한 역할을 한다는 기존 연구를 확장한다. 저자들은 miRBase에 등재된 miRNA 전구체들을 대상으로, 각 종의 최저 자유 에너지 구조와 서열 보존도를 동시에 고려한 ‘열역학·구조 합의 원리’를 제시한다. 이 원리는 “다수 종에서 관찰되는 열역학적·구조적 패턴이 가장 높은 가중치를 가진다”는 가정 하에, 각 전구체에 대해 서열‑구조‑열역학 3차원 공간에서 최적의 합의를 찾는다.

핵심적인 방법론은 다음과 같다. 첫째, 44종(동물·식물 포함)에서 2780개의 miRNA 전구체를 수집하고, RNAfold 등으로 가능한 서열‑구조 후보군을 생성한다. 둘째, 각 후보 구조에 대해 자유 에너지 ΔG와 루프·줄기·말단(ss‑RNA) 길이 비율을 정규화한다. 셋째, 다중 종에서 동일한 miRNA 패밀리(ortholog) 간의 ΔG와 구조적 파라미터 차이를 최소화하는 후보를 ‘합의 구조’로 선정한다. 넷째, 합의 구조를 기반으로 22 bp 성숙 부위를 예측하기 위해 서포트 벡터 머신(SVM) 기반 학습 모델을 훈련한다.

예측 성능은 79.4 % 정확도로, 기존의 전구체 전체 길이 예측(≈70 %)보다 현저히 향상되었다. 특히, 동물 miRNA에서는 성숙 부위가 루프 중심으로부터 줄기의 절반 지점에 위치하는 경향이 강하게 보존되었으며, 이를 ‘정규화 중심 거리’라고 명명한다. 반면, 식물 miRNA에서는 성숙 부위가 ss‑RNA 말단으로부터 일정 비율(≈0.2) 떨어진 위치에 존재함을 발견했다. 이는 Drosha/DCL1 절단 위치가 종마다 다를 수 있음을 시사한다.

또한, 저자들은 열역학·구조 합의에 기반한 서브그룹 분류를 수행했다. 기능적 정보를 투입하지 않고, 오직 보존된 열역학·구조 프로파일만으로 miRBase를 8개의 서브그룹으로 나누었으며, 각 서브그룹은 특정 조직 발현 패턴이나 질병 연관성 등과 통계적으로 연관됨을 사후 분석에서 확인했다. 이는 miRNA의 생물학적 기능이 열역학·구조적 특성과 깊게 연결되어 있음을 뒷받침한다.

이 논문의 주요 기여는 (1) 열역학과 구조 보존성을 동시에 고려한 새로운 합의 원리 제시, (2) 다종 간 일반화 가능한 성숙 miRNA 예측 모델 구축, (3) 동물·식물에서 서로 다른 절단 메커니즘을 제시하는 구조적 지표 도출, (4) 기능 정보를 배제한 보존 서브그룹 분류를 통해 miRNA 진화와 기능 사이의 내재적 연관성을 조명한 점이다. 다만, 실험적 검증이 제한적이며, 비정형 전구체(예: 비전형적인 헤어핀 구조)에서의 적용 가능성은 추가 연구가 필요하다.