[논문 리뷰] Heterogeneous path ensembles for conformational transitions in semi-atomistic models of adenylate kinase

Divesh Bhatt, Daniel M. Zuckerman|arXiv (Cornell University)|2009. 10. 08.

Protein Structure and Dynamics참고 문헌 54인용 수 5

한 줄 요약

이 연구는 반구조적 모델을 사용하여 아데닐레이트 키나제의 열린 상태와 닫힌 상태 사이의 구조적 전이를 탐색하기 위해 가중치가 부여된 앙상블(WE) 경로 샘플링 시뮬레이션을 적용한다. 이는 경로의 이질성을 드러내며, 뚜껑과 AMP 결합 도메인의 닫힘 순서에 따라 다릅니다. 평형 조건 하에서 정방향 및 역방향 전이 경로 간에 약간의 대칭성을 보이며, 이는 이론적 대칭 원칙을 검증하고 최소한의 계산 비용으로 효율적인 샘플링을 가능하게 한다.

ABSTRACT

We performed "weighted ensemble" path-sampling simulations of adenylate kinase, using several semi-atomistic protein models. Our study investigated both the biophysics of conformational transitions as well as the possibility of increasing model accuracy without sacrificing good sampling. Biophysically, the path ensembles show significant heterogeneity and the explicit possibility of two principle pathways in the Open-Closed transition. We recently showed, under certain conditions, a "symmetry of hetereogeneity" is expected between the forward and the reverse transitions: the fraction of transitions taking a specific pathway/channel will be the same in both the directions. Our path ensembles are analyzed in the light of the symmetry relation and its conditions. In the realm of modeling, we employed an all-atom backbone with various levels of residue interactions. Because reasonable path sampling required only a few weeks of single-processor computing time with these models, the addition of further chemical detail should be feasible.

연구 동기 및 목표

• 가중치가 부여된 앙상블 시뮬레이션에서 유도된 경로 앙상블을 사용하여 아데닐레이트 키나제의 구조적 전이의 생물물리학을 조사하기.
• 반구조적 표현 방식에서 모델 정확도를 향상시키되 샘플링 효율성을 희생하지 않는 것이 가능한지 평가하기.
• 평형 안정 상태에서 정방향 및 역방향 전이 간 '이질성의 대칭성' 이론 예측을 검증하기.
• 모델 내 잔기 상호작용 세부 정보 수준의 변화가 전이 경로에 미치는 민감도 평가하기.
• 화학적 세부 정보가 향상된 더 큰 생체분자에 대해 전체 경로 샘플링의 계산 가능성을 결정하기.

제안 방법

• 외부 힘 없이 아데닐레이트 키나제의 구조적 전이를 시뮬레이션하기 위해 가중치가 부여된 앙상블(WE) 경로 샘플링 방법 적용.
• 모든 원자 뼈대를 포함하고 잔기 및 잔기 상호작용 세부 정보 수준을 다양하게 설정한 반구조적 모델 사용.
• 정방향(열림 → 닫힘)과 역방향(닫힘 → 열림) 전이를 모두 시뮬레이션하여 비교 및 대칭성 검증 가능하게 하기.
• 경로 앙상블의 수렴을 보장하고 통계적 무게를 유지하기 위해 재샘플링 기법 활용.
• 도메인 닫힘 사건의 순서(뚜껑 대비 AMP 결합 도메인)를 추적하여 경로 이질성 분석.
• 정방향 및 역방향의 전이 속도와 경로 비율을 비교하여 대칭 조건 검증.

실험 결과

연구 질문

• RQ1 반구조적 모델링 하에서 아데닐레이트 키나제의 열린 상태 ↔ 닫힌 상태 전이에 대해 주요 구조적 전이 경로는 무엇인가?
• RQ2 잔기 상호작용 세부 정보 수준의 차이가 관측된 전이 메커니즘에 어느 정도 영향을 미치는가?
• RQ3 평형 통계역학에 따르면 예측된 바와 같이, 정방향 및 역방향 전이 간 경로 분포에 대칭성이 존재하는가?
• RQ4 반구조적 모델을 사용한 중간 크기 단백질에서 전체 경로 앙상블 생성에 대해 가중치가 부여된 앙상블 샘플링의 계산 효율성은 어떠한가?
• RQ5 아데닐레이트 키나제(214개 잔기)와 같은 작은 도메인을 초월한 더 큰 시스템에 대해 합리적인 계산 비용으로 경로 앙상블을 신뢰성 있게 생성할 수 있는가?

주요 결과

• 열린 상태 → 닫힌 상태 전이에서 뚜껑과 AMP 결합 도메인의 닫힘 순서에 따라 다릅니다.
• 모델 1 및 2에서 닫힌 상태 → 열린 상태 전이의 샘플링이 상당히 어려웠으며, 닫힌 상태-도메인-뚜껑-열린 상태 경로가 관측되지 않아, 분지 정의가 불완전하거나 안정 상태에 도달하지 못한 것으로 추정됨.
• 모델 3에서는 양방향에서 유사한 샘플링이 가능했으며, 경로 비율에 대해 약간의 대칭성을 보였고, 이는 이론적 '이질성의 대칭성'과 일치함.
• 모델 3에서 닫힌 상태 → 열린 상태 전이의 전이 속도는 2.5 × 10⁻⁶ per Monte Carlo step로 추정되었으며, 이는 브루트 포스 시뮬레이션으로 단일 전이를 위한 약 4년이 소요될 것임을 시사함.
• 가중치가 부여된 앙상블 시뮬레이션은 단일 프로세서 시간 약 2주 만에 50개의 수렴된 전이를 달성하여, 브루트 포스 방법 대비 상당한 효율성 향상을 보임.
• 이 연구는 반구조적 모델을 사용한 전체 경로 샘플링이 계산적으로 가능하며, 효소 기작의 고해상도 모델링으로의 길을 열어줌을 확인함.