 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 15 | 1/1 | 返回列表 |
| 查看: 7318 | 回復: 14 | ||||||||
| 【獎勵】 本帖被評價11次，作者sobereva增加金幣 8.8 個 | ||||||||
[資源]
使用Multiwfn結合VMD分析和繪制分子表面靜電勢分布
|
||||||||
|
使用Multiwfn結合VMD分析和繪制分子表面靜電勢分布 1 前言 分子表面靜電勢圖經常在文獻中出現，不同表面區(qū)域靜電勢大小通過不同顏色展現，使分子表面上靜電勢的分布一目了然。分子表面一般都用Bader定義的范德華表面，即電子密度為0.001 e/Bohr^3的等值面。只要提供所需的輸入數據，這種圖在許多程序中都可以作，例如Molekel。特別是在常用的Gaussview里作這種圖比較簡單快捷。本文介紹的通過Multiwfn的定量分子表面分析功能結合VMD和photoshop作分子表面靜電勢圖既不比使用Gaussview快，步驟也比它復雜，但是優(yōu)點十分顯著，就是可以在分子表面上顯示出靜電勢極值點位置，可以十分靈活地調節(jié)顯示效果，而且還可以通過相同的方法繪制出分子表面上靜電勢以外的實空間函數的分布，比如平均局部離子化能、局部電子親和能、Fukui函數等等。另外還可以順帶著獲得許多其它信息，比如原核與分子表面的距離，實空間函數在分子表面上分布的統計數據等等。總之，對于初學者、純粹圖省事者建議用Gaussview，本文是給那些想研究得更深入并且稍有動手能力的研究者讀的。只要了解每一步操作的意義，就會覺得其實根本不復雜，而且思想上獲得極大的解放。 本文所用Multiwfn版本為3.2(dev)（即3.2的開發(fā)版），可在http://multiwfn.codeplex.com上下載。VMD為1.9版，可在http://www.ks.uiuc.edu/Research/vmd/上免費下載。Photoshop為CS2版。本文將以一個簡單體系呋喃為例進行說明，波函數文件在B3LYP/6-31G**下由Gaussian產生。雖然本文主要是介紹作圖方法，但也會順帶著說一下定量分子表面分析的一些功能。 2 在Multiwfn中計算、統計、導出數據 啟動Multiwfn后依次輸入 furan.wfn 12 //定量分子表面分析 0 //開始計算 默認就是在電子密度為0.001的表面上計算靜電勢，并且格點間距為0.25Bohr。格點間距可以通過選項3來設定，格點間距越小，分子表面就會用越多的頂點來描述，定量統計值、極值點位置也會越精確，之后作出的分子表面填色圖的色彩過渡也越光滑，但是計算耗時將越長。如果你的目的僅僅是繪制分子表面靜電勢圖，可以忽略以下內容而直接跳到本節(jié)最后一段。 計算完畢后會看到分子表面上靜電勢極大、極小點的坐標和數值，并且輸出大量統計數據，比如分子表面上靜電勢的最大/最小值、平均值、方差、電荷平衡度、表面積等等，它們對于了解分子特征、建立QSPR/QSAR方程預測分子理化性質和生物活性等問題都十分有用，見Multiwfn手冊3.15.1節(jié)以及IJQC,85,676、JPCA,110,1005等文章的介紹。 此時會看到后處理菜單。點擊0進入圖形界面，并且將Ratio of atomic size設為4.0就可以清楚看到分子表面上靜電勢極大點（紅點）和極小點（藍點）的位置。如果點擊Minimal/Maximum label，就會顯示出極值點的編號，可以和命令行窗口顯示的極值點信息相對照，得知極值點上的具體數值。點擊右上角Return可以關閉圖形窗口。 此例中各個極值點靜電勢數值如下（由于格點精度有限，所以等價的極值點的數值不完全一致，這里取了平均） 極小點1（最小點）：-20.60 kcal/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 雖然是極大點，但是靜電勢數值為負，所以化學意義不大，可無視之 我們可以查看一下不同靜電勢區(qū)間內分子表面積，這對于了解分子表面靜電勢的定量分布很有益。做法是依次輸入 9 -25,22 //統計范圍。當前體系分子表面靜電勢范圍為-20.60~18.41kcal/mol，這里將范圍稍微擴大并取整數來定義范圍 15 //將-25~22kcal/mol均勻分為15個區(qū)間獲得表面積 3 //輸入的單位為kcal/mol 立刻屏幕上就輸出了不同靜電勢區(qū)間內的表面積，我們將Center這一列和Area這一列的數據分別從屏幕上拷下來（不會拷的讀者見手冊5.4節(jié)），并粘貼到諸如origin等作圖工具里并作成條形圖，就可以得到諸如以下圖像 靜電勢的定量分布在這張圖上一目了然。此分子不同靜電勢區(qū)間內的表面積分布還是相對比較均勻的。 如果選擇選項11，則程序會輸出每個原子附近的分子表面上的靜電勢統計數據，對于了解原子在此分子中的特征很有用。例如靜電勢平均值 Atom# All/Positive/Negative average 1 -10.08997 NaN -10.08997 2 -12.59141 NaN -12.59141 3 -12.64250 NaN -12.64250 4 -10.04843 NaN -10.04843 5 -13.54925 NaN -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），這也容易理解，畢竟其孤對電子對靜電勢有很大的負貢獻。由于其附近分子表面上沒有正值區(qū)域，所以正值區(qū)域的靜電勢平均值顯示的是NaN。呋喃中碳原子裸露在分子表面上的區(qū)域主要就是體現pi電子特征的區(qū)域，由于pi電子云使這部分區(qū)域靜電勢為負，所以碳原子附近靜電勢平均值也都為明顯負值，并且無正值區(qū)域。同時從靜電勢平均值上也體現出beta碳（2、3號）比alpha碳（1、4號）的pi電子云更富集因而靜電勢更負。這種討論分子表面上對應不同原子區(qū)域的定量數據的方法和分子表面極值點分析往往會得到共通的結論，但此方法可以得到更多的信息，比如alpha碳上沒有出現極值點（因此對它沒有任何描述），而通過分析它對應的分子表面上的局部區(qū)域的靜電勢平均值等數據就可以定量考察它的特征。這種分析方法是筆者獨創(chuàng)的，也僅有Multiwfn支持。另外，通過選項12還可以查看分子表面上對應于指定分子片段區(qū)域的定量性質。 選過選項9之后會問你是否輸出locsurf.pdb文件，通過此文件可以利用VMD查看不同原子對應的分子表面區(qū)域。由于這不是本文的重點，所以輸入n不讓程序輸出。 利用選項10可以查看分子表面相對于指定坐標點或指定原子核的最遠和最近距離，這對于討論分子間相互作用導致分子表面的穿透距離很有用。此功能也可以用于計算分子半徑或直徑，見《談談分子半徑的定義和計算方法》（http://hi.baidu.com/sobereva/item/3c9a79cac80a397588ad9e6e）的討論。 現在選擇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主窗口里(VMD Main)，它們的ID將分別是0、1、2。選擇Display-Depth Cueing將它關掉，否則圖像會有些朦朧。進入Graphics-Colors，選Display-Background-White將背景改為白色，并且在此界面的Color-Scale標簽頁里選擇BWR，使分子表面的色彩根據數值范圍由小到大以 藍-白-紅 的方式變化。選Display-Axes-Off不讓坐標軸顯示出來。 注：如果希望每次啟動VMD時都自動做如上操作，可以在VMD目錄下vmd.rc文件末尾添加以下4行 display depthcue off color scale method BWR color Display Background white axes location Off 進入Graphics-Representations，然后執(zhí)行以下3步來分別設定分子結構、分子表面極值點、分子表面頂點的顯示方式。（如果只是想簡單地看一下表面靜電勢分布，只要做第3步即可） (1)在Selected Molecule一欄里選擇furan.pdb，Drawing Method選Licorice，Bond Radius減小到0.1。 (2)將Selected Molecule一欄切換到surfanalysis.pdb，在控制臺輸入mol modstyle 0 1 VDW 0.06 （0和1是顯示方式編號和體系的ID，此命令代表用大小為0.06的VDW球顯示。由于GUI中VDW球最小只能設到0.1，故這里用命令行來實現）。然后在Selected Atoms里輸入carbon并回車，然后將Coloring Method選為ColorID，并且在右邊新出現的框里選Orange2。此時分子表面極大點就通過橙色小圓球顯示出來了。點擊Create Rep按鈕創(chuàng)建新顯示方式，在Selected Atoms里輸入oxygen并回車，然后將ColorID右邊的框設為Cyan，此時分子表面極小點就通過青色圓球顯示出來了。 (3)在Selected Molecule一欄里選擇vtx.pdb，Drawing Method選Points，Size設為25（設多大合決于視角的遠近，在當前視角下應當讓size恰好足夠大，使分子表面上的頂點緊密相連，不留明顯空隙），Coloring Method選Beta（根據pdb文件里B因子那一列的數據，此例即靜電勢數值進行填色），在Trajectory標簽頁里將Color Scale Data Range填上-22和22并點擊Set，代表色彩刻度設為-22~20kcal/mol（其實此例用默認的色彩刻度范圍就可以，這里只是為了取個整）�，F在分子表面填色圖就出現了。越藍的區(qū)域靜電勢越負，越紅的區(qū)域越正，白色區(qū)域的靜電勢數值在0附近。 之后給圖上加上色彩刻度軸。選Extensions-Visualization-Color Scale Bar，Color bar width設為0.08，Display title選on并且將Color bar title里寫上ESP (kcal/mol)，Minimum和Maximum scale value分別填-22和22，Number of axis labels輸入10，Color labels選Black，Label format選Decimal。然后點Draw Color Scale Bar按鈕，色彩刻度就出現在畫面中了，并且VMD Main窗口中多出了一個名為Color Scale Bar的一項。然后調整它的大小和位置，即雙擊VMD Main窗口中Color Scale Bar那一項當中的F標簽使之變?yōu)榧t色（即不讓色彩刻度軸在畫面中的位置凍結），而雙擊其它項目的F標簽使它們的F變?yōu)楹谏ㄗ屗鼈兊奈恢脙鼋Y�。�。然后激活VMD圖形窗口，按t鍵進入平移模式，然后拖動鼠標將色彩刻度軸放置到合適位置，并且用鼠標滾輪調整它的大小。調合適之后再按r鍵恢復旋轉視角模式，并且在VMD Main里將Color Scale Bar那一項的F重新雙擊成黑色，而其它三項的F重新雙擊為紅色。 當前的顯示效果如下 此圖還有許多地方值得進一步調整，也就是讓分子結構顯示出來（現在都被表面頂點掩蓋了）、給分子表面極值點標上具體數值、讓表面極值點更清楚地顯示，特別是藏在分子表面后方的極值點也顯示出來。為達到這些目的，需要利用photoshop。由于步驟比較細碎，筆者不把所有具體操作都敘述一遍，否則太羅嗦。只要會基本的photoshop操作的人都應該能理解應該怎么實現。 4 通過Photoshop改進作圖效果 將分子調整到一個合適的角度，然后在VMD main窗口里把所有條目的F標簽都雙擊成黑色來將它們固定住，以免隨后的操作過程中不慎旋轉了體系。 在VMD main窗口里面雙擊與furan.pdb和surfanalysis.pdb對應的條目的D標簽使其變紅，此時窗口內就只有分子表面和色彩刻度軸顯示了出來。然后按Alt+Printscreen鍵將窗口截圖，在photoshop里按Ctrl+N然后按OK，再用Ctrl+V把截的圖粘貼進去。此時圖像大小正好和VMD窗口大小完全一致。 在VMD main窗口里面只讓分子結構顯示出來，將背景改為藍色（只要不是白色就行，否則會和氫原子的白色連在一起），然后將窗口截圖并且粘貼進之前的ps窗口里成為新的圖層。選擇魔棒工具，Tolerance設0，Contiguous的對勾取消，然后點擊圖中藍色區(qū)域把背景區(qū)域選上，之后按delete鍵去除背景。之后將圖層的不透明度(Opacity)改為40%，這樣分子結構就會以半透明方式疊加到分子表面圖上了，而且疊加的位置完全精確。ps窗口里的圖像目前如下所示 讓VMD窗口里只顯示出表面極值點，然后將窗口截圖也粘貼進ps里成為新圖層，同樣將藍色的背景選中并刪掉。然后用矩形選框工具恰當地把圖像主體和色彩刻度軸圈上，構圖滿意后點Ctrl+Q將多余的區(qū)域裁掉，目前圖像如下 可見目前這些表面極值點無論是正面還是背面的都是以完全不透明方式展現，看不出前后。因此我們把那些明顯處在背面的極值點都用矩形選框圈上，注意圈的時候一直按著shift鍵使選區(qū)依次累加。然后在圖上點右鍵選Layer via Cut，這樣處在圖中正面（或者邊緣）的極值點與處在圖中背面的極值點就分別用兩個圖層來儲存了。將后者的圖層的不透明度設為50%。 最后，在圖上用文本工具標上一部分極值點的靜電勢數值，數值在Multiwfn當中已經輸出過了。如果極值點比較密，可以同時繪制箭頭避免混亂。如果箭頭或者文字疊加在分子表面上，為了讓其邊緣清楚好看，建議在圖層混合選項中設定外發(fā)光。最終圖像如下所示 在J. Phys. Org. Chem. 2013, 26, 473-483一文中，馬兜鈴酸的分子表面靜電勢圖也是通過類似方法繪制的，比此文的例子更復雜，在這里一起貼出 這是馬兜鈴酸的靜電勢定量分布圖，顯然沒有呋喃分布得均勻。正值區(qū)域、靜電勢接近0的區(qū)域占大部分分子表面，但是也有不小面積靜電勢非常負，這主要是羧基和氨基的氧的明顯負電荷導致的。 5 總結 本文介紹的繪制分子表面靜電勢圖的步驟比較瑣碎，需要很多手動操作，但這也是好處，使得我們可以很精細、隨意地調節(jié)顯示方式，而不受制于可視化程序所支持的選項。本文介紹的只是一般過程，建議大家多摸索以使顯示效果更好。利用本文相同的步驟，可以繪制各種各樣的實空間函數在分子表面的分布。例如繪制平均局部離子化能的圖，只要在Multiwfn的定量分子表面分析功能當中選選項2，然后再選2 Average local ionization energy，之后再選0啟動分子表面分析，隨后的步驟和前文都一樣。Fukui函數、雙描述符之類實空間函數并沒有出現在選項2給出的列表里，但是依然可以通過定量分子表面分析功能對它們在分子表面上的分布進行分析，并且隨后結合VMD繪制成填色圖，只不過操作過程稍微特殊一些，見3.2(dev)版手冊4.12.4節(jié)。 如果對分子表面分析有更多興趣，可以參看《使用Multiwfn的定量分子表面分析功能預測反應位點、分析分子間相互作用》（http://hi.baidu.com/sobereva/item/6601fff28c4fd9d643c36acc）。 6 其它：分析范德華表面穿透程度 將Multiwfn和VMD相結合，充分開動腦筋，還可以做很多有趣、有用的分析。例如，對水的二聚體中兩個水分別做定量分子表面分析，將得到的兩套表面頂點都放到VMD里同時顯示出來，就可以看到如下圖像 從圖中可以看到兩個水的范德華表面在形成二聚體后相互有明顯穿透，并且是靜電勢明顯為正（紅色）區(qū)域與靜電勢明顯為負（藍色）區(qū)域相互穿透，這表明水二聚體間的氫鍵作用本質主要是由于靜電吸引所導致的。 我們將視圖放大，并且不讓左邊的水顯示出來（因為它遮擋了表面頂點），然后找出兩個位置比較能反映穿透程度的表面頂點，按2鍵，然后分別點擊這兩個點，這兩個表面頂點間的距離就出現在圖上了，如下圖所示，其距離是1.12埃。這表明形成二聚體后，范德華表面總共被穿透了2*1.12埃，這叫做mutual penetration距離，其值越大通常說明弱相互作用越強。 注：Mutual penetration距離其實有不同具體計算方法，最簡單粗糙的是直接拿兩個原子的范德華半徑和減去原子間距。而本文這種計算方法無疑是最準確的，因為精確反映了范德華表面的實際形狀. |
回復此樓» 收錄本帖的淘帖專輯推薦
 量化軟件學習 | SOB 集錦 | 氣態(tài)分子化學動力學 | 量化軟件-Multiwfn |
| 和畢設有關的 | 路哥的淘貼 |
» 猜你喜歡
高效微納米儲能材料與產業(yè)化推進實驗室 蘇州 招收碩士研究生（含調劑）
已經有0人回復
-
求催化作用基礎第三版 （甄開吉) 課件
已經有3人回復
-
物理化學論文潤色/翻譯怎么收費?
已經有182人回復
-
26年材料/化學相關專業(yè)博士研究生招生 溫州大學化學與材料工程學院(26年秋季入學)
已經有0人回復
-
溫州大學新能源材料微觀溯源課題組招收2026級博士研究生（材料、化學方向）
已經有6人回復
-
材料與化工 323 英一+數二+物化，一志愿：哈工大 本人本科雙一流
已經有5人回復
-
26申博
已經有2人回復
-
求助電催化劑制備，怎么在碳布表面生長PdCo合金
已經有0人回復
-
求助電催化劑制備，怎么在碳布表面生長PdCo合金
已經有0人回復
-
304求調劑
已經有3人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 如何用Multiwfn軟件做出P DOS?
已經有4人回復
- NBO電荷比較親電活性大小
已經有8人回復
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結構及相對論效應的影響
已經有17人回復
- 靜電勢為什么有正負值？
已經有102人回復
- Multiwfn的初步學習，越詳細越好。
已經有6人回復
- 使用Multiwfn繪制NBO及相關軌道
已經有30人回復
- 用Multiwfn、Gaussview、Molekel、VMD觀看龍蝦、盆景、骨盆、大腦
已經有16人回復
- 引用了multiwfn的論文�，F在multiwfn軟件已經有19次引用率了呢
已經有13人回復
- 使用Multiwfn做拓撲分析以及計算孤對電子角度
已經有18人回復
- 回答一些關于Multiwfn的疑問以及未來Multiwfn的發(fā)展打算
已經有37人回復
- multiwfn輸出圖形只顯示一半
已經有3人回復
- 談談Molekel做靜電勢填色等值面圖的方法及誤區(qū)
已經有50人回復
3樓2013-07-28 10:27:28
|
本帖內容被屏蔽 |
4樓2013-07-28 14:50:23
5樓2013-07-29 06:37:24
6樓2013-07-29 08:52:48
8樓2013-08-12 08:38:47
10樓2014-10-19 22:09:33
11樓2015-04-05 10:36:44
12樓2015-12-16 15:23:26
13樓2015-12-16 23:17:56
14樓2015-12-16 23:18:49
15樓2020-06-14 22:32:03
簡單回復
安德2樓
2013-07-28 09:45
回復
五星好評 
windowtt7樓
2013-07-29 09:47
回復
五星好評 牛
zy543829樓
2013-08-31 16:34
回復
五星好評 
相關版塊跳轉
我要訂閱樓主 sobereva 的主題更新
| 15 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調劑 +3 | 睿08 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 290求調劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 293求調劑 +5 | 世界首富 2026-03-11 | 5/250 |
|
|
|
[考研] 278求調劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 考研化學學碩調劑，一志愿985 +3 | 張vvvv 2026-03-15 | 5/250 |
|
|
[考研] 機械專碩325，尋找調劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 070300化學學碩求調劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 0703 物理化學調劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 0856求調劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 289求調劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 0805，333求調劑 +3 | 112253525 2026-03-10 | 3/150 |
|
|
|
[考研] 337一志愿華南理工0805材料求調劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
|
[考研] 290求調劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
|
[考研] 工科，求調劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
|
[考研] 【考研調劑求收留】 +3 | Ceciilia 2026-03-11 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 295求調劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 化工學碩306求調劑 +9 | 42838695 2026-03-12 | 9/450 |
|
|
|
[考研] 270求調劑 085600材料與化工專碩 +3 | YXCT 2026-03-11 | 3/150 |
|
|
|
[考研] 298求調劑 +3 | Vv呀！ 2026-03-10 | 3/150 |
|