2010年5月31日 星期一

Docking

1. Prepare docking molecule  -->  d1.pdb  <-- download from http://www.drugbank.ca/

2. Prepare receptor molecule --> r1.pdb <-- download from http://www.ncbi.nlm.nih.gov/Structure/index.shtml

3. Automatic patch docking on http://bioinfo3d.cs.tau.ac.il/PatchDock/

4. Analysis by Pymol

2010年5月17日 星期一

SSD in Win7


1. hiberfil.sys 電腦休眠時的記憶體暫存檔案,將休眠功能取消後透過指令刪除
(1) 以系統管理員身分執行 [命令提示字元]
(2) 輸入 powercfg -h off 關閉休眠,按 Enter 鍵。假如要恢復的話請輸入 powercfg -h on 後按 Enter。
2. pagefile.sys 分頁檔,指定到其它硬碟
(1) 到 [控制台]、[系統及安全性]、[系統] 、 [進階系統設定]
(2) 出現 [系統內容] 視窗,切換至 [進階] 索引標籤
(3) 在 [效能] 群組中,點選 [設定] -->  [效能選項] --> [進階] 索引標籤
(4) 在 [虛擬記憶體] 群組中,按 [變更]。
3. 將IE的 cache 指定到 Ramdisk
4. 闗閉HD重組
5. 將系統 \tmp \temp 指定到 Ramdisk

新增簽名筆刷

1. 開啟簽名檔
2. 編輯--->定義筆刷預設集
3.按下確定-->新增在刷樣式清單



2010年5月13日 星期四

Phylogenetic analysis

Distance matrix methods:
1.    UPGMA --> 假設演化速率須一致才可用
2.    Neighbor-joining method --> 最廣泛應用,但會遺失序列特性, 全部轉為distance 數字畫樹
A.   Distance matrix 需經校正
校正方法:
l  Jukes-cantor (JC or JC69) : one rate of substitution, equal base frequencies
l  Kimura 2-Parameter (K2P or K80)  two types of substitution, equal base frequencies
l  Felsentein 1981(F81) : one rate of substitution, unequal base frequencies
l  Hasegawa, or Felsentein (HKY85 or F84) : two types of substitution, unequal base frequencies
l  Tamura and Nei (TN93) : tree types of substitution, unequal base frequencies
l  General time Reversible (GTR) : six types of substitution, unequal base frequencies
 à 若須精確計算(演化、流病、分子時鐘等可先跑 LRT (likelihood ratio test),決定選用的 model
à 一般採用 HKY85 , 鹼基頻率一致時採K2P即可
B.   Phylip 操作
                                               i.          Use [ClustalW method] run multiple sequence alignment à output format : phylip à rename output file as “infile”
                                             ii.          Use [dnadist.exe] run distance matrix table to get “outfile” à rename “outfile” as “infile”
                                           iii.          Run [neighbor.exe] to generate “outtree”
                                           iv.          View as [treeview]


Character state methods:
1.    Maximum parsimony method (最大簡約法) à 依據最簡約的演化步驟所得演化樹,以sequence 特性畫圖, 但只計算 information site (至少其中二組序列含相同突變才考量)(優點: 可得多種演化樹供選擇, 缺點: 不遵循minimal evolution不適用, HIV演化)
A.       演算法則:
                       i.              Exhaustive search : 一定可得最短演化路徑
                      ii.              Branch and bound : 以巧妙運算所得最短演化樹, 但可能非最短演化路徑
                    iii.              ….
Parsimony 方法可能會有 Hill-climbing algorithms, 所以輸入的序列必須選擇Random order, NJ法沒此問題。
會得到多種演化樹選擇, 可重新畫出 consensus tree
B.       Phylip 操作
                       i.              Run [dnapars.exe]  
l  S à choose “rearrange on one best tree”
l  J à Randomize ….--> 給一個奇數 à Mix 次數 1
l  Y à get “outtree” à rename it as “intree”
                      ii.              Run [consensus.exe] to generate “outtree”  
                    iii.              View as [treeview]
2.     Maximum likelihood method à 不假設演化為簡約方式,以所有序列計算最大likelihood, 故短序列 (<500bp)採用最佳。演算基於統計方法, 可得P value, 不須再進行 Bootstrap. (缺點: 演算慢)
A.       Phylip 操作
                       i.              Run [dnaml.exe] to generate “outfile” and “outtree”
                      ii.              View “outtree” with [treeview], “outfile” with [notepad]


Statistical test for phylogenetic tree
l  Bootstrap analysis à 演化樹經重新取樣後, topology 還是可重覆出現的比率,至少75%以上才為可信的單系群, 70%可勉強接受, 若值低於70%以下, 若不同演算法可得相同的結果, 也可佐證。
1.       Run “seqboot.exe” to generate 1000 rearrangement sequences file à rename “outfile” as “infile”
2.       In N-J method
u  Run “dnadist.exe” to get distance matrix file à rename outfileasinfile
   M: yes, 1000 data sets
u  Run [neighbor.exe] to generate 1000 “outtree” à rename “outtree” as “intree”
   M: yes, 1000 data sets
u  Run [consensus.exe] to generate consensus “outtree”
u  View “outfile” with [notepad]  
   發表時須使用 “original tree” 標上bootstrap 統計值, 不可直接使用 bootstrap 所得之演化樹
3.       In Parsimony methods
u  Run [dnapars.exe] 
   S: choose rearrange on one best tree
   J: Randomize .--> 給一個奇數 à Mix 次數 1
   M: D, 1000 data sets
   Y: getouttreerename it asintree
u  Run [consensus.exe] to generate “outtree” 
u  View “outfile” with [notepad] 
   發表時須使用 original tree 標上bootstrap 統計值, 不可直接使用 bootstrap 所得之演化樹
4.       In Likilihood methods
u  Run [dnaml.exe] 
u  Y: get “outtree” à rename it as“intree”
u  View “outtree” with [treeview], “outfile” with [notepad]
   發表時須使用 original tree 標上bootstrap 統計值, 不可直接使用 bootstrap 所得之演化樹


Guideline