摘要你要做共价抑制剂β-内酰胺、丙烯酰胺、环氧弹头的对接但 AutoDock Vina 本体不会成键商业 CovDock 又难调参数本文把 gnina 1.3 新增的共价对接从原理到命令一次讲清它不猜反应、只要你给出共价结合后的配体形式用 SMARTS 指认成键原子、OpenBabel 启发式摆位、再把残基-配体当作一个柔性残基只采样扭转角。你能在文中拿到 6 步流程表、6 个--covalent_*参数速查、README 官方示例吡唑-铁与真实 issue 案例doripenem-Ser的逐参数拆解外加 207 复合物基准成功率 36.2%–66.6%给出的选参策略和 8 条踩坑。照表填参数直接产出含共价键的复合物。关键字gnina共价对接分子对接SMARTS共价抑制剂AutoDock VinaCNN 打分OpenBabel相关教程与核心文献资源链接与本文关系gnina 项目仓库github.com/gnina/gnina本文命令与参数的官方出处GNINA 1.3 论文2025McNutt et al.J. Cheminformatics. 10.1186/s13321-025-00973-x共价对接流程与 207 复合物基准的原始描述GNINA 1.0 论文2021McNutt et al.J. Cheminformatics. 10.1186/s13321-021-00522-2CNN 打分与整体对接框架CNN 打分方法论文Ragoza et al.J. Chem. Inf. Model.2017. 10.1021/acs.jcim.6b007403D CNN 蛋白-配体打分的原始方法共价对接工具横评Scarpino et al.J. Chem. Inf. Model.2018. 10.1021/acs.jcim.8b00228GNINA 1.3 论文引用其共价基准gnina 官方 Colabcolab.research.google.com免安装在线试跑入口你项目里一旦出现弹头 残基侧链的组合普通对接就开始失效配体和受体之间是共价键而对接程序默认把两者当两个独立分子。gnina 1.3 把这件事收敛成了 6 个命令行参数下面从原理讲到命令。一、痛点Vina 默认不会成键共价抑制剂β-内酰胺打丝氨酸、丙烯酰胺打半胱氨酸、磺酰氟打赖氨酸的结合模式预测有两道坎成键几何弹头原子与残基侧链原子之间是真实化学键键长、键角都要合理不是距离近就行。反应问题成键常伴随化学修饰开环、离去基团脱落。程序要不要会预测反应市面上多数共价对接工具要么要求你预组装复合物要么内置反应规则。gnina 的选择很务实GNINA 1.3 论文原文不预设任何化学反应——它要求你提供共价结合后的配体结构含全部化学修饰例如环氧开环后的形式然后只负责把这颗已成键的配体摆对位置、采好构象。这和 CovDock 等主流工具的做法一致把反应化学留给你把几何和构象交给程序。二、共价对接的完整流程6 步走GNINA 1.3 论文给出的流程如下每一步都对应一个可调参数步骤做什么实现机制 / 参数1 指认原子对受体侧成键原子 配体侧成键原子--covalent_rec_atom链:编号:原子名 或 x,y,z--covalent_lig_atom_patternSMARTS标识单个原子2 确定摆位在受体原子上找合理的成键方向默认对受体原子先减氢用 OpenBabelGetNewBondVector启发式也可--covalent_lig_atom_position x,y,z手动指定3 搬运配体旋转 平移配体使成键原子到位、成键几何合理OpenBabelOBBuilder::Connect4 成键生成共价键--covalent_bond_order默认单键5 采样与优化残基-配体构建体作为一个柔性残基参与对接蒙特卡洛采样和能量最小化只处理内扭转角不做刚体平移/旋转沿用 gnina 标准 MC 最小化流程6 可选精修对整个残基-配体构建体做 UFF 优化进一步修匀成键几何会改变配体键长键角--covalent_optimize_lig三个容易被忽略的设计细节多个匹配原子全评估SMARTS 若匹配到多个配体原子所有配体原子 × 受体原子配对都会跑一遍输出 pose 数相应翻倍。手动锁位给了--covalent_lig_atom_position后再加--covalent_fix_lig_atom_position该原子在采样中不再移动——只定初始构型还是全程锁死是两种不同的实验设定。CNN 打分仍认得配体对接和打分时配体原子保持配体原子类型与受体原子分开编码CNN 网格输入不受成键影响。三、上手准备最小环境方式命令 / 入口说明预编译二进制推荐GitHub Releases 下载官方明确推荐WSL2 里也能跑DockerDocker Hubgnina镜像README 提供免编译源码编译cmake .. make需 CUDA ≥ 12.0 Boost/Eigen/protobuf 等cmake 会自动拉 libmolgrid 和 PyTorch版本门槛共价对接是gnina 1.32025新增功能跑gnina --version确认 ≥ 1.3。另一条 README 明示的硬坑OpenBabel3.1.1 及更早版本键级判断有错源码编译请用其指定的 OpenBabel 分支。四、命令实战README 官方示例拆解示例 1吡唑氮共价接到受体的铁原子上gnina -r rec.pdb.gz -l conformer.sdf.gz --autobox_ligand bindingsite.sdf.gz \ --covalent_rec_atom A:601:FE \ --covalent_lig_atom_pattern [$(n1nccc1)] \ -o output.sdf.gz逐参数看-r受体支持 .gz-l是共价形式的配体构象--autobox_ligand用结合位点分子定义搜索盒子——共价对接同样需要盒子约束非成键部分的构象空间A:601:FE A 链 601 号残基的 FE 原子金属配位场景SMARTS[$(n1nccc1)]是递归写法匹配吡唑环上成键的那个氮原子。示例 2手动摆位 UFF 精修gnina -r rec.pdb.gz -l conformer.sdf.gz --autobox_ligand bindingsite.sdf.gz \ --covalent_lig_atom_position -11.796,31.887,72.682 --covalent_optimize_lig \ --covalent_rec_atom A:601:FE \ --covalent_lig_atom_pattern [$(n1nccc1)] \ -o output.sdf.gz与示例 1 的差别只有两处不用 OpenBabel 启发式、改为手动给出成键原子的初始坐标成键后对残基-配体构建体做 UFF 优化。这是有晶体学先验信息时的推荐姿势——后文基准数据会说明它值多少成功率。真实案例issue #242doripenem 共价对接 β-内酰胺靶点PDB 5L2Fgnina -r 5l2f-Apo-Prep.pdb -l doripenem-prep.sdf --autobox_ligand doripenem-native.sdf \ --covalent_rec_atom A:80:OG \ --covalent_lig_atom_pattern [$(CO)] \ -o doripenem-cov.sdf配体侧指认 β-内酰胺羰基碳受体侧指认催化丝氨酸 OG。注意这条命令里的两个细节在踩坑一节展开第三段是原子名OG 不是残基名 SERSMARTS 匹配单个原子。把共价复合物整体导出输出 SDF 只有配体。要含共价键的蛋白-配体复合物加--out_flex cov_complex.pdb --full_flex_output柔性残基连同共价配体一起写出受体刚性部分用-r的 PDB 拼接即可。五、--covalent_*六个参数速查表参数必选作用--covalent_rec_atom是受体成键原子链:残基号:原子名或x,y,z--covalent_lig_atom_pattern是SMARTS 指认配体成键原子标识单个原子--covalent_lig_atom_position否手动给成键原子初始坐标替代启发式摆位--covalent_fix_lig_atom_position否配合上一项全程锁死该原子位置--covalent_bond_order否共价键键级默认 1--covalent_optimize_lig否成键后 UFF 优化残基-配体构建体六、207 复合物基准选参策略直接抄GNINA 1.3 论文用 Scarpino 等人的 207 个共价复合物重对接基准评估了这个功能结论对你选参数非常直接设定成功率生成构象 默认启发式摆位信息量最少36.2%实验构象 精确指定成键原子位置信息量最多66.6%对照不开共价对接CNN 打分27.5%对照不开共价对接Vina 打分15.8%三条选参结论先验信息很值钱从 36.2% 到 66.6% 之间就是摆位信息的价值。有共晶或同源模板时务必用--covalent_lig_atom_position给坐标甚至加--covalent_fix_lig_atom_position。共价设定本身就是最大杠杆最差的共价设定36.2%也远好于不开共价对接的最好成绩27.5%。只要化学上确定成键就该开--covalent_*。共价对接用 Vina 打分别用 CNN基准里 Vina 显著优于 CNN——原因是 CNN 训练集里没有共价复合物模型出了适用域。命令里给--cnn_scoring none用经验打分排序。论文原话点明这是把模型用在适用域之外的典型坑。整体上gnina 1.3 的共价对接与最先进工具相当但不超越换来的是极低的使用门槛。七、8 条实测踩坑--covalent_rec_atom第三段是原子名不是残基名。写A:80:SER必报Usage error: Could not find receptor atom丝氨酸成键原子要写A:80:OGissue #242 的真实报错。残基号用准备后受体的编号。加氢、去水、重编号之后编号可能挪位先在 PDB 里核对目标残基的链号和序号再填。SMARTS 只标识一个原子。递归写法[$(...)]最稳想精确到哪个原子就把环境写大如[$(n1nccc1)]限定吡唑环氮。多原子裸 SMARTS 会选错原子——issue #271 里磺酰氟弹头被连到了氟上就是 SMARTS 没锁到硫。配体必须是共价结合后的形式。弹头已反应完毕β-内酰胺开环、离去基团已脱gnina 不做反应预测。共价对接把 CNN 关掉--cnn_scoring none理由见第六节。多个 SMARTS 匹配会翻倍输出。pose 预算按匹配数放大--num_modes再用--min_rmsd_filter去冗余。导出复合物记得--out_flex--full_flex_output否则输出里只有配体。OpenBabel 版本≤ 3.1.1 键级判定有 bugREADME 明示源码编译用官方指定分支别用系统老版本。八、局限与下一步反应化学全靠你离去基团、开环、可逆弹头的可逆性都不在程序视野内错误的共价形式会被忠实地对接成错误的复合物。只处理一对成键原子双位点共价如双弹头不支持SMARTS 多匹配是逐个评估不是同时成键。受体仍以刚性为主成键残基并入柔性构建体但口袋其余部分默认刚性高度柔性口袋可搭配--flexres/--flexdistREADME 也提醒非必要不推荐。打分排序仍有误差即便用 Vina207 复合物基准的成功率说明共价复合物的构象排序比非共价更难建议产出 pose 后接 MM/GBSA 或短 MD 精排。下一步你该做什么挑一个自家体系有共晶最好先用默认启发式跑一遍 redock 算 RMSD再用实验坐标指定成键位置跑第二遍——两遍的差距就是你的体系里摆位信息的价值也顺手验证了整条流程。参考来源gnina GitHubhttps://github.com/gnina/gnina README共价对接示例与全部参数McNutt A, Li Y, Meli R, Aggarwal R, Koes D R.J. Cheminformatics, 2025. 10.1186/s13321-025-00973-x共价对接流程与基准McNutt A, et al.J. Cheminformatics, 2021. 10.1186/s13321-021-00522-2Ragoza M, et al.J. Chem. Inf. Model., 2017. 10.1021/acs.jcim.6b00740Scarpino A, Ferenczy G G, Keserű G M.J. Chem. Inf. Model., 2018. 10.1021/acs.jcim.8b00228gnina issue #242doripenem-5L2F 共价对接实例https://github.com/gnina/gnina/issues/242