[论文解读] Protein Folding: A Perspective From Statistical Physics

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

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），以极低的计算成本研究蛋白质折叠动力学，同时保留关键的物理相互作用。

提出的方法

• 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}$，与实验观测一致。
• 从预球状态到熔融球状态的转变伴随着氢键形成速率的突然增加，该现象在肌红蛋白的模拟中已被观察到。
• 推导出临界链长 $N_{\text{crit}} = 5.34 \times 10^7$，超过此长度后排斥体积变为正值，由于塌陷后氢键数量不足，稳定折叠将不再可能。
• 有效排斥体积 $\bar{c}$ 在 $\rho_{\max} \approx 6.85$ 处由负值（吸引）转变为正值（排斥），表明体系从紧密态向扩展态转变。
• 该模型识别出一种仅依赖于 $R_g$ 和 $N$ 的普遍弹性能量，为不同蛋白质家族中的标度律提供了物理解释。