프리온 AGAAAAGA 아밀로이드 섬유의 원자 수준 모델링 전략

초록

**
본 논문은 비결정성·불용성인 프리온 단백질의 AGAAAAGA 팔린드롬 영역(113‑120)의 아밀로이드 섬유 구조를, 기존 X‑ray·NMR이 불가능한 상황에서 컴퓨터 기반 모델링으로 규명한다. 최신 PDB 템플릿(3NHC)과 시뮬레이티드 어닐링·진화 알고리즘(SAEC)을 이용해 6개의 자유변수만으로 최적화하고, Amber 11로 최종 에너지 최소화를 수행해 원자 해상도의 3가지 모델을 제시한다.

**

상세 분석

**
이 연구는 프리온 단백질의 핵심 병리 부위인 AGAAAAGA(113‑120) 서열이 아밀로이드 섬유를 형성한다는 실험적 증거를 바탕으로, 전통적인 구조생물학적 방법이 적용될 수 없는 경우에 대한 대안을 제시한다. 저자들은 먼저 인간 M129 프리온 펩타이드(127‑132)의 고해상도 스테릭 짝퍼(PDB 3NHC)를 템플릿으로 선택하였다. 3NHC는 β‑시트 간 강한 반데르발스(vdw) 상호작용과 수소결합(H‑bond)으로 안정된 스테릭 짝퍼 구조를 가지고 있어, AGAAAAGA 서열을 동일한 스테릭 짝퍼 형태로 변형시키기에 적합하였다.

모델링 과정은 크게 네 단계로 구성된다. 첫째, SPDBV 툴을 이용해 3NHC의 AB 체인을 AGAAAAGA 서열에 맞게 돌연변이시킨 뒤, GH 체인을 복제해 4개의 체인(ABGH)으로 초기 구조를 만든다. 둘째, 초기 구조에서는 β‑시트 간 vdw 거리가 과도하게 멀어 실제 섬유 형태를 재현하지 못한다는 점을 확인한다. 셋째, A‑체인(ALA3)와 B‑체인(ALA4)의 좌표를 고정하고, G‑체인(ALA4)와 H‑체인(ALA3)의 좌표를 변수화하여 총 6개의 자유변수만 남긴 뒤, Lennard‑Jones(LJ) 퍼텐셜을 최소화한다. 이때 LJ식은 ε=σ=1의 무차원화된 형태로 전개되어, 전통적인 전역 최적화 문제와 동일하게 다룰 수 있다.

저자들은 SAEC(시뮬레이티드 어닐링·진화 컴퓨테이션) 알고리즘을 적용해 6변수 최적화를 수행하였다. SAEC는 기존의 시뮬레이티드 어닐링과 유전 알고리즘을 결합한 하이브리드 방식으로, 복잡한 에너지 지형에서도 효율적으로 전역 최소점을 탐색한다. 최적화 결과, G‑체인과 H‑체인의 CB 원자가 A‑체인·B‑체인과 적절한 vdw 접촉을 이루도록 위치가 조정되었으며, 이는 Fig.6‑8에서 확인할 수 있다.

마지막 단계에서는 Amber 11을 이용해 전체 시스템을 전자밀도 기반 힘장으로 재계산하고, 에너지 최소화와 짧은 동역학 최적화를 수행하였다. 이를 통해 전체 12개의 체인(AB‑GH, CD‑IJ, EK‑FL)으로 확장된 섬유 모델이 안정화되었으며, 최종 구조는 스테릭 짝퍼의 전형적인 특징(β‑시트 간 강한 vdw 및 H‑bond 네트워크)을 유지한다.

핵심적인 과학적 통찰은 다음과 같다. (1) 제한된 자유변수(6개)만으로도 복잡한 아밀로이드 섬유의 핵심 구조를 재현할 수 있다는 점은 고차원 최적화 문제를 저차원으로 축소하는 전략의 효용성을 보여준다. (2) SAEC와 같은 하이브리드 전역 최적화 기법은 비결정성 단백질 섬유와 같이 실험적으로 접근이 어려운 시스템에 적용 가능함을 입증한다. (3) 최종 모델은 향후 약물 설계에서 AGAAAAGA 부위에 결합하는 소분자 혹은 펩타이드 억제제의 구조적 기반을 제공한다.

**