[논문 리뷰] Flexibility of $\alpha$-helices: Results of a statistical analysis of database protein structures

Eldon Emberly, Ranjan Mukhopadhyay외 2명|arXiv (Cornell University)|2002. 09. 25.

Protein Structure and Dynamics참고 문헌 18인용 수 8

한 줄 요약

이 연구는 단백질 구조은 보관소(PDB)의 399개 단백질 구조를 대상으로 주성분 분석(PCA)을 수행하여 α-나선의 내재된 유연성을 정량화한다. 세 가지 주요 모드—두 개의 거의 디제너레이트된 굽힘 모드와 한 개의 비틀림 모드—를 규명하였으며, 이들의 진폭은 독립적인 가우시안 분포를 따르며, 준열적 평형과 일관된다. 수소 결합과 배타 부피를 고려한 단순한 스프링 모델이 이러한 모드를 성공적으로 재현하여, 저에너지의 집합 운동을 통합함으로써 단백질 구조의 정밀 조정과 설계를 향상시킬 수 있다.

ABSTRACT

$\alpha$-helices stand out as common and relatively invariant secondary structural elements of proteins. However, $\alpha$-helices are not rigid bodies and their deformations can be significant in protein function ({\it e.g.} coiled coils). To quantify the flexibility of $\alpha$-helices we have performed a structural principal-component analysis of helices of different lengths from a representative set of protein folds in the Protein Data Bank. We find three dominant modes of flexibility: two degenerate bend modes and one twist mode. The data are consistent with independent Gaussian distributions for each mode. The mode eigenvalues, which measure flexibility, follow simple scaling forms as a function of helix length. The dominant bend and twist modes and their harmonics are reproduced by a simple spring model, which incorporates hydrogen-bonding and excluded volume. As an application, we examine the amount of bend and twist in helices making up several coiled-coil proteins. Incorporation of $\alpha$-helix flexibility into structure refinement and design is discussed.

연구 동기 및 목표

• 고정체가 아닌, α-나선의 내재된 유연성을 정량화하기 위해 고려한다.
• α-나선 내의 집합적 변형이 열적 평형과 일관된 통계적 패턴을 따르는지 확인한다.
• 단백질 구조 모델링에 사용할 수 있는 주요 나선 유연성을 포괄하는 최소한의 물리 모델을 개발한다.
• 코일드 코일과 같은 생물학적으로 중요한 시스템에서 나선 유연성의 관련성을 평가한다.
• α-나선의 부드럽고 집합적인 모드를 통합함으로써 보다 정확한 단백질 구조의 정밀 조정과 설계를 가능하게 한다.

제안 방법

• SCOP에서 추출한 399개의 대표적 단백질 구조를 대상으로 Cα 위치를 이상적인 나선(3.6개 잔류물/회전, 1.5 Å/잔류물 상승)에 정렬하여 반복적으로 평균 나선 구조를 계산하였으며, 수렴 기준은 10−4 Å/잔류물 이내였다.
• 정렬된 나선들 간의 Cα 좌표 편차로부터 3L × 3L 구조 공분산 행렬을 구성하고, 이를 대각화하여 고유값과 고유벡터를 도출하였다.
• 에너지 함수의 수치적 이阶도를 사용하여 수소 결합 강도와 배타 부피를 두 가지 매개변수로 갖는 스프링 모델을 구축하여 주요 PCA 모드를 재현하였다.
• 모델의 정규 모드를 PCA에서 도출된 모드와 비교하여 검증하였으며, 저에너지 모드에서 뛰어난 일치를 보였다.
• 코일드 코일 단백질(루신 지퍼, Mnt 테트라머화 도메인, 피브리노겐)에 모델을 적용하여 관측된 변형에 기여하는 모드의 기여도를 평가하였다.

실험 결과

연구 질문

• RQ1 자연 상태의 단백질 구조에서 길이가 다른 α-나선의 주요 집합적 유연성 모드는 무엇인가?
• RQ2 이러한 유연성 모드의 진폭은 열적 평형과 일관된 통계적 분포를 따르는가?
• RQ3 수소 결합과 배타 부피를 고려한 단순한 스프링 모델이 실제 단백질 구조의 PCA에서 관측된 주요 모드를 재현할 수 있는가?
• RQ4 다양한 코일드 코일 단백질, 특히 긴 나선을 가진 단백질에서 굽힘 및 비틀림 모드의 기여도는 어떻게 달라지는가?
• RQ5 이러한 부드러운 모드를 통합함으로써 단백질 구조의 정밀 조정과 설계는 어느 정도 향상될 수 있는가?

주요 결과

• 세 가지 주요 유연성 모드—두 개의 거의 디제너레이트된 굽힘 모드와 한 개의 비틀림 모드—가 규명되었으며, 이들은 α-나선의 대부분의 구조적 분산을 설명한다.
• 이 모드들의 진폭은 상호 독립적인 가우시안 분포를 따르며, 특성 온도 T∗ ≈ 실온 수준의 준열적 평형과 일관된다.
• 주요 모드의 고유값은 나선 길이에 따라 단순한 거듭제곱 법칙에 따라 변화하며, 크기 의존성의 체계적 특성을 나타낸다.
• 수소 결합과 배타 부피만을 고려한 최소한의 스프링 모델이 저에너지 정규 모드를 포함한 고조화 모드까지 성공적으로 재현하였다.
• 코일드 코일 단백질에서, 짧은 나선(예: 루신 지퍼, Mnt)은 첫 번째 굽힘 및 비틀림 모드로 잘 기술되지만, 더 긴 나선(예: 피브리노겐)은 고차원 굽힘 조화 모드가 필요하다.
• 결과적으로, 단백질 모델링에 이러한 부드러운 모드를 통합할 경우, 전체 원자 유연성 없이도 에너지 최적화와 구조적 현실감을 향상시킬 수 있다.