24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 687  |  回复: 16
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

qzhaosdu

金虫 (著名写手)

[交流] 【求助】请教做膜模拟的朋友,有几个小问题

最近想做一个离子对于双层膜形貌影响的,但是有几个问题一直都还不是很 明白,

1.双层膜我已经建好,膜上下的水层也已经加上。然后就是加盒子了吗?那盒子的大小怎么选呢?因为加了盒子以后我就不需要加水了,所以我觉的盒子大小挺关键的,要是选不对的话2个盒子之间的膜靠的太远,周期性边界条件就没有意义了。

2.如果所要研究的分子本来就是离子型的,比如说SDS(12烷基硫酸钠),是不是先把十二烷基硫酸根离子作为一个分子建模,然后加了再用Na+替换水呢?因为实验上用的SDS是有Na+包含的(这个我不是很清楚,应该不会只是硫酸根吧?)那么这样研究离子对于膜的影响是不是有意义呢?还需不需要加别的盐呢?

[ Last edited by lei0736 on 2009-11-25 at 12:18 ]
回复此楼
Anewday,anewhour,anewminute,anewpeople.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)

引用回帖:
Originally posted by bay__gulf at 2009-6-13 19:45:
感觉做cgmd 和dpd 的关键是相互作用的形式太少,不能准确说明问题。
而我们模拟就是为了预测,如果仅仅为了得到好的数据,而不顾实际才是得不偿失。

我大概就知道这么多了,对于胶束这些也没有做过。
现在只 ...

还是要谢谢专家。。
Anewday,anewhour,anewminute,anewpeople.
16楼2009-06-13 19:49:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
xuefei06(金币+3,VIP+0):3Q! 6-13 18:34
yuhuobuku(金币+2,VIP+0):欢迎参加讨论 6-13 19:12
1 盒子的大小应该略大于体系的尺寸,然后用NpT系综,使达到合适的体积。一般的做法是限制膜平面2个方向的比例,垂直于膜的方向自由变化。

靠的太远PBC没有意义了,太近了又会引起模拟体系的不稳定。
最好大于1-2A,以肉眼看不到明显重叠为准。

2 SDS 本身是离子,而就算长程力的ewald 算法要求体系为电中性,体系必须加对离子抗衡电荷。这也是现实的情形。如果要做特定离子环境下的模拟,就要根据比例添加电荷,如果只是做水溶液,这样就可以了。
2楼2009-06-13 17:34:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)

引用回帖:
Originally posted by bay__gulf at 2009-6-13 17:34:
1 盒子的大小应该略大于体系的尺寸,然后用NpT系综,使达到合适的体积。一般的做法是限制膜平面2个方向的比例,垂直于膜的方向自由变化。

靠的太远PBC没有意义了,太近了又会引起模拟体系的不稳定。
最好大于 ...

多谢专家解答了,还有几点不清楚的啊:
1.限制膜平面2个方向的比例,这个怎么做到呢?只能用NPT系综模拟吗?可不可以先优化cell,然后再用NVT模拟呢?
2.如果只是加Na+,也就是中和SDS的负离子,这样算不算研究离子对于膜的影响呢?文献上是这么模拟的。
3.我主要是想研究离子对于膜整体形貌的破坏,不知道这个要模拟多少个ns?

谢谢专家!
Anewday,anewhour,anewminute,anewpeople.
3楼2009-06-13 18:37:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yuhuobuku(金币+6,VIP+0):欢迎参加讨论 6-13 19:12
1.限制膜平面2个方向的比例,这个怎么做到呢?
不同的软件有不同的方法,比如在NAMD中可以用 useConstantRatio 设定。
可以参考
Christian Kandt, Walter L. Ash, D. Peter Tieleman
Setting up and running molecular dynamics simulations of membrane proteins
Methods 41 (2007) 475–488


只能用NPT系综模拟吗?可不可以先优化cell,然后再用NVT模拟呢?
动力学体系包含的自由度太多,只用EM 方法很难达到好的预平衡效果。
需要使用长时间(数个ns)的动力学模拟才能达到平衡。
这一点同样参考上一篇文献



2.如果只是加Na+,也就是中和SDS的负离子,这样算不算研究离子对于膜的影响呢?
SDS 的膜,要么加对离子(counter ion) ,有么做质子化的,但需要体系电荷为0.
MD 方法用了太多的近似,它不是什么性质都能模拟,也不能太过迷信它的结果。


3.我主要是想研究离子对于膜整体形貌的破坏,不知道这个要模拟多少个ns?
J. Phys. Chem. B 2009, 113, 2378 做D2O对DPPC膜的形态破坏,用了50 ns
J. Phys. Chem. B 2009, 113, 1983 做乙醇对POPC膜的形态破坏,用了100 ns
4楼2009-06-13 18:51:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见