UNG（尿嘧啶-DNA 糖苷酶） 分子动力学模拟（隐式溶剂）
UniProt P12295 ｜ 起始结构 PDB 3UF7 ｜ OpenMM 8.6.1 · AMBER14 + GBSA-OBC2 · 300 K ｜ 截至 2026-10-17
5.00 ns模拟时长
3739主链原子数（含氢）
5.56 Å主链 RMSD 均值
1.18 ÅCA RMSF 均值
轨迹动画 每帧 5.0 ps，共 1000 帧
若视频不自动播放，点一下播放键；源文件 数据/md/UNG_5ns/md.mp4
左：轨迹回放（结构按残基编号着色）｜右：上=主链 RMSD，下=各残基 CA 涨落
这条蛋白为什么值得看动力学 防污染体系的核心：切除既往扩增产物中的尿嘧啶
UNG 为单域蛋白（UDG 催化域 1–229），D64 为催化质子受体。本次模拟为 apo 态（无 DNA 底物、无 Mg²⁺），只观察该域在 300 K 下的整体柔性。
柔性最大的残基（RMSF 前 8）：
| 残基 | RMSF (Å) |
|---|
| HIS234 | 7.51 |
| HIS237 | 7.30 |
| HIS236 | 7.12 |
| LEU230 | 7.09 |
| HIS233 | 7.04 |
| HIS235 | 6.52 |
| GLU229 | 6.18 |
| GLU231 | 5.98 |
结果解读 自动按本次轨迹生成，不预设结论
- 回转半径由 19.0 Å 降到 17.7 Å（6.9%），结构在模拟中明显压实。这是隐式溶剂（GB）模型的已知偏差——缺少显式水与离子，疏水塌缩被放大；不能据此判断酶在真实缓冲体系中会收缩。
- 最后 1 ns 的 RMSD 均值 6.03 Å，与前 0.5 ns 仅差 0.01 Å，已进入平台段；但 5 ns 仍不足以证明充分平衡，只能说明这一段没有持续单向漂移。
- 柔性最大的区域集中在 N/C 两端（HIS234、HIS237、HIS236、LEU230）；末端涨落在多数酶里属正常，不应解读为活性位点不稳。
- 本体系未包含 DNA 底物、Mg²⁺ 与结晶水：UNG 的催化构象需要底物结合才被组织起来，无底物模拟只能看 apo 态的呼吸运动。
方法与边界 必须一起讲
这是隐式溶剂 MD，不是全原子显式水模拟。
- 溶剂用 GBSA-OBC2 连续介质近似，省去水分子与离子，速度快 10–50 倍，但不能用来判断结合水、离子桥、界面水网络；
- 时长 5.00 ns，只能看局部构象涨落与早期形变，不能据此判断热稳定性排序或货架期；
- 起始结构为晶体结构（含缺失段由 pdbfixer 补齐），预测段与补齐段的运动学解释力弱；
- 单条轨迹、无重复，结论属“观察”，不属“统计显著”。
销售怎么用 只说观察，不说保证
- 讲构象不是静止的：酶在 300 K 下持续涨落，活性位点附近的柔性决定它“容不容易被温度/离子条件带偏”；
- 指到具体残基：RMSF 高的环区往往是突变设计、耐热改造的着力点；
- 边界话术：客户问“能不能证明更耐热”，回答用实测（半衰期、Tm、加速老化），不用 MD 动画代替。
