2. 体系建模与初始构型工程

分子动力学模拟的可靠性，从初始构型确定的那一刻就已经被大幅锁定。拓扑文件描述“谁和谁成键、原子带多少电荷”，坐标文件描述“每个原子在空间中的位置”，二者共同构成模拟的起点。本节围绕不同化学体系的构建策略、溶剂化与离子添加、以及周期性边界条件（PBC）设置展开，重点说明每一步的输入、操作与判断依据。

一、不同体系的构建策略

蛋白质体系。最常见路径是从 PDB 结构出发。第一步是补全缺失原子与缺失残基：晶体结构中柔性 loop 常因电子密度不足而缺失，需用建模工具补全；第二步是加氢并确定质子化态，这一步必须结合 pH 判断。例如 Asp、Glu 在 pH 7 去质子化带 −1，Lys、Arg 带 +1，His 的 pKa 约 6.0，在 pH 7 通常为中性，但活性位点中的 His 需按局部环境单独判断。第三步是处理二硫键、金属配位和辅因子，这些不能靠自动加氢解决，必须显式指定成键关系。

核酸体系。DNA/RNA 的构建难点在于骨架磷酸基团带负电、碱基堆积与配对氢键的几何约束强。从 PDB 提取时要注意 5′/3′ 端、缺失碱基和修饰核苷酸。核酸对力场版本敏感，建议整套使用同一核酸力场参数，不要混用不同来源的碱基参数。初始构型中若出现碱基对严重扭曲，应先做约束优化再进入自由模拟。

膜体系。膜蛋白或纯脂质双层的构建通常采用“预组装膜 + 插入蛋白”的策略。关键参数是脂质种类、上下叶不对称性、膜面积与蛋白疏水区长度的匹配。若蛋白疏水带厚度与膜厚不匹配，会出现局部膜变形甚至蛋白倾斜。构建后应检查脂质头基朝向、水相是否贯穿膜内部，以及是否存在脂质与蛋白的原子重叠。

聚合物与无机材料。聚合物需要先确定链长分布、立构规整度和端基，再通过单链构建或熔融弛豫得到无规缠结构型。无机材料（如 SiO₂、金属氧化物、黏土矿物）通常从晶体结构出发切面、扩胞，再按表面化学修饰加羟基或配体。这类体系的力场往往需要自定义参数，电荷分配方式会显著影响表面吸附行为。

二、溶剂化与离子添加

溶剂化的核心是水模型选择与盒子形状、尺寸的确定。常见水模型包括 TIP3P、SPC/E、TIP4P/2005 等，选择时需与溶质力场保持兼容。盒子形状上，球形或十二面体盒子在相同截断半径下可比立方盒子减少约 30% 的水分子数，从而降低计算量。盒子尺寸的通用判据是：溶质原子到盒子边界的最近距离不小于截断半径，通常取 1.0～1.2 nm；若使用 PME 静电，还需保证盒子最小边长大于两倍截断半径。

离子添加有两个目的：中和体系净电荷，以及达到目标离子强度。操作上先加反离子中和，再按目标浓度加入盐。浓度换算可用近似式：

N_ion ≈ C × N_A × V_box
示例：V_box = 343 nm³（约 7 nm 立方盒）
目标 0.15 mol/L NaCl
N ≈ 0.15 × 6.022e23 × 343e-24 ≈ 31 对离子
实际添加时需扣除已加入的中和离子数

添加位置应优先替换远离溶质的高能水分子，避免在蛋白表面直接放置离子造成局部结构破坏。添加后必须检查体系净电荷是否为零，并确认没有离子重叠。

三、周期性边界条件设置

PBC 的本质是让盒子在三个方向上无限复制，从而消除表面效应。设置时需关注三点：一是盒子向量与截断半径的关系，二是最小镜像约定下溶质与自身镜像的距离，三是压力耦合是否允许盒子形状变化。对于膜体系，通常只允许 z 方向独立变化以维持膜面积；对于溶液体系，可采用各向同性耦合。若模拟中观察到溶质与自身镜像发生非物理相互作用，说明盒子过小，应扩大盒子或减小截断半径。

四、构建后的检查清单

无论哪种体系，构建完成后都应执行以下检查，这些检查直接决定后续模拟是否可信：

• 原子重叠：任意两原子间距小于 0.1 nm 即视为严重重叠，需能量最小化或手动调整。
• 净电荷：体系总电荷必须为零，否则 PME 下会出现非物理漂移。
• 盒子边界距离：溶质到边界最近距离应大于截断半径，通常 ≥ 1.0 nm。
• 密度合理性：水盒子密度应在 0.98～1.02 g/cm³ 附近，膜体系需检查脂质面积。
• 拓扑一致性：原子数、残基编号、二硫键、质子化态与坐标文件一一对应。

完成上述步骤后，体系才具备进入能量最小化和平衡阶段的条件。初始构型的质量不会在后续模拟中被“自动修复”，它只会被放大或掩盖。因此，把构建当作模拟的第一道科学判断，而不是机械的预处理步骤，是本节最核心的工程态度。