단백질 구조에서 잔기 변동성에 대한 서열 구성 및 환경 효과

초록

본 연구는 단백질‑단백질 복합체를 대상으로, 잔기의 질량 분포와 주변 환경이 잔기 변동성에 미치는 영향을 통계‑열역학적으로 분석한다. 질량 비균일성을 반영한 1‑node‑per‑residue 탄성 네트워크 모델을 구축하고, 재정규화된 잔기‑잔기 포텐셜을 적용하였다. 결과는 표면 잔기가 코어 잔기보다 크게 진동하며, 변동성 순위에 따라 아미노산을 고변동·중변동·저변동 세 그룹으로 구분할 수 있음을 보여준다. 특히 루프와 같은 불규칙 2차 구조에서 고변동 잔기의 비율이 높고, 저변동 잔기가 부족한 경우 구조적 불안정성이 증가한다. 인터페이스 잔기는 전체 표면보다 평균적으로 더 강직하며, Gly, Ala, Ser, Cys, Leu, Trp이 안정적인 결합 패치를 형성한다는 점도 확인되었다.

상세 분석

본 논문은 기존의 균일 질량 가정에 기반한 탄성 네트워크 모델(ENM)을 확장하여, 각 잔기의 실제 질량을 반영한 비균일 질량 ENM을 제안한다. 이를 위해 원자 수준의 상호작용을 잔기 수준으로 평균화하고, 잔기‑잔기 포텐셜을 재정규화함으로써 비선형적인 힘 상수를 도입하였다. 모델은 라플라시안 행렬의 고유값·고유벡터를 이용해 각 잔기의 평균 제곱 변위(⟨ΔR²⟩)를 계산하고, 이를 전체 평균으로 정규화하여 상대 변동성 지표를 얻는다.

분석 결과, 두 가지 주요 요인이 변동성에 결정적인 역할을 한다는 것이 밝혀졌다. 첫째, 질량 분포 자체가 변동성을 억제하거나 증폭한다. 무거운 잔기(예: Trp, Tyr)는 관성 효과로 인해 변동이 작고, 가벼운 잔기(Gly, Ala)는 반대로 크게 진동한다. 둘째, 잔기가 위치한 환경, 즉 표면·코어·인터페이스 여부가 변동성에 큰 차이를 만든다. 표면 잔기는 용매와 직접 접촉하면서 제약이 적어 평균 변동이 코어보다 약 1.5배 높았다. 흥미롭게도, 단백질‑단백질 인터페이스는 전체 표면보다 평균 변동이 낮아, 결합 부위가 상대적으로 강직함을 시사한다.

아미노산별 변동성 순위는 Gly > Ala > Ser > Cys > Leu > Trp > … > Ile > Val > Phe > Tyr > Met > His > Glu > Asp > Gln > Asn > Lys > Arg > His > Lys > Arg > … 와 같이 고변동·중변동·저변동 세 그룹으로 구분된다. 특히 루프와 같은 비정형 2차 구조에서는 고변동 잔기의 비율이 평균보다 20 % 이상 높았으며, 이는 루프가 구조적 불안정성의 온상임을 뒷받침한다. 또한, ‘카멜레온 서열’이라 불리는 구조 전이 가능성이 높은 서열에서도 고변동 잔기가 과다하게 존재함을 확인하였다.

이러한 결과는 단백질 설계와 변형 예측에 실질적인 함의를 제공한다. 예를 들어, 변동성이 큰 잔기를 목표로 변형하거나, 저변동 잔기를 보강함으로써 단백질의 열적 안정성을 향상시킬 수 있다. 또한, 인터페이스에서 저변동 잔기가 풍부하게 배치되는 경향은 결합 친화도와 특이성을 높이는 설계 원칙으로 활용될 수 있다.