Amber Mmgbsa 没有提供足够的上下文信息来进行准确的翻译。然而，如果我们假设 "Amber" 和 "Mmgbsa" 是特定的术语或名称，那么它们在中文中仍然可以保持原样，即 "Amber Mmgbsa"。如果您能提供更多的上下文信息或解释，这将有助于提供更准确的翻译。

Amber MM/GBSA — 单结构工作流概述 本技能描述了使用AmberTools / Amber计算单个静态结构的MM/GBSA（分子力学/广义Born表面面积）结合自由能的完整、生产就绪的工作流 —— 无需分子动力学采样。它设计用于两个主要用例： 快速姿势评分 —— 评估晶体结构或顶级对接姿势的结合能，而无需运行MD模拟。 机制分解 —— 获取标准分解（VDWAALS / EEL / EGB / ESURF / ΔG_bind），以解释哪些相互作用驱动或反对结合。 范围说明：本工作流计算单结构的类似于焓的MM/GBSA估计。它不包括构象熵、水介导相互作用或ensemble平均。对于多个配体的量化排名或与实验ΔG值比较的结果，强烈推荐使用完整的MD轨迹 + MMPBSA.py ensemble平均。

前提 软件工具 版本 目的 tleap AmberTools 20+ 系统构建和拓扑 antechamber AmberTools 20+ 配体参数化 parmchk2 AmberTools 20+ 缺失的GAFF2参数 cpptraj AmberTools 20+ 轨迹提取 MMPBSA.py AmberTools 20+ MM/GBSA计算 parmed 任意最近的拓扑检查和合并

设置环境： export AMBERHOME=/path/to/amber24 # 调整到您的安装 export PATH=$AMBERHOME/bin:$PATH PY=$AMBERHOME/miniconda/bin/python # 首选Python

输入文件 文件 描述 receptor.pdb 受体（蛋白质 ± DNA ± 金属）具有正确的残基名称 ligand.pdb 配体具有3D坐标、正确的键序和已知的净电荷

推荐的目录布局 project_mmgbsa/ ├── input/ │ ├── receptor.pdb │ └── ligand.pdb ├── prep/ # 清理结构 ├── build/ # 拓扑和坐标 ├── traj/ # 单帧轨迹 ├── mmgbsa/ # MMPBSA.py输入和输出 └── logs/ # 所有日志文件

步骤协议 步骤0 —— 验证输入结构 在任何计算之前，确认： 受体：标准Amber残基名称（蛋白质：标准20个氨基酸；DNA：DC/DA/DG/DT；水：WAT/HOH）。 移除缓冲盐和无参数的小分子。 配体：3D坐标存在；键序化学合理；净电荷已知。 复合体几何：受体和配体必须已经处于正确的结合姿势（晶体结构或顶级对接姿势）。 它们不能在单独的文件中随机对接而没有参考复合体。

步骤1 —— 清理受体PDB 重命名晶体水为Amber兼容的WAT： mkdir -p prep logs awk ' BEGIN{OFS=""} /^ATOM|^HETATM/ { res=substr($0,18,3) if (res=="HOH") { print substr($0,1,17),"WAT",substr($0,21) } else { print $0 } next } {print} ' input/receptor.pdb > prep/receptor_clean.pdb 是否保留晶体水？通常除非特定的水分子被结构验证为关键桥接分子，否则不保留晶体水。 如果保留，记录此决定并在所有比较系统中一致应用。

步骤2 —— 参数化配体（GAFF2） 2a. 生成GAFF2 mol2文件，带有AM1-BCC电荷 cd prep antechamber \ -i ../input/ligand.pdb \ -fi pdb \ -o ligand.mol2 \ -fo mol2 \ -c bcc \ -nc 0 \ -at gaff2 \ -j 4 \

../logs/antechamber.log 2>&1

# 验证mol2文件是否生成 ls -lh ligand.mol2 关键：-nc值必须与配体的真实形式电荷相匹配。 常见示例：中性有机分子 → -nc 0 羧酸盐 → -nc -1 磷酸盐 → -nc -2或-3 金属复合物 → 单独确认电荷

2b. 检查缺失的力场参数 parmchk2 \ -i ligand.mol2 \ -f mol2 \ -o ligand.frcmod \

../logs/parmchk2.log 2>&1

# 如果"frcmod"文件为空或缺失的参数是关键的， # 可以： # a) 在PDB中重新绘制问题子结构 # b) 手动将参数添加到frcmod文件中

步骤3 —— 构建受体拓扑 仅蛋白质 mkdir -p build cat > build/tleap_rec.in <<'EOF' source leaprc.protein.ff14SB source leaprc.water.tip3p rec = loadpdb ../prep/receptor_clean.pdb desc rec saveamberparm rec receptor.prmtop receptor.inpcrd quit EOF tleap -f build/tleap_rec.in > logs/tleap_rec.log 2>&1

蛋白质 + DNA cat > build/tleap_rec.in <<'EOF' source leaprc.protein.ff14SB source leaprc.DNA.OL15 source leaprc.water.tip3p rec = loadpdb ../prep/receptor_clean.pdb desc rec saveamberparm rec receptor.prmtop receptor.inpcrd quit EOF tleap -f build/tleap_rec.in > logs/tleap_rec.log 2>&1

蛋白质 + DNA + 金属离子或特殊基团 cat > build/tleap_rec.in <<'EOF' source leaprc.protein.ff14SB source leaprc.DNA.OL15 loadoff $AMBERHOME/dat/leap/lib/terminal_monophosphate.lib source leaprc.water.tip3p rec = loadpdb ../prep/receptor_clean.pdb desc rec saveamberparm rec receptor.prmtop receptor.inpcrd quit EOF tleap -f build/tleap_rec.in > logs/tleap_rec.log 2>&1

验证成功： ls -lh build/receptor.prmtop build/receptor.inpcrd