别再只数氢键个数！HBCalculator 可视化氢键密度与强度空间分布

3月前浏览165

大家好，我是刘建川，14 年深耕 GROMACS 分子模拟领域，始终专注分子模拟实战。我将拆解分子模拟核心原理，帮你理清技术脉络、避开选择误区，真正做到学以致用。

在分子动力学模拟中，氢键的分析几乎是一项标配任务。但长期以来，主流工具（如VMD内置命令、MDAnalysis等）主要提供氢键的数量统计，却无法直接回答一个关键问题：氢键到底分布在空间的哪些位置？密度和强度如何变化？这一局限，让许多涉及界面、纳米限域、生物分子识别的研究难以深入。2024年发表在 Journal of Chemical Information and Modeling 上的一项工作，带来了一款专门解决这一痛点的工具——HBCalculator。

论文链接：https://pubs.acs.org/doi/10.1021/acs.jcim.4c00054

一、原理：从“数个数”到“看分布”

HBCalculator 的核心创新在于：它不仅判断氢键是否形成，还能将每个氢键的位置信息投影到一维或二维坐标轴上，并通过空间分箱（binning）统计每个区域内的氢键数量或强度。如论文中的图1 所示，代码包含六个清晰的流程步骤：

加载文件 → 设置盒子与分箱 → 输出初始化 → 氢键计算循环 → 投影与分箱 → 数据精炼与输出。整个工具基于 Tcl 脚本语言编写，运行在 VMD 环境中，充分利用了 VMD 强大的原子选择与轨迹处理能力。

二、功能：1D/2D分布 + 强度可加权

HBCalculator 的核心功能可归纳为三个层次：

1、一维分布计算

沿 x, y, z任一方向，输出氢键数密度或强度的曲线分布。

2、二维分布计算

在 xy, yz, xz 平面上生成氢键分布的彩色等高线图，直观揭示空间异质性。

3、氢键强度的加权处理

支持三种模式：

mode = 1：仅统计氢键数量密度（单位：Å⁻³）

mode = 2：使用预定义的供体-受体距离加权函数

mode = 3：允许用户自定义加权函数，灵活性极高

论文中明确指出：“用户可以获取氢键中的供体-受体距离、氢-受体距离以及角度等变量，用于自定义强度计算。”

三、使用方法：一条命令完成全部计算

HBCalculator 的使用方式非常直接，完全依赖一条位置严格的命令：

hbc file_1 file_2 projection firstFrame lastFrame step cell atomSelection cutoff angle mode

论文中给出了一个真实示例：

hbc gra-water.psf main.dcd xy 551 600 1 {0 0 0} {36 81 71} {36 81 1} {water} 4.0 30.0 2

这条命令的意思是：加载 gra-water.psf 和 main.dcd在 xy 平面上做二维投影，分析第 551 到 600 帧（共 50 帧），盒子范围：x∈[0,36],y∈[0,81],z∈[0,71]，分箱数：36×81×1（即在z 方向不做分割）仅分析水分子，氢键判据：距离 < 4.0 Å，角度 < 30°，使用 mode 2 计算强度分布。（如果不清楚氢键定义，可以查看我们的推文：刘建川・GROMACS 实战｜氢键核心原理与氢键计算全解析）。

所有参数都必须给出，没有默认值，这一点在论文表1中有清晰说明。