作者，Evil Genius
我们2026的课程，这就算结束了，包括两部分。
针对空间转录组的在2026年空间转录组系列课程。
针对分子对接分子动力学的在2026生化小课--CADD计算机辅助药物设计。
每个组学都需要学习很多的内容，大家如果选择学习，就一定要学会，不学就好好放松一下，分析的时候该找公司找公司，只能说自己学会，性价比最高，学会了之后积累经验即可，请公司做，那就只能一直请公司。
虽然现在都有Ai了，但是Ai的过程与结果还是需要人来判断，内部逻辑参数如果不理解，最终的结果还是不行，改都不知道怎么改。
这一篇我们来补充关于分子动力学模拟的过程中添加辅酶、金属离子等内容，课程中我们学习了如何添加配体分子，我们来添加一些其他必需的部分。
对于 GROMACS 分子动力学（MD），如果蛋白体系中存在辅酶、金属离子、底物、配体等非标准组分，不能简单地把它们当成普通水分子加入 solv.gro。正确做法是：结构 → 参数 → 拓扑 → 组装 → 溶剂/离子 → 能量最小化 → 平衡 → 生产模拟。
1. 整体流程
例如一个蛋白中有 Zn²⁺ + NAD⁺ + 底物：
蛋白 Protein
│
┌─────────────────┼─────────────────┐
│ │ │
Zn²⁺ NAD⁺ Ligand
│ │ │
└─────────────────┴─────────────────┘
│
▼
建立复合物结构
│
▼
为非标准组分准备参数
│
┌──────────────┼──────────────┐
▼ ▼ ▼
金属参数 辅酶参数 配体参数
│ │ │
└──────────────┼──────────────┘
▼
topology.top
│
▼
加水 TIP3P/SPC
│
▼
添加 Na⁺ / Cl⁻
│
▼
能量最小化
│
▼
NVT
│
▼
NPT
│
▼
Production MD
2. 最重要的一点：先判断它是什么
不同非标准组分，处理方式完全不同。
| 组分 | 典型例子 | 参数化难度 | 是否需要单独拓扑 |
|---|---|---|---|
| 普通离子 | 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
需要定义：
bond
angle
dihedral
atom type
charge
这种方法可以维持特定配位结构，但它引入了人为约束。
方法 C：更高级的金属模型
例如针对特定金属中心采用：
bonded model
nonbonded model
dummy atom model
12-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 10
7. 一个很关键的问题：辅酶属于“配体”吗？
从 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.mdp
9. 拓扑文件的关系
最终：
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 topology
topol.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/PBSA
11. 特别推荐：保留实验结构中的辅酶和金属
如果蛋白是从 PDB 晶体结构得到的，例如：
Protein + NAD + Mg
通常不要简单地：
remove HETATM
然后重新“猜”辅酶和金属位置。
因为晶体结构已经给出了一个实验观察到的：
Protein
↓
Cofactor
↓
Metal
↓
Water
空间关系。
比较合理的是：
实验 PDB
│
├── Protein
├── Cofactor
├── Metal
└── Key waters
│
▼
参数化/检查
│
▼
GROMACS
12. 重点掌握这张图
实验结构 / 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 会有明显区别。
