주변 표면 정보 엔트로피를 활용한 단백질‑펩타이드 결합 예측

초록

본 연구는 비상호작용 표면(NIS)의 무극·전하 비율 변동성을 정량화한 ‘주변 표면 정보(PSI) 엔트로피’를 도입한다. HADDOCK3 기반 도킹과 명시적 용매 MD 시뮬레이션을 통해, 친화적인 결합 파트너는 제한된 수의 N‑state(매크로스테이트)로 수렴하며 낮은 PSI 엔트로피를 보인다. WW, PDZ, MDM2 도메인 및 36개의 실험적 WW 구조를 분석한 결과, 이러한 저엔트로피 N‑state는 방법론에 독립적으로 존재함이 확인되었다. PSI 엔트로피는 전역적인 결합 제약을 포착하는 새로운 열역학‑정보학 지표로 제시된다.

상세 분석

이 논문은 기존의 결합 자유에너지(ΔG) 중심 접근법이 비국소적인 표면 효과를 간과한다는 점을 지적하고, 비상호작용 표면(NIS)의 무극(apolar)·전하(charged) 비율을 통계적으로 기술하는 새로운 엔트로피 개념, PSI SΨ를 제안한다. NIS는 인터페이스에 직접 관여하지 않지만, 전체 용매 접근성 및 전기적 환경을 조절함으로써 결합 친화도에 간접적인 영향을 미친다. 저자들은 각 복합체에 대해 잔류물별 용매 접근성(RSA)을 계산하고, RSA≥0.05인 잔류물을 A(무극), C(전하), P(극성) 세 클래스로 분류한다. 인터페이스 잔류물은 제외하고 남은 (n_A, n_C, n_P) 튜플을 매크로스테이트 라벨로 사용한다.

도킹은 HADDOCK3의 3단계 파이프라인(강체, 반유연, 명시적 용매 MD)으로 수행했으며, 에너지 가중치가 단계마다 다르게 적용된다. 각 단계에서 생성된 수천 개의 모델을 확률적으로 샘플링해 마이크로스테이트 Ω를 정의하고, 동일한 (N_a, N_c) 좌표에 속하는 마이크로스테이트 수 g(N_i)를 매크로스테이트의 중복도라 한다. Shannon 엔트로피 S′Ψ=−∑_i (g_i/Ω) log₂(g_i/Ω) 로 정의한 뒤, 접촉 통계에서 얻은 총 가중 접촉 질량 M과 고유 접촉 쌍 수 Q를 이용해 정규화 상수 K=Q/M을 도입, 최종 PSI 엔트로피 SΨ=K·S′Ψ 로 표현한다.

핵심 가설은 ‘친화적인 결합은 N‑state의 변동성을 감소시켜 낮은 SΨ를 보인다’는 것이다. 이를 검증하기 위해 동일 수용체(WW, PDZ, MDM2)에 대해 정식 펩타이드와 비정식 펩타이드를 동일 파라미터로 도킹했다. 결과는 모든 경우에 걸쳐 친화적인 복합체가 소수의 매크로스테이트(예: N≈20~30)로 집중하고, 비친화적 복합체는 수백 개의 매크로스테이트에 분산되는 패턴을 보였다. 또한, 30 ns 명시적 용매 MD 후에도 N‑state 분포는 크게 변하지 않아, PSI 엔트로피가 동역학적 안정성을 반영함을 확인했다.

실험적 검증으로는 36개의 고해상도 WW‑펩타이드 복합체 구조를 메타‑앙상블링해 동일한 N‑state 분석을 수행했다. 여기서도 특정 (N_a, N_c) 조합이 과잉 나타나 저엔트로피 매크로스테이트가 존재함을 발견했으며, 이는 진화적 선택압이 특정 NIS ‘지문’을 보존한다는 가설을 뒷받침한다.

마지막으로 저자들은 SΨ를 이용한 분류 기준을 제시한다. Regime I(낮은 SΨ, 낮은 Q/M)는 고친화, 전형적인 ‘핵심’ 결합을 의미하고, Regime II(높은 SΨ, 높은 Q/M)는 다중특이성 혹은 허브 단백질과 같은 광범위한 상호작용을, Regime III(중간 SΨ, 낮은 Q/M)는 강한 인터페이스 접촉을 유지하면서도 주변 표면에서 다중 N‑state를 허용하는 ‘유연한 친화’를 나타낸다. 이러한 구분은 약물 설계 시 목표 결합부위의 고정성 vs. 유연성을 판단하는 새로운 지표로 활용될 수 있다.

전반적으로 PSI 엔트로피는 기존 에너지 기반 스코어링에 비해 비국소적인 구조적 정보를 정량화함으로써, 결합 친화도 예측의 정확도를 향상시키고, 진화적 보존 패턴을 해석하는 데도 유용함을 입증한다.