24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1365  |  回复: 7

lyfxie8607

木虫 (著名写手)

[求助] 300金币再开一帖,欢迎了解的朋友赐教

本人已经发过一个有关ABF模拟的帖子,在此首先感谢虫友yuxiaolin,本人受益匪浅。
但还有新的疑问,请各位了解的朋友不吝赐教,欢迎指正。
原帖内容如下:
所用软件:NAMD
模拟方法:ABF simulation
模拟对象:离子通道(膜蛋白) 单离子穿越通道模拟PMF(Potential means of force)
问题:投文章,审稿人给的意见如下 It is particularly striking that the PMF is asymmetric (different end point values between cytosolic and extracellular sides although - under periodic boundary conditions - the final ion activities are identical). This hints at a massive hysteresis or nonergodicity problem that needs to be clarified and remedied.
大概意思是生成的PMF曲线不对称。一般情况下PMF给出的数据曲线为对称的(说对称也不是很准确,应该是曲线的起点和终点的数值应该相近,红圈标出的位置应该都在0附近,如图一)。而我做的两条曲线(红线和蓝线)的起点和终点的差距比较大,如图二中的红圈部分所示。审稿人提示说是hysteresis(磁滞)或nonergodicity(非遍历性)的问题。

新的疑问:
个人理解的,觉得可能是这个原因导致的PMF起始端和末端不一致。

审稿人认为曲线两端应该是“对称”的,因为我设置了周期边界条件(periodic boundary conditions),所以可以认为离子在其穿越方向是连续跨膜,既超出周期边界的话离子会回到其原点,所以曲线两端是一致的。

可是我仔细查看了体系和模拟的条件发现个问题,模拟的离子在z轴方向跨通道的距离(82 Å)远小于设置的周期边界z轴的距离(130 Å),所以离子不可能超出周期边界,而使曲线两端的值一致。因为没有达到是离子相当于回到原点的条件,所以PMF不是对称的(如图三)。我是自己这样理解的,但是怕审稿人是这方面的“专家”,说错了反而伤到自己。
请问各位虫友,这个观点是否站得住脚?

集中的焦点在于周期边界条件,因为ABF模拟会给离子在z轴方向施加一个力,如果z轴的周期边界很大是不是就不算是连续跨通道模型了呢?

示意图.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hp2633819

新虫 (初入文坛)

★ ★ ★
zh1987hs: 金币+3, 谢谢 2013-04-08 20:20:07
同意楼主,两面环境不同正离子和膜的静电作用应该还挺强,如果初始和结束位置离膜与蛋白距离有差距的话pmf还是会差不少的,当然这也取决与simulation的离子强度
2楼2013-04-05 12:23:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuxiaolin

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
chaizhm: 金币+5, 谢谢~ 2013-04-06 09:02:50
lyfxie8607: 金币+150, ★★★很有帮助 2013-04-06 09:36:18
图形很说明问题。起始端和终止端的环境不等价,即使是加上周期性边界条件。
PMF应该是不对称的。否则,能量曲线在两端也应该是对称的。
G=E-TS+PV
你应该是做NPT模拟,在起始端和终端,S 和 V变化应该很小。
所以两端的自由能差deltaG~deltaE
从你的结果上看,G和Econsistent。
3楼2013-04-06 04:22:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wlalnut

至尊木虫 (职业作家)

睛朗

【答案】应助回帖

感谢参与,应助指数 +1
很专业呀,不太明白,但仍很专注。祝你成功。
4楼2013-04-06 07:06:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
5楼2013-04-06 08:25:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
6楼2013-04-06 09:03:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyfxie8607

木虫 (著名写手)

引用回帖:
3楼: Originally posted by yuxiaolin at 2013-04-06 04:22:42
图形很说明问题。起始端和终止端的环境不等价,即使是加上周期性边界条件。
PMF应该是不对称的。否则,能量曲线在两端也应该是对称的。
G=E-TS+PV
你应该是做NPT模拟,在起始端和终端,S 和 V变化应该很小。
所 ...

最后一句话似乎没说完?
您大概的意思就是,我的这种解释站得住脚吧?
7楼2013-04-06 09:37:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

david6307

银虫 (小有名气)

楼主,我想问一下,你这个能量值是怎么得到的,我也是用namd进行的模拟。
8楼2017-03-01 20:07:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lyfxie8607 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 295求调剂 +8 FZAC123 2026-04-03 8/400 2026-04-05 17:46 by 蓝云思雨
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 341求调剂 +3 洛多罗 2026-04-02 4/200 2026-04-04 21:36 by 智能智慧
[考研] 材料调剂 +12 一样YWY 2026-04-02 13/650 2026-04-04 20:49 by 蓝云思雨
[考研] 277求调剂 +4 12A3 2026-04-02 5/250 2026-04-04 20:28 by 蓝云思雨
[考研] 306求调剂 +3 hyb上名工 2026-04-02 3/150 2026-04-04 18:12 by 热情沙漠
[考研] 292求调剂 +21 是妍子也是研子 2026-03-30 22/1100 2026-04-03 21:44 by qlm5820
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 338求调剂 +4 zzz,,r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 266求调剂 +3 08电气工程 2026-04-03 3/150 2026-04-03 14:05 by 1753564080
[考研] 重庆大学材料与化工085600,初试370+,求求调剂建议 +8 shzhou_ 2026-04-01 9/450 2026-04-03 09:31 by 蓝云思雨
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-04-02 3/150 2026-04-02 15:06 by cal0306
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[考研] 08工科求调剂290分 +5 1314捧花 2026-04-02 8/400 2026-04-02 13:16 by 乔哒哒哒
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 285求调剂 +7 AZMK 2026-03-30 13/650 2026-04-01 17:00 by 七度不信任
[考研] 调剂申请 +8 张张张张zy 2026-03-31 9/450 2026-04-01 08:29 by zjbkx
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
信息提示
请填处理意见