4. 溶剂化模型与 PB/GB 模型归纳总结
做 MM/GBSA 自由能计算，溶剂化模型是绕不开的坎。说白了，你模拟的生物分子泡在水里，水分子怎么处理，直接影响最终的能量值。我入行那会儿，觉得用水盒子（显式溶剂）才真实，结果一跑计算，噪声大得离谱，后来才明白——不同问题要用不同模型。
这一章，我就带你捋清楚溶剂化模型的大家族，重点把 PB（Poisson-Boltzmann）和 GB（Generalized Born）这两个隐式模型的来龙去脉、优缺点、选型依据说透。最后我会画张图，帮你把知识点串起来。
4.1 溶剂化模型全家福
溶剂化模型分两大派：显式和隐式。
| 模型类型 | 代表 | 精度 | 速度 | 适用范围 |
|---|---|---|---|---|
| 显式溶剂 | TIP3P、SPC/E | 高（含微观结构） | 慢（需大量水分子） | 分子动力学、水结构研究 |
| 隐式溶剂（连续介质） | PB、GB | 中（忽略水分子细节） | 快（无溶剂粒子） | 自由能计算、结合预测 |
显式溶剂模型把水分子一个个建出来，最真实，但计算量巨大。我曾经跑过一个 10 ns 的 MD，光水分子就占了 6 万个原子，服务器直接报警。后来我改用隐式模型，同样的体系几小时出结果，虽然丢了一些水桥细节，但能量趋势完全对得上。
隐式模型的核心思想：把水看作均匀连续介质，溶剂对溶质的效应用一个「溶剂化能」Gsolv 来近似。而这个 Gsolv，通常分成极性和非极性两项：
G_solv = G_polar + G_nonpolar非极性项简单，一般用表面积拟合；极性的静电贡献，就是 PB 和 GB 模型的主战场了。
4.2 PB 模型——精确但昂贵
PB 模型基于经典静电学，直接求解泊松-玻尔兹曼方程：
∇·[ε(r) ∇φ(r)] − κ² ε(r) φ(r) = −4π ρ(r)这个方程描述了带电溶质在含盐电解质溶液中的电势分布。求解它，你要给体系划分网格、迭代数值计算。你想想看，对于一个包含几千个原子的蛋白质，网格动辄几百万格点。
PB 的优点很明显：它从第一性原理出发，对任意形状的溶质都能给出相当准确的静电势。我在做蛋白-配体结合自由能计算时，PB 的结果和实验值往往更接近，特别是有盐桥或强电荷相互作用的地方。
但缺点也突出：
- 慢：一次 PB 计算可能花几分钟到几小时，而且难以并行加速。
- 依赖参数：介电常数、盐浓度、网格间距、边界条件……每个参数都可能影响结果。我见过有人把网格间距从 0.5 Å 换成 0.25 Å，能量差出 5 kcal/mol。
- 对体系形状敏感：内部空腔、凹槽处理不好，会出现人工高能量。
4.3 GB 模型——快糙猛
GB 模型是 PB 的近似。它不求解偏微分方程，而是用一个解析表达式直接估算静电溶剂化能。核心是所谓「广义 Born 公式」：
ΔG_polar ≈ −(1/2) · (1/ε_in − 1/ε_out) · Σ_i Σ_j q_i q_j / f_GB这里的 fGB 是一个基于有效 Born 半径的函数。每个原子的有效半径怎么算？GB 模型内部有一套算法（比如利用分子体积积分），不同的参数化方案（igb=1, 2, 5, 8…）给出不同结果。
GB 的好处：
- 极快：通常在毫秒级，比 PB 快两到三个数量级。
- 容易并行、梯度清晰：很适合做分子动力学模拟中的隐式溶剂。
- 数值稳定：没有网格和边界问题。
但代价是精度不稳定。我遇到过一种体系，GB 给出的溶剂化能和 PB 差 3 kcal/mol，根源是配体表面有一个深沟槽，GB 无法准确算出里面的水排斥效应。
| 对比维度 | PB | GB |
|---|---|---|
| 理论基础 | 解 Poisson-Boltzmann 方程 | 解析近似 Born 公式 |
| 计算耗时 | 几分钟～几小时 | 毫秒～秒 |
| 对复杂形状处理 | 优秀 | 较差（依赖 Born 半径） |
| 参数敏感度 | 高（网格、边界） | 中（igb 变体、盐浓度） |
| 标准 MM/GBSA 常用 | 少用（太慢） | 常用（igb=5/8） |
| 梯度（用于MD） | 困难（数值梯度） | 光滑解析梯度 |
4.4 如何选？我的经验法则
你在实际项目中，选 PB 还是 GB？我一般这样判断：
- 排筛和粗筛：用 GB 最快，igb=5（OBC 模型）加 saltcon=0.1，能快速跑几百个配体的结合能排序。我曾在两周内用 GB 筛掉 80% 的候选分子，精度足够。
- 精修和终局：对最后几个候选再做 PB 精细计算。这时慢一点可以接受，但你得小心网格参数。我建议在 PB 计算前，先用
pbsa(1)或AMBER Tools的pbsa.py做一轮收敛性测试。 - 带电配体 vs 中性配体：当配体带强电荷，或者结合界面有盐桥，PB 明显更可靠。中性疏水分子用 GB 完全够用。
一句话总结：想要速度且趋势对，用 GB；想要绝对数值且舍得等，用 PB。大多数论文里的 MM/GBSA 其实就是 GB，因为 PB 计算成本实在太高了。
4.5 避坑指南（我曾踩过的坑）
我曾经有一次，用 GB 算出一个结合能 -45 kcal/mol，对照实验值只有 -8 kcal/mol。查了两天，发现是 igb 参数设错了——默认 igb=1 对某些力场不兼容。后来我核对规范，发现 igb=2 或 igb=5 才是推荐值。
还有一个坑：PB 计算中，如果你设置 inp=2（内部/外部介电常数） 为 1/80，但对配体内部空腔不处理，能量会严重高估。一定要在输入文件里开空腔修正：
&pb
istrng=0.1,
inp=2,
cavity=1,
solve_order=2,
/另一个常见问题是盐浓度单位。PB 和 GB 的 saltcon 单位是 mol/L，但有些老代码里是 mM，换算搞错的话，静电项能差出几十。我习惯写注释：
# saltcon 单位 M（mol/L）
saltcon=0.15 # 模拟生理浓度4.6 溶剂化模型决策图谱
下面这张 SVG 图帮你把整个知识体系串起来。对照着看，哪个环节该用哪种模型，一目了然。