24小时热门版块排行榜    

查看: 1694  |  回复: 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的回帖

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的回帖
查看全部 8 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 帮忙看看fileCode +4 wwncly 2026-08-10 8/400 2026-08-10 22:07 by demonstreets
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 为什么网上很多人说本周 12号出结果 +3 瞬息宇宙 2026-08-10 3/150 2026-08-10 18:27 by Tide man
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见