大分子结构核磁预测：二级化学位移与实验化学位移关联

本项目基于 120 个残基的蛋白结构预测 ¹³Cα 与 ¹³Cβ 二级化学位移，并与 BMRB 实验值对比。使用 SHIFTS/SHIFTX2 等经验模型，共 460 个被归属核。通过 DSSP 定义二级结构（螺旋、折叠、无规），分析不同结构元素的位移分布。

案例分析

预测值与实验值的散点图显示 ¹³Cα R²=0.982，¹³Cβ R²=0.971，RMSD 为 0.94 ppm，说明模型对 460 个核的化学位移预测精度高。二级结构元素统计显示螺旋区 ¹³Cα 化学位移整体正偏移，折叠区负偏移，与典型二级化学位移规律一致。94% 残基的 DSSP 二级结构归属与基于化学位移的判定一致。

模型细节、验证与工程意义

位移预测用 SHIFTX2 网页服务器 / SHIFTS 软件，输入为 PDB 结构。二级结构由 DSSP 2.0 判定。460 个核从 BMRB 27418 条目中提取，已排除未归属核与化学交换峰。验证时 ¹³Cα R²=0.982、¹³Cβ R²=0.971，RMSD 0.94 ppm 接近 SHIFTX2 方法本身精度上限，说明结构质量高。

总结教训

第一个教训是化学位移预测对结构精度敏感。若晶体结构中含未解析 loop，RMSD 会上升 0.2–0.3 ppm。第二个教训是 ¹³Cβ 比 ¹³Cα 更易受侧链构象影响，预测偏差略大。第三个教训是氢键与环电流效应是二级化学位移的主要来源，经验模型已较好包含这些效应。

项目经验

1. 蛋白 NMR 位移预测基于 120 残基结构。
2. 预测核：460 个被归属的 ¹³Cα/¹³Cβ。
3. ¹³Cα R² = 0.982，¹³Cβ R² = 0.971。
4. RMSD = 0.94 ppm。
5. 二级结构由 DSSP 定义：helix、strand、coil。
6. 螺旋区 ¹³Cα 正偏移，折叠区负偏移。
7. DSSP 与位移法二级结构一致性 94%。
8. 结构精度对 RMSD 影响显著，未解析 loop 会增大约 0.2–0.3 ppm。
9. ¹³Cβ 比 ¹³Cα 更易受侧链构象影响。
10. 氢键与环电流是二级位移的主要来源。
11. 经验模型（SHIFTX2/SHIFTS）已较好涵盖主要物理效应。
12. 二级结构 DSSP 2.0 自动判定。
13. 460 个核来自 BMRB 27418，已排除未归属核。
14. 本项目方法可迁移到结构验证、无序区域识别。

补充经验

1. 可考虑加入分子动力学系综平均，提升柔性区域的位移预测。
2. pH 与质子化状态会影响羧基/胺基 nearby 核的位移。
3. 柔性 loop 区域建议做系综平均，可进一步降低 RMSD。
4. 金属结合位点附近核的化学位移需单独参数化。

复现参数速查（案例图硬参数）

更多相关仿真经验，可参考本站的 MD 仿真案例栏目。