计算化学公社

 找回密码 Forget password
 注册 Register

求助:SMD拉伸后如何选择关键帧来做PMF计算?

查看数: 160 | 评论数: 2 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-12-29 12:23

正文摘要:

蛋白质复合物体系是通过比较有无抑制剂情况下拉动drug来判断需要能量差。把drug拉动约3.75nm从结合位置拉出跨膜蛋白,mdp设置如下: nsteps      = 3750000    ; 7.5ns pull_coord1 ...

回复 Reply

FrancisCho 发表于 Post on 2025-12-30 09:12:44
1、推荐在做SMD的过程中k常数可以设置得稍微大一些,这样移动距离更好地与模拟时间线性吻合;
2、推荐模拟窗口按照已发表文献设置,我做的方向一般采样要几纳秒;缺少数据可以将histo.xvg文件拖进Origin,全选绘数据图看看哪部分没有重叠;
3、我一般是使用gmx trjconv输出质心(我做的是离子穿过路径的PMF)的坐标,通过Python代码提取每一帧对应的z坐标及其对应的帧数,先通过gnuplot查看一下以帧数为横轴,z坐标为纵轴的曲线,然后再用Python设定每0.1 nm(我自己做的体系已发表文献常用窗口间隔)来提取z坐标对应的帧数,然后输出成对应的con*.gro文件,然后复制进US脚本内。
student0618 发表于 Post on 2025-12-29 12:40:35
本帖最后由 student0618 于 2025-12-29 12:45 编辑

Umbrella sampling的好处是哪个点采样不足直接补窗口就可以,不用全部重跑。

先看Histogram重叠,哪点缺了就补那个窗口。


btw 有些论点是protein-ligand体系的话,直接拉伸的umbrella sampling不足以采样"natural unbinding pathway",可能要考虑一下相关的审稿意见。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-24 04:22 , Processed in 0.166328 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list