| 查看: 3039 | 回复: 12 | |||||
[求助]
高斯计算,溶液中蒙脱石与黄曲霉毒素的相互作用已有1人参与
|
|||||
|
想用高斯计算,在水中,蒙脱石对黄曲霉毒素的吸附作用,即二者间的相互作用。 一、研究物质: 1、蒙脱石晶体,是由两层硅氧四面体,包夹一层铝氧八面体,形成的层状结构铝硅酸盐矿物。八面体中Al常被Mg取代,因而整个蒙脱石带负电。阳离子置换所引起的电荷不平衡主要由层间阳离子(Ca离子)来补偿,层间的阳离子常以水化状态出现。 图为单层蒙脱石的结构,中间铝氧八面体有一个Al被Mg取代,由MS绘制,作为高斯输入文件 图为蒙脱石晶体结构示意图; 2、黄曲霉素结构如图,在数据库http://webbook.nist.gov/cgi/cboo ... O6&Units=SI查找得到: 二、经验假设二者之间的作用: 1、黄曲霉毒素通过Ca离子,和蒙脱石之间形成“阳离子键桥”; 2、黄曲霉毒素和蒙脱石矿物,通过静电作用、氢键等弱作用力发生; 3、黄曲霉毒素,通过以上作用,尤其是通过作用1,甚至进入了蒙脱石阳离子层间。 以上假设有实验现象支持:阳离子溶出量随吸附量增加而增加,吸附量较大时蒙脱石XRD图发生变化。 小弟即想通过高斯计算,对这三个假设进行分析; 三、小弟目前情况: 1、初步计算计划: 方法、机组和溶剂模型 : b3lyp/6-311g scrf=(cpcm, solvent=water); (下述A为黄曲霉素,B为单层蒙脱石晶体) (1)在A上,分别连接Ca离子,为A-Ca,用上述方法和模型优化并计算能量; (2) 优化并计算B的能量; (3)将A-Ca和B放到一起计算; 作用能=E(A-CaB)-E(A)-E(B),并作BSSE矫正。 (4)然后,分析A和B之间的相互作用,键能键长等。 2、其他 小弟想做初步的计算,所以没有加极化和弥散。不过现在看起来算的挺慢的,第(3)部花了两天还没有算完(g09 A for linux,64GB内存,AMD32核,计算时加了scf=qc),也不知道出来会是什么结果。 四、存在的疑问: 小弟这学期断断续续接触计算,我们学校也没有人做量化计算的,基础很差,一路摸索,不会只能在论坛上求助,承蒙虫友照顾,但还想让大家帮我把这个工作给提提批判批判,帮小弟解答一下其中疑问: 1、蒙脱石结构(在第(3)步计算前),是应当在的MS中优化后在高斯中计算好,还是直接在gaussian中优化好呢? 2、合适的初始构型: (1)小弟在A(黄曲霉素)的两个羰基上,用gview直接和Ca连接(A-Ca),这样做是否合适?下图为A-Ca优化前后(优化后O-Ca之间的较长的键是静电作用吗?) (2) 此外,将优化后的A-Ca,放入B(蒙脱石)中,是否有什么讲究(如角度和距离等,该如何设置才具有代表性,还是应该多试试看)? (3) 我这样的做法,是不是根本就不合适呢? 3、体系电荷和自旋多重度如何考虑呢? 之前看了一些关于体系带电和自旋多重度的设置(单电子数+1=自旋多重度),但是还是不太理解,如A-Ca,应该怎么设置呢? B(单层蒙脱石,一个Al被Mg取代,体系带1个负电荷),在计算时应该设置多少呢(我设置的是“-1 1”)?把A-Ca和B放到一起,应该怎么设置呢?(单电子数目我不太会看啊) 4、溶剂模型的选择 这样的设置是否就足够了呢? scrf=(cpcm, solvent=water)。 专家ChemAndy曾在一篇计算有机物之间相互作用能的帖子中建议,“如果溶剂分子与溶质分子之间有H键,你就摆些溶剂分子在H键位点上,然后依然做CPM。”但是我感觉着不太适合这个情况,因为体系中有离子。 5、最关键的问题:预期的结果是否存在可能性? 预计的情况是: 用双层蒙脱石结构(而不是现在的单层),层间含有水合Ca离子。A(黄曲霉素)和蒙脱石在水中,进入层间的极性水化阳离子层,通过某些作用被吸附(静电作用,阳离子键桥等),而这些作用可以通过gaussian的计算得到直观的表征。 不知道这样的预计,能否实现(因为蒙脱石包括其层间是极性区域,而黄曲霉素是非极性的),如何实现(就是上述的几个问题)。 6、各位老师和前辈,请多给一点批评和建议,衷心感谢。 QQ截图20140328145844.png |
» 收录本帖的淘帖专辑推荐
矿物材料与应用 | 计算 | 量子化学和计算化学 |
» 猜你喜欢
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
2026博士申请—化学能源方向,本人熟练分子动力学模拟的使用
已经有4人回复
物理化学论文润色/翻译怎么收费?
已经有90人回复
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生
已经有0人回复
量子点电致发光器件(QLED):蓝光量子点(QDs)和纳米氧化锌电子传输材料
已经有0人回复
书籍求助:汽车市场营销理论与实务(电子版)——章小平
已经有0人回复
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有26人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生(还有1个名额)
已经有0人回复
需要一个圆偏振光诱导化学反应,请问使用圆偏振光作为光源,应该需要哪些设备搭建
已经有0人回复
华北电力大学(北京)(第一性原理计算)博士招生——学博,专博各1人
已经有40人回复
» 本主题相关商家推荐: (我也要在这里推广)

ChemiAndy
木虫 (正式写手)
- QC强帖: 3
- 应助: 52 (初中生)
- 金币: 2573.3
- 散金: 18
- 红花: 104
- 帖子: 590
- 在线: 553.8小时
- 虫号: 836907
- 注册: 2009-08-30
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★ ★ ★
zhou2009: 2014-03-29 05:38:30
zhou2009: 2014-03-29 07:05:06
zhou2009: 2014-03-29 15:28:41
nima0332(zhou2009代发): 金币+8 2014-03-29 15:30:47
zhou2009: 2014-03-29 05:38:30
zhou2009: 2014-03-29 07:05:06
zhou2009: 2014-03-29 15:28:41
nima0332(zhou2009代发): 金币+8 2014-03-29 15:30:47
|
不敢称专家,受邀请说点自己的想法。 首先你的问题的本质是通过实验现象得出了一些A-Ca(2+)-B离子键桥(Ion bridge)作用的假说(hypothesis),然后想用计算进行验证。其实这里面包含多个问题,第一,该假说是否正确?即三者是否有这样的ion bridge作用。一般计算验证者很少质疑实验的推论,而只是根据实验推论进行搭构型进行计算,因为实验是我们研究的基础。因此在给出结论时要记住这一点。第二,如果假说是正确的,那么实际的作用构形是怎样的的?有几种可能?三者在何位点相互作用?三者比例多大? 这些都不知道。因此本质上这个问题属于构形搜索,而且不涉及化学变化。构形搜索一般有两种方法,一是根据化学直觉和已知的实验证据自己搭构型,当然,这要求尽可能穷举出所有最可能的构形进行计算。楼主显然属于这种情况,而且陷在穷举构型的泥沼里不知道如何突围。而且这种穷举法还有个弊端,即只能比较已经举出的构型进行比较研究,并不能确认不存在其它可能,或者肯定就是这些构型。因此,这种方法只适用于简单体系,而不适用于复杂体系。复杂体系的构型搜索一般要用分子动力学或者蒙特卡洛模拟方法,按照真实体系的比例搭好盒子,然后进行采样,多数情况下还需要加速采样技术,然后从结果中看看哪种构型的概率最大,意味着最可能的构型,通常还能找到一些意想不到的可能构型,它们的势能稍高一些,实际存在的时间稍短一些,但也是可能的。 我不知道楼主的项目是什么性质的。如果是你的博士论文,那么毫无疑问你得去搞分子动力学的作用机理搜索。如果你只是一个硕士论文,或者短期项目,我不建议另起炉炤,中途换炮。因为这个问题的复杂性,使得它不适合作为一个硕士论文的题目。你可以继续量子化学计算,把目标定义在构型预研究(Preliminary study)。尽最大努力例举出A-Ca(2+), Ca(2+)-B的可能作用构形,分别计算作用能,得出最大可能性。然后在此基础上给出A-Ca(2+)-B的最可能作用构型,并给出作用能大小的数量级估算。这样的话,纯粹是验证ion-bridge作用机理,但是为进一步的作用机理搜索提供一个基础。 所以,我的建议是: (1) 研究A-Ca(2+)作用,提出1~2个Ca2+与A作用的可能构型,并比较作用能大小。 (2) 研究Ca(2+)-B蒙脱石作用。建议搭建一个更小的蒙脱石cluster即可,重点是表现出Ca2+与蒙脱石作用位点的化学环境,以及化学作用本质。比如,Al-Mg替代效应,离子-氢键作用,还是离子-离子作用等。我想在蒙脱石表面寻找几个可能的离子作用位点还是不难的。 (3) 最后把最可能的集中A-Ca(2+)和Ca(2+)-B连起来,优化一下构型,算一下结合能。 当然,所有的计算最好使用隐式溶剂模型PCM。在剩余的可能氢键位点上加显式的水分子。加水分子可能对稳定分子团簇有贡献。当然,你也可以比较加水与不加水分子到底有什么差别。 关于多重度,既然你的体系中没有自由基,中间体之类的含未成对电子的化学种类,因此,全部为1。电荷呢,看有没有离子。该多大就设为多大。 关于你所说的预期的结果是否存在的问题,正如我上面说的,你的题目应该是对“黄曲霉素、钙离子、蒙脱石之间离子键桥作用假说的预研究计算验证”。你的结果将提出最可能的离子键桥作用构型,和估算的作用能。但该假说是否正确,超出简单量子化学计算所能回答的问题范围。 PS:你的黄曲霉素的构型有误。注意右下角的氧-甲基。这里有正确构型:http://en.wikipedia.org/wiki/Aflatoxin_B1 |
» 本帖已获得的红花(最新10朵)
4楼2014-03-29 00:57:38
qchem
铁杆木虫 (著名写手)
- 应助: 280 (大学生)
- 金币: 6871.6
- 散金: 1411
- 红花: 43
- 帖子: 2586
- 在线: 668.9小时
- 虫号: 183011
- 注册: 2006-02-12
- 专业: 理论和计算化学
5楼2014-03-29 07:45:16
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 553.8
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+8 2014-03-29 15:32:13
nima0332: 金币+19, ★★★★★最佳答案, 渐渐有了一点入门的感觉!谢谢您! 2014-04-01 10:45:27
感谢参与,应助指数 +1
zhou2009: 金币+8 2014-03-29 15:32:13
nima0332: 金币+19, ★★★★★最佳答案, 渐渐有了一点入门的感觉!谢谢您! 2014-04-01 10:45:27
|
楼主一个短消息把我喊来了 发现楼上专家说的好极了 我补充几句吧 第一 你明显是把所谓的蒙脱石作为一个cluster来处理了 问题是 你需要对边界进行饱和吗?如果需要,用什么原子饱和?如果不需要,有什么理由不需要?你切出来的这个蒙脱石cluster是按照哪个晶面切出来的?还是自己随便切割出来的? 抱歉我不做cluster,我只能问这个问题供楼主思考。 第二 你的黄曲霉素结构的初始结构简直是太差了。在这个结构中我举几个例子,看3a,6a-dihydrofuro[2,3-b]furan这个残基结构吧,两颗中心C居然能让它呈平面结构,这简直是不可原谅的错误。另一个例子,你的甲氧基直接画错了! 可以相信数据库里的分子拓扑结构但是绝不能把拓扑结构直接拿来计算。自己去搭建模型 磨刀不误砍柴工。我觉得对你来说,先把《Exploring Chem with Elect Struct Method》读完,再读一本linux入门书,你再来咱们版继续讨论吧!真心建议 |
8楼2014-03-29 11:30:18

9楼2014-03-29 12:32:42

10楼2014-03-29 12:36:07

2楼2014-03-28 21:36:46

3楼2014-03-28 22:53:28

6楼2014-03-29 09:52:44

7楼2014-03-29 10:57:46













回复此楼

nima0332