Prediction of the self-association of the nuclear inhibitor of DNA binding and differentiation proteins

기술 요약: 핵 내 DNA 결합 억제 단백질의 자기 결합(Self-Association) 예측

문제 제로 (Problem Statement)
DNA 결합 및 분화 억제(Inhibitor of DNA binding and differentiation, ID) 단백질(ID1–ID4)은 세포 주기 진행과 분화에 결정적인 역할을 하는 나선-루프-나선(helix-loop-helix, HLH) 전사 조절 인자이다. 이들이 E-단백질과의 헤테로 이량체 형성(heterodimerization)을 통해 우성 음성 조절자(dominant-negative regulators)로서 작용한다는 사실은 잘 알려져 있으나, 이들의 자기 결합(homodimerization)에 관한 생물물리학적 거동은 여전히 미지의 영역으로 남아 있다. 악성 종양의 맥락에서 ID 단백질이 과발현되는 경우가 많으므로, 단량체(monomer), 동종 이량체(homodimer), 그리고 이종 이량체(heterodimer) 사이의 균형은 단백질의 안정성, 핵 내 국소화 및 기능에 결정적인 영향을 미칠 수 있다. 특히, ID 단백질의 자기 결합을 결정하는 구조적 결정 요인과 생리적 환경에서의 복합체 열역학적 안정성을 규명하는 것은 병리적 상호작용을 차단하기 위한 잠재적 치료 표적을 식별하는 데 필수적이다.

방법론 (Methodology)
본 연구는 ID 단백질의 자기 결합을 모델링하고 분석하기 위해 다단계 계산 생물물리학 및 생물정보학 파이프라인을 채택하였다:

1. 구조 모델링 (Structural Modeling): 실험적 구조가 부재한 이들 고유 무질서 단백질(intrinsically disordered proteins)을 위해 AlphaFold3, I-TASSER, Phyre2를 사용하여 ID1–ID4의 3D 모델을 생성하였다. 생성된 모델은 MolProbity(Ramachandran plot, rotamer analysis, clash scores)를 통해 엄격하게 검증되었으며, GalaxyRefine 서버를 통해 정교화되었다. 입체 화학적 품질이 가장 우수한 AlphaFold3 유래 모델을 주요 후보로 선정하였다.
2. 단량체 분자 동역학 (Monomer Molecular Dynamics, MD): 정교화된 단량체 구조는 CHARMM36m 포스 필드와 GROMACS를 사용하여 명시적 용매(TIP3P) 환경에서 100 ns 전원자(all-atom) MD 시뮬레이션을 거쳤다. 궤적(trajectories)은 구조적 역동성을 특징짓기 위해 RMSD, RMSF, 회전 반경(Radius of Gyration, Rg), 용매 접근 가능 표면적(Solvent Accessible Surface Area, SASA), 그리고 분자 내 수소 결합을 분석하는 데 사용되었다.
3. 클러스터링 (Clustering): 대표적인 구조적 형태를 식별하기 위해, 대규모의 유연한 시스템에서 구조적 유사성을 더 잘 포착할 수 있는 전통적인 RMSD 대신 공통 접촉 분율(Fraction of Common Contacts, FCC) 지표를 사용하여 궤적 클러스터링을 수행하였다.
4. 단백질-단백질 도킹 (Protein-Protein Docking): 동종 이량체(ID1:ID1, ID2:ID2 등) 및 이종 이량체(예: ID1:ID4) 상호작용을 조사하기 위해, 가장 밀도가 높은 클러스터로부터 추출된 대표적 구조들을 HADDOCK을 사용하여 도킹하였다. HLH 도메인을 활성 부위로 정의하였으며, 복합체는 HADDOCK 점수와 예측 결합 친화도(PRODIGY)에 따라 순위가 매겨졌다.
5. 복합체 안정성 및 결합 자유 에너지 (Complex Stability and Binding Free Energy): 상위 점수를 받은 복합체(ID1–ID1 및 ID1–ID4)는 안정성을 검증하기 위해 추가적인 100 ns MD 시뮬레이션을 거쳤다. 결합 자유 에너지($\Delta G_{binding}$)는 분자 역학/포아송-볼츠만 표면적(MM/PBSA) 방법을 사용하여 계산되었으며, 에너지 기여도를 반데르발스, 정전기적, 극성 용매화 및 비극성 용매화 항으로 분해하였다.

주요 결과 (Key Results)

• 구조 모델링: AlphaFold3는 rotamer 비율 및 Ramachandran 선호 영역 측면에서 I-TASSER 및 Phyre2보다 일관되게 우수한 초기 모델 품질을 보였다. GalaxyRefine을 통한 정교화 후, 고품질 구조(예: ID3는 100% Ramachandran 선호 잔기 달성)를 확보하여 후속 분석에 적합한 상태가 되었다.
• 단량체 역동성: MD 시뮬레이션은 ID1–ID4의 고유 무질서 특성을 확인해주었다. 모든 단량체는 높은 RMSD 값(26–30 Å)과 상당한 유연성(RMSF), 특히 N-말단 및 C-말단에서의 유연성을 보였다. 그러나 모든 단백질은 자발적인 소수성 압축(hydrophobic compaction, Rg 및 SASA 감소)을 경험하였고 일시적인 분자 내 수소 결합을 형성하였는데, 이는 용액 내에서 확장된 형태에서 부분적으로 구조화된 앙상블로의 전이를 시사한다.
• 이량체화 안정성: 도킹 및 후속 MD 시뮬레이션 결과, 이량체 형성이 ID 단백질을 유의미하게 안정화하는 것으로 나타났다.
  • RMSD 감소: ID1–ID1 동종 이량체(6.84 Å)와 ID1–ID4 이종 이량체(9.93 Å)의 RMSD는 각각의 단량체보다 약 3배 낮았으며, 이는 이량체 형성이 구조적 엔트로피를 억제함을 나타낸다.
  • 유연성: 두 이량체 모두 단량체에 비해 현저히 낮은 RMSF(평균 ~3.8–4.2 Å)를 보였으며, 10 Å를 초과하는 잔기가 없었다. 이는 HLH 계면이 상호 간의 구조적 구속을 가함을 확인시켜 준다.
  • 조밀함: 단량체와 달리, 이량체는 100 ns 시뮬레이션 동안 안정적이고 조밀한 4차 구조(일정한 Rg)를 유지하였다.
• 결합 에너지 (MM/PBSA):
  • ID1–ID1 동종 이량체: 주로 반데르발스 상호작용($\Delta E_{VDW} = -43.70$ kcal/mol)에 의해 안정화되었으며, 이는 효율적인 소수성 패킹을 의미한다. 총 $\Delta G_{binding}$은 $-28.07$ kcal/mol이었다.
  • ID1–ID4 이종 이량체: 더 강한 전체 결합 친화도($\Delta G_{binding} = -41.55$ kcal/mol)를 나타냈다. 이 복합체는 상당한 탈용매화 비용($\Delta E_{GB} = +418.97$ kcal/mol)에도 불구하고, 주로 정전기적 상호작용($\Delta E_{EL} = -433.58$ kcal/mol)에 의해 구동되었다.
  • 비교: ID1–ID4 복합체가 ID1–ID1 동종 이량체보다 에너지적으로 더 유리한 것으로 나타났으며, 이는 시뮬레이션된 조건 하에서 이종 결합(heterotypic association)을 선호함을 시사한다.

의의 및 주장 (Significance and Claims)
저자들은 본 연구가 ID 단백질의 자기 결합에 대한 구조적 토대를 성공적으로 지도화하였으며, 이들 단백질이 HLH 도메인을 통해 안정적인 동종 이량체 및 이종 이량체를 형성할 수 있음을 입증했다고 주장한다. 핵심적인 발견은 ID1–ID4 이종 이량체가 ID1–ID1 동종 이량체보다 더 강한 결합 친화도와 열역학적 안정성을 가진다는 점이다.

본 논문은 ID1–ID4 복합체의 형성이 핵 내 자유 ID1 및 ID4 단량체의 농도를 낮추는 메커니즘으로 작용하여 이들의 생물학적 활성을 조절할 수 있다고 상정한다. ID 단백질이 암 진행에서 수행하는 역할과 이들의 잠재적 치료 표적임을 고려할 때, 이러한 특정 상호작용 패턴과 결합 핫스팟(binding hotspots)의 식별은 ID 단백질의 활성이 어떻게 조절되는지 이해하기 위한 분자적 근거를 제공한다. 저자들은 이러한 계산적 통찰이 ID1–ID4 복합체의 생물학적 관련성에 대한 향후 실험적 조사와 이러한 병리적 상호작용을 표적으로 하는 억제제 개발을 위한 기초를 마련한다고 결론짓는다.