Quantum chemical modelling of sequence encoded energetics and electrostatic landscapes governing early peptide association in α- synuclein and amyloid β

기술 요약: 초기 펩타이드 결합 과정에서의 서열 인코딩된 에너지 및 정전기적 지형의 양자 화학적 모델링

문제 제기
단백질 미스폴딩(misfolding)과 응집은 $\alpha$-시뉴클레인($\alpha$-synuclein) 및 아밀로이드 $\beta$(A$\beta$)와 같은 펩타이드의 초기 분자 결합에 의해 구동되는 파킨슨병(PD) 및 알츠하이머병(AD)과 같은 신경 퇴행성 질환의 핵심 요소이다. 실험 연구를 통해 성숙한 섬유(fibril)의 키네틱스와 구조는 규명되었으나, 과도기적이고 불균일한 초기 결합 단계 동안 개별 분자 구성이 기여하는 에너지적 기여도는 여전히 명확히 해결되지 않은 상태로 남아 있다. 특정 선형 잔기 순서(서열)가 전체 아미노산 조성과는 별개로, 초기 결합 사건을 지배하는 분자적 특성을 어떻게 형성하는지에 대한 이해에는 중요한 격차가 존재한다. 본 연구는 응집 경향이 있는 펩타이드 세그먼트와 조성이 동일한 스크램블드(scrambled) 대조군을 비교함으로써 서열 배열의 기여도를 분리해내고, 이를 통해 잔기 순서가 에너지, 구조적 조직 및 정전기적 지형에 미치는 영향을 결정하는 것을 목표로 한다.

방법론
본 연구는 단량체(monomer) 및 동종 이량체(homodimer) 상태를 모델링하기 위해 준경험적 양자 화학(semi-empirical quantum chemistry)과 밀도 범함수 이론(DFT)을 통합한 일관된 계산 워크플フロー를 채택한다.

• 시스템: 응집 경향이 있는 펩타이드 세그먼트($\alpha$-시뉴클레인: VTGVTAVAQK; A$\beta$: LVFFAEDVGS)를 조성 기반의 스크램블드 대조군(AGATQTVVVK 및 EFLVGSDVAF)과 비교하였다.
• 기하 구조 최적화: 단량체 및 동종 이량체 구조는 xTB에 구현된 GFN2-xTB(Geometry, Frequency, Noncovalent extended Tight Binding, second generation) 방법을 사용하여 최적화되었다. 이량체 구조는 최종 최적화 전 GFN 포스 필드를 사용하여 먼저 완화(relaxation)되었다.
• 에너지론: 결합 에너지($\Delta E_{assoc}$)는 최적화된 이량체의 수렴된 GFN2-xTB 전자 에너지와 독립적으로 최적화된 단량체 에너지의 두 배 사이의 차이로 계산되었다.
• 전자 구조 분석: 최적화된 단량체에 대해 PBE 교환-상관 범함수(exchange-correlation functional)와 실공간 그리드(0.20 Å 간격) 상의 투영 평면파(Projector Augmented Wave, PAW) 형식론을 적용한 GPAW를 사용하여 DFT 계산을 수행하였다.
• 기술자(Descriptors):
  • 정전기학: 분자 정전기 퍼텐셜(ESP)을 전자 밀도 표면(0.001 e Bohr$^{-3}$)에 매핑하고, 상대적 지형을 시각화하기 위해 평균 중심화(mean-centering)를 통해 분석하였다.
  • 전하 분포: 원자 해상 전하 분포를 정량화하기 위해 Hirshfeld 인구 분석을 사용하였으며, 여기에는 통계적 지표(평균, 중앙값, 표준 편차) 및 잔기별 합계가 포함된다.
  • 구조/서열: 기술자로는 서열 전하 장식(SCD), 국소 친수성(hydropathy), 소수성 실행 길이(hydrophobic run length), 회전 반경($R_g$), 말단 간 거리(end-to-end distance) 및 비구형성(asphericity)이 포함된다.
  • 인터페이스: 중원자 거리($\le$ 4.5 Å 및 5.0 Å)와 특정 N/O/S 접촉($\le$ 3.5 Å)을 통해 분자 간 접촉을 정량화하였다.

주요 결과
연구 결과, 잔기 재배열은 전체 아미노산 조성은 보존됨에도 불구하고 에너지 및 구조 영역에서 펩타이드 특이적인 뚜렷한 반응을 생성함을 보여준다.

1. 단량체 에너지론: 잔기 순서는 고립된 단량체의 전자 에너지를 유의미하게 변화시켰다. A$\beta$의 경우, 스크램블드 대조군이 응집 경향이 있는 서열보다 약 45.33 kJ/mol 더 낮았다. 반대로, $\alpha$-시뉴클레인의 경우 응집 경향이 있는 서열이 스크램블드 대조군보다 약 151.12 kJ/mol 더 낮았다.
2. 결합 에너지 vs 인터페이스 패킹:
   • A$\beta$: 응집 경향이 있는 서열(115.617 kJ/mol)과 스크램블드 대조군(232.552 kJ/mol) 모두 양(+)의 결합 에너지를 나타냈으며, 이는 최적화된 이량체가 독립적으로 최적화된 단량체보다 전자 에너지가 높음을 의미한다. 그러나 응집 경향이 있는 서열은 스크램블드 대조군보다 적은 수의 중원자 접촉(40 대 81)과 낮은 접촉 밀도를 보였다. 이는 역설적으로 응집 경향 서열이 더 낮은 양의 결합 에너지를 가짐을 나타낸다.
   • $\alpha$-시뉴클레인: 스크램블드 대조군은 응집 경향이 있는 서열(-172.819 kJ/mol)보다 더 유리한(더 음의 값인) 결합 에너지(-494.371 kJ/mol)를 나타냈다. 그러나 응집 경향 서열은 더 광범위한 인터페이스를 형성하였으며(166개 중원자 접촉, 대 105개), 더 높은 밀도의 근접 N/O/S 접촉(15 대 5)을 보였다.
   • 결론: 결합 에너지와 인터페이스 패킹(접촉 수/밀도)은 선형적으로 변하지 않으며, 결합의 상보적이지 않은, 즉 서로 보완적인 척도임을 나타낸다.
3. 정전기적 지형 (ESP): 잔기 재배열은 ESP의 공간적 분포를 재구성하였다. A$\beta$와 $\alpha$-시뉴클레인 쌍은 유사한 전역 전하 상태를 공유했으나, 양(+) 및 음(-)의 퍼텐셜 영역의 구체적인 공간적 조직은 응집 경향 서열과 스크램블드 서열 사이에서 크게 변화했다. 이는 서열 순서가 정전기적 특성의 3D 제시 방식을 결정함을 시사한다.
4. Hirshfeld 전하 분포: 원자 전하 분석 결과, 잔기 재배열은 전역 전하 크기를 균일하게 이동시키기보다는 주로 원자 전하 인구(population)를 재분배하였다. $\alpha$-시뉴클레인의 경우 양/음성 원자의 비율은 거의 일정하게 유지되었으나, 전하 값의 확률 분포가 이동하였다. A$\beta$에서는 양/음성 원자의 균형에 대한 완만한 변화와 함께 재분배가 관찰되었다.
5. 서열 및 구조 기술자:
   • A$\beta$: 재배열은 서열 전하 장식(SCD), 소수성 실행 길이(5 잔기 대 3 잔기), 그리고 분자 압축성을 변화시켰다. 응집 경향 펩타이드는 스크램블드 대조군(약 11.0 Å)에 비해 약 10.1 Å의 $R_g$를 나타냈다.
   • $\alpha$-시뉴클레인: 재배열은 주로 국소 친수성(최대 3-잔기 평균 친수성이 대조군에서 ~2.6에서 ~4.2로 증가)을 변화시켰으며, SCD, $R_g$, 비구형성은 유사하게 유지되었다.

의의 및 주장
본 논문은 잔기 순서가 아미노산 조성과는 독립적으로, 초기 펩타이드 결합의 근본적인 결정 요인임을 확립한다. 본 연구는 동일한 조성이라도 서열 의존적인 분자 조직을 통해 뚜렷한 에너지적, 구조적, 정전기적 특성을 생성할 수 있음을 입증한다.

주요 주장은 다음과 같다:

• 결정 요인으로서의 서열: 잔기의 배열은 단순히 아미노산의 조성을 넘어, 분자 간 상호작용의 균형을 지배하며, 결과적으로 발생하는 결합의 열역학 및 기하학적 인터페이스를 결정한다.
• 펩타이드 특이적 반응: 잔기 재배열의 분자적 결과는 보편적이지 않다. $\alpha$-시뉴클레인과 A$\beta$는 서열 변화가 에너지적 유리함과 인터페이스 패킹의 서로 다른 조합으로 이어지는 대조적인 반응을 보인다.
• 상보적 기술자: 본 연구는 결합 에너지와 인터페이스 접촉 지표가 분자 결합의 서로 다른 측면을 포착한다는 점을 강조한다. 두 펩타이드 시스템에서 나타난 상이한 경향처럼, 더 많은 인터페이스가 반드시 더 유리한 결합 에너지를 의미하는 것은 아니다.
• 정전기적 재구성: 동일한 아미노산 조성이라도 전하 패턴과 3D 구조의 재구성을 통해 뚜렷하게 다른 정전기적 환경을 생성할 수 있으며, 이는 서열 특이적 결합 행동의 분자적 기초를 제공한다.

이 연구는 서열 인코딩된 분자 결정 요인을 식별하기 위한 정량적 프레임워크를 제공하며, 초기 응집의 "코드"가 잔기의 구체적인 선형 순서에 있으며, 이것이 결과적인 에너지 및 정전기적 지형을 결정한다는 점을 시사한다.