pymol 基本操作

1、精選文檔簡介&安裝Pymol是一個開放源碼，由使用者贊助的分子三維結(jié)構(gòu)顯示軟件，由Warren Lyford DeLano編寫，并且由DeLano Scientific LLC負(fù)責(zé)商業(yè)發(fā)行。Pymol被用來創(chuàng)作高品質(zhì)的分子（特殊是生物大分子如蛋白質(zhì)）三維結(jié)構(gòu)。據(jù)軟件作者宣稱，在全部正式發(fā)表的科學(xué)論文中的蛋白質(zhì)結(jié)構(gòu)圖像中，有四分之一是使用Pymol來制作的。Pymol名字的來源：“Py”表示該軟件基于python這個計算機語言，“Mol”則是英文分子（molucule）的縮寫，表示該軟件用來顯示分子結(jié)構(gòu)。由于試驗需要，本人正在學(xué)習(xí)該軟件，在這里把學(xué)習(xí)過程記錄下來，期望對有需要的伴侶有所掛念。今日

2、先來說說安裝吧。自2006年8月1日起，DeLano Scientific 對事先編譯好的PyMOL執(zhí)行程序（包括beta版）實行限定下載的措施。目前，只有付費用戶可以取得。不過源代碼目前還是可以免費下載，供使用者編譯。假如你和我一樣，不想為此花錢的話： 1. 假如你是Windows用戶，首先下載Pymol的源代碼。然后安裝CygWin，并且確保正確安裝以下模塊： C+ (gcc or g+ package name) Python OpenGL PNG 然后在源代碼名目里面依次運行： 2. 假如你是Linux用戶，首先確保以下東東已安裝： Python Pmw OpenGL driver（我

3、用的是NVdia） libpng Subversion client（下載源代碼需要） 然后下載Pymol的源代碼$ mkdir pymol-src$ svn co pymol-src然后進入源代碼名目# cd pymol-src開頭依次編譯# python setup.py install# python setup2.py install拷貝執(zhí)行腳本到某個$PATH，安裝就搞定了# cp ./pymol /usr/bin假如運行時得到錯誤信息ImportError: No module named Pmw，那么你應(yīng)當(dāng)運行# python setup2.py install pmw假如你在使

4、用Gentoo，請確保編譯python時添加了tcl/tk支持，否則運行是會提示錯誤ImportError: No module named _tkinter# USE=tcl tk emerge python好了，下面我們就可以進入Pymol的世界了。 基本的鼠標(biāo)操作里主要介紹一下Pymol的基本操作，包括窗口菜單、加載文件、圖像的基本鼠標(biāo)操作等等。 當(dāng)你打開Pymol后，你將會看到如下圖所示的界面：該界面分為2窗口，上面的外部GUI窗口（External GUI）和下面的Viewer Window。Viewer Window又分為左右兩塊，左邊用來顯示結(jié)構(gòu)圖像的（Viewer），右邊則是一

5、個內(nèi)部GUI窗口（Internal GUI)。Viewer自身包含一個命令行（如圖中左下方的PyMOL提示符），可以用來輸入Pymol命令；在Inernal GUI中則可以選定一些特定的對象并完成一些操作。External GUI則包含一個標(biāo)準(zhǔn)菜單、一個輸出區(qū)、一個命令行輸入?yún)^(qū)以及右邊的一些常用命令按鈕。請留意，標(biāo)準(zhǔn)的“復(fù)制、剪切和粘貼”操作只能在External GUI中完成，并且必需使用“CtrlC、CtrlX以及CtrlV”來完成，這也是這個所謂的外部GUI的最重要的優(yōu)點。加載文件，有二種方法：1. 在External GUI中選擇File Open 2. 使用命令行： load 例如我

6、們現(xiàn)在從上下載了一個離子通道蛋白的pdb文件（PENTAMERIC LIGAND GATED ION CHANNEL FROM ERWINIA CHRYSANTHEMI），名字為2vl0.pdb，然后用pymol打開它： load 2vl0該蛋白質(zhì)的結(jié)構(gòu)就被顯示出來啦，如下圖：基本的圖像操作:是不是覺得上面的那個三維結(jié)構(gòu)圖看起來亂七八糟的阿，那是由于蛋白質(zhì)分子都是由成千上萬個原子組成的，而Pymol打開pdb文件時是默認(rèn)把全部的原子都顯示在那個小小的Viewer窗口里面的，當(dāng)然看起來就很亂了。這時候就需要我們對這個圖像進行一些操作，來得到秀麗的清楚的蛋白質(zhì)三維結(jié)構(gòu)圖。先說

7、說鼠標(biāo)吧。 任意旋轉(zhuǎn)圖像： 對準(zhǔn)圖像的任意處點住鼠標(biāo)左鍵然后移動鼠標(biāo)。 放大/縮小圖像： 對準(zhǔn)圖像的任意處點住鼠標(biāo)右鍵然后移動鼠標(biāo)：向上是縮小，向下則是放大。 移動圖像： 對準(zhǔn)圖像的任意處點住鼠標(biāo)中鍵或者滾輪，然后移動鼠標(biāo)。 設(shè)定圖像旋轉(zhuǎn)中心： CtrlShift鼠標(biāo)中鍵或滾輪。 移動剪切平面： Shift鼠標(biāo)右鍵。鼠標(biāo)上下移動：調(diào)整前剪切平面（離你近的）；鼠標(biāo)左右移動：調(diào)整后剪切平面（離你遠(yuǎn)的）。 最終一項“移動剪切平面”有點不簡潔理解，需要多試幾次。協(xié)作下面的示意圖你會發(fā)覺Pymol的這項設(shè)定其實很便利。今日沒時間了，明天還要出遠(yuǎn)門，就學(xué)到這里吧，用下面這個圖作為結(jié)束，其實就是用carto

8、on的形式顯示了上面的那個蛋白質(zhì)，不過還比較難看。By wei luPyMOL用法（教程二）基礎(chǔ)Pymol命令這里主要介紹一下Pymol的一些基本命令操作。就像Linux一樣，要想更好的操作Pymol，把握一些常用的命令是必不行少的。 Pymol是區(qū)分大小寫的，不過目前為止Pymol還是只用小寫，所以記住，全部的命令都是使用小寫字母的。當(dāng)你開頭用Pymol來完成一個項目時，你或許想會讓Pymol自動保存你全部輸入過的命令，以便利日后你再次讀取并修改。這個可以通過創(chuàng)建一個log文件來達(dá)到，該文件的后綴名應(yīng)為.pml，記住，Pymol像Linux一樣支持Tab鍵命令補全：Pymol log_ope

9、n log-file-name.pml假如你想終止記錄，只需要鍵入：Pymol log_close好了，現(xiàn)在載入pdb文件（連續(xù)前用的pdb文件）：Pymol load 2vlo.pdb現(xiàn)在Pymol就創(chuàng)建了一個叫2vlo的對象，你可以在內(nèi)部GUI窗口里面觀察這個項目的名字。但是你也可以自己定義該項目的名字（如test）：Pymol load 2vlo.pdb, test下面說說如何來操作你新建的對象。首先：Pymol show representationPymol hide representation其中representation可以為：cartoon, ribbon, dots, s

10、pheres, surface和mesh。使用這2個命令可以讓Pymol以不同的方式顯示蛋白質(zhì)結(jié)構(gòu)。例如當(dāng)我們鍵入：Pymol hide linesPymol show ribbon我們將得到如下結(jié)果：或許你已經(jīng)留意到結(jié)構(gòu)中有2個一模一樣的蛋白質(zhì)分子，只是方向不同而已，那么如何讓Pymol只顯示當(dāng)中的一個分子呢？首先輸入如下命令：Pymol label all, chains這個命令的作用是讓Pymol給蛋白質(zhì)結(jié)構(gòu)中的“鏈”編號，你會發(fā)覺，第一個分子由“鏈”AE組成，其次個則由FJ組成。好了，假如我們想把一個蛋白質(zhì)分子去掉，那么只要把“鏈”AE或者FJ去掉即可：Pymol hide ribbo

11、n, chain f+g+h+i+j上面的東東還可以這樣完成：Pymol select test, chain f+g+h+i+jPymol hide ribbon, test上面的第一句命令的作用是選擇“鏈”FJ，并命名為test，然后在其次句命令中隱蔽它。這樣做的好處是，一旦你選擇并命名了某個目標(biāo)，你可以在后面隨時對它進行各種操作。并且你在右邊的把握面板里面也可以看到你選定的目標(biāo)，并可以對其進行操作。比如你可以：Pymol hide everything, testPymol show cartoon, test這樣你會得到：說到這里就提到了Pymol的一個比較重要的東東，就是選擇并命名目

12、標(biāo)，它的基本語法就是：Pymol select selection-name, selection-expression其中名字可以由字母A/aZ/z，數(shù)字09已經(jīng)下劃線_組成，但是要避開使用：! # $ % & * ( ) | ? /假如你要刪除你選定的目標(biāo)或者整個對象，你可以：Pymol delete selection-namePymol delete object-name下面講講如何給對象以及目標(biāo)轉(zhuǎn)變顏色。預(yù)定義的顏色名字可以在外部GUI窗口的Settings Colors中找到：Pymol color color-namePymol color color-name, select

13、ion-expression比如我們可以：Pymol color red, ss hPymol color yellow, ss sPymol color green, ss l+其中“ss”代表secondary structure，“h”代表Helix，“s”代表Beta sheet，l+代表Loop和所以其他結(jié)構(gòu)。這3句的作用分別是把全部的Helix變成紅色；把全部的Beta sheet變成黃色；把全部的Loop以及其他部分變成綠色，于是我們得到：Pymol可以同時打開多個pdb文件：Pymol load object-name-1.pdbPymol load object-name-2

14、.pdb假如你想臨時關(guān)閉/打開某個對象，可以這樣：Pymol disable object-name-1Pymol enable object-name-1你也可以用disable命令去除最終一個選擇的目標(biāo)上消滅的粉紅色的小點，但是該命令并不會使你選定的目標(biāo)不行見。Pymol disable selection-name使用鼠標(biāo)通常是轉(zhuǎn)變圖像視角的最便利的方法，不過命令如zoom，orient等等有時候使用起來也是很有用的，它們供應(yīng)了另一種轉(zhuǎn)變圖像視角的方法。放大選定目標(biāo)：Pymol zoom selection-name定向選定目標(biāo)，可以使選定目標(biāo)最大的尺寸水平顯示，次大的尺寸豎直顯示：Py

15、mol orient selection-name你也可以用view命令保存你目前的視角，留意，該命令只保存視角，并不保存你的對象顯示方式：Pymol view key, action其中“key”是你任憑給當(dāng)前視角定的名字，“action”可以為：store或者recall。假如不加任何“action”，則默認(rèn)為recall：Pymol view v1, storePymol view v1, recallPymol view v1說了這么多，最終說說如何保存文件吧。Pymol有3個層面的保存方式，下面來分別說說。1. 使用log_open命令把你全部使用過的命令記錄為一個文本文檔： Pym

16、ol log_open script-file-name這樣以后當(dāng)你再次調(diào)用該文檔時，Pymol將執(zhí)行上面的全部命令：Pymol script-file-name不過留意，假如你想記錄當(dāng)前視角，則必需使用get_view命令。你可以選擇外部GUI窗口中的File Append/Resume/Close Log來分別暫停記錄/恢復(fù)記錄/停止記錄該文檔。你可以隨時編輯該文檔。在linux下，該文檔的默認(rèn)保存名目為當(dāng)前用戶的home名目。2. 假如你想下次打開Pymol時直接回到當(dāng)前所在的狀態(tài)，那么你可以選擇外部GUI窗口里面的File Save Session，創(chuàng)建一個會話文件（.pse）。 該會

17、話文件和上面提到的文檔文件的區(qū)分在于，首先文檔文件可以編輯，但會話文件不行以；記錄文檔文件前必需先運行l(wèi)og_open命令，而會話文件可以隨時創(chuàng)建；最終文檔文件以文檔形式運行（），而打開會話文件則必需選擇外部GUI窗口中的File Open。什么時候需要創(chuàng)建會話文件呢？比如當(dāng)你在某時有多個選擇時，你可以保存當(dāng)前狀態(tài)，然后一一嘗試這些選擇，不滿足時只需要重新打開該會話文件即可。也就是說創(chuàng)建會話文件起到了“undo”的作用，這正是Pymol所缺少的。期望開發(fā)者能趕快加入該功能，那么這個會話文件好像就沒什么大用了，呵呵。3. 假如你覺得當(dāng)前顯示窗口里面顯示的結(jié)構(gòu)圖像已經(jīng)滿足你的要求了，你可以把它保存

18、為圖片。在這之前你可以使用ray命令來優(yōu)化你的圖像，它可以使你的圖像具有三維的反射及陰影特效： Pymol rayPymol pngyour_path/image_name最終就用該命令導(dǎo)出的圖片結(jié)束這次筆記吧。Pymol命令的語法與目標(biāo)選擇的表達(dá)上次介紹一些Pymol的基本命令。現(xiàn)在來具體說說Pymol命令的語法，還有在選擇操作目標(biāo)應(yīng)當(dāng)假如表達(dá)。個人覺得這部分內(nèi)容對學(xué)習(xí)Pymol來說是至關(guān)重要的。從上次講的一些例子中不難看出，Pymol的命令都是由關(guān)鍵詞（keyword）加上一些變量（argument）組成，格式如下：Pymol keyword argument其中關(guān)鍵詞（keyword）當(dāng)

19、然是必需的，而變量則不是必需的，比如退出命令quit就不需要附加變量：Pymol quit當(dāng)然更多的命令通常是需要加變量的，比如放大命令zoom，但是你會發(fā)覺即使你不加任何變量該命令也可以被執(zhí)行，這是由于Pymol的很多命令有一個默認(rèn)變量，下面兩個命令的作用是一樣的，其中的目標(biāo)選擇all就是zoom的默認(rèn)變量：Pymol zoomPymol zoom all還有些命令可以帶多個參數(shù)，比如加色命令color，它的用法如下：Pymol color color-namePymol color color-name, selection-expression第一個color雖然只帶一個變量color-

20、name，但其實它包含了其次個默認(rèn)變量all，所以它的作用是把整個結(jié)構(gòu)變成color-name的顏色。其次個color帶兩個變量，和第一個的區(qū)分就是把默認(rèn)的目標(biāo)選擇變量all變成了selection-expression，也就是說只有被這個變量選中的部分才會被變成color-name定義的顏色。要留意的是，假如一個命令帶多個變量，則這些變量之間必需用逗號,隔開。通過這個例子，大家可以發(fā)覺，有些變量本身是很簡潔的，比如color-name，就是一個顏色名字而已，沒什么簡單的。另一些則不一樣，比如selection-expression，它可以很簡潔，也可以格外的簡單。這個東東，我稱之為選擇表達(dá)，

21、對Pymol命令的使用格外重要，所以下面要具體的講一下。選擇表達(dá)（selection-expression）表示的實際就是一些被選中的部分，它們可以是一些個原子，一些個Helix，一些個Beta sheet，或者它們的混合物。你可以給你的選擇表達(dá)起個名字，以便可以多次使用它們。名字可以由大小寫字母，數(shù)字以及下劃線_組成，但是因避開使用下列符號：! # $ % & * ( ) | ? /選擇表達(dá)由所謂的selector加上identifier組成，其中selector定義了某類屬性，而identifier則在該類屬性下需要被選擇的部分。如下例：Pymol select test, name c+

22、o+n+ca其中name就是一個selector，它表示在pdb文件中描述的原子的名字；c+o+n+ca則是對應(yīng)的indentifier，它表示我們要選擇pdb文件中名字叫ca+cb的原子（ca代表alphacarbon，cb代表betacarbon）。整個語句的作用就是選擇上訴的原子并命名為test，這樣我們可以在后面連續(xù)使用它。下表列出了大多數(shù)的selector：Selector簡寫Identifier及例子symbole.chemical-symbol-list周期表中的元素符號Pymol select polar, symbol o+nnamen.atom-name-listpdb文件

23、中的原子名字Pymol select carbons, name ca+cb+cg+cdresnr.residue-name-list氨基酸的名字Pymol select aas, resn asp+glu+asn+glnresii.residue-identifier-listpdb文件中基團的編號Pymol select mults10, resi 1+10+100residue-identifier-rangePymol select nterm, resi 1-10altaltalternate-conformation-identifier-list一些單字母的列表，選擇具有2種構(gòu)型

24、的氨基酸Pymol select altconf, alt a+bchainc.chain-identifier-list一些單字母或數(shù)字的列表Pymol select firstch, chain asegis.segment-identifier-list一些字母（最多位）的列表Pymol select ligand, segi ligflagf.flag-nummer一個整數(shù)（）Pymol select f1, flag 0numeric_typent.type-nummer一個整數(shù)Pymol select type1, nt. 5text_typett.type-string一些字母（

25、最多位）的列表Pymol select subset, tt. HA+HCididexternal-index-number一個整數(shù)Pymol select idno, id 23ernal-index-number一個整數(shù)Pymol select intid, index 23sssssecondary-structure-type代表該類結(jié)構(gòu)的單字母Pymol select allstrs, ss h+s+l+下表是另一些Selector，有關(guān)比較的：Selector簡寫Identifier及例子bbcomparison-operator b-factor-valu

26、e一個實數(shù)，用來比較b-factorPymol select fuzzy, b 12qqcomparison-operator occupancy-value一個實數(shù)，用來比較occupancyPymol select lowcharges, q 0.5formal_parison-operator formal charge-value一個整數(shù)，用來比較formal chargePymol select doubles, fc. = -1partial_parison-operator partial charge-value一個實數(shù)，用來比較

27、partial chargePymol select hicharges, pc. -1另外有一些Selector是不需要Identifier的，它們被列在下表中：Selector簡寫描述all*全部當(dāng)前被Pymol加載的原子nonenone什么也不選hydroh.全部當(dāng)前被Pymol加載的氫原子hetatmhet全部從蛋白質(zhì)數(shù)據(jù)庫HETATM記錄中加載的原子visiblev.全部在被“可見”的顯示的對象中的原子presentpr.全部的具有定義坐標(biāo)的原子在Identifier中用到的原子以及氨基酸的命名規(guī)章可以在下面的網(wǎng)址中找到：/docs.html在選

28、擇表達(dá)中，selector還可以協(xié)作規(guī)律操作子（logical operator）使用，這樣我們可以表達(dá)更加簡單的選擇。這些操作子被列于下表中：Operator簡寫效果與例子not s1! s1選擇原子但不包括s1中的Pymol select sidechains, ! bbs1 and s2s1 & s2選擇既在s1又在s2中的原子Pymol select far_bb, bb & farfrm_tens1 or s2s1 | s2選擇s1或者s2中的原子（也就是包含全部的s1和s2原子）Pymol select all_prot, bb | sidechains1 in s2s1 in s

29、2選擇s1中的那些原子，其identifiers (name, resi, resn, chain, segi)全部符合s2中對應(yīng)的原子Pymol select same_atom, pept in prots1 like s2s1 l. s2選擇s1中的那些原子，其identifiers (name, resi)符合s2中對應(yīng)的原子Pymol select similar_atom, pept like prots1 gap Xs1 gap X選擇那些原子，其van der Waals半徑至少和s1的van der Waals半徑相差XPymol select farfrm_ten, res

30、i 10 gap 5s1 around Xs1 a. X選擇以s1中任何原子為中心，X為半徑，所包括的全部原子Pymol select near_ten, resi 10 around 5s1 expand Xs1 e. X選擇以s1中任何原子為中心，X為半徑，然后把s1擴展至該新的范圍所包含的全部原子Pymol select near_ten_x, near10 expand 3s1 within X of s2s1 w. X of s2選擇以s2為中心，X為半徑，并包含在s1中的原子Pymol select bbnearten, bb w. 4 of resi 10byres s1br.

31、s1把選擇擴展到全部residuePymol select complete_res, br. bbnear10byobject s1bo. s1把選擇擴展到全部objectPymol select near_obj, bo. near_resneighbor s1nbr. s1選擇直接和s1相連的原子Pymol select vicinos, nbr. resi 10這些規(guī)律選擇還可以組合使用。比如你想選擇chain b，但是不選擇其中的residue 88：Pymol select chain b and (not resi 88)在使用多重規(guī)律選擇時，為了讓Pymol正確處理挨次，請使用

32、括號，這樣最里層括號里面的內(nèi)容將會被最先處理，以此類推。好了，目標(biāo)選擇就先說到這里。其實關(guān)于目標(biāo)選擇還有所謂的“宏”可以用，可以簡化表達(dá)式，預(yù)備下次說說。by Wei L - PyMOL用法（教程三）Pymol的選擇宏上次具體講了如何在Pymol中怎么用selection-expression選取目標(biāo)，其實在某些狀況下，還可以用Pymol供應(yīng)的宏來選擇操作目標(biāo)。使用這個選擇宏往往可以是一個原本很簡單的表達(dá)變得簡潔緊湊。例如我們想選擇2vlo這個pdb文件中的chain a中的第100個基團的炭原子，假如用selection-expression來表達(dá)的話是這

33、樣：Pymol select chain a and resi 100 and ca假如用宏的，可以這樣：Pymol select a/100/ca是不是覺得簡潔了很多。好了，下面就來具體講講這個宏吧。由于這個宏是用來選擇目標(biāo)的，所以我稱之為選擇宏，它用斜杠/來定義Identifier，并且它使用上次介紹過的規(guī)律操作子and。一個完整的，按挨次的選擇宏的表達(dá)如下：/object-name/segi-identifier/chain-identifier/resi-identifier/name-identifier之所以說選擇宏是有挨次的，是由于Pymol就是靠挨次推斷每個斜杠后面的東東都是什

34、么東東。假如再細(xì)分一下的話，其實這個選擇宏有2種寫法，一個是帶開頭的斜杠，另一個是不帶開頭斜杠。區(qū)分是：假如不以斜杠開頭，那么Pymol則認(rèn)為你的表達(dá)式的最終一項就是選擇宏的末尾的最終一項，也就是name-identifier。例如：Pymol show lines, a/100/caPymol show lines, 100/ca如過以斜杠開頭，那么Pymol就認(rèn)為你是從選擇宏的表達(dá)式的頂端開頭的，也就是從/object-name開頭的。例如：Pymol zoom /2vl0/a/100/caPymol zoom /2vl0/a/100細(xì)心的讀者確定發(fā)覺了上面的例子中有兩道斜杠中間什么內(nèi)容也

35、沒有，不會是寫錯了吧？當(dāng)然不是，其實在這種狀況下Pymol會默認(rèn)選擇這個兩道斜杠中被省略的Identifier列表中的全部元素，也就是說被省略的部分被Pymol當(dāng)作了一個通配符。例如上例中我要選擇全部的segment，所以我就把它給省略不寫了，呵呵，便利吧。在舉些例子來說明一下：Pymol color green, a/142/斜杠后面的name-identifier被省略了，所以第142號基團的素有原子都會變成綠色。Pymol shwo cartoon, a/a斜杠后面的resi-identifier以及最終斜杠后面的name-identifier被省略了，所以整個a鏈將以cartoon的方

36、式被顯示。Pymol zoom /2vl0/b2個斜杠間的segi-identifier被省略，所以全部的b鏈將被放大。最終總結(jié)一下，Pymol的選擇宏必需至少包含一個斜杠/，以此來和Pymol的select-expression區(qū)分；并且不能包含空格，由于Pymol是把宏作為一個詞來讀取的；還有就是其實Pymol在執(zhí)行宏的時候首先是把它翻譯成正常的select-expression，然后再執(zhí)行的。關(guān)于cartooncartoon經(jīng)常被用來顯示一個蛋白質(zhì)的總體結(jié)構(gòu)，看起來也很秀麗。這次就來說說它的具體用法。不久前本人剛搞定了一個Glucosyltransferase的結(jié)構(gòu)，所以下面全部的例子都

37、用來它來說明。cartoon的命令格式如下：Pymol cartoon type, (selection)總結(jié)一下cartoon的顯示類型：automatic：默認(rèn)的顯示方式looptube: 比loop粗點putty: 這個比較好玩，依據(jù)R-factor來顯示，越高越粗ovalrectanglearrow：和rectangle幾乎一樣，就是多了個箭頭dumbbell：在oval的基礎(chǔ)上，在helix的邊緣加上一個cylinderskip：隱蔽，該圖中隱蔽了6120號氨基酸。下面說說如何設(shè)置cartoon的一些具體細(xì)節(jié)。比較下面的2幅圖：你會發(fā)覺第一張圖中sheets是平的，而當(dāng)中的那個氨基酸

38、的支鏈并沒連接在sheet上，這是由于為了顯示的秀麗，把sheet人為的抹平了。而其次張圖中的sheets則表達(dá)了蛋白質(zhì)的真實走向，所以氨基酸的支鏈也顯示正常。也就是說，假如你想表達(dá)某個局部的具體細(xì)節(jié)的時候，最好接受其次張圖中的顯示方式。2張圖對應(yīng)的命令分別是：Pymol set cartoon_flat_sheets, 1Pymol set cartoon_flat_sheets, 0類似的命令對應(yīng)于loop，就不舉例子了：Pymol set cartoon_smooth_loops, 1Pymol set cartoon_smooth_loops, 0下面再說說cartoon尺寸。Heli

39、x的厚度和寬度：Pymol set cartoon_oval_width, 0.2Pymol set cartoon_oval_length, 1.5sheet的厚度和寬度：Pymol set cartoon_rect_width, 0.5Pymol set cartoon_rect_length, 1.5loop的半徑：Pymol set cartoon_loop_radius, 0.2假如你設(shè)置了cartoon的顯示風(fēng)格為fancyPymol set cartoon_fancy_helices, 1Pymol set cartoon_fancy_sheets, 1這樣你得到的helix的邊

40、上會帶有一個很細(xì)的cylinder，也就是上面幾張圖中的顯示方式。此時設(shè)置helix的厚度，寬度，以及這個cylinder的半徑分別是：Pymol set cartoon_dumbbell_width, 0.1Pymol set cartoon_dumbbell_length, 2Pymol set cartoon_dumbbell, 0.2依此類推，還可以設(shè)置和putty，tube等等顯示類型相關(guān)的尺寸，就不一一類舉了。最終再加幾個還用的著的命令吧：上色：Pymol set cartoon_color, green竟然還可以refine，呵呵，逗號后面可以接數(shù)字，好像120都可以，數(shù)字越大優(yōu)

41、化的越大，感覺的確能變秀麗點：Pymol set cartoon_refine, 20設(shè)置透亮：Pymol set cartoon_transparency, 0.5關(guān)于cartoon還有些命令，感覺不怎么常用，有些我也不知道是干什么的。有愛好再爭辯吧。PyMOL用法（教程四）關(guān)于label在顯示一個蛋白結(jié)構(gòu)的某個細(xì)節(jié)的時候，經(jīng)常會需要給某些重要的氨基酸打上標(biāo)簽，這就需要用到label命令。label的命令格式如下：Pymol label selection, expressionselection當(dāng)然就是你要加標(biāo)簽的對象，expression就是標(biāo)簽的內(nèi)容，可選的有：name, resn,

42、resi, chain等等。你也可以組合使用它們。expression也可以是你自定義的一段內(nèi)容，這時候只要把內(nèi)容用引號包含起來就行：Pymol label selection, user-defined expression在下面這個例子中，我想把glucosyltransferase中UDP-Glucose的binding pocket標(biāo)注出來：首先說明一下，該pdb文件中A鏈?zhǔn)堑鞍踪|(zhì)，B鏈?zhǔn)荱DP-Glucose。Pymol load glucosyltransferase.pdb, tmpPymol extract upg, chain bPymol extract pro, chai

43、n aPymol delete tmpPymol select near, pro within 4.5 of upgPymol hide allPymol show sticks, upgPymol show lines, nearPymol label near, (%s/%s) % (resn, resi)#(%s/%s): 設(shè)定顯示格式。 上面的圖看起來有點亂，由于默認(rèn)Pymol在每個原子上都打上了標(biāo)簽。要想看起來順眼點，需要一點加工。在這之前，讓我們先看一下關(guān)于label的其他設(shè)置：投影模式，可選值0（無投影），1（object有投影到label上，但是label本身無投影），2（o

44、bject有投影到label上，label也有投影），3（object不投影到label上，label本身有投影）:Pymol set label_shadow_mode, 3文字顏色：Pymol set label_color, color-name, selection標(biāo)簽文字的輪廓的顏色，這樣就讓在例如白色背景上加白色標(biāo)簽成為了可能：Pymol set label_outline_color, color-name, selection字體，pymol內(nèi)置了12種字體，編號從516。15號和16號字體是unicode的：Pymol set label_font_id, 5字體大小，假如為

45、正值，則單位就是正常的px。你也可以用負(fù)值，則單位是：Pymol set label_size, -0.5Pymol set label_size, 4設(shè)置label位置，用下列命令可以設(shè)置label離默認(rèn)位置的三維偏移值，在需要給spheres加標(biāo)簽的時候有用：Pymol set label_position, (x,y,z)最終說說怎樣用單個字母標(biāo)注氨基酸。首先在$HOME/.pymolrc中加入：# start $HOME/.pymolrc modificationsingle =VAL:V, ILE:I, LEU:L, GLU:E, GLN:Q, ASP:D, ASN:N, HIS:H

46、, TRP:W, PHE:F, TYR:Y, ARG:R, LYS:K, SER:S, THR:T, MET:M, ALA:A, GLY:G, PRO:P, CYS:C# end modification用法，用singleresn代替resn：Pymol label n. CG and i. 230+246, singleresn下面是改進過的圖片：是不是看起來好多了，下面是具體的代碼，其中關(guān)于電子密度圖的設(shè)置請見下一節(jié)：Pymol select near, resi 139+229+230+233+246+499+519+520Pymol as cartoon, proPymol 鼠標(biāo)操作：顯示near的sidechainPymol set_color grey1, 224,224,224Pymol set cartoon_color, grey1Pymol set cartoon_transparency, 0.3Pymol set cartoon_fancy_helices, 1Pymol label n. CB and near, (%s%s) % (singleresn, resi)Pymol 進入Editing模