刘建川・GROMACS 实战 | 如何处理二硫键

大家好，我是刘建川，14 年深耕 GROMACS 分子模拟领域，始终专注分子模拟实战。在使用pdb2gmx为蛋白质、多肽等生物体系生成拓扑文件时，体系中的二硫键无法被程序默认识别，必须单独进行配置。本文将介绍二硫键的基础作用，同时详细讲解 GROMACS 中二硫键的添加方法与核验技巧。

一、二硫键基础介绍

二硫键是一种稳定性较强的共价键，用于连接同一肽链或不同肽链上两个半胱氨酸残基的巯基。它是蛋白质二级结构与三级结构的重要组成部分，能够牢牢稳固肽链的空间构象。一般而言，蛋白质体系内的二硫键数量越多，其结构抵抗外界影响的能力也就越强。

二、pdb2gmx 交互式添加二硫键

想要在生成拓扑文件时构建二硫键，需要在pdb2gmx命令中添加 -ss 参数，开启交互式选择功能，完整执行命令如下：

gmx pdb2gmx -f O.pdb -ff amber99s b -ignh -o O.gro -ss

运行命令后，程序会自动检索体系内潜在的二硫键位点并弹出交互提示。若确认当前两个残基需要构建二硫键，输入 y 并回车即可。如果体系包含多个二硫键，程序会逐个进行询问，按照实际需求依次选择即可。

三、核验二硫键配置是否正确

仅含单个二硫键的体系，通常可以直接生成成功，无需额外检查。若体系存在多个二硫键，或是中途手动跳过了部分位点，建议人工核验二硫键的添加位置是否准确。

首先使用 GaussView、VMD 等常用可视化软件，打开pdb2gmx输出的.gro结构文件，定位形成二硫键的两个硫原子，记下对应的原子编号。随后打开配套的top或itp拓扑文件，查看[ bonds ]区块，确认上述两个硫原子编号已正常录入，即可判定二硫键配置无误。

按照这套流程完成二硫键的配置与核验后，就能保障蛋白质、多肽类生物体系初始结构的完整性，让后续的分子动力学模拟可以稳定开展。