在进行分子动力学（MD）模拟前，正确处理蛋白结构中的二硫键（disulfide bond）是确保模拟稳定性和结果可靠性的关键步骤。如果二硫键未被正确识别和建模，会导致力场参数错误、结构不稳定甚至模拟崩溃。AMBER 等主流力场（如 ff14SB、ff19SB）不会完全自动识别 PDB 文件中的二硫键， 如果遇到pdb4amber未能正确识别二硫键，必须通过人工干预显式定义共价连接。以下是针对AMBER 力场体系的处理流程。✅ 一、核心原则在 Amber（特别是 tleap）中，如果自动工具（如 pdb4amber）未能正确识别二硫键，你可以显式地手动定义二硫键。原因：1. AMBER 不会自动识别 PDB 中的 S–S 距离（2.2 Å）并创建共价键。2. 某些 PDB 文件会包含 SSBOND 记录，明确标出哪些半胱氨酸（CYS）残基之间形成了二硫键；有些 PDB 结构可能缺失 SSBOND 行，但实际存在二硫键，需手动确认 SG–SG 距离必须作为共价键显式建模。否则：两个半胱氨酸（CYS）被视为独立残基；SG 原子无化学键约束 → 自由旋转或断裂；蛋白三级结构在模拟初期迅速崩塌；GB/PB 模型中静电/范德华参数错乱，自由能计算（如 MM/GBSA）结果严重失真。✅ 二、标准处理流程方法一：首先尝试使用pdb4amber + tleap处理体系若在后续的MD阶段失败（体系崩溃或者难以稳定），可使用下面方法二。本例以PDB结构6xie.pdb为例，去除结构中的溶剂和配体分子，将蛋白保存为6xie_protein.pdb。从PDB文件可以看到，有3个二硫键，存在于以下6个CYS之间：... SSBOND 1 CYS B 223 CYS B 255 1555 1555 2.05 SSBOND 2 CYS B 323 CYS B 358 1555 1555 2.05 SSBOND 3 CYS B 375 CYS B 378 1555 1555 2.04 ...PDB结构中二硫键记录清楚，我们可以直接使用pdb4amber(--nohyd 或者 -y删除氢原子，让tleap根据力场自动加载;--no-conect 去掉硫原子连接，让自动处理)：pdb4amber -i 6xie_protein.pdb -o cleaned_6xie_protein.pdb --nohyd --no-conect⚠️ 输入 PDB 中的半胱氨酸必须保留为CYS（不是CYX），SG 原子名必须为SG，且 SG–SG 距离 2.2 Å。AmberTools 提供了 pdb4amber 工具可自动检测二硫键（基于 SG 原子间距），将参与二硫键的 CYS 残基重命名为 CYX，输出pdb文件对残基从1开始编号。可继续进行下一步tleap，例如：## leap.in ## 加载力场，原子半径参数 source leaprc.protein.ff19SB # 推荐蛋白力场（替代旧 ff19SB） source leaprc.water.tip3p # 水模型 set def Amber分子动力学模拟6: 蛋白结构准备中二硫键的处理
在进行分子动力学（MD）模拟前，正确处理蛋白结构中的二硫键（disulfide bond）是确保模拟稳定性和结果可靠性的关键步骤。如果二硫键未被正确识别和建模，会导致力场参数错误、结构不稳定甚至模拟崩溃。AMBER 等主流力场（如 ff14SB、ff19SB）不会完全自动识别 PDB 文件中的二硫键， 如果遇到pdb4amber未能正确识别二硫键，必须通过人工干预显式定义共价连接。以下是针对AMBER 力场体系的处理流程。✅ 一、核心原则在 Amber（特别是 tleap）中，如果自动工具（如 pdb4amber）未能正确识别二硫键，你可以显式地手动定义二硫键。原因：1. AMBER 不会自动识别 PDB 中的 S–S 距离（2.2 Å）并创建共价键。2. 某些 PDB 文件会包含 SSBOND 记录，明确标出哪些半胱氨酸（CYS）残基之间形成了二硫键；有些 PDB 结构可能缺失 SSBOND 行，但实际存在二硫键，需手动确认 SG–SG 距离必须作为共价键显式建模。否则：两个半胱氨酸（CYS）被视为独立残基；SG 原子无化学键约束 → 自由旋转或断裂；蛋白三级结构在模拟初期迅速崩塌；GB/PB 模型中静电/范德华参数错乱，自由能计算（如 MM/GBSA）结果严重失真。✅ 二、标准处理流程方法一：首先尝试使用pdb4amber + tleap处理体系若在后续的MD阶段失败（体系崩溃或者难以稳定），可使用下面方法二。本例以PDB结构6xie.pdb为例，去除结构中的溶剂和配体分子，将蛋白保存为6xie_protein.pdb。从PDB文件可以看到，有3个二硫键，存在于以下6个CYS之间：... SSBOND 1 CYS B 223 CYS B 255 1555 1555 2.05 SSBOND 2 CYS B 323 CYS B 358 1555 1555 2.05 SSBOND 3 CYS B 375 CYS B 378 1555 1555 2.04 ...PDB结构中二硫键记录清楚，我们可以直接使用pdb4amber(--nohyd 或者 -y删除氢原子，让tleap根据力场自动加载;--no-conect 去掉硫原子连接，让自动处理)：pdb4amber -i 6xie_protein.pdb -o cleaned_6xie_protein.pdb --nohyd --no-conect⚠️ 输入 PDB 中的半胱氨酸必须保留为CYS（不是CYX），SG 原子名必须为SG，且 SG–SG 距离 2.2 Å。AmberTools 提供了 pdb4amber 工具可自动检测二硫键（基于 SG 原子间距），将参与二硫键的 CYS 残基重命名为 CYX，输出pdb文件对残基从1开始编号。可继续进行下一步tleap，例如：## leap.in ## 加载力场，原子半径参数 source leaprc.protein.ff19SB # 推荐蛋白力场（替代旧 ff19SB） source leaprc.water.tip3p # 水模型 set def