 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1539 | 回復(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ù)此樓» 猜你喜歡
354求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料專業(yè)求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有3人回復(fù)
-
085600材料與化工
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有6人回復(fù)
-
286求調(diào)劑
已經(jīng)有9人回復(fù)
-
焦慮
已經(jīng)有12人回復(fù)
-
344求調(diào)劑
已經(jīng)有6人回復(fù)
-
266求調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助該怎樣回復(fù)國外教授
已經(jīng)有7人回復(fù)
- 【求助】2012考研南京大學(xué)356分求調(diào)劑
已經(jīng)有48人回復(fù)
- 新人求助，環(huán)境專業(yè)去哪里比較好
已經(jīng)有14人回復(fù)
- XRD表征求助
已經(jīng)有13人回復(fù)
- 新手求助關(guān)于中科院的問題一二··
已經(jīng)有7人回復(fù)
- 求助苯環(huán)上的醛基保護方法，反應(yīng)的催化劑最好容易得到。
已經(jīng)有4人回復(fù)
- 求助：倫敦國王學(xué)院、荷蘭萊頓大學(xué)，難以選擇哪個申請CSC
已經(jīng)有9人回復(fù)
- MTT法求助--藥物溶解問題
已經(jīng)有10人回復(fù)
- 求2011年國基標(biāo)書一份，愿傾所有金幣【無人應(yīng)助，求助結(jié)束】
已經(jīng)有5人回復(fù)
- 求助翻譯一段話
已經(jīng)有1人回復(fù)
- 新人求助 關(guān)于公派讀博
已經(jīng)有7人回復(fù)
- 求助，誰錯了？
已經(jīng)有5人回復(fù)
- 【求助】硫酸蒽酮法測多糖出現(xiàn)渾濁
已經(jīng)有13人回復(fù)
- 【求助】求助Easy Modeller 的使用
已經(jīng)有11人回復(fù)
- 【求助】求endnote 的一個manuscript templates
已經(jīng)有6人回復(fù)
- 【求助】中國藥品生物制品檢定所的考研真題
已經(jīng)有5人回復(fù)
2樓2012-02-29 09:20:25
|
對于gap_penalties_1d明白了，在modeller網(wǎng)站上查到這樣一句話 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 ，算是懂一點了： 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ā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +6 | llllkkkhh 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 9/450 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 293求調(diào)劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +4 | 芬達46 2026-03-12 | 4/200 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|