中科科翼 专注模拟

结合自由能的定义

结合自由能（Binding Free Energy）是指在一定温度和压力下，两个或多个分子（如配体与受体）结合形成复合物时，自由能的变化量。它是衡量分子间相互作用强度和稳定性的重要热力学参数，负的结合自由能值表示结合过程是自发的，且绝对值越大，结合越稳定。结合自由能在药物研发、蛋白质 - 配体相互作用、酶催化等领域具有关键作用，例如用于评估药物分子与靶点的结合能力，筛选潜在的药物候选物等。

· 基本原理：通过计算机模拟分子体系在不同状态下的能量变化，利用统计力学或热力学原理，计算出结合自由能。这些方法通常基于分子动力学模拟或其他模拟技术，通过改变体系的状态（如将配体从溶液中转移到受体的结合位点），并计算相应的能量变化来得到结合自由能。

计算模拟方法：

自由能微扰（FEP）：如上述原理所述，是一种较为精确但计算量非常大的方法，适用于对结合自由能计算精度要求较高的情况，尤其是对于小分子配体与蛋白质受体的结合研究。

热力学积分（TI）：与 FEP 类似，也是一种精确的计算方法，但在某些情况下计算效率可能略高于 FEP，同样适用于研究分子间的结合自由能。

分子力学/ 泊松 - 玻尔兹曼表面积（MM/PBSA）：计算相对简单，计算量较小，适用于大规模的分子对接和虚拟筛选后的结合自由能估算，但精度相对 FEP 和 TI 较低。

线性相互作用能（LIE）：计算速度较快，基于经验公式，在一定程度上能够快速估算结合自由能，适用于初步筛选和快速评估，但准确性依赖于经验参数的选择和体系的适用性。

常用软件：

GROMACS：广泛用于生物分子模拟，包括结合自由能的计算。例如，通过 g_mmpbsa 等工具可以进行 MM/PBSA 计算，也可以通过一些插件或用户自定义脚本实现 FEP 等方法的计算。具有高效的并行计算能力，能够处理大规模的分子体系。

AMBER：包含一系列用于分子动力学模拟、自由能计算等的程序和工具。其力场在生物分子模拟中被广泛认可和使用，提供了丰富的功能和教程，方便用户进行结合自由能的计算，如通过 pmx 等模块进行 FEP 计算。

Schrödinger Suite：包含多个模块，如 Prime、FEP + 等，可用于蛋白质结构预测、分子对接、自由能计算等。FEP + 模块提供了高效的自由能微扰计算功能，适用于高精度的结合自由能计算，尤其在药物研发领域应用广泛。

NAMD：适用于生物大分子体系的分子动力学模拟，结合一些后处理工具和脚本，可以进行 MM/PBSA 等结合自由能计算，具有良好的可视化界面和用户支持。