使用Multiwfn结合VMD分析和绘制分子表面静电势分布Word下载.doc
- 文档编号:3798566
- 上传时间:2023-05-02
- 格式:DOC
- 页数:26
- 大小:1.01MB
使用Multiwfn结合VMD分析和绘制分子表面静电势分布Word下载.doc
《使用Multiwfn结合VMD分析和绘制分子表面静电势分布Word下载.doc》由会员分享,可在线阅读,更多相关《使用Multiwfn结合VMD分析和绘制分子表面静电势分布Word下载.doc(26页珍藏版)》请在冰点文库上搜索。
虽然本文主要是介绍作图方法,但也会顺带着说一下定量分子表面分析的一些功能。
2在Multiwfn中计算、统计、导出数据
启动Multiwfn后依次输入
furan.wfn
12
//定量分子表面分析
0
//开始计算
默认就是在电子密度为0.001的表面上计算静电势,并且格点间距为0.25Bohr。
格点间距可以通过选项3来设定,格点间距越小,分子表面就会用越多的顶点来描述,定量统计值、极值点位置也会越精确,之后作出的分子表面填色图的色彩过渡也越光滑,但是计算耗时将越长。
如果你的目的仅仅是绘制分子表面静电势图,可以忽略以下内容而直接跳到本节最后一段。
计算完毕后会看到分子表面上静电势极大、极小点的坐标和数值,并且输出大量统计数据,比如分子表面上静电势的最大/最小值、平均值、方差、电荷平衡度、表面积等等,它们对于了解分子特征、建立QSPR/QSAR方程预测分子理化性质和生物活性等问题都十分有用,见Multiwfn手册3.15.1节以及IJQC,85,676、JPCA,110,1005等文章的介绍。
此时会看到后处理菜单。
点击0进入图形界面,并且将Ratioofatomicsize设为4.0就可以清楚看到分子表面上静电势极大点(红点)和极小点(蓝点)的位置。
如果点击Minimal/Maximumlabel,就会显示出极值点的编号,可以和命令行窗口显示的极值点信息相对照,得知极值点上的具体数值。
点击右上角Return可以关闭图形窗口。
此例中各个极值点静电势数值如下(由于格点精度有限,所以等价的极值点的数值不完全一致,这里取了平均)
极小点1(最小点):
-20.60kcal/mol
体现氧的孤对电子对静电势的负贡献
极小点2、3:
-15.23kcal/mol
在两个beta碳(邻位碳)正上方,体现pi电子对静电势的负贡献
极大点1、4(最大点):
18.41kcal/mol
体现alpha位的氢原子所带的正电
极大点5、6:
15.21kcal/mol
体现beta位的氢原子所带的正电
极大点2、3:
-9.63kcal/mol
虽然是极大点,但是静电势数值为负,所以化学意义不大,可无视之
我们可以查看一下不同静电势区间内分子表面积,这对于了解分子表面静电势的定量分布很有益。
做法是依次输入
9
-25,22
//统计范围。
当前体系分子表面静电势范围为-20.60~18.41kcal/mol,这里将范围稍微扩大并取整数来定义范围
15
//将-25~22kcal/mol均匀分为15个区间获得表面积
3
//输入的单位为kcal/mol
立刻屏幕上就输出了不同静电势区间内的表面积,我们将Center这一列和Area这一列的数据分别从屏幕上拷下来(不会拷的读者见手册5.4节),并粘贴到诸如origin等作图工具里并作成条形图,就可以得到诸如以下图像
静电势的定量分布在这张图上一目了然。
此分子不同静电势区间内的表面积分布还是相对比较均匀的。
如果选择选项11,则程序会输出每个原子附近的分子表面上的静电势统计数据,对于了解原子在此分子中的特征很有用。
例如静电势平均值
Atom#
All/Positive/Negativeaverage
1
-10.08997
NaN
-10.08997
2
-12.59141
-12.59141
3
-12.64250
-12.64250
4
-10.04843
-10.04843
5
-13.54925
-13.54925
6
6.87357
9.80155
-3.13475
7
4.84667
8.19922
-4.56994
8
4.84069
8.21673
-4.50535
9
6.88875
9.78376
-3.16678
氧原子(5号)附近的分子表面静电势平均值最负(-13.55kcal/mol),这也容易理解,毕竟其孤对电子对静电势有很大的负贡献。
由于其附近分子表面上没有正值区域,所以正值区域的静电势平均值显示的是NaN。
呋喃中碳原子裸露在分子表面上的区域主要就是体现pi电子特征的区域,由于pi电子云使这部分区域静电势为负,所以碳原子附近静电势平均值也都为明显负值,并且无正值区域。
同时从静电势平均值上也体现出beta碳(2、3号)比alpha碳(1、4号)的pi电子云更富集因而静电势更负。
这种讨论分子表面上对应不同原子区域的定量数据的方法和分子表面极值点分析往往会得到共通的结论,但此方法可以得到更多的信息,比如alpha碳上没有出现极值点(因此对它没有任何描述),而通过分析它对应的分子表面上的局部区域的静电势平均值等数据就可以定量考察它的特征。
这种分析方法是笔者独创的,也仅有Multiwfn支持。
另外,通过选项12还可以查看分子表面上对应于指定分子片段区域的定量性质。
选过选项9之后会问你是否输出locsurf.pdb文件,通过此文件可以利用VMD查看不同原子对应的分子表面区域。
由于这不是本文的重点,所以输入n不让程序输出。
利用选项10可以查看分子表面相对于指定坐标点或指定原子核的最远和最近距离,这对于讨论分子间相互作用导致分子表面的穿透距离很有用。
此功能也可以用于计算分子半径或直径,见《谈谈分子半径的定义和计算方法》(
现在选择2将分子表面极值点导出到当前目录下surfanalysis.pdb文件中,此文件中碳和氧原子分别对应分子表面静电势极大点和极小点,pdb文件的B因子那一列的数据是静电势数值(kcal/mol)。
然后再选择6将所有分子表面顶点导出到当前目录下vtx.pdb文件中,B因子是每个顶点的静电势数值。
另外,如果你还没有当前分子的几何结构文件,应再选择5然后输入furan.pdb把当前体系的坐标导出到当前目录下furan.pdb中。
3在VMD中作图
由于Multiwfn所用的图形库的限制,Multiwfn自身无法直接产生不同颜色填充的分子表面图,所以需要借助VMD来实现。
下面的操作步骤对于当前体系比较适合,对于其它体系,应举一反三,根据实际情况进行略微调整。
启动VMD,然后把furan.pdb、surfanalysis.pdb和vtx.pdb按顺序依次拖动到VMD主窗口里(VMDMain),它们的ID将分别是0、1、2。
选择Display-DepthCueing将它关掉,否则图像会有些朦胧。
进入Graphics-Colors,选Display-Background-White将背景改为白色,并且在此界面的Color-Scale标签页里选择BWR,使分子表面的色彩根据数值范围由小到大以蓝-白-红的方式变化。
选Display-Axes-Off不让坐标轴显示出来。
注:
如果希望每次启动VMD时都自动做如上操作,可以在VMD目录下vmd.rc文件末尾添加以下4行
displaydepthcueoff
colorscalemethodBWR
colorDisplayBackgroundwhite
axeslocationOff
进入Graphics-Representations,然后执行以下3步来分别设定分子结构、分子表面极值点、分子表面顶点的显示方式。
(如果只是想简单地看一下表面静电势分布,只要做第3步即可)
(1)在SelectedMolecule一栏里选择furan.pdb,DrawingMethod选Licorice,BondRadius减小到0.1。
(2)将SelectedMolecule一栏切换到surfanalysis.pdb,在控制台输入molmodstyle01VDW0.06(0和1是显示方式编号和体系的ID,此命令代表用大小为0.06的VDW球显示。
由于GUI中VDW球最小只能设到0.1,故这里用命令行来实现)。
然后在SelectedAtoms里输入carbon并回车,然后将ColoringMethod选为ColorID,并且在右边新出现的框里选Orange2。
此时分子表面极大点就通过橙色小圆球显示出来了。
点击CreateRep按钮创建新显示方式,在SelectedAtoms里输入oxygen并回车,然后将ColorID右边的框设为Cyan,此时分子表面极小点就通过青色圆球显示出来了。
(3)在SelectedMolecule一栏里选择vtx.pdb,DrawingMethod选Points,Size设为25(设多大合决于视角的远近,在当前视角下应当让size恰好足够大,使分子表面上的顶点紧密相连,不留明显空隙),ColoringMethod选Beta(根据pdb文件里B因子那一列的数据,此例即静电势数值进行填色),在Trajectory标签页里将ColorScaleDataRange填上-22和22并点击Set,代表色彩刻度设为-22~20kcal/mol(其实此例用默认的色彩刻度范围就可以,这里只是为了取个整)。
现在分子表面填色图就出现了。
越蓝的区域静电势越负,越红的区域越正,白色区域的静电势数值在0附近。
之后给图上加上色彩刻度轴。
选Extensions-Visualization-ColorScaleBar,Colorbarwidth设为0.08,Displaytitle选on并且将Colorbartitle里写上ESP(kcal/mol),Minimum和Maximumscalevalue分别填-22和22,Numberofaxislabels输入10,Colorlabels选Black,Labelformat选Decimal。
然后点DrawColorScaleBar按钮,色彩刻度就出现在画面中了,并且VMDMain窗口中多出了一个名为ColorScaleBar的一项。
然后调整它的大小和位置,即双击VMDMain窗口中ColorScaleBar那一项当中的F标签使之变为红色(即不让色彩刻度轴在画面中的位置冻结),而双击其它项目的F标签使它们的F变为黑色(让它们的位置冻结住)。
然后激活VMD图形窗口,按t键进入平移模式,然后拖动鼠标将色彩刻度轴放置到合适位置,并且用鼠标滚轮调整它的大小。
调合适之后再按r键恢复旋转视角模式,并且在VMDMain里将ColorScaleBar那一项的F重新双击成黑色,而其它三项的F重新双击为红色。
当前的显示效果如下
此图还有许多地方值得进一步调整,也就是让分子结构显示出来(现在都被表面顶点掩盖了)、给分子表面极值点标上具体数值、让表面极值点更清楚地显示,特别是藏在分子表面后方的极值点也显示出来。
为达到这些目的,需要利用photoshop。
由于步骤比较细碎,笔者不把所有具体操作都叙述一遍,否则太罗嗦。
只要会基本的photoshop操作的人都应该能理解应该怎么实现。
4通过Photoshop改进作图效果
将分子调整到一个合适的角度,然后在VMDmain窗口里把所有条目的F标签都双击成黑色来将它们固定住,以免随后的操作过程中不慎旋转了体系。
在VMDmain窗口里面双击与furan.pdb和surfanalysis.pdb对应的条目的D标签使其变红,此时窗口内就只有分子表面和色彩刻度轴显示了出来。
然后按Alt+Printscreen键将窗口截图,在photoshop里按Ctrl+N然后按OK,再用Ctrl+V把截的图粘贴进去。
此时图像大小正好和VMD窗口大小完全一致。
在VMDmain窗口里面只让分子结构显示出来,将背景改为蓝色(只要不是白色就行,否则会和氢原子的白色连在一起),然后将窗口截图并且粘贴进之前的ps窗口里成为新的图层。
选择魔棒工具,Tolerance设0,Contiguous的对勾取消,然后点击图中蓝色区域把背景区域选上,之后按delete键去除背景。
之后将图层的不透明度(Opacity)改为40%,这样分子结构就会以半透明方式叠加到分子表面图上了,而且叠加的位置完全精确。
ps窗口里的图像目前如下所示
让VMD窗口里只显示出表面极值点,然后将窗口截图也粘贴进ps里成为新图层,同样将蓝色的背景选中并删掉。
然后用矩形选框工具恰当地把图像主体和色彩刻度轴圈上,构图满意后点Ctrl+Q将多余的区域裁掉,目前图像如下
可见目前这些表面极值点无论是正面还是背面的都是以完全不透明方式展现,看不出前后。
因此我们把那些明显处在背面的极值点都用矩形选框圈上,注意圈的时候一直按着shift键使选区依次累加。
然后在图上点右键选LayerviaCut,这样处在图中正面(或者边缘)的极值点与处在图中背面的极值点就分别用两个图层来储存了。
将后者的图层的不透明度设为50%。
最后,在图上用文本工具标上一部分极值点的静电势数值,数值在Multiwfn当中已经输出过了。
如果极值点比较密,可以同时绘制箭头避免混乱。
如果箭头或者文字叠加在分子表面上,为了让其边缘清楚好看,建议在图层混合选项中设定外发光。
最终图像如下所示
在J.Phys.Org.Chem.2013,26,473-483一文中,马兜铃酸的分子表面静电势图也是通过类似方法绘制的,比此文的例子更复杂,在这里一起贴出
这是马兜铃酸的静电势定量分布图,显然没有呋喃分布得均匀。
正值区域、静电势接近0的区域占大部分分子表面,但是也有不小面积静电势非常负,这主要是羧基和氨基的氧的明显负电荷导致的。
5总结
本文介绍的绘制分子表面静电势图的步骤比较琐碎,需要很多手动操作,但这也是好处,使得我们可以很精细、随意地调节显示方式,而不受制于可视化程序所支持的选项。
本文介绍的只是一般过程,建议大家多摸索以使显示效果更好。
利用本文相同的步骤,可以绘制各种各样的实空间函数在分子表面的分布。
例如绘制平均局部离子化能的图,只要在Multiwfn的定量分子表面分析功能当中选选项2,然后再选2Averagelocalionizationenergy,之后再选0启动分子表面分析,随后的步骤和前文都一样。
Fukui函数、双描述符之类实空间函数并没有出现在选项2给出的列表里,但是依然可以通过定量分子表面分析功能对它们在分子表面上的分布进行分析,并且随后结合VMD绘制成填色图,只不过操作过程稍微特殊一些,见3.2(dev)版手册4.12.4节。
如果对分子表面分析有更多兴趣,可以参看《使用Multiwfn的定量分子表面分析功能预测反应位点、分析分子间相互作用》(
6其它:
分析范德华表面穿透程度
将Multiwfn和VMD相结合,充分开动脑筋,还可以做很多有趣、有用的分析。
例如,对水的二聚体中两个水分别做定量分子表面分析,将得到的两套表面顶点都放到VMD里同时显示出来,就可以看到如下图像
从图中可以看到两个水的范德华表面在形成二聚体后相互有明显穿透,并且是静电势明显为正(红色)区域与静电势明显为负(蓝色)区域相互穿透,这表明水二聚体间的氢键作用本质主要是由于静电吸引所导致的。
我们将视图放大,并且不让左边的水显示出来(因为它遮挡了表面顶点),然后找出两个位置比较能反映穿透程度的表面顶点,按2键,然后分别点击这两个点,这两个表面顶点间的距离就出现在图上了,如下图所示,其距离是1.12埃。
这表明形成二聚体后,范德华表面总共被穿透了2*1.12埃,这叫做mutualpenetration距离,其值越大通常说明弱相互作用越强。
Mutualpenetration距离其实有不同具体计算方法,最简单粗糙的是直接拿两个原子的范德华半径和减去原子间距。
而本文这种计算方法无疑是最准确的,因为精确反映了范德华表面的实际形状.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 使用 Multiwfn 结合 VMD 分析 绘制 分子 表面 静电 分布