 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (5838)
- >導(dǎo)師招生 (1158)
- >考研 (1101)
- >文獻(xiàn)求助 (424)
- >蟲(chóng)友互識(shí) (412)
- >休閑灌水 (185)
- >考博 (136)
- >碩博家園 (128)
- >招聘信息布告欄 (125)
- >論文投稿 (106)
- >基金申請(qǐng) (95)
- >博后之家 (87)
- >綠色求助(高懸賞) (43)
- >SciFinder/Reaxys (36)
- >論文道賀祈福 (32)
- >找工作 (32)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1535 | 回復(fù): 11 | ||
[求助]
modeller multiple templates求助
|
|
from modeller import * log.verbose() env = environ() env.io.atom_files_directory = './:../atom_files/' aln = alignment(env) for (code, chain) in (('2mdh', 'A'), ('1bdm', 'A'), ('1b8p', 'A')): mdl = model(env, file=code, model_segment=('FIRST:'+chain, 'LAST:'+chain)) aln.append_model(mdl, atom_files=code, align_codes=code+chain) for (weights, write_fit, whole) in (((1., 0., 0., 0., 1., 0.), False, True), ((1., 0.5, 1., 1., 1., 0.), False, True), ((1., 1., 1., 1., 1., 0.), True, False)): aln.salign(rms_cutoff=3.5, normalize_pp_scores=False, rr_file='$(LIB)/as1.sim.mat', overhang=30, gap_penalties_1d=(-450, -50), gap_penalties_3d=(0, 3), gap_gap_score=0, gap_residue_score=0, dendrogram_file='fm00495.tree', alignment_type='tree', # If 'progresive', the tree is not # computed and all structues will be # aligned sequentially to the first feature_weights=weights, # For a multiple sequence alignment only # the first feature needs to be non-zero improve_alignment=True, fit=True, write_fit=write_fit, write_whole_pdb=whole, output='ALIGNMENT QUALITY') aln.write(file='fm00495.pap', alignment_format='PAP') aln.write(file='fm00495.ali', alignment_format='PIR') aln.salign(rms_cutoff=1.0, normalize_pp_scores=False, rr_file='$(LIB)/as1.sim.mat', overhang=30, gap_penalties_1d=(-450, -50), gap_penalties_3d=(0, 3), gap_gap_score=0, gap_residue_score=0, dendrogram_file='1is3A.tree', alignment_type='progressive', feature_weights=[0]*6, improve_alignment=False, fit=False, write_fit=True, write_whole_pdb=False, output='QUALITY') 這一段中，粗體字部分看不懂，期待得到解答，先謝謝了 |
回復(fù)此樓» 猜你喜歡
277調(diào)劑
已經(jīng)有6人回復(fù)
-
341求調(diào)劑
已經(jīng)有7人回復(fù)
-
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
293求調(diào)劑
已經(jīng)有12人回復(fù)
-
070300化學(xué)319求調(diào)劑
已經(jīng)有4人回復(fù)
-
268求調(diào)劑
已經(jīng)有9人回復(fù)
-
301求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)336分求調(diào)劑
已經(jīng)有5人回復(fù)
-
301求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料專(zhuān)碩306英一數(shù)二
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助該怎樣回復(fù)國(guó)外教授
已經(jīng)有7人回復(fù)
- 【求助】2012考研南京大學(xué)356分求調(diào)劑
已經(jīng)有48人回復(fù)
- 新人求助，環(huán)境專(zhuān)業(yè)去哪里比較好
已經(jīng)有14人回復(fù)
- XRD表征求助
已經(jīng)有13人回復(fù)
- 新手求助關(guān)于中科院的問(wèn)題一二··
已經(jīng)有7人回復(fù)
- 求助苯環(huán)上的醛基保護(hù)方法，反應(yīng)的催化劑最好容易得到。
已經(jīng)有4人回復(fù)
- 求助：倫敦國(guó)王學(xué)院、荷蘭萊頓大學(xué)，難以選擇哪個(gè)申請(qǐng)CSC
已經(jīng)有9人回復(fù)
- MTT法求助--藥物溶解問(wèn)題
已經(jīng)有10人回復(fù)
- 求2011年國(guó)基標(biāo)書(shū)一份，愿傾所有金幣【無(wú)人應(yīng)助，求助結(jié)束】
已經(jīng)有5人回復(fù)
- 求助翻譯一段話(huà)
已經(jīng)有1人回復(fù)
- 新人求助 關(guān)于公派讀博
已經(jīng)有7人回復(fù)
- 求助，誰(shuí)錯(cuò)了？
已經(jīng)有5人回復(fù)
- 【求助】硫酸蒽酮法測(cè)多糖出現(xiàn)渾濁
已經(jīng)有13人回復(fù)
- 【求助】求助Easy Modeller 的使用
已經(jīng)有11人回復(fù)
- 【求助】求endnote 的一個(gè)manuscript templates
已經(jīng)有6人回復(fù)
- 【求助】中國(guó)藥品生物制品檢定所的考研真題
已經(jīng)有5人回復(fù)
2樓2012-02-29 09:20:25
|
對(duì)于gap_penalties_1d明白了，在modeller網(wǎng)站上查到這樣一句話(huà) The gap initiation and extension penalties are specified by gap_penalties_1d. The default values of -900 -50 for the 'as1.sim.mat' similarity matrix were found to be optimal for pairwise alignments of sequences that share from 30% to 45% sequence identity (RS and AŠ, in preparation). |
3樓2012-02-29 10:39:13
4樓2012-02-29 15:02:49
5樓2012-03-03 16:36:22
6樓2012-03-03 16:40:30
7樓2012-03-03 16:42:06
8樓2012-03-03 16:53:56
|
gap_penalties_3d ，算是懂一點(diǎn)了： This optimal alignment is obtained by an iterative procedure as follows (the same as in alignment.align3d()). The average structures for both sub-alignments are calculated for all sub-alignment positions with at least one defined selected atom. This calculation is straightforward because the structures within the two sub-alignments are already superposed with each other (see below). Then, the distance matrix for dynamic programming with affine gap penalties is calculated as the matrix of Euclidean distances between the two averages. The dynamic programming results into a new alignment, dependent also on the gap initiation and extension penalties gap_penalties_3d (a reasonable setting is (0, 3)). gap_penalties_3d[0] is a gap creation penalty (usually 0), and gap_penalties_3d[1] is a gap extension penalty, say 3. When the gap initiation penalty is 0, pairs of positions are identified as equivalent when they have their selected atoms at most 2 times gap_penalties_3d[1] angstroms apart in the current superposition, as described for the alignment.align3d() command. The new alignment is then used to generate the new superposition of the two averages, and the iteration of the distance matrix calculation, alignment and superposition is repeated until there are no changes in the number of equivalent positions and in the rotation matrix relating the two averages. |
9樓2012-03-03 20:33:20
10樓2012-03-03 20:36:05
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 nufang19a 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +7 | zjl的號(hào) 2026-03-16 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | 大口吃飯 身體健 2026-03-13 | 4/200 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽(tīng)風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專(zhuān)業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | 柯淮然 2026-03-10 | 8/400 |
|