[논문 리뷰] Enumerating Designing Sequences in the HP Model

Anders Irbäck, Carl Troein|arXiv (Cornell University)|2002. 01. 02.

Protein Structure and Dynamics참고 문헌 10인용 수 7

한 줄 요약

이 논문은 최적화된 구조적 공간 압축 및 서열 탐색 알고리즘을 사용하여 2D HP 격자 모델에서 설계 서열의 정확한 수를 N=25까지 확장한다. 설계된 서열이 무작위 서열보다 높은 수의 수소친화성/비수소친화성 덩어리 수를 유지함을 입증하며, 이는 음성 수소친화성 상관관계를 확인하는 바이며, 이 특성은 최대 밀도 구조로 제한된 경우 사라지며, 단백질 접힘 모델링에서 구조적 집합 선택의 중요성을 부각시킨다.

ABSTRACT

The hydrophobic/polar HP model on the square lattice has been widely used to investigate basics of protein folding. In the cases where all designing sequences (sequences with unique ground states) were enumerated without restrictions on the number of contacts, the upper limit on the chain length N has been 18-20 because of the rapid exponential growth of the numbers of conformations and sequences. We show how a few optimizations push this limit by about 5 units. Based on these calculations, we study the statistical distribution of hydrophobicity along designing sequences. We find that the average number of hydrophobic and polar clumps along the chains is larger for designing sequences than for random ones, which is in agreement with earlier findings for N up to 18 and with results for real enzymes. We also show that this deviation from randomness disappears if the calculations are restricted to maximally compact structures.

연구 동기 및 목표

• 2D HP 모델에서 설계 서열의 정확한 수를 N=18–20을 초월해 상한을 확장하기 위해.
• 서열 길이 증가에 따라 설계 서열의 통계적 특성, 특히 수소친화성 상관관계가 어떻게 변화하는지 조사하기 위해.
• 최대 밀도 구조로 제한된 구조에 영향을 미치는 설계 서열의 통계적 특성에 미치는 영향을 평가하기 위해.
• HP 모델 내에서 설계된 서열이 실제 단백질과 유사한 협동성 접힘 전이를 보일 수 있는지 평가하기 위해.

제안 방법

• 구조적 공간을 압축하기 위해 접촉 집합 표현을 사용하여 고유한 구조의 수를 줄이되, 필수 구조적 정보를 손상시키지 않음.
• 계산 시간을 기존의 난이도 높은 탐색 방식 대비 약 천 배 감소시킨 최적화된 서열 공간 탐색 알고리즘을 구현함.
• 열역학적 분석을 위해 고안된 안정성이 높은 서열을 생성하기 위해 몬테카를로 기반 서열 설계 알고리즘을 사용함.
• 접힘 전이 및 비열 행동을 연구하기 위해 우산 샘플링 기반 몬테카를로 시뮬레이션을 수행함.
• 전체 구조적 공간과 최대 밀도 구조 제약 조건 하에서의 설계 서열 간 통계적 특성(특히 덩어리 수)을 비교함.
• 블록 변수 통계를 통해 수소친화성 상관관계를 분석하여 서열 규칙성을 정량화함.

실험 결과

연구 질문

• RQ1 N=18에서 N=25로 서열 길이를 연장했을 때, 설계 서열의 음성 수소친화성 상관관계 특성이 유지되는가?
• RQ2 전체 구조적 공간과 최대 밀도 구조 서브셋 간에 설계 서열의 통계적 특성(예: 덩어리 수)은 어떻게 다를까?
• RQ3 최대 밀도 구조로 제한된 구조적 집합으로의 제약이 실제 단백질 서열과 비교해 설계 서열의 통계적 특성에 얼마나 큰 영향을 미치는가?
• RQ4 HP 모델 내에서 설계된 서열이 실제 단백질과 유사한 협동성 접힘 전이를 보일 수 있는가?
• RQ5 모델 내에서 최적화된 안정성 서열의 열역학적 특성은 일반적인 서열과 어떻게 비교되는가?

주요 결과

• 최적화된 알고리즘을 사용하여 2D HP 모델에서 설계 서열의 정확한 수를 N=18–20에서 N=25까지 성공적으로 확장함.
• N=25에서 설계 서열은 무작위 서열보다 수소친화성 및 비수소친화성 덩어리의 평균 수가 유의미하게 높으며, 이는 음성 수소친화성 상관관계를 확인함.
• N을 25로 연장해도 덩어리 통계에서의 랜덤성 이탈은 강건하게 유지되며, 이는 소형 시스템 크기의 오류가 아니라는 것을 시사함.
• 최대 밀도 구조로 제한된 경우, 설계 서열의 통계적 특성(특히 덩어리 수)은 무작위 서열과 구분되지 않으며, 생물학적 관련성이 상실됨을 나타냄.
• 안정성 최적화된 서열의 접힘 전이에서는 날카로운 비열 피크와 저에너지 상태와 옹호 상태의 공존을 보이며, 일반적인 서열에 비해 협동성이 향상됨을 나타냄.
• 결과는 구조적 집합 선택이 설계 서열의 통계적 특성에 결정적인 영향을 미치며, 최대 밀도 서브셋은 실제 단백질 서열과의 일치도가 낮은 결과를 도출함을 시사함.