[论文解读] CSAW: a dynamical model of protein folding

Kerson Huang|arXiv (Cornell University)|Jan 12, 2006

Protein Structure and Dynamics参考文献 4被引用 4

一句话总结

CSAW（条件化自避行走）是一种基于随机过程与蒙特卡洛方法的动力学模型，通过结合自避行走与基于能量的接受规则来模拟蛋白质折叠，以描述疏水效应和氢键作用。该模型再现了关键折叠动力学，包括快速的疏水塌陷形成熔融球状中间体，以及随后缓慢转变为稳定的天然构象，其半径回转量随时间的变化趋势与实验数据定性一致。

ABSTRACT

CSAW (conditioned self-avoiding walk) is a model of protein folding that combines SAW (self-avoiding walk) with Monte-Carlo. It simulates the Brownian motion of a chain molecule in the presence of interactions, both among chain residues, and with the environment. In a first model that includes the hydrophobic effect and hydrogen bonding, a chain of 30 residues folds into a native state with stable secondary and tertiary structures. The process starts with a rapid collapse into an intermediate "molten globule", which slowly decays into the native state afer a relatively long quiescent period. The behavior of the radius of gyration mimics experimental data.

研究动机与目标

• 开发一种计算高效、物理机制透明的蛋白质折叠模型，以捕捉关键的生物物理作用力。
• 研究疏水相互作用与氢键如何协同驱动三级结构与二级结构的形成。
• 模拟折叠的动态路径，包括熔融球状中间体，并与实验观测结果进行比较。
• 实现对结构与能量特征（如主链、侧链、氢键）的逐步引入，以系统研究其作用机制。
• 提供一个灵活的框架，用于研究折叠机制，同时避免全原子分子动力学模拟的计算开销。

提出的方法

• 采用由扭转角定义的、由'曲柄'连接而成的主链模型来表示蛋白质链。
• 通过硬球排斥实现条件化自避行走（SAW），以强制残基之间空间不重叠。
• 采用蒙特卡洛更新与梅特罗普利斯接受准则，模拟热涨落与能量最小化过程。
• 引入构象能量函数：E = -g₁K₁ - g₂K₂，其中K₁表示疏水残基接触数，K₂表示氢键数量。
• 将疏水性与氢键作用设为可调参数（g₁, g₂），并以kBT = 1进行归一化处理。
• 通过几何与角度约束建模氢键：O–H距离 ≤ 2 Å，且呈反平行排列。

实验结果

研究问题

• RQ1 仅靠疏水相互作用如何驱动蛋白质链形成明确的三级结构？
• RQ2 氢键在折叠过程中对稳定α-螺旋等二级结构起到何种作用？
• RQ3 该模型能否再现实验观测到的熔融球状中间体状态？
• RQ4 半径回转量如何随时间演变？其动力学是否与实验折叠动力学一致？
• RQ5 从熔融球状状态向天然构象的转变是否表现出成核控制的首阶相变特征？

主要发现

• 仅疏水相互作用即可形成明确的三级结构，尽管存在显著的构象波动。
• 引入氢键后，三级结构更加稳定，并诱导出α-螺旋的形成。
• 链经历快速塌陷形成熔融球状中间体，随后缓慢过渡至天然构象，与实验折叠动力学一致。
• 半径回转量的演化在塌陷阶段表现出短暂扩张，可能暗示类似雷诺稀化（Reynolds dilatancy）的行为。
• 熔融球状中间体表现出与最终天然构象相似的致密性与结构特征（如α-螺旋与环结构），表明其是合理的折叠中间体。
• 从熔融球状状态到天然构象的转变在不同运行中并非确定性过程，表明其为成核驱动的随机过程。