MOE分子模拟

结合案例图分析，第一张案例图是配体周围空间的静电势映射。红色区域表示负静电势（亲核/富电子），蓝色区域表示正静电势（亲电/缺电子），颜色标尺从−1.0到+1.0 kJ/mol。可以看到配体中羰基氧附近是深蓝色区域（强负电），而胺基氢附近是红色区域（正电）。这种极性分布决定了配体如何与蛋白口袋的带电残基互补。图中配体原子被圈出，说明关注的是配体表面的静电互补性，而不是整个分子的总电荷。

第二张案例图是MOE生成的2D interaction map。图中是配体骨架，碳、氮、氧原子分别用不同颜色圆球表示。周围标注了与配体相互作用的蛋白残基：Asp 46、Ala 61、Lys 118、Phe 80。相互作用类型用虚线区分：绿色虚线是氢键（H-bond），橙色虚线是π-π堆积（π-stack），灰色虚线是疏水相互作用（hydrophobic）。例如，配体的一个羰基氧与Asp 46侧链形成氢键，一个芳环与Phe 80形成π-stack，一个烷基部分与Ala 61产生疏水接触。这些相互作用共同贡献了结合自由能。

从静电势面和2D相互作用图的组合来看，这是一个MOE分子模拟项目——用的是Chemical Computing Group的MOE软件包，主要做配体-受体相互作用的可视化和能量分析。静电势面展示了配体表面的极性分布，2D interaction map展示了具体残基贡献的相互作用类型和距离。

MOE的工作流程与AutoDock不同。我先用MOE的Protein Builder给受体加氢、修补缺失环，然后用Site Finder识别潜在结合位点，再用Dock模块做柔性对接。MOE的优势在于后处理可视化，特别是2D interaction map能直观展示关键残基贡献。能量分解方面，我用MOE的GBVI/WSA dG打分函数估算结合自由能，它综合考虑了范德华、静电、氢键、去溶剂化等项。

我在我这个项目里学到两个教训。第一个是pH环境：配体和可电离残基的质子化状态必须根据实验pH设置，我在pH=7.4下给Asp和Glu去质子化、Lys和Arg质子化，否则静电作用会完全错误。第二个是侧链柔性：MOE对接时如果只让配体柔性而受体刚性，可能会错过一些由诱导契合带来的结合模式。我后来对活性位点5 Å范围内的残基做了部分柔性处理，得到的新构象与晶体结构RMSD从1.8 Å降到1.2 Å。另外ESP可视化的等值面数值范围要手动调，默认范围有时会压掉关键特征。

这张图对客户的价值是：定量给出了配体-受体之间的关键相互作用类型和残基，可以用于指导下一步的SAR改造。关于MOE分子模拟的更多实战技巧，我整理在MOE分子模拟栏目里。