作者Evil Genius每个组学都需要学习很多的内容大家如果选择学习就一定要学会不学就好好放松一下分析的时候该找公司找公司只能说自己学会性价比最高学会了之后积累经验即可请公司做那就只能一直请公司。虽然现在都有Ai了但是Ai的过程与结果还是需要人来判断内部逻辑参数如果不理解最终的结果还是不行改都不知道怎么改。这一篇我们来补充关于分子动力学模拟的过程中添加辅酶、金属离子等内容课程中我们学习了如何添加配体分子我们来添加一些其他必需的部分。对于GROMACS 分子动力学MD如果蛋白体系中存在辅酶、金属离子、底物、配体等非标准组分不能简单地把它们当成普通水分子加入solv.gro。正确做法是结构 → 参数 → 拓扑 → 组装 → 溶剂/离子 → 能量最小化 → 平衡 → 生产模拟。1. 整体流程例如一个蛋白中有Zn²⁺ NAD⁺ 底物蛋白 Protein │ ┌─────────────────┼─────────────────┐ │ │ │ Zn²⁺ NAD⁺ Ligand │ │ │ └─────────────────┴─────────────────┘ │ ▼ 建立复合物结构 │ ▼ 为非标准组分准备参数 │ ┌──────────────┼──────────────┐ ▼ ▼ ▼ 金属参数 辅酶参数 配体参数 │ │ │ └──────────────┼──────────────┘ ▼ topology.top │ ▼ 加水 TIP3P/SPC │ ▼ 添加 Na⁺ / Cl⁻ │ ▼ 能量最小化 │ ▼ NVT │ ▼ NPT │ ▼ Production MD2. 最重要的一点先判断它是什么不同非标准组分处理方式完全不同。组分典型例子参数化难度是否需要单独拓扑普通离子Na⁺、Cl⁻、K⁺★GROMACS 力场通常已有金属离子Mg²⁺、Zn²⁺、Ca²⁺、Fe²⁺★★–★★★★★通常需要谨慎处理辅酶NAD⁺、FAD、FMN、CoA★★★★通常需要参数小分子配体ATP、底物、抑制剂★★★通常需要参数糖/脂质UDP-glucose、膜脂★★★★通常需要参数血红素Heme★★★★★金属配位尤其重要金属簇Fe-S cluster★★★★★需要专门参数/模型所以Na⁺、Cl⁻ 和 NAD⁺、Zn²⁺并不是同一种处理方式。3. 最简单的情况Na⁺ / Cl⁻如果只是为了中和体系电荷设置 150 mM NaCl模拟生理盐环境GROMACS 本身就可以处理。例如gmx genion \ -s ions.tpr \ -o solv_ions.gro \ -p topol.top \ -pname NA \ -nname CL \ -neutral或者gmx genion \ -s ions.tpr \ -o solv_ions.gro \ -p topol.top \ -pname NA \ -nname CL \ -conc 0.15这类离子通常不需要你自己制作.itp。4. 如果是 Mg²⁺、Zn²⁺、Ca²⁺怎么办这里就复杂很多。例如Protein │ ├── ASP ├── GLU ├── HIS │ └── Zn²⁺不能简单认为Zn²⁺ Protein ↓ 加入一个 ZN 原子 ↓ MD因为金属和蛋白之间可能存在配位键静电作用配位几何约束水分子配位氨基酸侧链配位例如 Zn²⁺可能形成HIS \ Zn²⁺ / \ CYS HIS如果参数处理不合理MD过程中可能出现初始结构 Zn / | \ HIS CYS HIS ↓ MD Zn \ HIS 其他配位键逐渐消失因此金属离子不能只看“有没有这个离子”还要考虑你希望模拟什么样的金属配位模型。5. 金属参数有几种思路大体可以分成方法 A非键模型把金属作为普通带电粒子Zn²⁺ electrostatic interaction van der Waals interaction优点简单容易建立适合某些体系缺点不能很好描述方向性的配位键配位结构可能不稳定方法 B键合模型人为建立Zn—N Zn—O Zn—S例如HIS-N │ │ Zn²⁺ │ │ CYS-S需要定义bondangledihedralatom typecharge这种方法可以维持特定配位结构但它引入了人为约束。方法 C更高级的金属模型例如针对特定金属中心采用bonded modelnonbonded modeldummy atom model12-6-4 型模型QM/MM对于复杂金属酶QM/MM 往往比简单经典力场更合适尤其当金属中心涉及反应、价态变化或明显的配位化学变化时。6. 辅酶NAD⁺、FAD、ATP 等这部分和普通金属离子不同。例如Protein │ ├──── NAD │ └──── substrate你需要准备NAD │ ├── structure ├── atom types ├── charges ├── bonds ├── angles └── dihedrals最终可能得到NAD.itp然后在topol.top里面#include NAD.itp最后[ molecules ] Protein 1 NAD 1 SOL 30000 NA 10 CL 107. 一个很关键的问题辅酶属于“配体”吗从 GROMACS 拓扑角度看通常可以把辅酶作为独立的分子类型处理。例如Protein NAD ATP Mg2 Water Na Cl-它们可以分别存在于体系中。但从生物化学角度辅酶是否应该独立参数化、是否需要与蛋白形成特殊键合模型要根据实际结构和力场来决定。例如 ATP-Mg²⁺Protein │ │ ATP ─── Mg²⁺这里 Mg²⁺可能与 ATP 的磷酸基团发生明显配位。因此不能简单地把ATP.itp Mg.itp扔进去就结束。8. 实际工作中推荐这样组织文件例如做一个蛋白 NAD⁺ Mg²⁺ ligand可以组织成MD/ │ ├── protein.pdb ├── NAD.pdb ├── ligand.pdb ├── Mg.pdb │ ├── NAD.itp ├── ligand.itp ├── Mg.itp │ ├── topol.top ├── solv.gro ├── ions.tpr ├── solv_ions.gro │ ├── em.mdp ├── nvt.mdp ├── npt.mdp └── md.mdp9. 拓扑文件的关系最终topol.top │ ┌────────────────┼────────────────┐ │ │ │ ▼ ▼ ▼ Protein topology NAD.itp ligand.itp │ │ │ │ ▼ ▼ │ NAD Ligand │ ▼ Mg.itp │ ▼ Mg²⁺然后[ molecules ] Protein 1 NAD 1 LIG 1 MG 1 SOL 30000 NA 8 CL 8这里最重要的是坐标文件中的分子数量必须和[ molecules ]中一致。10. 推荐的实际操作顺序建议把流程固定成① 准备蛋白protein.pdb去除不需要的水缓冲液杂质crystallization additives但不要一上来把辅酶、金属、底物全部删掉。② 判断 PDB 中的辅酶/金属例如Protein NAD MG HOH确认NAD 是否是真正的辅酶 MG 是否是功能性 Mg²⁺ 水分子是否参与金属配位③ 给非标准组分参数化例如NAD ↓ NAD.itp NAD.gro以及Ligand ↓ LIG.itp LIG.gro④ 组装复合物最终protein NAD Mg ligand ↓ complex.pdb⑤ 建立 GROMACS topologytopol.top包含Protein NAD LIG MG⑥ 加水gmx solvate得到solv.gro⑦ 加离子gmx genion加入Na Cl-⑧ 能量最小化EM ↓ 检查 Potential Energy Max Force⑨ NVT控制Temperature⑩ NPT控制Pressure Density⑪ Production MD例如100 ns 200 ns 500 ns然后分析RMSD RMSF Rg SASA H-bond Distance Interaction MM/PBSA11. 特别推荐保留实验结构中的辅酶和金属如果蛋白是从 PDB 晶体结构得到的例如Protein NAD Mg通常不要简单地remove HETATM然后重新“猜”辅酶和金属位置。因为晶体结构已经给出了一个实验观察到的Protein ↓ Cofactor ↓ Metal ↓ Water空间关系。比较合理的是实验 PDB │ ├── Protein ├── Cofactor ├── Metal └── Key waters │ ▼ 参数化/检查 │ ▼ GROMACS12. 重点掌握这张图实验结构 / AlphaFold / Docking │ ▼ ┌──────────────┐ │ Protein │ │ Cofactor │ │ Metal ion │ │ Ligand │ └──────┬───────┘ │ ▼ 参数化 PARAMETER │ ┌──────────────┼──────────────┐ ▼ ▼ ▼ Protein FF Cofactor Ligand │ │ │ │ NAD/FAD/ATP │ │ │ │ │ ▼ │ │ Metal model │ │ │ │ └──────────────┼──────────────┘ ▼ topology.top │ ▼ Solvate │ ▼ Na / Cl- / ions │ ▼ EM │ ▼ NVT │ ▼ NPT │ ▼ Production MD │ ┌──────────────┼──────────────┐ ▼ ▼ ▼ RMSD RMSF H-bond/Distance │ │ │ └──────────────┼──────────────┘ ▼ MM/PBSA / MM/GBSA可以把整个问题归纳成一句话辅酶、配体、金属离子真正进入 GROMACS 的关键不是“怎么加进去”而是先解决“它们的力场参数和拓扑如何描述”然后才能正确加入体系进行 EM → NVT → NPT → MD。NAD⁺/FAD、ATP-Mg²⁺、Zn²⁺金属酶这三类最好分别处理尤其是“蛋白 NAD/FAD 金属离子”参数化方法和普通小分子 LIG 会有明显区别。生活很好有你更好。