 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3577 | 回復(fù): 11 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
cnlics木蟲 (小有名氣)
|
[交流]
【原創(chuàng)】從頭建模和從頭計(jì)算 已有6人參與
|
||||
|
“從頭”這個(gè)名詞，學(xué)化學(xué)的人更熟悉，量子化學(xué)、量子力學(xué)中，利用薛定諤方程解析波函數(shù)，進(jìn)而計(jì)算原子、分子的能量狀態(tài)分布，是十分常見行為，這對(duì)于研究蛋白質(zhì)結(jié)構(gòu)的人來說，實(shí)在是太有誘惑力了，因?yàn)楹茱@然地，如果能計(jì)算出蛋白質(zhì)結(jié)構(gòu)，從實(shí)用性這一點(diǎn)來說，誰(shuí)計(jì)算得多，意味著誰(shuí)就能發(fā)現(xiàn)、擁有更多的藥靶。但受限于目前的計(jì)算能力，對(duì)蛋白質(zhì)或者較大的肽段進(jìn)行波函數(shù)計(jì)算需要的時(shí)間非常長(zhǎng)，在實(shí)際工作中要做大量簡(jiǎn)化，引入分子動(dòng)力學(xué)方法，使得計(jì)算蛋白質(zhì)結(jié)構(gòu)成為可能，隨著計(jì)算機(jī)計(jì)算的進(jìn)步，蛋白質(zhì)結(jié)構(gòu)計(jì)算方法也會(huì)一直進(jìn)化著。 “從頭”這個(gè)詞對(duì)應(yīng)的拉丁語(yǔ)翻譯為ab initio，當(dāng)然，也可以是更為粗略的理解，比如理解為de novo。對(duì)于精細(xì)的化學(xué)計(jì)算方法，學(xué)生物的一般會(huì)感覺頭疼的，米氏方程已經(jīng)是生化中相當(dāng)復(fù)雜的公式了，多變量方程中，系數(shù)的擬合、參數(shù)求解會(huì)難到一大批學(xué)生物的。我把以前整理的東西發(fā)出來，可能會(huì)加上一些說明，希望純粹學(xué)生物的人能理解，更希望能利用現(xiàn)成的一些工具實(shí)現(xiàn)自己的目的，也因?yàn)榇耍泳筒环旁凇胺肿幽M”目錄下了，希望能起到拋磚引玉的作用。 內(nèi)容提要： 第一部分 蛋白質(zhì)結(jié)構(gòu)與分子建模 第二部分 蛋白質(zhì)結(jié)構(gòu)建�；A(chǔ) 第三部分 從頭預(yù)測(cè)方法 第四部分 Rosetta方法 第五部分 從頭計(jì)算 [ Last edited by cnlics on 2010-9-16 at 15:23 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 分子生物實(shí)驗(yàn)及蛋白純化結(jié)晶相關(guān)鏈接 | 蛋白質(zhì)生物學(xué)實(shí)驗(yàn)經(jīng)驗(yàn) |
» 猜你喜歡
焦慮
已經(jīng)有10人回復(fù)
-
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
-
297求調(diào)劑
已經(jīng)有3人回復(fù)
-
302求調(diào)劑
已經(jīng)有8人回復(fù)
-
環(huán)境工程調(diào)劑
已經(jīng)有7人回復(fù)
-
268求調(diào)劑
已經(jīng)有6人回復(fù)
-
311求調(diào)劑
已經(jīng)有11人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有4人回復(fù)
-
301求調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 第一性原理、分子模擬和量子化學(xué)計(jì)算之間的區(qū)別是什么
已經(jīng)有23人回復(fù)
- 【討論】《量子化學(xué)——基本原理和從頭計(jì)算法》關(guān)于GTO積分的一處小錯(cuò)誤
已經(jīng)有3人回復(fù)
- 【討論】DFT方法屬不屬于從頭算？
已經(jīng)有13人回復(fù)
- 【分享】量子化學(xué)——基本原理和從頭計(jì)算法（上、中、下冊(cè)）.pdf
已經(jīng)有246人回復(fù)
cnlics
木蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 3014.2
- 紅花: 4
- 帖子: 270
- 在線: 422.4小時(shí)
- 蟲號(hào): 795158
- 注冊(cè): 2009-06-16
- 性別: GG
- 專業(yè): 當(dāng)代宗教
|
第二部分 蛋白質(zhì)結(jié)構(gòu)建�；A(chǔ) 2.1建模一般過程見下圖  2.2結(jié)構(gòu)建模方法： •從頭建模 •比較或同源建模 •折疊識(shí)別或者穿針引線（threading）方法 2.3力場(chǎng)決定分子折疊機(jī)理 力場(chǎng)是一種描述分子構(gòu)象或者說原子坐標(biāo)決定分子能量的數(shù)學(xué)函數(shù)。力場(chǎng)是分子內(nèi)部化學(xué)鍵的拉伸，彎曲，扭曲以及范德華力、靜電相互作用等作用的總和。 常用力場(chǎng)有： CHARMM (Harvard) http://yuri.harvard.edu/ GROMOS96 (Groningen/ETHZ) http://www.igc.ethz.ch/gromos AMBER (Scripps) http://amber.scripps.edu SYBYL Tripos Inc. DISCOVER MSI Inc. …….. 2.4非鍵相互作用（即，不是通過化學(xué)鍵而產(chǎn)生的作用） •靜電相互作用。計(jì)算方法是：  •氫鍵作用。計(jì)算方法：  •疏水作用。計(jì)算方法為HINT •范德華力：利用Lennard-Jones函數(shù)計(jì)算  2.5簡(jiǎn)化的結(jié)構(gòu)模型方法 •格點(diǎn)模型 •離散狀態(tài)非格點(diǎn)（off-lattice）模型 •通過局部結(jié)構(gòu)預(yù)測(cè)，減少構(gòu)象搜索空間 Bonneau R. et al.Annu. Rev. Biophys. Biomol. Struct. 30:173–89 (2001 ) 2.5.1格點(diǎn)模型：簡(jiǎn)化方法為，用格點(diǎn)表示肽鏈 優(yōu)點(diǎn):分析、計(jì)算簡(jiǎn)單 缺點(diǎn): 難以表示精細(xì)的幾何取向（例如鏈扭曲）；骨架結(jié)構(gòu)精確性不會(huì)超過格點(diǎn)間距的一半。  2.5.2離散狀態(tài)格點(diǎn)模型 降低復(fù)雜度的方法：僅允許特定的支鏈構(gòu)象且肽鏈的化學(xué)鍵轉(zhuǎn)動(dòng)受限（即將支鏈限定為單一的構(gòu)象；肽鏈骨架限定為特殊的φ/ψ角度組合等） 缺點(diǎn): ω角總是平面角（0或者180°） 2.5.3通過局部構(gòu)象預(yù)測(cè)，限定總體構(gòu)象搜索空間 降低復(fù)雜度的方法：利用局部構(gòu)象的偏向性（biases）或者序列基元（motif）降低復(fù)雜度 缺點(diǎn)：局部構(gòu)象的偏向性（biases）程度、強(qiáng)度非常依賴序列；不同的序列中，序列基元非常傾向于采取單一的局部構(gòu)象 2.5打分函數(shù) 基于溶劑化的打分方法：總結(jié)已知蛋白質(zhì)中各位點(diǎn)的溶劑化程度；弄清各氨基酸在各位點(diǎn)出現(xiàn)的頻率。 配對(duì)相互作用：某倆殘基有多大可能互相靠近 二級(jí)結(jié)構(gòu)安排：打分評(píng)估二級(jí)結(jié)構(gòu)的元件之間匹配的程度。 [ Last edited by cnlics on 2010-9-18 at 19:42 ] |
3樓2010-09-14 00:55:43
cnlics
木蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 3014.2
- 紅花: 4
- 帖子: 270
- 在線: 422.4小時(shí)
- 蟲號(hào): 795158
- 注冊(cè): 2009-06-16
- 性別: GG
- 專業(yè): 當(dāng)代宗教
|
第一部分 蛋白質(zhì)結(jié)構(gòu)與分子建模 分子建模對(duì)象可以分為大、小分子兩類： •大分子:蛋白質(zhì)，酶，聚合物，DNA等 •小分子：水分子，普通有機(jī)、無機(jī)化合物，藥物，除草劑，殺蟲劑等 分子建模內(nèi)容包括： •顯示分子 •研究構(gòu)象、結(jié)構(gòu) •研究分子性質(zhì) •預(yù)測(cè)新分子的性質(zhì) •分子比較 •研究分子相互作用 •研究反應(yīng)機(jī)理 •預(yù)測(cè)反應(yīng)機(jī)理 1.1 蛋白質(zhì)的結(jié)構(gòu)層次和結(jié)構(gòu)元件 •活性位點(diǎn)，模塊，結(jié)構(gòu)核心，折疊 •結(jié)構(gòu)域，結(jié)構(gòu)基元 •家族，超家族 •結(jié)構(gòu)元件 •結(jié)構(gòu)的類型 •1-4級(jí)結(jié)構(gòu) 1.2 蛋白質(zhì)結(jié)構(gòu)分類規(guī)律：拉氏圖（Ramchandran plot）  1.3結(jié)構(gòu)相關(guān)的假設(shè)與實(shí)踐 結(jié)構(gòu)假設(shè): •蛋白質(zhì)的肽鏈會(huì)折疊形成唯一的3D結(jié)構(gòu) •處于自由能最低狀態(tài) •不同的序列，可以有相同的折疊方式 •結(jié)構(gòu)形成存在驅(qū)動(dòng)力 實(shí)踐經(jīng)驗(yàn): •蛋白質(zhì)天然結(jié)構(gòu)與未折疊狀態(tài)之間的能量差別非常小 •蛋白質(zhì)的天然結(jié)構(gòu)往往不是能量最低狀態(tài)，但與其他可能的結(jié)構(gòu)相比，最接近最低自由能狀態(tài) [ Last edited by cnlics on 2010-9-18 at 18:40 ] |
2樓2010-09-14 00:41:44
cnlics
木蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 3014.2
- 紅花: 4
- 帖子: 270
- 在線: 422.4小時(shí)
- 蟲號(hào): 795158
- 注冊(cè): 2009-06-16
- 性別: GG
- 專業(yè): 當(dāng)代宗教
|
Part III: ab initio prediction methods 3.1 Scoring functions Molecular Dynamics simulations(MD) Monte Carlo (MC) simulations Pathway Models Combined Hierarchical Approach Genetic Algorithms more……. 3.2 Ab initio prediction:Using pathway models Pathway models combine the scoring function and the search. HMMSTR-CM: a fragment library (knowledgebased potentials ) + a set of nucleation/propagation-based rules(for building a protein contact maps) 3.3 Ab initio prediction: TOUCHSTONE ----- threading based tertiary restraints SICHO (SIde CHain Only) model Prediction of tertiary restraints:side chain contact(PROSPECTOR); consensus contacts Structure selection with an atomic potential:Monte Carlo simulations; Kihara D. et al .PNAS , 98 (18) :10125–10130(2001) 3.4 Ab initio prediction: Combined Hierarchical Approach highly simplified tetrahedral lattice model:all-atom models combined allatom knowledge-based scoring function:three smaller subsets consensus-based distance geometry procedure Samudrala R. et al.PROTEINS: Structure, Function, and Genetics Suppl 3:194–198 (1999) 3.5 Ab initio prediction: more…. Distance geometry-based Ramachandran Plots-based Rosetta Huang ES et al. J. Mol. Biol. (1999) 290, 267-281. Bernasconi A. et al.ERCIM News No.43 (2000 ) |
4樓2010-09-14 00:59:20
cnlics
木蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 3014.2
- 紅花: 4
- 帖子: 270
- 在線: 422.4小時(shí)
- 蟲號(hào): 795158
- 注冊(cè): 2009-06-16
- 性別: GG
- 專業(yè): 當(dāng)代宗教
|
Part IV: Rosetta 4.1 What’s CASP：Critical Assessment of Structure Prediction (homology modeling, threading, and ab-initio) Every two years in Dec, community-wide blind test of prediction methods Experimentalists announce some protein sequences that they are going to resolve structurally CASP put these sequence on web for pediction with deadline Computational biologists submit their predictions CASP evaluates the predictions according to the results resolved by experimentalists http://PredictionCenter.llnl.gov/caspi/ ( i = 1,2,3,4,5. ) http://www.forcasp.org/ 4.2 Rosetta at CASP:an example taken from CASP5 4.3 Rosetta : the method Model: Narrow the search with local structure Prediction Scoring function(Solvation-based & Pair interactions) Method outline (two steps) Get tiny pieces:sequence profile alignment Put them together:Monte-Carlo method; Bayesian scoring function Chivian D. et al.PROTEINS: Structure, Function, and Genetics 53:524–533 (2003) 4.3.1 Get tiny pieces:Construction of I-sites library Assumption Distribution of conformations sampled for a given nine residue segment of the chain is reasonably well approximated by the distribution of structures adopted the sequence(and closely related sequences) in known protein structures. Method Fragment libraries for each three and nine residue segment of the chain are extracted from PDB using sequence profile alignment 4.3.2 Get tiny pieces: construction procedure Construct profiles (PSI-BLAST like) for each solved structure Collect each possible segments of fixed length (len = 3, 9, 15) Perform k-means clustering of segments Check each cluster for a “coherent” structure (in terms of dihedral angles Prune incoherent structures Iteratively refine remaining clusters by removing structurally different segments, redefining cluster membership, etc. 4.4.1 Put them together: Procedure For representative proteins, backbones were assembled from a library of 1000 different 5-residue fragments. 4.4.2 Put them together: Monte Carlo Search the resulting conformational space with Monte-Carlo method Bayesian scoring function:Chose the most likely structure given the sequence: 4.4.3 Put them together: Scoring Function 4.5 Using Rosetta: Comparative modeling Detection of the best parent for each putative domain: Blast or PSI-Blast parents or Pcons parents Sequence alignment to that parent: K*SYNC (kitchen sink) Modeling of structurally variable regions:match with DSSP assigned secondary structure Optimization to increase the physical reasonableness of the final model:fragment replacement and random angle perturbations Reassemble the complete chain when domains were parsed and processed individually:evaluated by a coarse energy function 4.6 Using Rosetta: De Novo structure predictions Fragment libraries for each three and nine residue segment Monte Carlo procedure with energy function favoring compact structures, buried hydrophobic residues, and paired beta strands low free energy models : MC Minimization procedure to relieve backbone atomic clashes MC minimize an all-atom energy function Bonneau R. et al.J. Mol. Biol. (2002) 322, 65–78 4.7 Using Rosetta: Automated Method for Full Chain Structure Prediction Robetta: de novo, comparative, or mixed models Secondary structure prediction from the JUFO-3D method |
5樓2010-09-14 01:20:03
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +5 | zbzihdhd 2026-03-15 | 6/300 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|