提示词：利用gromacs version 2025.2   力场 AMBER99SB-ILDN 研究多肽和LBP（脂多糖结合蛋白）结合自由能，已经有二者的复合物（binder_l37_s368566.pdb），给出具体步骤，以及参数文件
# 检查非标准原子、缺失残基、链标识
```请在此输入代码块语言
grep "^ATOM" binder_l26_s130261.pdb | awk '{print $4,$5,$6}' | sort | uniq -c
```
# 生成拓扑（选择 AMBER99SB-ILDN 和 TIP3P 水模型）
```请在此输入代码块语言
gmx pdb2gmx -f binder_l26_s130261.pdb \
            -o complex.gro \
            -p topol.top \
            -i posre.itp \
            -ff amber99sb-ildn \
            -water tip3p \
            -ignh
```
# 建立十二面体盒子，边缘距离 1.0 nm
```请在此输入代码块语言
gmx editconf -f complex.gro -o box.gro -c -d 1.0 -bt dodecahedron
```
# 填充溶剂（SPC/E 或 TIP3P，与力场匹配）
```请在此输入代码块语言
gmx solvate -cp box.gro -cs spc216.gro -o solv.gro -p topol.top
```
# 创建 ions.mdp（仅用于生成 tpr 文件）：
```请在此输入代码块语言
; ions.mdp - only for grompp/genion
integrator  = steep
nsteps      = 1000
cutoff-scheme   = Verlet
coulombtype     = PME
rcoulomb        = 1.2
rvdw            = 1.2
pbc             = xyz
```
# 执行离子化：
```请在此输入代码块语言
gmx grompp -f ions.mdp -c solv.gro -p topol.top -o ions.tpr -maxwarn 2
```
# 添加离子：中性化 + 150 mM NaCl
```请在此输入代码块语言
gmx genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -nname CL \
           -neutral -conc 0.15 -seed 42
```
# 能量最小化 em.mdp
```请在此输入代码块语言
; Run Control
integrator      = steep
nsteps          = 50000         ; 最大步数
emtol           = 500.0         ; 最大力收敛标准 (kJ/mol/nm)
emstep          = 0.001         ; 初始步长 (nm)

; Output
nstlog          = 100
nstenergy       = 100

; Neighbor searching
cutoff-scheme   = Verlet
nstlist         = 20
ns_type         = grid
verlet-buffer-tolerance = 0.005

; Electrostatics
coulombtype     = PME
pme_order       = 4
fourierspacing  = 0.12
rcoulomb        = 1.2

; van der Waals
vdwtype         = Cut-off
rvdw            = 1.2
rvdw-switch     = 1.0
DispCorr        = EnerPres     ; 对能量和压力进行长程校正

; Bonds
constraints     = none         ; EM 阶段不约束
```
# 执行能量最小化：
```请在此输入代码块语言
gmx grompp -f em.mdp -c solv_ions.gro -p topol.top -o em.tpr -maxwarn 1
gmx mdrun -deffnm em -v -ntmpi 1 -ntomp 16
```
# NVT 平衡（控温，控体积）nvt.mdp
```请在此输入代码块语言
; Run Control
integrator      = md
dt              = 0.002         ; 2 fs 步长
nsteps          = 500000        ; 1 ns
nstxout         = 5000
nstvout         = 5000
nstfout         = 5000
nstlog          = 1000
nstenergy       = 1000
nstcalcenergy   = 100

; Constraints
constraints             = h-bonds
constraint_algorithm    = LINCS
lincs_iter              = 1
lincs_order             = 4

; Neighbor searching
cutoff-scheme   = Verlet
nstlist         = 20
coulombtype     = PME
pme_order       = 4
fourierspacing  = 0.12
rcoulomb        = 1.2
rvdw            = 1.2
DispCorr        = EnerPres

; Temperature coupling
tcoupl          = V-rescale
tc-grps         = Protein Non-Protein
tau_t           = 0.1 0.1
ref_t           = 300 300

; Pressure coupling
pcoupl          = no

; PBC
pbc             = xyz

; Velocity generation
gen_vel         = yes
gen_temp        = 300
gen_seed        = -1

; COM motion removal
comm-mode       = Linear
comm-grps       = Protein Non-Protein

; Position restraints (需在 topol.top 中引入 #include "posre.itp")
define          = -DPOSRES
```
# 执行 NVT 平衡：
```请在此输入代码块语言
gmx grompp -f nvt.mdp -c em.gro -r em.gro -p topol.top -o nvt.tpr -maxwarn 1
gmx mdrun -deffnm nvt -v -ntmpi 1 -ntomp 16
```
# NPT 平衡（控温，控压）npt.mdp
```请在此输入代码块语言
; Run Control
integrator      = md
dt              = 0.002
nsteps          = 500000        ; 1 ns
nstxout         = 5000
nstvout         = 5000
nstfout         = 5000
nstlog          = 1000
nstenergy       = 1000

; Constraints
constraints             = h-bonds
constraint_algorithm    = LINCS
lincs_iter              = 1
lincs_order             = 4

; Neighbor searching
cutoff-scheme   = Verlet
nstlist         = 20
coulombtype     = PME
pme_order       = 4
fourierspacing  = 0.12
rcoulomb        = 1.2
rvdw            = 1.2
DispCorr        = EnerPres

; Temperature coupling
tcoupl          = V-rescale
tc-grps         = Protein Non-Protein
tau_t           = 0.1 0.1
ref_t           = 300 300

; Pressure coupling (GROMACS 2020+ 推荐 C-rescale，Parrinello-Rahman 也可)
pcoupl                  = C-rescale
pcoupltype              = isotropic
tau_p                   = 2.0
ref_p                   = 1.0
compressibility         = 4.5e-5

; PBC
pbc             = xyz

; Continue from NVT
continuation    = yes
gen_vel         = no

; COM motion removal
comm-mode       = Linear
comm-grps       = Protein Non-Protein

; Position restraints (逐步放开，若需要可在此阶段去掉)
define          = -DPOSRES
```
# 执行 NPT 平衡：
```请在此输入代码块语言
gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -t nvt.cpt -p topol.top -o npt.tpr -maxwarn 1
gmx mdrun -deffnm npt -v -ntmpi 1 -ntomp 16
```
# 生产 MD 模拟 md.mdp （100 ns）
```请在此输入代码块语言
; Run Control
integrator      = md
dt              = 0.002
nsteps          = 50000000      ; 100 ns (建议至少 200-500 ns，多肽-蛋白需充分采样)
nstxout         = 50000         ; 全精度坐标输出频率 (100 ps)
nstvout         = 50000
nstfout         = 50000
nstlog          = 5000
nstenergy       = 5000
nstcalcenergy   = 100

; Output compression
nstxout-compressed  = 5000      ; XTC 轨迹每 10 ps 写一次
compressed-x-grps   = System

; Constraints
constraints             = h-bonds
constraint_algorithm    = LINCS
lincs_iter              = 1
lincs_order             = 4

; Neighbor searching
cutoff-scheme   = Verlet
nstlist         = 20
coulombtype     = PME
pme_order       = 4
fourierspacing  = 0.12
rcoulomb        = 1.2
rvdw            = 1.2
DispCorr        = EnerPres

; Temperature coupling
tcoupl          = V-rescale
tc-grps         = Protein Non-Protein
tau_t           = 0.1 0.1
ref_t           = 300 300

; Pressure coupling
pcoupl                  = C-rescale
pcoupltype              = isotropic
tau_p                   = 2.0
ref_p                   = 1.0
compressibility         = 4.5e-5

; PBC
pbc             = xyz

; Continue from NPT
continuation    = yes
gen_vel         = no

; COM motion removal
comm-mode       = Linear
comm-grps       = Protein Non-Protein
```
# 生产 MD 模拟 md.mdp （50 ns）
```请在此输入代码块语言
; Run Control
integrator      = md
dt              = 0.002
nsteps          = 25000000      ; 50 ns (25000000 × 0.002 fs = 50000 ps)
nstxout         = 50000
nstvout         = 50000
nstfout         = 50000
nstlog          = 5000
nstenergy       = 5000
nstcalcenergy   = 100
nstxout-compressed  = 5000
compressed-x-grps   = System

; ... 其余参数保持不变 ...
```
# 1. 创建 screen 并运行（带续跑保护）
```请在此输入代码块语言
screen -S md_prod1
# 3. 重连查看
screen -r md_prod1
# 执行 MD 模拟：
gmx grompp -f md.mdp -c npt.gro -t npt.cpt -p topol.top -o md.tpr -maxwarn 1
gmx mdrun -deffnm md -v -ntmpi 1 -ntomp 16 -gpu_id 0 2>&1
gmx mdrun -deffnm md -v -ntmpi 1 -ntomp 16 -gpu_id 0 2>&1
gmx mdrun -deffnm md -cpi md.cpt -append -v -ntmpi 1 -ntomp 16 -gpu_id 0 -maxh 47.5  | tee -a md_run.log

# 1. 创建 screen 并运行（带续跑保护）
screen -S cluster 
gmx trjconv -f md.xtc -s md.tpr -o cluster.xtc -pbc cluster -n index.ndx
```
# 
#**仅仅为MD流程示例，实际场景还需要根据蛋白场景改变。**


提示词：利用gromacs version 2025.2   力场 AMBER99SB-ILDN 研究多肽和LBP（脂多糖结合蛋白）结合自由能，已经有二者的复合物（binder_l37_s368566.pdb），给出具体步骤，以及参数文件

MD gromacs

利用gromacs version 2025.2   力场 AMBER99SB-ILDN 研究多肽和LBP（脂多糖结合蛋白）结合自由能……

编程语言

算法

人工智能

文章

问答

视频

教程

学习中心

腾讯云实验室

直播

竞赛

腾讯云代码分析专区

腾讯iOA零信任安全管理系统专区

腾讯云架构师技术同盟交流圈

腾讯云数据库专区

腾讯云智能顾问专区

腾讯云原生专区

腾讯混元专区

腾讯云TCE专区

腾讯云Lighthouse专区

腾讯云HAI专区

腾讯云Edgeone专区

腾讯云存储专区

腾讯云智能专区

腾讯轻联专区 

腾讯云开发专区

TAPD专区

腾讯轻量云游戏服专区

EdgeOne AI 安全实战专区

腾讯云最具价值专家

腾讯云架构师技术同盟

腾讯云创作之星

腾讯云开发者先锋

腾讯云代码助手

云原生构建

TAPD 敏捷项目管理

Cloud Studio

SDK中心

API中心

命令行工具

涵盖代码开发、场景应用、自动测试全流程，助你从零构建专属AI助手

一站式MCP教程库，解锁AI应用新玩法

聚焦“写作效率、视觉美观与运行性能”三方面进行全面升级，为您提供更高效、稳定的创作环境

社区富文本&Markdown编辑器全新改版上线，欢迎大家体验!

诚挚邀请您参与本次调研，分享您的真实使用感受与建议。您的反馈至关重要，感谢您的支持与参与！

社区新版编辑器体验调研

MD gromacs-腾讯云开发者社区-腾讯云

MD gromacs

MD gromacs

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐