分子动力学--非标残基的处理二（配体）

作者，Evil Genius

这一篇我们继续分子动力学，关于配体分子的处理，那么分析之前，我们需要了解很多的背景知识，每一个组学想要学通都需要花费大量的时间。

问题一：分子动力学模拟为什么不能直接用autodock对接的结果进行模拟

简单来说，这是两套 “操作系统” ，底层逻辑截然不同：

• AutoDock 的输出（如PDBQT格式）：这是为快速筛选而设计的简化结构。它缺失了MD模拟所必需的力场参数，例如原子的电荷、键、键角等“专属说明书”，其原子类型也与MD力场不匹配。如果强行使用，常会导致模拟时配体从结合位点脱离这一经典问题，或者在pdb2gmx等步骤直接报错终止。

• 分子动力学模拟的输入：它需要一套能描绘整个分子物理运动的完整“说明书”（力场参数文件），涵盖了每个原子的电荷、质量，以及原子间的键、键角、二面角等所有信息。

问题2：配体分子文件用Gaussian处理的核心目的是什么

通过量子化学计算，从分子的三维坐标出发，获得其能量、稳定结构、电子性质等关键理化信息。

具体而言：

几何结构优化 (Opt)

找到分子在势能面上的最稳定构型（能量极小点）。

为什么重要？ 实验测得的或建模软件（如Open Babel）给出的初始结构往往带有不合理的键长、键角，优化后才能代表真实分子。

单点能计算 (SP)

在给定结构下计算体系的总能量及其电子相关性质（如前线分子轨道、电子密度、原子电荷等）。

为什么重要？ 能量是所有稳定性和反应趋势的根源，后续分析反应能垒、结合能等都依赖于它。

振动频率分析 (Freq)

验证优化得到的结构是否为真正的极小点（无虚频），同时获得红外光谱、拉曼光谱、热力学量（熵、焓、自由能）。

为什么重要？ 没有频率分析，你无法确认“优化”是否成功，也无法预测分子的光谱特征。

预测反应路径与过渡态 (IRC, QST2等)

寻找反应物和产物之间的能垒（活化能）以及过渡态结构。

为什么重要？ 解释化学反应如何发生，对比不同反应路径的快慢。

模拟外部环境影响

通过隐式溶剂模型（如 SMD、PCM）或更高级的显式溶剂模型，计算分子在溶液、蛋白质环境中的能量与性质。

为什么重要？ 大多数实际化学反应发生在溶液或生物环境中，气相模拟往往不准确。

计算激发态性质 (TD-DFT等)

预测紫外-可见吸收光谱、荧光光谱、激发态能量等。

为什么重要？ 解释光物理、光化学过程（如染料、发光材料、光催化剂）。

问题3、分子对接的时候我们已经添加了电荷，为什么分子动力学需要重新分配电荷？

分子对接时添加的电荷，确实不能直接用于Gaussian计算。

虽然都叫“电荷”，但两者在计算方法、物理意义和用途上有着本质的不同。理解这个区别，就能明白为什么需要Gaussian重新计算。

为什么对接用的“电荷”不能用？

➡️ 目的不同：精准模型 vs. 快速筛选

分子对接的首要任务是速度，它需要在短时间内处理上万个小分子，因此采用了大量经验性近似，其电荷就是这些近似中的一环。而Gaussian的目标是精度，它通过求解量子力学方程来得到最准确的电子结构。

➡️ “电荷”的定义不同：净电荷 vs. 原子电荷

Gaussian需要的“电荷” (Net Charge)：是指整个分子的净电荷数，比如是中性（0）、正一价（+1）还是负一价（-1）。它不关心电荷在原子间如何分布，只关注总和。

分子对接中的“电荷” (Partial Atomic Charge)：是指电荷在各个原子上的具体分配，比如一个水分子中氧原子带负电，两个氢原子带正电，这就是“原子电荷”。分子对接时添加的正是这种。

同时，Gaussian的输入文件（.gjf）还需要你指定“自旋多重度”。例如稳定有机分子通常是“0 1”，即中性分子（0）且为单重态（1），这代表所有电子都已配对，没有未成对电子。

电荷 和 自旋多重度 在 .gjf 文件中是不可省略的核心信息，写在分子坐标部分的最前面。

➡️ 精度要求不同：粗略估算 vs. 高精度拟合

以AutoDock最常用的Gasteiger电荷为例，它是基于原子电负性粗略估算的。在一些研究中甚至发现，使用更精确的量子化学（QC）电荷，对AutoDock的打分函数没有带来明显的统计改进。这说明AutoDock的力场本身就不需要特别精确的电荷，够用就行。

但Gaussian的计算需要原子电荷能精确反映分子的真实静电势（Electrostatic Potential, ESP）分布。通常，在获得Gaussian的高精度计算结果后，会采用RESP（约束静电势拟合） 方法或CHELPG等方式，拟合出一套能与HF/6-31G*等高水平理论计算相匹配的原子电荷，用于后续更精确的分子动力学模拟（如使用AMBER或CHARMM力场）。

简而言之：对接用的电荷（如Gasteiger）是“经验估算”，而Gaussian计算出的电荷是基于量子力学的“高精度拟合”。因此，直接用在对接中算好的部分电荷，是无法作为Gaussian的输入文件进行计算的。

问题4、那么如何对分子对接的结果进行分子动力学模拟？

要把 AutoDock 的结果“翻译”成 MD 软件能懂的语言，主要是两个步骤：

1. 准备蛋白的拓扑文件

首先要用 PyMOL 或 AutoDock Tools (ADT) 等软件将蛋白从对接的 PDBQT 文件还原为标准的 PDB 格式。之后，利用 MD 软件自带的工具（如GROMACS的pdb2gmx）或通过 CHARMM-GUI 平台，为该蛋白指定一个力场（如Amber、CHARMM）并生成模拟所需的拓扑文件和结构文件。

2. 准备配体的拓扑文件

GROMACS 无法自动处理配体，这是最需要“手动”操作的步骤。我们需要借助第三方工具为配体生成拓扑文件：

* Sobtop (强烈推荐)：这是由国内著名的计算化学家卢天（Sobereva）老师开发的免费工具。你只需将配体文件（mol2或pdb）拖入即可在交互式界面中完成操作，非常方便。

* ATB (Automated Topology Builder)：一个免费的在线服务器，上传配体的PDB或SMILES，它会返回对应的ITP/GRO文件，可能需要手动注册。

* CGenFF (CHARMM General Force Field)：服务于CHARMM力场的在线服务器，与CHARMM-GUI无缝衔接。

* PRODRG：一个免费的在线服务器，但准确性相对较低，生成的参数可能需要谨慎。

* ACPYPE (AnteChamber PYthon Parser interfacE)：一个功能强大的免费命令行工具，它充当桥梁，调用 Ambertools 中的 antechamber 来生成GROMACS格式的拓扑文件。

💎 总结

最直接的核心信息就是：AutoDock的输出文件不能直接用。但只要能遵循“提取分子、生成拓扑、合并提交”这一核心准则，并把为配体创建准确的拓扑文件作为关键步骤来处理，你就可以顺利地将对接的优秀结果延续到更精确的分子动力学模拟中了。

在文章分子动力学--非标残基的处理一（配体）演示了在win系统下的Gaussian计算，那么接下来来到我们的重点，用linux系统计算，首先我们需要了解一下昨天设置的参数

首先第一点：在GaussView的Calculate -> Gaussian Calculation Setup窗口中，这部分定义了Gaussian程序要对你的分子执行什么计算任务。这是所有计算的起点。

参数二、方法 (Method Tab) – 选择理论方法

这部分定义了如何近似求解描述你分子体系的薛定谔方程。

如何选择？（快速决策指南）

常规基态计算（优化、频率、能量） → 选择 Ground State + DFT（如 B3LYP）+ 适当基组（如 6-31G(d)）。

预测中小分子（<50 原子）的紫外-可见光谱 → TD-SCF (TD-DFT)，泛函推荐 ωB97X-D 或 CAM-B3LYP（含长程校正）。

预测大分子（>100 原子）的紫外-可见光谱 → ZINDO，快速定性。

TD-DFT 计算遇到严重问题（如大量虚激发或收敛失败） → 尝试 TDA。

需要接近实验精度（如高精度光谱归属） → 小分子用 EOM-CCSD，中等分子用 SAC-CI。

极早期方法（如 CIS） → 已过时，不推荐常规使用。

参数3、基组 (Basis Set Tab) – 定义计算精度

实战选择指南

• 常规初筛、大体系估算：STO-3G, 3-21G, 6-31G。速度最快，但精度有限。

• 有机小分子结构优化与频率分析：6-31G(d,p) / def2-SVP。常规DFT计算的黄金起点。

• 高精度计算：6-311+G(d,p) / cc-pVTZ / def2-TZVP。对能量和结构精度要求高的体系，计算量适中。

• 弱相互作用分析：aug-cc-pVTZ / 6-31+G(d)。包含弥散函数是必要的。

• 过渡金属、镧系锕系元素：必须使用赝势基组，如 LANL2DZ, SDD, def2-TZVP+ECP。有效处理内层电子相关性和相对论效应。

• 激发态计算：aug-cc-pVTZ / 6-31+G(d,p)。弥散函数对描述激发态电子云扩张至关重要。

参数4、溶剂模型 (Solvation Model) – 模拟真实环境

Solvation标签页设置，这部分内容让计算考虑溶剂环境的影响，通常在与溶液反应或生物体系模拟相关时是必须打开的。

参数5、附加计算与属性 (Additional Calculation & Properties)

这部分在Gaussian Calculation Setup窗口的Properties标签页中，用于控制程序输出除基础信息（如能量和轨道）外的更多分子属性。

参数6、资源命令 (Link 0 Tab) – 提交资源与系统控制

这部分在Gaussian Calculation Setup窗口的Link 0标签页中，定义了计算所需的硬件资源。

总结与关键建议

参数真的是非常多，理解起来还是有一定的难度的。

生活很好，有你更好。