SBDS 단백질 구조 예측과 변이 분석을 통한 샤와흐만 보디안 다이아몬드 증후군 이해

초록

본 연구는 인간 SBDS 단백질의 3차원 구조를 동종 모델링으로 예측하고, 질환 관련 미스센스 변이를 구조적 관점에서 해석한다. 20여 종의 병원성 변이를 모델에 매핑하여, 핵심 도메인(특히 N‑말단 FYSH 도메인과 C‑말단 RNA 결합 부위)의 안정성 및 상호작용 파괴가 기능 손실을 초래함을 제시한다. 이러한 구조‑기능 연결 고리는 향후 치료제 설계와 진단 마커 개발에 활용될 수 있다.

상세 요약

본 논문은 SBDS 단백질이 현재 실험적으로 확인된 3차원 구조가 부재함에도 불구하고, 생물정보학적 접근을 통해 구조적 단서를 확보하고자 한다. 먼저, 인간 SBDS 서열(Q9Y3A5)을 이용해 BLAST와 PSI‑BLAST를 수행해 고유한 동종 서열을 탐색했으며, 가장 높은 동질성을 보인 아키아균의 YvrA 단백질(약 30% 동일성)을 템플릿으로 삼아 MODELLER 기반 동종 모델링을 진행하였다. 모델의 품질 평가는 DOPE 스코어와 Verify3D, Ramachandran 플롯을 통해 전반적으로 92% 이상의 잔기가 허용 구역에 위치함을 확인하였다.

구조는 크게 세 부분으로 구분된다. 첫 번째는 N‑말단 FYSH 도메인(1‑90 residues)으로, α‑헬릭스와 β‑시트가 교차하는 복합 구조를 이루며, 리보솜 조립 및 RNA 결합에 관여한다는 선행 연구와 일치한다. 두 번째는 중앙의 베타‑프린즈(91‑170 residues)로, 단백질-단백질 상호작용 플랫폼 역할을 한다. 세 번째는 C‑말단 RNA‑결합 도메인(171‑250 residues)으로, 풍부한 양전하를 가진 표면이 tRNA와 같은 작은 RNA와 결합한다는 가설을 뒷받침한다.

이후, 20여 개의 질환 관련 미스센스 변이를 모델에 매핑하였다. 변이들은 크게 두 군으로 나뉜다. 첫 번째 군은 FYSH 도메인 내부에 위치한 고보존 잔기(Lys33, Asp45, Gly58 등) 변이로, 이들 변이는 α‑헬릭스의 코어를 교란하거나 수소 결합 네트워크를 파괴해 전반적인 도메인 안정성을 저하시킨다. 두 번째 군은 중앙 베타‑프린즈와 C‑말단 도메인의 표면에 위치한 잔기(예: Arg124, Phe176, Lys210) 변이로, 이들 변이는 전기적 특성 변화를 일으켜 RNA 혹은 단백질 파트너와의 결합 친화도를 감소시킨다. 특히, Arg124His와 Lys210Glu 변이는 전하 반전으로 인해 기존에 형성되던 염성 상호작용을 소실시키며, 구조역학 시뮬레이션에서 RMSD 증가와 전반적인 유연성 상승을 보였다.

분자 동역학(MD) 시뮬레이션(100 ns, GROMACS) 결과, 변이 단백질은 야생형 대비 평균 1.8 Å의 RMSD 상승과 루프 영역의 B‑factor 증가를 나타냈다. 이는 변이가 국소적인 구조 붕괴를 초래해 전체적인 단백질 안정성을 저해함을 의미한다. 또한, 표면 전하 분포 분석(Pymol APBS)에서는 변이 전후 전하 패턴이 크게 변해, RNA 결합 부위의 전기적 친화도가 감소함을 확인하였다.

이러한 구조‑기능 연관성 분석을 통해, SBDS의 기능 손실은 단순히 단백질 발현 감소가 아니라, 핵심 도메인의 구조적 파괴와 결합 부위의 전하 변형에 기인한다는 결론에 도달한다. 이는 향후 작은 분자 안정화제 혹은 보조 단백질을 설계해 변이로 인한 구조적 결함을 보완하는 전략적 기반을 제공한다.