[논문 리뷰] Equivalence Classes of Optimal Structures in HP Protein Models Including Side Chains

Martin Mann, Rolf Backofen외 1명|ArXiv.org|2009. 10. 20.

Protein Structure and Dynamics참고 문헌 13인용 수 4

한 줄 요약

이 논문은 HP 단백질 모델에서 최적의 구조에 대해 동치 관계를 도입하여, 동일한 수소화성 잔기 배치를 가진 구조를 그룹화함으로써 대표 구조 수를 크게 감소시킨다. CPSP 방법을 확장하여 각 동치 클래스당 하나의 최적 대표 구조(핵심 다중성)를 열거함으로써, 구조적 다양성은 최대 세 계수 감소하고, 이는 국소 탐색 시뮬레이션의 효율적인 초기화를 가능하게 하며, 총 다중성보다도 핵심 다중성이 더 안정적인 구조적 안정성 측정치임을 제안한다.

ABSTRACT

Lattice protein models, as the Hydrophobic-Polar (HP) model, are a common abstraction to enable exhaustive studies on structure, function, or evolution of proteins. A main issue is the high number of optimal structures, resulting from the hydrophobicity-based energy function applied. We introduce an equivalence relation on protein structures that correlates to the energy function. We discuss the efficient enumeration of optimal representatives of the corresponding equivalence classes and the application of the results.

연구 동기 및 목표

• 특히 측쇄 포함 모델에서 발생하는 HP 격자 단백질 모델의 최적 구조의 조합 폭발 문제를 해결하기 위해.
• 국소 탐색 초기화에 사용하기 위한 대표 최적 구조 수를 줄여 시뮬레이션 효율성을 향상시키기 위해.
• 핵심 다중성—수소화성 핵심 구성에 대응하는 대표 구조 수—을 제안하여 총 다중성보다 더 의미 있는 구조적 안정성 측정치가 되도록 하기 위해.
• 제약 프로그래밍을 활용해 최적의 동치 클래스를 효율적으로 열거함으로써 CPSP 프레임워크를 확장하기 위해.

제안 방법

• 두 구조가 격자 내 동일한 위치에 수소화성(H) 잔기를 가진 경우 동치로 간주하는 동치 관계를 정의한다.
• 최적의 H-핵심 구성에 집중하여, 각 동치 클래스당 하나의 대표 구조를 열거하는 데 목적이 있는 CPSP 접근법을 확장한다.
• 제약 프로그래밍을 사용하여 주어진 시퀀스에 대해 모든 최적의 H-핵심 포장(최대 접촉 구성)을 생성한다.
• 자기 회피와 에너지 최소화를 유지하면서, 각 최적의 H-핵심에 대해 나머지 극성(P) 잔기를 스레드한다.
• 3D 입방격자 및 면심격자에서 백본 전용 및 측쇄 포함 HP 모델 모두에 이 방법을 적용한다.
• HPrep 소프트웨어를 사용하여 테스트 시퀀스 전반에 걸쳐 핵심 다중성과 총 다중성을 비교한다.

실험 결과

연구 질문

• RQ1 접힘과 관련된 구조적 다양성을 유지하면서도 HP 모델의 최적 단백질 구조 수를 어떻게 줄일 수 있는가?
• RQ2 수소화성 핵심 배치를 기반으로 한 동치 관계가 최적 구조의 더 효율적이고 의미 있는 표현을 가능하게 할 수 있는가?
• RQ3 핵심 다중성은 구조적 안정성과 어느 정도 상관이 있으며, 총 다중성보다 더 나은 안정성 지표가 될 수 있는가?
• RQ4 모든 최적 구조를 전부 열거하는 것에 비해, 각 동치 클래스당 하나의 대표 구조를 열거하는 데 얼마나 효율적인가?

주요 결과

• 백본 전용 모델에서는 평균 핵심 다중성이 4,856이었으며, 평균 다중성 124,800에 비해 두 계수 감소한 결과를 보였다.
• 측쇄 포함 모델의 경우 평균 핵심 다중성은 약 1,550이었고, 평균 다중성은 최소 10^6으로 추정되어, 최소 세 계수 감소한 것으로 나타났다.
• 모든 테스트 시퀀스에서 백본 전용 모델의 경우 핵심 다중성 대표 수가 10^6 미만이었고, 총 다중성이 10^6 이하인 경우는 66%에 그쳤으며, 측쇄 포함 모델의 경우 이 비율은 4.08%에 불과했다.
• 핵심 다중성 대표를 열거하는 데 평균 실행 시간은 시퀀스당 2초였으며, 총 다중성 열거보다 6배 빠른 속도를 기록했다.
• 핵심 다중성 측정치는 제약이 없는 P-단량체 배열로 인한 다중성의 영향을 효과적으로 걸러내어, 더 안정적이고 의미 있는 구조적 내성 지표로 기능한다.
• 이 방법은 수소화성 붕괴를 반영하는 고품질의 소규모 다양성 있는 시작 구조 집합을 제공함으로써 국소 탐색 시뮬레이션의 효율적인 초기화를 가능하게 한다.