[论文解读] Folding Mechanism of Small Proteins
本研究提出了一种新颖的蒙特卡洛折叠模拟方法，采用单一优化的原子级势函数，将四种小蛋白（betanova、1fsd、HP-36 和 protein A）折叠为天然态类似结构。该方法揭示折叠过程受热力学稳定性和动力学因素共同调控，早期折叠轨迹会收敛至特定的坍塌态，从而决定折叠路径是横向还是对角线地接近天然态，这一动态特性无法仅通过自由能面捕捉。
Extensive Monte Carlo folding simulations for four proteins of various structural classes are carried out, using a single atomistic potential. In all cases, collapse occurs at a very early stage, and proteins fold into their native-like conformations at appropriate temperatures. The results demonstrate that the folding mechanism is controlled not only by thermodynamic factors but also by kinetic factors: The way a protein folds into its native structure, is also determined by the convergence point of early folding trajectories, which cannot be obtained by the free energy surface.
研究动机与目标
- 开发一种统一的模拟方法，仅使用单一原子级势函数即可将多种小蛋白折叠至其天然态。
- 在传统自由能面分析之外，研究折叠过程中动力学与热力学因素的相互作用。
- 理解早期折叠轨迹如何影响最终折叠路径及向天然结构的收敛。
- 考察结构类别（如β-折叠与α-螺旋束）在决定折叠机制与协同性中的作用。
- 探索低温折叠模拟中玻璃态行为与非遍历动力学的出现。
提出的方法
- 采用包含静电、范德华力及多体项的连续原子级相互作用的联合残基（UNRES）力场。
- 使用构象空间退火方法同时识别低能量极小点并优化四种蛋白的力场参数。
- 应用马尔可夫链蒙特卡洛动力学，通过小角度扰动模拟从非天然构象出发的折叠过程，持续时间达 10⁹ MCS。
- 以与天然结构的RMSD、回转半径（Rg）以及天然接触分数（Q 和 ρ）作为关键可观测量追踪折叠进程。
- 将 ρ 定义为天然态模拟中天然接触概率的加权平均值，以捕捉涨落效应。
- 在 (ρ, Rg) 空间分析轨迹收敛性，以区分动力学折叠路径（横向或对角线）。
实验结果
研究问题
- RQ1折叠轨迹的早期收敛点在 (ρ, Rg) 空间中如何影响整体折叠路径？
- RQ2与仅依赖热力学稳定性相比，动力学因素在多大程度上决定了折叠机制？
- RQ3单一优化的原子级势函数能否成功将多种不同结构类别的蛋白折叠为天然态类似构象？
- RQ4二级结构组成（β-折叠与α-螺旋）在塑造折叠动力学与协同性方面发挥何种作用？
- RQ5温度如何影响折叠模拟中玻璃态非遍历行为的出现？
主要发现
- 所有四种蛋白均成功折叠为天然态类似构象，RMSD 值低至 0.78 Å（betanova）和 1.07 Å（1fsd），表明结构精度极高。
- 快速坍塌在约 10⁴ MCS 内发生，随后为较慢的折叠过程，在 (ρ, Rg) 空间中观察到明显的不同折叠路径。
- 对于 betanova 和 1fsd，早期轨迹分别收敛至 (ρ, Rg) ≈ (0.35, 8.5 Å) 和 (0.3, 9 Å)，随后沿横向路径向天然态移动。
- 对于 HP-36 和 protein A，轨迹分别收敛至 (ρ, Rg) ≈ (0.2, 11 Å) 和 (0.25, 12 Å)，表现出更多对角线折叠路径。
- 在低温下（betanova 的 T ≤ 50，HP-36 的 T ≤ 60）观察到玻璃态非遍历行为，该行为在 T ≥ 70–90 时消失。
- 在蛋白 A 中观察到一级相变样坍塌转变，发生在 T ≈ 120 附近，从 (Q, Rg) ≈ (0.15, 18 Å) 转变为 (0.15, 12 Å)。
更好的研究，从现在开始
从阅读论文到最终审阅，大幅缩短您的研究时间。
无需绑定信用卡
本解读由 AI 生成，并经人工编辑审核。