在做完单点突变后（无论你是用 PyMOL、FoldX 还是 Rosetta 得到的突变体 PDB 文件），如果手头没有 GPU 显卡，利用廉价的纯 CPU 云服务器（如 8 核或 16 核的按量付费实例）跑完前期的能量最小化（EM）、**NVT（恒温）和NPT（恒压）**平衡是一个非常务实的选择。
由于没有 GPU 加速，在纯 CPU 环境下运行 GROMACS 必须格外注意线程分配和编译优化，否则计算效率会极低。以下是完整的实操工作流。
一、 云端纯 CPU 环境准备与优化
很多初学者习惯直接用 sudo apt install gromacs 安装，但这往往会得到一个没有针对你 CPU 指令集（如 AVX2 或 AVX_512）优化的版本，速度会慢几倍。
1. 推荐使用 Conda 快速安装优化版
在云服务器上，最简单且兼顾性能的方法是通过 Miniconda 安装由 conda-forge 编译的高性能版本。它会自动匹配你主机的 CPU 指令集：
# 创建并激活一个独立环境
conda create -n gmx_cpu -c conda-forge gromacs -y
conda activate gmx_cpu
验证安装及 CPU 支持情况：
gmx -version
在输出中注意查看 SIMD instructions 一栏，确保其显示为 AVX2 或 AVX_512，这意味着 GROMACS 可以全力调用你 CPU 的矢量计算单元。
二、 体系构建与准备（以突变体 mutant.pdb 为例）
确保你的 mutant.pdb 已经去除了非标准的杂原子、结晶水（如有必要），并且残基编号连续。
1. 产生拓扑文件
选择适合你体系的力场。如果是普通可溶性蛋白，推荐使用经典的 Amber99sb-ildn 或 Amber14sb。
gmx pdb2gmx -f mutant.pdb -o mutant_processed.gro -water spce -ignh
提示：-ignh 会忽略原 PDB 中的氢原子，由力场根据 pH=7 自动加氢，避免突变工具加氢不规范导致的电荷或结构冲突。
2. 定义模拟盒子与加水
这里定义一个距离蛋白边缘 1.0 nm 的立方体盒子：
gmx editconf -f mutant_processed.gro -o mutant_newbox.gro -c -d 1.0 -bt cubic
gmx solvate -cp mutant_newbox.gro -cs spc216.gro -o mutant_solv.gro -p topol.top
3. 双离子的加入（中和体系）
突变往往会改变体系的总电荷，必须加入钠离子或氯离子使体系呈电中性。
先准备一个极简的 ions.mdp（只需包含基本的参数即可，用于让 grompp 顺利编译生成 tpr）：
; ions.mdp
cutoff-scheme = Verlet
然后运行编译与离子替换：
gmx grompp -f ions.mdp -c mutant_solv.gro -p topol.top -o ions.tpr
gmx genion -s ions.tpr -o mutant_solv_ions.gro -p topol.top -pname NA -nname CL -neutral
执行 genion 时，终端会提示选择溶剂组，通常输入 13（或代表 SOL 的数字）回车即可。
三、 能量最小化（EM）
突变后，突变位点处的侧链极有可能与周围残基产生严重的原子碰撞（VdW 冲突）。直接跑平衡会瞬间崩掉，必须先进行严格的能量最小化。
新建 em.mdp 文件：
; em.mdp
define = -DFLEXIBLE ; 允许水分子内部柔性，有利于消除严重碰撞
integrator = steer ; 陡峭下降法
emtol = 1000.0 ; 停止收敛的标准 (kJ/mol/nm)
emstep = 0.01 ; 初始步长
nsteps = 50000 ; 最大步数
nstlist = 1 ; 纯CPU下更新邻区表可以稍微频繁点
cutoff-scheme = Verlet
ns-type = grid
coulombtype = PME
rcoulomb = 1.0
rvdw = 1.0
pbc = xyz
运行编译与计算：
gmx grompp -f em.mdp -c mutant_solv_ions.gro -p topol.top -o em.tpr
gmx mdrun -v -deffnm em
运行完毕后，检查 em.edr 确保势能（Potential）为负值，且最大力 Maximum force 小于 1000。
四、 NVT 恒温平衡（100 ps）
在无显卡服务器上，控制好温度耦合是防止局部结构融化的关键。
新建 nvt.mdp 文件：
; nvt.mdp
define = -DPOSRES ; 对蛋白重原子进行位置限制
integrator = md
dt = 0.002 ; 2 fs
nsteps = 50000 ; 50000步 = 100 ps
nstxout-compressed = 5000 ; 减少磁盘I/O写入，纯CPU机器的磁盘性能通常一般
cutoff-scheme = Verlet
ns-type = grid
coulombtype = PME
rcoulomb = 1.0
rvdw = 1.0
tcoupl = V-rescale ; CPU模拟推荐使用 V-rescale 控温，极其稳定
tc-grps = Protein Non-Protein
tau-t = 0.1 0.1
ref-t = 300 300
pcoupl = no ; NVT不控压
constraints = h-bonds ; 限制氢原子键
运行编译与计算：
由于是通过 SSH 连接云服务器，为防止网络波动导致连接中断、进程死掉，强烈建议使用 nohup 或 screen/tmux 挂在后台运行。
gmx grompp -f nvt.mdp -c em.gro -r em.gro -p topol.top -o nvt.tpr
# 使用 nohup 后台运行，并将输出重定向到 log
nohup gmx mdrun -v -deffnm nvt > nvt.log 2>&1 &
你可以通过 tail -f nvt.log 查看实时进度。
五、 NPT 恒压平衡（100 ps）
引入压力耦合，让溶剂盒子的密度达到真实状态（约 1000 g/L）。
新建 npt.mdp 文件：
; npt.mdp
define = -DPOSRES
integrator = md
dt = 0.002
nsteps = 50000 ; 100 ps
nstxout-compressed = 5000
cutoff-scheme = Verlet
ns-type = grid
coulombtype = PME
rcoulomb = 1.0
rvdw = 1.0
tcoupl = V-rescale
tc-grps = Protein Non-Protein
tau-t = 0.1 0.1
ref-t = 300 300
pcoupl = Berendsen ; 平衡阶段推荐使用Berendsen控压，比Parrinello-Rahman更稳健
pcoupltype = isotropic
tau-p = 2.0
ref-p = 1.0
compressibility = 4.5e-5
refcoord-scaling = com
constraints = h-bonds
运行编译与计算：
gmx grompp -f npt.mdp -c nvt.gro -t nvt.cpt -r nvt.gro -p topol.top -o npt.tpr
nohup gmx mdrun -v -deffnm npt > npt.log 2>&1 &
六、 纯 CPU 核心算力压榨技巧（核心避坑指南）
在无显卡环境下，gmx mdrun 的性能完全取决于你对 CPU 核心的合理分配。如果设置不当，多个核心之间频繁通信反而会导致速度断崖式下跌。
1. 合理设置 -nt、-ntmpi 与 -ntomp
-nt(总线程数)：应当等于你租用的云服务器的物理 CPU 核心数（如果是超线程，可以尝试物理核心数或逻辑核心数）。- 单机多核优化法则：在没有显卡时，通常不需要使用 MPI 分区。直接使用 OpenMP 并行往往效率更高。
例如，如果你租用的是一台 8核 的云服务器，运行 mdrun 时显式指定线程：
gmx mdrun -v -deffnm npt -nt 8 -ntmpi 1 -ntomp 8
这会强制 GROMACS 在单一进程内启用 8 个 OpenMP 线程，避免了局域网/进程间 MPI 通信的开销。
2. 避免频繁的文件写入（降低 I/O 损耗）
很多低配云服务器的云硬盘读写 IOPS 有限制。如果在 .mdp 文件中把 nstxout、nstvout、nstenergy 设置得太小（比如 100 步一写），CPU 就会频繁等待硬盘写入，计算速度慢得像乌龟。
- 在
nvt.mdp和npt.mdp中，将不必要的坐标和速度输出关闭，仅使用压缩轨迹nstxout-compressed = 5000（每 10 ps 记录一次）。
3. 如何判断平衡已经完成？
完成 NPT 后，利用 gmx energy 工具提取密度（Density）和压力（Pressure）：
gmx energy -f npt.edr -o density.xvg
在提示中输入对应 Density 的编号。用绘图工具或命令行查看数据，若密度在 1000 kg/m³ 附近趋于一条水平线（通常在 980 - 1020 之间波动），则说明你的突变体结构在水环境中已经基本平衡完毕，可以开展后续的生产跑（MD Production）了。
