中大唯信-唯信计算|Protocols | 计算配体结合自由能的新工具FMAP,Nature( 二 )


作者在程序中提供了设置漏斗形偏置势所需的参数的更详细说明 , 成功案例研究中使用的FM参数如图1所示 。
中大唯信-唯信计算|Protocols | 计算配体结合自由能的新工具FMAP,Nature
文章图片
图1FM参数设置示例
【中大唯信-唯信计算|Protocols | 计算配体结合自由能的新工具FMAP,Nature】图片源自NatureProtocols.
与目前可用的结合自由能计算手段相比 , FM具有以下优点:
1)不需要关于配体结合模式的先验知识 。 研究者可以随机选择相对于蛋白质结合位点的配体的起始位置 , 并让配体找到其到达结合位点的途径;
2)所研究的配体之间不需要结构相似性 , 例如在FEP的情况下;
3)配体、蛋白质和溶剂分子是明确的 , 并且在FM中具有完全的柔性 。 因此 , 在配体结合过程中可以细致研究和评估水分子和蛋白质的运动 。 这代表了相对于其他方法的重大改进 , 因为溶剂化和去溶剂化作用、构象选择和诱导拟合作用通常在正确识别配体结合模式中起决定性作用;
4)整个配体结合过程是可重复的;
5)FM对非收敛结果能够自我诊断 。 在CV选择不当的情况下 , 计算不会收敛 , 并且可以轻松检测到采样问题(例如磁滞现象) 。
当然 , FM方法也有一定的局限或者难点所在 , 主要体现在:
1)CV的选择 。 在FM中 , CV的选择是至关重要的一步 , 因为错误的CV选择会导致结果无法收敛 。 CV必须描述配体结合过程中系统的最慢自由度 , 从而能够区分模拟过程中访问到的系统的不同状态;
2)维度问题 。 FMAP十分灵活 , 可以适应使用任何类型的CV 。 但是 , CV的数量应尽可能少 , 因为重建自由能的计算成本随CV的数量呈指数增长;
3)自由能收敛问题 。 在错误的CV选择或忽略重要的慢自由度的情况下 , 自由能计算就不会收敛;
4)准确度 。 FM已被证明能够以小于1.5kcalmol的精度预测实验结合自由能 。 但是 , 应该考虑到FM是基于MD模拟的 。 因此 , 计算结果取决于蛋白质和配体所用力场参数的准确性以及系统的设置(例如氨基酸-配体质子化状态);
5)计算成本 。 FM模拟代码的最新版本和计算能力的指数级增长使得如今有可能在合理的时间内在本地计算机和高性能计算集群(HPC)上执行计算 。 值得一提的是在原文中 , 作者特别针对CV的选择和FM计算收敛检查给了详细的建议和步骤描述 , 本文由于篇幅所限 , 不在此一一列出 。
FMAP利用GUI进行交互以准备输入文件 , 以进行FM模拟和结果分析 。 GUI通过以下三个阶段组成的分步过程来帮助研究者进行快速地上手使用:
1)预处理阶段 。 在靶标蛋白上设置漏斗形偏置势;定义FM参数并准备用于FM模拟的输入文件(通过GUI建议标准值并显示帮助消息);
2)模拟阶段 。 用户使用其首选的MD模拟软件(包括Gromacs、Amber和NAMD)启动FM , 同时提供了命令和故障排除信息的列表;
3)后处理阶段 。 FM结果和数据检索工具的分析 。
FMAPGUI(https://github.com/limresgrp/FMAP_v1/tree/master/funnel_gui)以VMD程序插件形式体现 , 模拟代码在PLUMED-v.2(https://www.plumed.org/doc-master/user-doc/html/_installation.html)中实现 。 图2和图3展示了FMAPGUI在预处理阶段和后处理阶段的软件设置界面 。 注:原文对软件每个设置选项以及使用过程中可能遇到的问题都进行了详细的描述和给出了合理的建议 , 具体请参看原文 。
中大唯信-唯信计算|Protocols | 计算配体结合自由能的新工具FMAP,Nature
文章图片
图2FMAPGUI预处理界面(Funnel.tcl)
图片源自NatureProtocols.
中大唯信-唯信计算|Protocols | 计算配体结合自由能的新工具FMAP,Nature
文章图片
图3FMAPGUI后处理界面(FFS.tcl)


推荐阅读