[논문 리뷰] Protein Folding: A Perspective From Statistical Physics

Jinzhi Lei, Kerson Huang|arXiv (Cornell University)|2010. 02. 26.

Protein Structure and Dynamics참고 문헌 50인용 수 6

한 줄 요약

이 논문은 통계역학 기반의 군집화된 접근법인 조건부 자가피避보행(CSAW) 모델을 소개한다. 이 모델은 수용액 내에서 고분자로서의 폴리펩타이드 사슬을 확률적 구조로 간주하며, 수소결합과 친수성 수축을 조건부 제약 조건을 통해 통합한다. 모델은 편평한 상태, 프리글로불, 녹은 글로불 상태에 대응하는 반경의 궤도($R_g \sim N^\nu$)에 대한 보편적 스케일링 법칙을 예측하며, 각각의 지수는 $\nu = 3/5$, $3/7$, $2/5$이다. 또한, 안정적인 접힘을 위한 이론적 한계($N_{\text{crit}} \approx 5.34 \times 10^7$)를 규명하여, 이 이상의 크기에서는 수축 후 수소결합의 부족으로 인해 안정적인 접힘이 불가능하다고 밝힌다.

ABSTRACT

In this paper, we introduce an approach to the protein folding problem from the point of view of statistical physics. Protein folding is a stochastic process by which a polypeptide folds into its characteristic and functional 3D structure from random coil. The process involves an intricate interplay between global geometry and local structure, and each protein seems to present special problems. We introduce CSAW (conditioned self-avoiding walk), a model of protein folding that combines the features of self-avoiding walk (SAW) and the Monte Carlo method. In this model, the unfolded protein chain is treated as a random coil described by SAW. Folding is induced by hydrophobic forces and other interactions, such as hydrogen bonding, which can be taken into account by imposing conditions on SAW. Conceptually, the mathematical basis is a generalized Langevin equation. To illustrate the flexibility and capabilities of the model, we consider several examples, including helix formation, elastic properties, and the transition in the folding of myoglobin. From the CSAW simulation and physical arguments, we find a universal elastic energy for proteins, which depends only on the radius of gyration $R_{g}$ and the residue number $N$. The elastic energy gives rise to scaling laws $R_{g}\sim N^ν$ in different regions with exponents $ν=3/5,3/7,2/5$, consistent with the observed unfolded stage, pre-globule, and molten globule, respectively. These results indicate that CSAW can serve as a theoretical laboratory to study universal principles in protein folding.

연구 동기 및 목표

• 서열에 의존하지 않는 보편적인 물리적 원리를 반영하면서도 계산 효율성이 높은 모델을 개발하는 것.
• 접힘 경로를 이끄는 친수성 수축과 수소결합 간의 상호작용을 조사하는 것.
• 접힘 중간체인 편평한 상태, 프리글로불, 녹은 글로불 상태에서 반경의 궤도($R_g$)에 대한 스케일링 법칙을 유도하는 것.
• 안정적인 접힘이 열역학적으로 불가능한 단백질 크기의 이론적 상한선을 결정하는 것.
• 최소한의 계산 비용으로도 핵심 물리적 상호작용을 유지하면서 단백질 접힘 동역학을 연구할 수 있는 이론적 실험실을 구축하는 것.

제안 방법

• CSAW 모델은 수용액 내에서의 자가피避보행(SAW)으로서, 배제 체적 효과를 반영한 무작위 고리 상태를 나타내는 페프타이드 사슬을 처리한다.
• 친수성 상호작용은 SAW에 조건부 제약 조건을 도입하여 밀도 있는 구조를 선호하도록 하여 친수성 수축을 시뮬레이션한다.
• 수소결합은 친수성과 별개의 에너지 항으로서 펩타이드 뼈대 상호작용을 기반으로 잠재력에 통합된다.
• 모델은 열역학적 분위기 내에서 사슬의 확률적 동역학을 기술하기 위해 일반화된 랑주방정식 프레임워크를 사용한다.
• 몬테카를로 시뮬레이션을 통해 배제 체적, 친수성, 수소결합의 복합 영향을 받는 구조를 샘플링한다.
• 효과적인 배제 체적 $\bar{c}$ 는 물리적으로 가능한 사슬 구성의 집합 평균을 통해 계산되며, 이는 임계 사슬 길이 유도에 기여한다.

실험 결과

연구 질문

• RQ1 단백질 접힘의 다양한 단계에서 반경의 궤도($R_g$)에 대한 보편적 스케일링 법칙은 무엇이며, 이는 편평한 사슬, 프리글로불, 녹은 글로불 상태와 어떻게 관련이 있는가?
• RQ2 친수성 수축과 수소결합 간의 균형이 접힘 중간체의 조밀함과 안정성에 어떻게 영향을 미치는가?
• RQ3 CSAW와 같은 단순화된 군집화 모델이 최소한의 계산 비용으로도 단백질 접힘의 동역학 및 열역학적 특징을 재현할 수 있는가?
• RQ4 친수성 수축과 수소결합 형성의 제약 조건을 고려할 때, 안정적으로 접을 수 있는 단백질의 이론적 최대 크기는 얼마인가?
• RQ5 접힘 모델의 물리적 기반에서 친수성과 수소결합 효과의 분리가 어느 정도 정당화될 수 있는가?

주요 결과

• CSAW 모델은 20잔기 폴리알라닌 사슬을 안정적인 알파 나선으로 접히게 하여, 서열만으로도 천연 구조를 재현할 수 있음을 입증한다.
• 모델은 반경의 궤도에 대해 세 가지 명확한 스케일링 영역을 예측한다: 편평한 상태에서는 $R_g \sim N^{3/5}$, 프리글로불 상태에서는 $R_g \sim N^{3/7}$, 녹은 글로불 상태에서는 $R_g \sim N^{2/5}$이며, 이는 실험 관측과 일치한다.
• 프리글로불에서 녹은 글로불로의 전이가 myoglobin의 시뮬레이션에서 관찰된 것처럼 수소결합 형성 속도의 급격한 증가로 특징지어진다.
• 임계 사슬 길이 $N_{\text{crit}} = 5.34 \times 10^7$ 가 도출되었으며, 이 이상의 크기에서는 배제 체적이 양이 되어 수축 후 수소결합의 부족으로 안정적인 접힘이 불가능하다.
• 효과적인 배제 체적 $\bar{c}$ 는 $\rho_{\max} \approx 6.85$ 에서 음성(끌림)에서 양성(반발력)으로 전이되며, 이는 조밀한 상태에서 편평한 상태로의 전이를 나타낸다.
• 모델은 $R_g$ 와 $N$ 에만 의존하는 보편적인 탄성 에너지를 규명하여, 다양한 단백질 가족에 걸쳐 스케일링 법칙의 물리적 기반을 제공한다.