欢迎关注我的博客Blockbuster-drug 的CSDN 博客主页专栏推荐《多肽性质预测模型实践》《开源蛋白结构预测》《蛋白生成》《开源多肽设计模型部署》《Amber分子动力学系列》摘要本文围绕 Amber 分子动力学模拟的六大应用场景展开逐一拆解蛋白–配体结合自由能、构象稳定性、蛋白–蛋白相互作用、膜蛋白、核酸与 IDP 等场景的核心问题、模拟要点、常见坑位与代码写法并补充增强采样、位置约束与 SHAKE 三个关键时机的判断依据帮助读者在开跑前选对方法、定好收敛指标。关键字AMBER、分子动力学、应用场景、MM/GBSA、结合自由能、增强采样、位置约束、SHAKE前面几篇把单一体系蛋白-配体 13.1/13.2、蛋白-蛋白 9、蛋白-核酸 10、膜蛋白 11的完整操作走了一遍。这篇换一个视角先选场景再谈操作。实际项目里最常见的翻车不是某条命令写错而是场景和方法不匹配——用 10 ns 轨迹算 MM/GBSA 就敢给 ΔΔG 排序、拿隐式溶剂跑 PPI、用 ntp1 各向同性压浴跑膜体系。本文把 Amber 的六大应用场景逐个拆开这个场景回答什么问题、模拟要点是什么、坑在哪、代码怎么写让你在开跑之前就知道终点长什么样。相关教程与核心文献官方教程教程内容与本文关系AmberTools 官方教程索引全部官方教程入口各场景章节的命令出处Tutorial 3MM-PBSA雌激素受体–雷洛昔芬复合物 MM-PBSA Nmode 熵§二、§四自由能计算的官方参考Amber 官方模型页蛋白/核酸力场与推荐水模型§三力场-水模型配对的权威依据Amber24 手册 §21.6sander 输入 §24.3GaMDntr/ntc/ntf 与 igamd 参数的原始定义§六/§七/§八参数表的出处核心文献文献为什么值得先读Tian C 等,JCTC16, 528 (2020), DOI 10.1021/acs.jctc.9b00591ff19SB 蛋白力场原始论文——§三力场选择的依据Case D 等,JCIM63, 6183 (2023), DOI 10.1021/acs.jcim.3c01153AmberTools 2023 综述——工具全貌与推荐用法Zgarbová M 等,JCTC11, 5723 (2016), DOI 10.1021/acs.jctc.5b00716OL15 DNA 力场 β 二面角精修——§五核酸场景OL15 ε/ζ OL1χOL4β OL1仅 DNARNA 配 OL3Dickson C 等,JCTC18, 1726 (2022), DOI 10.1021/acs.jctc.1c01217Lipid21 膜力场——§五膜蛋白场景Miao Y 等,JCTC11, 3584 (2015), DOI 10.1021/acs.jctc.5b00436GaMD 增强采样原始论文——§四/§五 IDP 场景Wang J, Miao Y,JCTC18, 1275 (2022), DOI 10.1021/acs.jctc.1c00974PPI-GaMD2 μs 捕获 barnase–barstar 解离/再结合——§四 PPI 采样上限的参照Robustelli P 等,PNAS115 (2018), DOI 10.1073/pnas.1800690115a99SB-disp 力场折叠无序蛋白兼顾——§五 IDP 场景Love O, Winkler L, Cheatham TE,JCTC20, 625 (2023), DOI 10.1021/acs.jctc.3c01164dsDNA 推荐 OL21OPC 的 vdW 参数扫描——§五核酸场景OL21 推荐依据一、先选场景再动手Amber 的应用场景可以按你要回答什么问题分成六类场景回答什么问题核心方法时长量级系列姊妹篇1. 蛋白–配体结合自由能这批分子谁结合更强ΔΔG 多少MM/GBSA 粗排 → TI/FEP 精算MM/GBSA 5 ns × 6-8 replicatesTI 每窗口 5 nsAmber分子动力学模拟13.1: MD模拟要点汇总一/Amber分子动力学模拟13.2: MD模拟要点汇总二2. 蛋白构象稳定性突变/温度/配体是否破坏折叠RMSD/RMSF/Rg/DSSP 收敛分析100 ns 起柔性区 μs本文 §三3. 蛋白–蛋白相互作用界面哪些残基是热点结合多强界面接触/氢键分析 MM/GBSA 分解100 ns–1 μsAmber分子动力学模拟9: 蛋白-蛋白相互作用模拟操作及亲和力计算示例4. 膜蛋白跨膜信号传导、通道门控Lipid21 半各向异性压浴500 ns 起Amber分子动力学模拟11: 膜蛋白动力学模拟操作5. 核酸/核酸–药物DNA/RNA 构象、小分子嵌入模式OL15/bsc1 离子环境控制100 ns–μsAmber分子动力学模拟10: 蛋白-核酸复合物模拟操作6. IDP/相分离无序蛋白的构象系综a99SB-disp / GaMD / REMDμs 级或增强采样本文 §五选场景的顺序是问题 → 方法 → 成本预估 → 验证指标。最后一步最重要——开跑前先定好什么结果算收敛否则跑完只会在轨迹里找自己想看的东西。二、场景一蛋白–配体结合自由能2.1 场景说明药物设计里最常用的 Amber 应用对虚拟筛选命中的分子、SAR 系列类似物做亲和力排序或解释为什么这个分子比那个强。结构准备、加氢、力场配对、电荷计算在Amber分子动力学模拟13.1: MD模拟要点汇总一与Amber分子动力学模拟13.2: MD模拟要点汇总二已详述这里只讲方法分层MM/GBSA / MM/PBSA快粗一条复合物轨迹 三个拓扑就能出 ΔG适合 10–100 个分子的相对排序。精度预期 ±2 kcal/mol只能排序别当绝对结合能。TI / FEP慢准alchemical 路径逐窗口消耦合配体ΔΔG 精度可达 1 kcal/mol 以内适合 lead optimization 后期关键决策。成本是 MM/GBSA 的 10–100 倍λ 窗口 ×2 相 × 重复次数叠加。2.2 模拟要点要点说明系综平衡与生产都用NPT300 K, 1 atm。NVT 只用于加热阶段生产时长MM/GBSA5 ns × 6-8 replicates不同初始速度详见Amber分子动力学模拟13.2: MD模拟要点汇总二§5.1单条 100 ns 轨迹并不优于多 replicate 短轨迹Hou 2011JCIM51:69 Pushkaran Arabi 2025IJBM306:141408 bootstrap 共识力场组合蛋白 ff14SB/ff19SB 配体 GAFF2 AM1-BCC 电荷水模型跟着蛋白力场走ff19SB→OPCff14SB→TIP3P 均可别把 OPC 配 ff14SB 还当成升级轨迹格式NetCDFioutfm1, ntxo2帧间隔 10–50 psMM/GBSA 参数igb5GBn配 mbondi2 半径最常用igb8GBn2需 mbondi3 半径切换前确认拓扑半径集TI 设置每窗口 ≥5 ns12–21 个 λvdW 消耦合窗口在两端加密软核势ifsc1两相溶液复合物都要跑与Amber分子动力学模拟13.2: MD模拟要点汇总二§5.1 口径一致2.3 避免踩坑❌ 踩坑✅ 正确做法用 10 ns 轨迹 MM/GBSA 就报 ΔΔG5 ns × 6-8 replicates 起详见Amber分子动力学模拟13.2: MD模拟要点汇总二§5.1报 ΔΔG 时给均值 ± SDΔΔG 差距小于 2 倍 SD 时明说不可分辨只跑复合物轨迹做 MM/GBSA 单轨迹近似就下结论单轨迹近似仅适合结构高度相似的类似物系列跨骨架比较用独立轨迹法把 MM/GBSA 绝对值与实验 ΔG 对标只对比相对排序或做线性回归校正TI 忘开软核势直接消耦合ifsc1scmask1/scmask2否则窗口两端能量发散体系净电荷不为零就开 0.15 M 盐先加抗衡离子中和再addIonsRand加盐2.4 代码实例MM/GBSA 输入文件mmpbsa.ingeneral startframe101, endframe1000, interval10, ! 跳过前 0.5 ns100 帧interval10 → 每 10 帧取 1 帧 0.05 ns/快照约 20 帧/ns verbose2, keep_files0, / gb igb5, saltcon0.15, /# 溶剂化拓扑作 -sp三个干拓扑用 cpptraj parmstrip 从 -sp 派生不要混用来源不同的拓扑 MMPBSA.py -O -i mmpbsa.in -o FINAL_RESULTS.dat \ -sp complex_solv.prmtop \ -cp complex_dry.prmtop -rp receptor.prmtop -lp ligand.prmtop \ -y prod_1.nc # 多 replicate逐条跑后取 N 个 ΔG 的均值 ± SDMMPBSA.py 的 -y 一次只接一条轨迹的帧TI 的 λ 窗口批量提交完整双拓扑准备见Amber分子动力学模拟9: 蛋白-蛋白相互作用模拟操作及亲和力计算示例§5.2# 均匀 11 窗口起步模板vdW 消耦合大基团突变时在 λ≈0/1 两端加密到 12–21 窗口与 §2.2 口径一致 for L in 0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00; do sed s/LAMBDA/${L}/ ti.in.template ti_${L}.in pmemd.cuda -O -i ti_${L}.in -p complex_solv.prmtop \ -c eq.rst7 -o ti_${L}.out -r ti_${L}.rst7 -x ti_${L}.nc done # ti.in.template 里icfe1, ifsc1, clambdaLAMBDA, timask1/scmask1配体1 # 注意本例两位小数定点命名ti_0.10…ti_1.00字典序恰好正确若改用整数索引命名ti_0, ti_1, …, ti_10则 ti_10 会排在 ti_2 前——汇总时用 ti_${L} 显式循环或 sort -t_ -k2 -g三、场景二蛋白水溶液构象稳定性3.1 场景说明评估野生型 vs 突变体、不同温度/pH 下蛋白是否保持天然构象或验证同源建模/Alphafold 结构在溶剂里是否稳定。典型问题SOD1 A4V 突变为什么致病——是整体去折叠还是局部 loop 松动这类场景的核心不是跑出 ΔG而是把收敛的结构分析做扎实。3.2 模拟要点要点说明系综全程 NPT。研究热稳定性可加 350/400 K 对照组时长小蛋白150 残基100 ns 起步含柔性 loop 或 IDR 片段 500 ns–1 μs力场-水模型配对使用ff19SBOPC官方推荐组合或 ff14SBTIP3P传统配对与 §2.2 口径一致初始结构优先高分辨率 X-ray2.0 ÅNMR ensemble 取代表结构或分别跑末端处理缺失末端用 tleap 补全或加 ACE/NME 封端避免假电荷末端续跑长模拟分段跑每 100–500 ns 重启一次irest1, ntx5保留速度3.3 避免踩坑❌ 踩坑✅ 正确做法RMSD 还在爬升就停止模拟开始分析主链 RMSD 进入平台波动 1–2 Å Rg 稳定才算采样充分只看 RMSD 一条曲线RMSF定位柔性区 Rg整体展开 DSSP二级结构占比 关键疏水核心距离四件套一起看100 ns 没看到构象变化就断言稳定μs 级构象切换如 SH3 domain常规 MD 捕捉不到改 GaMD/REMD§五续跑用.inpcrd无速度用.rst7/.ncrstntx5, irest1且ig-1换随机种子忘开iwrap1长模拟坐标漂出盒子iwrap1把分子绕回主盒分析/可视化都需要3.4 代码实例生产段输入prod.in5 ns NPT × 6-8 replicates5 ns production NPT (×6-8 replicates 共识时长) cntrl imin0, irest1, ntx5, nstlim2500000, dt0.002, ntt3, temp0300.0, gamma_ln2.0, ig-1, ntb2, ntp1, barostat2, pres01.0, taup2.0, ntc2, ntf2, cut10.0, ntxo2, ioutfm1, ntpr2500, ntwx2500, ntwr500000, iwrap1, /换算关系要门儿清nstlim × dt是总时长5 ns 对应nstlim2500000 × dt0.002100 ns 对应nstlim50000000 × dt0.002ntwx × dt是帧间隔2500 × 2 fs 5 ps/帧——5 ns 产 1000 帧100 ns 产 20 000 帧。写过ntwx1000 是每 1 ps这种注释的都是没算过——1000 步 × 2 fs 2 ps。cpptraj 四件套分析parm complex.prmtop trajin prod.nc rms first :1-243CA out rmsd_ca.dat # 主链 RMSD对比首帧 atomicfluct :1-243CA out rmsf.dat byres # 残基 RMSFbyres 按残基聚合 gyrate :1-243 out rg.dat mass # 回转半径 secstruct :1-243 out dssp.dat sumout dssp_sum.dat # DSSP 二级结构占比 run四、场景三蛋白–蛋白相互作用PPI4.1 场景说明研究复合物界面稳定性、识别热点残基、解释突变对结合的影响或给 PPI 抑制剂设计提供界面细节p53–MDM2 这类界面抑制剂项目就是典型。PPI 与蛋白-配体的本质差异界面大、重排慢、界面水常参与介导结合所以对采样时长和溶剂处理的要求都高一档。4.2 模拟要点要点说明初始结构优先实验复合物PDB对接起点HADDOCK/ZDOCK必须先验证界面合理性系综NPT显式溶剂——界面水参与氢键网络隐式溶剂在此场景不可用时长界面稳定化 ≥100 ns结合/解离事件常需 μs 级 → 用 GaMD/PPI-GaMDWang Miao 2022 用 6×2 μs PPI-GaMD 捕获了 barnase–barstar 的重复解离/再结合体系规模常 10 万原子必须 GPUpmemd.cuda分析界面接触数4.5 Å、界面氢键/盐桥、ΔSASA、界面残基 RMSF、质心距离自由能MM/GBSA per-residue 分解找热点可用但 PPI 界面大、熵贡献复杂精度低于蛋白-配体场景——结论要谨慎4.3 避免踩坑❌ 踩坑✅ 正确做法对接姿势没验证就开 500 ns 生产先检查界面互补性/关键残基接触短模拟20 ns观察界面是否保持只看全蛋白 RMSD 判断稳定用界面残基 RMSD 接触数时间序列全蛋白 RMSD 平台可能掩盖界面重排丢掉界面结晶水建体系时保留介导氢键的界面水solvatebox前处理轨迹只存溶质原子省空间PPI 分析常需要界面水ntwprt慎用MM/GBSA 分解出的热点直接当实验事实与丙氨酸扫描实验对齐验证后再用hbond不给距离/角度判据就数界面氢键hbond :1-108 :109-195 out hb.dat dist 3.5 angle 120——两 mask 才是界面间氢键默认 3.0 Å/135° 偏严4.4 代码实例cpptraj 界面分析全套parm complex.prmtop trajin prod.nc rms first :1-195CA out rmsd_all.dat # 界面接触A 链 1-108 vs B 链 109-1954.5 Å 内残基对 nativecontacts :1-108 :109-195 byresidue distance 4.5 \ out nc.dat writecontacts contacts.pdb resout nc_residues.dat # 质心距离监测解离趋势mass 即按质量加权质心 distance COMdist :1-108 :109-195 out com_dist.dat mass # SASAΔSASA SASA_A SASA_B − SASA_AB molsurf :1-108 out sasa_A.dat molsurf :109-195 out sasa_B.dat molsurf :1-195 out sasa_AB.dat # 界面氢键两个 mask 只数 A↔B 之间的氢键intramol 是错的那是链内氢键 hbond interHB out hb.dat :1-108 :109-195 dist 3.5 angle 120 runGaMD 增强ctrl 段关键开关参数含义见 Amber24 手册 §24.3 GaMD 一节官方教程暂无 GaMD 独立篇cntrl imin0, irest0, ntx1, dt0.002, nstlim5000000, ! 生产 10 nsnstlim 必须给GaMD 全程含 ntcmdnteb 统计段 ! PPI 体系用 igamd16PPI-GaMD 单加势或 igamd17dual-boost 双加势 ! 通用 dual-boost 用 igamd3势能 二面角双加势更激进用 igamd5势能 非键双加势 igamd3, iE1, iEP1, iED1, ntcmdprep200000, ntcmd1000000, ntebprep200000, nteb1000000, ntave50000, sigma0P6.0, sigma0D6.0, ! 各参数须为 ntave 的整数倍 ntt3, temp0300.0, gamma_ln2.0, ntb2, ntp1, pres01.0, taup2.0, ntxo2, ioutfm1, ntpr5000, ntwx5000, /五、场景四~六速览膜蛋白、核酸、IDP这三个场景各有一整篇姊妹篇Amber分子动力学模拟11: 膜蛋白动力学模拟操作/Amber分子动力学模拟10: 蛋白-核酸复合物模拟操作这里只给场景层判断 最高频的坑。5.1 膜蛋白Amber分子动力学模拟11: 膜蛋白动力学模拟操作说明GPCR 激活、离子通道门控、转运蛋白构象循环蛋白必须嵌在脂质双层里才有意义。要点力场 ff14SB Lipid21每条磷脂拆头基两条尾部三残基建体系用 PACKMOL-MemgenAmberTools 自带或 CHARMM-GUI 转换生产段压浴用ntp2半各向异性——膜平面与膜法向分开缩放膜面积才能松弛。避坑ntp1各向同性压浴会人为压制膜面积涨落八面体盒子 ntp2直接报 Nonisotropic scaling on nonorthorhombic unit cellsLipid21 不含胆固醇cpptraj 分析时脂质残基名是:PC,:PA,:OL不是:POPC。# PACKMOL-Memgen 建 POPC 膜关键参数2026.3.25 版实测 packmol-memgen --pdb receptor.pdb \ --lipids POPC --distxy_fix 100 --dist_wat 22 \ --salt --saltcon 0.15 --notprotonate # 注意--salt 是开关、浓度另给 --saltcon只给 --saltcon 会静默降级5.2 核酸与核酸–药物Amber分子动力学模拟10: 蛋白-核酸复合物模拟操作说明DNA/RNA 构象动力学、转录因子–DNA 识别、抗癌药嵌入模式如阿霉素。要点DNA 推荐leaprc.DNA.OL21Amber24 当前推荐或leaprc.DNA.OL15/leaprc.DNA.bsc1RNA 推荐leaprc.RNA.OL3没有 leaprc.RNA.OL15 这个文件——OL15 ε/ζ OL1χOL4β OL1 仅 DNARNA 对应的是 OL3 命名。磷酸骨架带来强负电12 bp 双链 24 个核苷酸其中 22 个带磷酸——5′ 端核苷酸无磷酸净电荷 ≈ −22公式 −2×(链长−1)与Amber分子动力学模拟10: 蛋白-核酸复合物模拟操作官方教程口径一致离子环境直接决定结构稳定——先中和再加 0.15 M 盐RNA 体系考虑 Mg²⁺。避坑整链 HETATM 的核酸结构不要按杂原子清理会把整条 DNA 删掉糖环 puckerDNA C2′-endo / RNA C3′-endo是判断模拟可信度的硬指标跑完用 cpptraj 查 pucker 分布老教程里leaprc.ff99SB 自拼 frcmod 的核酸参数已过时DNA 直接用 OL21/OL15RNA 用 OL3。5.3 IDP / 相分离本文收尾场景说明FUS、TDP-43 这类内在无序蛋白没有单一天然态目标是构象系综而非单结构——Rg 分布、瞬态接触图谱才是输出物。要点力场换a99SB-dispRobustelli 2018为折叠无序蛋白同时训练标准蛋白力场会把 IDP 压得过度紧致采样靠 GaMD 或温度副本交换 REMDμs 级总采样是入场券。避坑用 ff14SB/ff19SB 跑 IDP 再报紧凑构象为主是力场伪影单条长轨迹不等于系综采样IDP 结论必须来自多条轨迹的合并分布REMD 的温度阶梯要覆盖目标温度并保证相邻副本交换率 20–40%。六、什么时候该上增强采样——场景信号 注意事项常规 MD 撞上下面这些信号时说明体系的时间尺度超出单轨迹能力该换增强采样了具体操作在系列增强采样篇展开Amber分子动力学模拟18.0: Amber增强采样介绍总览 Amber分子动力学模拟18.1: GaMD-高斯加速分子动力学/18.2: T-REMD 温度副本交换分子动力学/18.3: SMD拉伸分子动力学/18.4: Umbrella Sampling伞形采样/18.5: CpHMD 恒pH分子动力学五个分篇。六个该换的信号#信号对应方法增强采样分篇为什么它合适1常规 MD 里构象从不翻转——二面角/loop 取向整条轨迹锁死在一个盆地GaMDigamd1/3双加势详见Amber分子动力学模拟18.1: GaMD-高斯加速分子动力学高斯加势抹平底能垒无需预设反应坐标2蛋白-配体结合/解离事件要看的k_on/k_off、解离路径LiGaMD / PPI-GaMDigamd10/11/16/17详见Amber分子动力学模拟18.1: GaMD-高斯加速分子动力学只对配体或界面加势靶标不动事件加速几个数量级3IDP / 柔性肽需要系综分布Rg 分布、瞬态接触单条轨迹永远偏GaMD 重加权或T-REMD详见Amber分子动力学模拟18.2: T-REMD 温度副本交换分子动力学GaMD 单卡可跑T-REMD 用温度换遍历但要 N 张卡4问题本身是沿一条已知坐标的 PMF孔道通透、去折叠路径、结合模式差异Umbrella Sampling WHAM详见Amber分子动力学模拟18.4: Umbrella Sampling伞形采样沿反应坐标开窗每窗采样要求低积分即得 PMF5需要考察质子化态随 pH 变化催化残基、可滴定口袋CpHMD详见Amber分子动力学模拟18.5: CpHMD 恒pH分子动力学普通 MD 的质子化是冻结的CpHMD 让 MC 在 λ∈[0,1] 间切换6想拖动体系走一条路径看力学响应/构象过渡SMD详见Amber分子动力学模拟18.3: SMD拉伸分子动力学外力牵引定性看路径与中间态不直接给自由能配 WHAM 才行上手前五条注意事项细节在 18 篇各分篇先确认常规 MD 真的不够把没看到变化与看不到变化分开——前者是时长问题后者才是采样方法问题。判断依据多 replicate 短轨迹见 §二仍锁死同一盆地才有换增强采样的正当性。GaMD 只在 pmemd 系pmemd / pmemd.cuda / pmemd.cuda.MPIsander 没有 GaMD——Amber24 手册 §24.3 原话 GaMD is not available in Sander。GaMD 生产必须重加权PyReweighting二阶 cumulant直接拿加势轨迹算性质是错的。T-REMD 的成本是 N 倍温度阶梯副本数常 8–16每个副本一份完整体系相邻副本交换率目标 20–40%阶梯设计错了交换率过低等于白烧卡。多副本通信要求 MPI 形态的可执行文件multipmemd -ng每组 ≥2 rank。dt 与加势的兼容GaMD/LiGaMD 软核相关实现实测 dt0.002 可能爆温TEMP 飙到 10⁴ K降到 dt0.001 更稳先短跑验证再上生产。增强采样 ≠ 万能它加速遍历不修正力场——力场伪影如 ff14SB 压缩 IDP在增强采样下只会更快地收敛到错误的分布。先把力场选对§五 IDP 换 a99SB-disp再谈采样。七、什么时候需要位置约束ntr1——场景、操作与注意事项位置约束positional restraintntr1用简谐势把指定原子拴在参考坐标上E k·Δx²k restraint_wtΔx 是偏离参考位置的距离。它只该出现在平衡阶段——生产的目的是自由采样约束不放开等于没跑。四个需要约束的场景#场景约束对象力常数kcal/mol/Ų出处1加热阶段heat蛋白重原子/骨架CA,C,N,O或!:WAT !H5–10通用协议防热冲击把结构震散2NPT 密度平衡equil1同上1–2阶梯递减Amber分子动力学模拟13.1: MD模拟要点汇总一§五10→5→1→0 或 10→2→1→0 阶梯释放3核酸体系平衡整条核酸:1-2425 → 0.5两段Amber分子动力学模拟10: 蛋白-核酸复合物模拟操作官方教程25 高位起步0.5 收尾4配体结合位/共价键邻近骨架 配体CA,C,N,O:LIG1–5操作mdin 命令行各一半缺一不可# equil1.in —— NPT 平衡骨架配体弱约束 cntrl imin0, nstlim50000, dt0.002, ntt3, temp0300.0, gamma_ln2.0, ntb2, ntp1, barostat2, pres01.0, taup2.0, ntr1, restraint_wt1.0, restraintmaskCA,C,N,O | :LIG !H, /# pmemd 系ntr1 必须显式给 -ref参考坐标rst7 格式sander 不给时默认读 -c 同名文件 pmemd.cuda -O -i equil1.in -p complex.prmtop \ -c heat.rst7 -ref heat.rst7 \ -o equil1.out -r equil1.rst7 -x equil1.nc五条注意事项生产段 ntr0约束是给平衡用的。生产还拴着RMSD漂亮是假的——构象根本没采样。例外只有 targeted MD-ref 配合 tgtrmsd这类特殊研究。pmemd 的 -ref 是硬要求ntr1 不给-ref直接 OPEN 报错Amber分子动力学模拟11: 膜蛋白动力学模拟操作的 CHARMM-GUI 案例cp step5_input.rst7 xxx.refc补上即好参考坐标用上一阶段的 rst7不是晶体 PDB——已经过最小化几何合理。restraintmask 用平衡阶段的重排编号tleapcombine会重排残基号写 mask 前用desc complex核对Amber分子动力学模拟10: 蛋白-核酸复合物模拟操作的坑mask 字符串上限 256 字符。力常数阶梯释放不要一步放开10→5→1→0或 10→2→1→0每档 50–100 ps从 10 直接到 0 会看到 RMSD 跳变——那是约束势能瞬间消失的回弹不是物理构象变化Amber分子动力学模拟13.1: MD模拟要点汇总一。想冻结原子用 ibelly不是加大 k真要钉死某几个原子金属表面、QM/MM 活性区外用ibelly1, bellymask...只让 mask 内原子动——ntr 加大 k 只是近似冻结且力常数大易与 SHAKE 打架。注意 ibelly 与 igb0隐式溶剂互斥。八、SHAKEntc/ntf——什么时候开、怎么配对SHAKE 是什么把含氢键的键长冻结约束求解不是力去掉体系里最快的运动C–H/O–H 伸缩周期 ~10 fs从而允许 2 fs 步长。手册原话The SHAKE option should be used for most MD calculations——默认该开。ntc/ntf 取值与配对Amber24 手册 §21.6ntc含义ntf 配对用在1不约束默认ntf1最小化特殊动力学见注意事项 32约束含氢键ntf2MD 生产标配TIP3P 水走三点专用算法手册明示 NTFNTC23约束所有键含重原子键ntf3少用sander 并行/QM-MM 不支持操作cntrl imin0, dt0.002, nstlim2500000, ! 2 fs 步长SHAKE 开启的前提 ntc2, ntf2, ! 约束含氢键 跳过这些键的力计算 ... /四条注意事项最小化不开 SHAKEntc1SHAKE 基于动力学最小化器看不见它短 min 为了去掉 bad contacts 除外手册原文。系列命令链里 min1/min2 都是 ntc1。ntf 跟 ntc 走NTFNTCntc2 配 ntf1 不报错但浪费——含氢键还在算力却没被用ntc2/ntf2 才是省力配对Amber分子动力学模拟13.2: MD模拟要点汇总二§5.2 均此口径。例外场景 ntf1GaMD/LiGaMD 含软核的实现要求 ntf1力全算且此时 dt0.002 实测可能爆温降到dt0.001见 §六注意事项 4。升到 4 fs 要 HMRdt0.004 必须先做氢质量重分配parmedHMassRepartition把 H 质量提到 ~3 a.m.u.、从重原子挪过来只开 SHAKE 不做 HMR 就上 4 fs第一步就崩。HMR 后注意水扩散性质轻微偏快。九、跨场景通用避坑清单系综顺序最小化两步先固定溶质再全放开→ NVT 加热0→300 K50–100 ps→ NPT 平衡1–2 ns密度到 ~1.0 g/cm³→ 生产。跳过密度平衡直接生产前几十 ns 都在补密度。续跑三件套irest1, ntx5, ig-1——读速度、换种子。只给坐标不给速度会看到能量跳变。格式统一 NetCDFioutfm1轨迹ntxo2重启。ASCII 只留给调试。帧间隔自己算ntwx × dt。写进输入文件注释前先乘一遍别抄别人错的。iwrap1长模拟不绕盒坐标漂移会让分析和可视化全部错位。先中和再加盐净电荷不为零的体系谈 0.15 M 生理盐浓度没有意义。独立重复MM/GBSA 排序结论至少 6–8 条不同初始速度的轨迹ig不同见 §二机制/稳定性类结论至少 3 条。单轨迹的稳定性不构成统计证据。收敛才下结论能量、RMSD、分析指标都看后半段平台的均值 ± 波动用 block average 估误差。跑完 100 ns本身不是结论。十、结语六大场景共用一套 Amber 引擎差别在力场组合、系综设置、采样时长、分析指标这四个旋钮怎么拧。方法选型的顺序永远是先定要回答的问题排序热点系综再选方法MM/GBSA / TI / GaMD再估成本条数 × 时长 × 体系大小最后定验证指标——指标不收敛数字再漂亮也不能写进报告。MD 不是跑完就行是采样充分、指标收敛才有效。参考来源Tian C, Kasavajhala K, Belfon KAA, et al. ff19SB: Amino-Acid-Specific Protein Backbone Parameters Trained against Quantum Mechanics Energy Surfaces in Solution.JCTC16, 528–552 (2020). DOI 10.1021/acs.jctc.9b00591Case DA, Aktulga HM, Belyanaya K, et al. AmberTools.JCIM63, 6183–6191 (2023). DOI 10.1021/acs.jcim.3c01153Zgarbová M, Šponer J, Otyepka M, et al. Refinement of the Sugar–Phosphate Backbone Torsion Beta for AMBER Force Fields.JCTC11, 5723–5736 (2016). DOI 10.1021/acs.jctc.5b00716Dickson CJ, Walker RC, Gould IR. Lipid21: Complex Lipid Membrane Simulations with AMBER.JCTC18, 1726–1736 (2022). DOI 10.1021/acs.jctc.1c01217Miao Y, Feher VA, McCammon JA. Gaussian Accelerated Molecular Dynamics: Unconstrained Enhanced Sampling for Biomolecular Simulations.JCTC11, 3584–3595 (2015). DOI 10.1021/acs.jctc.5b00436Wang J, Miao Y. Protein–Protein Interaction-Gaussian Accelerated Molecular Dynamics (PPI-GaMD).JCTC18, 1275–1285 (2022). DOI 10.1021/acs.jctc.1c00974Robustelli P, Piana S, Shaw DE. Developing a molecular dynamics force field for both folded and disordered proteins.PNAS115 (2018). DOI 10.1073/pnas.1800690115Love O, Winkler L, Cheatham TE. van der Waals Parameter Scanning with Amber Nucleic Acid Force Fields: Revisiting Means to Better Capture the RNA/DNA Structure through MD.JCTC20, 625–643 (2023). DOI 10.1021/acs.jctc.3c01164AMBER 官方教程索引Amber Tutorials