前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >rosetta-fixbb

rosetta-fixbb

作者头像
DrugSci
发布2021-02-04 15:04:03
6850
发布2021-02-04 15:04:03
举报
文章被收录于专栏:FindKeyFindKey

purpose:

use fixbb application to invoke the Rosetta packer to optimize the side-chain conformations of an input structure.

basic command:

fixbb.static.linuxgccrelease -in:file:s 1l2y.pdb -in:file:fullatom -resfile resfile.txt -nstruct 5 >log.txt eplain:

-in:file:s 1l2y.pdb : the input PDB file

-in:file:fullatom: the input file is a full-atom respresentation (as opposed to a centroid-only model)--(全原子模型区别于质心模型)

-resfile resfile.txt : controlling the packer(rosetta中的文件,用于控制packer).

-nstruct 5: tells Rosetta to repeat the protocol five times, producing five independent outputs. (重复5次,得到5个结构,但由于rosetta使用随机算法,所以实际数目需要提升)

resfile:

代码语言:javascript
复制
NATAA
start

```

input:

代码语言:javascript
复制
1l2y.pdb

output:

代码语言:javascript
复制
1l2y_0001.pdb
1l2y_0002.pdb
1l2y_0003.pdb
1l2y_0004.pdb
1l2y_0005.pdb
log.txt
score.sc

step:

(1)pymol->open 1l2y.pdb(observe)

(2)output:

score.sc

代码语言:javascript
复制
SEQUENCE:
SCORE: total_score dslf_fa13    fa_atr    fa_dun   fa_elec fa_intra_rep fa_intra_sol_xover4              fa_rep              fa_sol hbond_bb_sc hbond_lr_bb    hbond_sc hbond_sr_bb lk_ball_wtd       omega     p_aa_pp pro_close rama_prepro         ref yhh_planarity description
SCORE:     -13.104     0.000   -75.943    18.637   -17.270        0.173               3.270              12.601              57.560      -1.336      -1.005       0.000      -6.434      -2.802       2.415      -6.345     3.741       1.959      -2.324         0.000 1l2y_0001
SCORE:     -13.104     0.000   -75.943    18.637   -17.270        0.173               3.270              12.601              57.560      -1.336      -1.005       0.000      -6.434      -2.802       2.415      -6.345     3.741       1.959      -2.324         0.000 1l2y_0002
SCORE:     -13.104     0.000   -75.943    18.637   -17.270        0.173               3.270              12.601              57.560      -1.336      -1.005       0.000      -6.434      -2.802       2.415      -6.345     3.741       1.959      -2.324         0.000 1l2y_0003
SCORE:     -13.104     0.000   -75.943    18.637   -17.270        0.173               3.270              12.601              57.560      -1.336      -1.005       0.000      -6.434      -2.802       2.415      -6.345     3.741       1.959      -2.324         0.000 1l2y_0004
SCORE:     -13.104     0.000   -75.943    18.637   -17.270        0.173               3.270              12.601              57.560      -1.336      -1.005       0.000      -6.434      -2.802       2.415      -6.345     3.741       1.959      -2.324         0.000 1l2y_0005

pdb:

better command:

代码语言:javascript
复制
fixbb.static.linuxgccrelease -in:file:s 1l2y.pdb -in:file:fullatom -resfile resfile.txt -nstruct 1 -ex1 -ex2 -ex3 -ex4 -extrachi_cutoff 0 -out:prefix run2_ -out:file:scorefile run2.sc >log2.txt

``` "-ex1", "-ex2", "-ex3", and "-ex4" :activate extra rotamers for the first, second, third, and fourth side-chain chi angles(控制侧链旋转角度).

-extrachi_cutoff: flag determines the number of residue neighbours (nearby residues) that a resiude must have in order to receive the additional rotamers(如果你想使所有的残基侧链进行旋转,将它设置为0).

-out:prefix run2_ :write PDB files with the prefix "run2_" in the filename(以run2_作为前缀)

-out:file:scorefile run2.sc:directs the final scores to a new file(输出得分文件)

run2.sc

代码语言:javascript
复制
SEQUENCE:
SCORE: total_score dslf_fa13    fa_atr    fa_dun   fa_elec fa_intra_rep fa_intra_sol_xover4              fa_rep              fa_sol hbond_bb_sc hbond_lr_bb    hbond_sc hbond_sr_bb lk_ball_wtd       omega     p_aa_pp pro_close rama_prepro         ref yhh_planarity description
SCORE:     -21.078     0.000   -80.039    20.857   -24.261        0.176               3.379              11.104              60.743      -0.800      -1.005      -1.813      -6.434      -2.428       2.415      -6.345     3.741       1.959      -2.324         0.000 run2_1l2y_0001
SCORE:     -21.078     0.000   -80.039    20.857   -24.261        0.176               3.379              11.104              60.743      -0.800      -1.005      -1.813      -6.434      -2.428       2.415      -6.345     3.741       1.959      -2.324         0.000 run2_1l2y_0002
SCORE:     -20.799     0.000   -77.560    19.762   -19.161        0.165               3.106               9.948              55.829      -0.800      -1.005      -0.917      -6.434      -3.178       2.415      -6.345     3.741       1.959      -2.324         0.000 run2_1l2y_0003
SCORE:     -21.078     0.000   -80.039    20.857   -24.261        0.176               3.379              11.104              60.743      -0.800      -1.005      -1.813      -6.434      -2.428       2.415      -6.345     3.741       1.959      -2.324         0.000 run2_1l2y_0004
SCORE:     -20.798     0.000   -77.364    19.993   -18.947        0.163               3.085               9.608              55.796      -0.800      -1.005      -1.031      -6.434      -3.307       2.415      -6.345     3.741       1.959      -2.324         0.000 run2_1l2y_0005

reference:

http://www.ccp14.ac.uk/ccp/web-mirrors/garlic/garlic/commands/dihedrals.html

本文参与 腾讯云自媒体分享计划,分享自微信公众号。
原始发表:2020-07-17,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 FindKey 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档