Rosetta가 풀지 못하는 네 가지 작은 퍼즐

초록

Rosetta는 대형 단백질·RNA 구조 예측에 성공했지만, 20‑아미노산 이하의 초소형 시스템에서는 실험 구조와 차이를 보인다. 저자들은 Trp‑cage, 콘오톡신, 세린 프로테아제 억제제 루프, UUCG RNA 테트라루프 네 가지 사례를 통해 샘플링은 충분히 수행되었음에도 불구하고 현재 Rosetta 전원 원자 에너지 함수의 근사와 누락된 항목이 원자 수준 정확도의 모델을 구분하지 못한다는 점을 입증한다. 이러한 “퍼즐”은 Rosetta 에너지 함수의 근본적 개선을 위한 벤치마크로 활용될 수 있다.

상세 분석

본 논문은 Rosetta가 가장 기본적인 구조 예측 과제에서도 한계를 드러낸다는 점을 네 가지 초소형 모델을 통해 정량적으로 보여준다. 첫 번째 사례인 20‑잔류 Trp‑cage는 기존 Rosetta 프로토콜로 수천 개의 구조를 생성했음에도 불구하고, 최저 에너지 모델이 실험 구조와 RMSD ≈ 2 Å 이상 차이 나는 것이 관찰되었다. 이는 샘플링이 충분히 이루어졌음에도 불구하고, 전원 원자 Lennard‑Jones, 전하‑전하 상호작용, 수소 결합 방향성 등 현재 에너지 함수가 작은 루프와 고밀도 패키징을 정확히 평가하지 못함을 의미한다.

두 번째 퍼즐인 4‑잔류 디설파이드 콘오톡신은 이론적으로는 단일 이황화 결합과 짧은 β‑시트만으로 정의될 수 있다. Rosetta는 이 결합을 올바르게 형성하지만, 주변 잔기의 회전 자유도와 용매 효과를 과소평가해 비실험적 구조가 에너지적으로 우세하게 된다. 특히, 용매 접근성(FA)와 탈수 효과를 반영하는 term이 부재하거나 과도하게 보정된 것이 원인으로 지적된다.

세 번째 사례인 세린 프로테아제 억제제의 반응 루프는 6‑잔류 길이의 고유동성 구간으로, 루프의 입체적 제약과 전하 분포가 핵심이다. Rosetta는 루프의 전반적인 토폴로지를 재현하지만, 정확한 전하‑전하 상호작용과 물 분자 매개 결합을 무시함으로써 실험 구조와 에너지 격차가 발생한다. 이는 현재의 implicit solvent 모델이 작은 루프의 전자 구름을 충분히 묘사하지 못한다는 증거다.

마지막으로 UUCG RNA 테트라루프는 RNA 2차 구조의 대표적인 안정화 모티프이며, 4‑잔류만으로도 특유의 비틀림과 스택을 형성한다. Rosetta의 RNA 에너지 함수는 스택 상호작용과 비인접 염기 간 전기적 상호작용을 단순화하고, Mg²⁺와 같은 이온 효과를 무시한다. 결과적으로, 가장 낮은 에너지 모델은 비실험적 비틀림을 보이며, RMSD ≈ 1.8 Å 수준의 오차를 남긴다.

공통적으로 저자들은 “샘플링은 충분히 수행되었다”는 가정을 검증하기 위해 대규모 Monte‑Carlo 및 fragment‑assembly 시뮬레이션을 10⁴ ~ 10⁵ 회 수행했으며, 생성된 구조 집합 내에서 실험 구조와 거의 동일한 RMSD를 갖는 모델이 존재함을 확인했다. 따라서 근본적인 병목은 에너지 함수의 정밀도 부족이며, 특히 작은 시스템에서의 비극성·극성 상호작용, 용매 효과, 이온 매개 결합을 정밀하게 모델링할 필요가 있다.

이러한 분석은 Rosetta 개발자들에게 전원 원자 에너지 함수의 재구축, 더 정교한 implicit/explicit solvent 모델 도입, 그리고 작은 루프·RNA 모티프에 특화된 파라미터 튜닝이 필요함을 강력히 시사한다.