拉伸动力学模拟

我做 拉伸动力学模拟 最初是客户想知道一个蛋白域被 AFM 针头拽开要多少力。我设了个 0.1 nm/ns 的拉速跑了一发，给出 200 pN，文献里却是 80 pN。排查发现拉速差了两个数量级，外力把体系拽进了非平衡态。那次让我把 拉伸动力学模拟 的”拉速依赖”当成铁律来对待。

一、为什么拉伸要先讲拉速

拉伸动力学模拟 里，拉速直接决定测到的是平衡解离力还是动力学阻力。拉得越快，测到的力越高，因为键来不及按平衡路径断。我做拉力扫描时一般从 0.01 到 0.5 nm/ns 取几档，外推到低速极限才接近实验 AFM/光镊的量级。曾经一个 DNA 双链解链，高速给 120 pN、低速给 15 pN，差出一个量级，不扫描根本不知道取哪个。

二、核心原理：SMD 与 PMF

拉伸动力学模拟 常用受约束动力学（SMD）：把反应坐标（如两端距离）以恒定速度拉动，记录约束力。单条轨迹给的是力-伸长曲线，有涨落；要得自由能剖面（PMF），得用伞形采样（US）或 ABF 沿反应坐标分段加偏置势，再用 WHAM 整合。我做拉伸动力学模拟时，定性看趋势用 SMD，定量报解离能必上 US+WHAM，否则数字没有统计意义。

三、关键技术要点：力场与溶剂

拉伸断的是非共价键（氢键、疏水），力场对这类相互作用的精度直接决定力值。我用 CHARMM36m 或 Amber 系列，并开显式水（拉伸在水相才有意义）。曾经一个肽-表面粘附，隐式水给的解离力比显式水高 40%，因为没算脱水合代价。另一个坑是弹簧常数，SMD 的耦合弹簧太硬会引入额外刚度，我一般取软弹簧让系统弛豫跟得上。

四、实操流程：从建系到 PMF

我的步骤：建拉伸体系（两端选锚点，一端固定一端拉）→ 平衡 → SMD 拉力扫描定大致力窗 → 沿坐标布伞形窗口 → 各窗独立 MD 采样 → WHAM 得 PMF → 取势垒为解离能。做 拉伸动力学模拟 最容易被省的是平衡和窗口重叠检查，但 WHAM 要求相邻窗口数据重叠 1~2 Å，缺了整合就发散，我一定画图确认。

实算一个蛋白域 unfolding：拉速 0.01、0.1、0.5 nm/ns 三档，解离力分别是 85、140、210 pN，外推到 0 nm/ns 附近约 70 pN，和光镊实测 75 pN 吻合。如果不扫描，只取常用的 0.1 nm/ns，会报 140 pN，高估近一倍，客户据此设计的界面强度会系统性偏强。所以拉伸的拉速扫描不是锦上添花，是结果可信的前提——你报的每一个力值，都得先说清楚是在哪个拉速下测的。

五、常见坑：反应坐标选错

拉伸动力学模拟 若反应坐标没抓住真正的断键路径，会得到假平台。我见过把端点距离当坐标，结果中间发生构象翻转绕过了真实断点，PMF 峰偏低。换成因子的内禀反应坐标或先短 MD 观察断键顺序再定坐标才稳。另一个是忽略温度耦合，拉伸产热大， thermostat 不力会让局部升温软化。

六、复盘：拉伸是”测机械强度的探针”

回过头看，拉伸动力学模拟 最值钱的是给出解离力/粘附能随序列、界面的变化规律，帮设计更韧或更脆的材料与生物界面。我接项目时默认先要拉伸场景（蛋白 unfolding、界面剥离还是高分子拽断），再定协议。把拉速扫描、显式水、PMF 三件事做对，拉伸结果才敢和 AFM 对标。被证明有用的，是愿意为拉速多跑三档的耐心。下次有人只要一个拉力值，我会先反问——你这个拉速，对应实验里的哪种加载。