| 5 | 1/1 | 返回列表 |
| 查看: 3579 | 回復(fù): 11 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||||
cnlics木蟲(chóng) (小有名氣)
|
[交流]
【原創(chuàng)】從頭建模和從頭計(jì)算 已有6人參與
|
||||
|
“從頭”這個(gè)名詞,學(xué)化學(xué)的人更熟悉,量子化學(xué)、量子力學(xué)中,利用薛定諤方程解析波函數(shù),進(jìn)而計(jì)算原子、分子的能量狀態(tài)分布,是十分常見(jiàn)行為,這對(duì)于研究蛋白質(zhì)結(jié)構(gòu)的人來(lái)說(shuō),實(shí)在是太有誘惑力了,因?yàn)楹茱@然地,如果能計(jì)算出蛋白質(zhì)結(jié)構(gòu),從實(shí)用性這一點(diǎn)來(lái)說(shuō),誰(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ì)感覺(jué)頭疼的,米氏方程已經(jīng)是生化中相當(dāng)復(fù)雜的公式了,多變量方程中,系數(shù)的擬合、參數(shù)求解會(huì)難到一大批學(xué)生物的。我把以前整理的東西發(fā)出來(lái),可能會(huì)加上一些說(shuō)明,希望純粹學(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 ] |
分子生物實(shí)驗(yàn)及蛋白純化結(jié)晶相關(guān)鏈接 | 蛋白質(zhì)生物學(xué)實(shí)驗(yàn)經(jīng)驗(yàn) |
木蟲(chó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 ) |
木蟲(chóng) (小有名氣)
木蟲(chóng) (小有名氣)
|
第二部分 蛋白質(zhì)結(jié)構(gòu)建;A(chǔ) 2.1建模一般過(guò)程見(jiàn)下圖 ![]() 2.2結(jié)構(gòu)建模方法: •從頭建模 •比較或同源建模 •折疊識(shí)別或者穿針引線(threading)方法 2.3力場(chǎng)決定分子折疊機(jī)理 力場(chǎng)是一種描述分子構(gòu)象或者說(shuō)原子坐標(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非鍵相互作用(即,不是通過(guò)化學(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)模型 •通過(guò)局部結(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ì)超過(guò)格點(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通過(guò)局部構(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 ] |
木蟲(chó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 AssumptionDistribution 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. MethodFragment 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 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +3 | Ymlll 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開(kāi)始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-12 | 5/250 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[論文投稿]
投稿問(wèn)題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|