版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

sijifengsd

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 273.7
散金: 50
红花: 2
帖子: 94
在线: 12.1小时
虫号: 794988
注册: 2009-06-15
性别: GG
专业: 物理化学

[交流] 【经验】用Modeller补全残基学习心得已有8人参与

最近专注于分子动力学的模拟研究。发现好多的蛋白质晶体结构或多或少的有一部分残基丢失了，要做动力学模拟，就需要把丢失的残基给补上。看文献发现比较公认的、又是免费的软件当属Modeller了。所以本人花了两个周的时间研究了这个软件的使用方法。现把个人的心得体会奉献出来，希望给需要的同学们提供点帮助：
首先大家可以到这个网站申请个序列号，下载个免费的软件:
http://salilab.org/modeller/registration.html
安装成功就可用了。
这个软件要用到一些脚本文件，补全残基可以带着配体，水，也可以不带。
一没有配体、水分子：
第一步处理大分子pdb 文件，删除杂原子和水分子等。
第二步建立一个get_ali.py 脚本文件

# get sequece of the 2FMX PDB file, and write to an alignment file.

log.verbose()
env = environ()

env.io.atom_files_directories = ['.', '../atom_file']

code = '2FMX'
m = model(env, file=code)
aln = alignment(env)

aln.append_model(m, align_codes=code)
aln.write(code+'.seq')

第三步：运行 modpv7 get_ali.py 产生一个sequence 文件：2FMX.seq

>P1;2FMX
structureX:2FMX:  13 :A:+556 :B:MOL_ID  1; MOLECULE  GTP-BINDING PROTEIN SAR1B; CHAIN  A, B; FRAGMENT  RESIDUES 10-198; SYNONYM  SAR1, GTBPB; ENGINEERED  YES:MOL_ID  1; ORGANISM_SCIENTIFIC  CRICETULUS GRISEUS; ORGANISM_COMMON  CHINESE HAMSTER; ORGANISM_TAXID  10029; EXPRESSION_SYSTEM  ESCHERICHIA COLI; EXPRESSION_SYSTEM_TAXID  562; EXPRESSION_SYSTEM_VECTOR_TYPE  PLASMID; EXPRESSION_SYSTEM_PLASMID  PET11D: 1.82:-1.00
SSVLQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEELTIAGMTFTTFDLGRVWKNYLPAI
NGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEERLREMFGLYGQTTGKGSVSLKEL
NARPLEVFMCSVLKRQGYGEGFRWMAQYID
第四步：用产生的序列文件产生一个alignment.ali脚本文件：把上面的序列复制到脚本中

>P1;2FMX
structureX:2FMX.pdb: 13 :A:4872 :A:undefined:undefined:-1.00:-1.00
SSVL---QFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEELTIAGMTFTTFDLGRVWKNYLPAI
NGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEERLREMFGLYGQTTGKGSVSLKEL
NARPLEVFMCSVLKRQGYGEGFRWMAQYID---

>P1;2FMX_fill
sequence:2FMX_fill:    :A  :    :A ::: 0.00: 0.00
SSVLHJKQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEELTIAGMTFTTFDLGRVWKNYLPAI
NGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEERLREMFGLYGQTTGKGSVSLKEL
NARPLEVFMCSVLKRQGYGEGFRWMAQYIDLPO
注意：上面的一个序列是pdb源文件的残基，而下面是加上缺失残基的全序列。13 :A:4872 A:表示在源文件中开始和结束的残基序号。
第五步建立一个select.py文件，在这个文件中我们只优化补全的那部分残基。

from modeller import *
from modeller.automodel import * # Load the automodel class

log.verbose()
env = environ()

# directories for input atom files
env.io.atom_files_directories = ['.', '../atom_file']
class MyModel(loopmodel):
def select_atoms(self):
      return selection(self.residue_range('1:A', '9:A'),
                     self.residue_range('68:A', '73:A'),

a = MyModel(env,
         alnfile = 'alignment.ali',
         knowns = '2FMX',
         sequence = '2FMX_fill')

a.starting_model = 1
a.ending_model = 2

a.make()

注释: 上面的'1:A', '9:A'，'68:A', '73:A' 是加入的那部分残基在从1开始中的序列编号；下面的红字表示可以改动的地方。class MyModel(automodel) 中的automodel也可以换成loopmodel,这样对已loop区的优化更好。
第六步运行 mod9v7 select.py 结果提取pdb文件即可。
二包含配体和水分子
第一步大分子就不要处理了。
第二步建立一个get_ali.py 脚本文件

# get sequece of the 2FMX PDB file, and write to an alignment file.

log.verbose()
env = environ()

env.io.atom_files_directories = ['.', '../atom_file']

code = '2FMX'
m = model(env, file=code)
aln = alignment(env)

aln.append_model(m, align_codes=code)
aln.write(code+'.seq')

第三步：运行 modpv7 get_ali.py 产生一个sequence 文件：2FMX.seq

>P1;2FMX
structureX:2FMX:  13 :A:+556 :B:MOL_ID  1; MOLECULE  GTP-BINDING PROTEIN SAR1B; CHAIN  A, B; FRAGMENT  RESIDUES 10-198; SYNONYM  SAR1, GTBPB; ENGINEERED  YES:MOL_ID  1; ORGANISM_SCIENTIFIC  CRICETULUS GRISEUS; ORGANISM_COMMON  CHINESE HAMSTER; ORGANISM_TAXID  10029; EXPRESSION_SYSTEM  ESCHERICHIA COLI; EXPRESSION_SYSTEM_TAXID  562; EXPRESSION_SYSTEM_VECTOR_TYPE  PLASMID; EXPRESSION_SYSTEM_PLASMID  PET11D: 1.82:-1.00
SSVLQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEELTIAGMTFTTFDLGRVWKNYLPAI
NGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEERLREMFGLYGQTTGKGSVSLKEL
NARPLEVFMCSVLKRQGYGEGFRWMAQYID
第四步：用产生的序列文件产生一个alignment.ali脚本文件：把上面的序列复制到脚本中

>P1;2FMX
structureX:2FMX.pdb: 13 :A:4890 :A:undefined:undefined:-1.00:-1.00
SSVL---QFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEELTIAGMTFTTFDLGRVWKNYLPAI
NGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEERLREMFGLYGQTTGKGSVSLKEL
NARPLEVFMCSVLKRQGYGEGFRWMAQYID---.wwwwwwwwwwww

>P1;2FMX_fill
sequence:2FMX_fill:    :A  :    :A ::: 0.00: 0.00
SSVLHJKQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEELTIAGMTFTTFDLGRVWKNYLPAI
NGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEERLREMFGLYGQTTGKGSVSLKEL
NARPLEVFMCSVLKRQGYGEGFRWMAQYIDLPO. wwwwwwwwwwww
注意：到这一步和上面有区别了，因为有配体和水分子就要在序列的后面加上”.(表示配体，有几个配体杂原子加几个点)” 和”w(表示水分子，有几个加几个w)”。如果水分子和配体时单独成链的要加/例如：/. Wwwwwwwwwwww,或者 . /wwwwwwwwwwww
第五步建立一个select.py文件，在这个文件中我们只优化补全的那部分残基。

from modeller import *
from modeller.automodel import * # Load the automodel class

log.verbose()
env = environ()
env.io.water = True
env.io.hetatm = True

# directories for input atom files
env.io.atom_files_directories = ['.', '../atom_file']
class MyModel(loopmodel):
def select_atoms(self):
      return selection(self.residue_range('1:A', '9:A'),
                     self.residue_range('68:A', '73:A'),

a = MyModel(env,
         alnfile = 'alignment.ali',
         knowns = '2FMX',
         sequence = '2FMX_fill')

a.starting_model = 1
a.ending_model = 2

a.make()

注释:这里加入了env.io.water = True  env.io.hetatm = True还识别水分子和配体。
第六步运行 mod9v7 select.py 结果提取pdb文件即可
三复杂点的比如两条链的

第一步下载pdb
第二步建立一个get_ali.py 脚本文件

# get sequece of the 2FMX PDB file, and write to an alignment file.

log.verbose()
env = environ()
env.io.water = True
env.io.hetatm = True

env.io.atom_files_directories = ['.', '../atom_file']

code = '2FMX'
m = model(env, file=code)
aln = alignment(env)

aln.append_model(m, align_codes=code)
aln.write(code+'.seq')
注：这里直接用env.io.water = True  nv.io.hetatm = True 会直接把配体、水给你表示出来

第三步：运行 modpv7 get_ali.py 产生一个sequence 文件：2FMX.seq
>P1;2FMX
structureX:2FMX:  13 :A:+556 :B:MOL_ID  1; MOLECULE  GTP-BINDING PROTEIN SAR1B; CHAIN  A, B; FRAGMENT  RESIDUES 10-198; SYNONYM  SAR1, GTBPB; ENGINEERED  YES:MOL_ID  1; ORGANISM_SCIENTIFIC  CRICETULUS GRISEUS; ORGANISM_COMMON  CHINESE HAMSTER; ORGANISM_TAXID  10029; EXPRESSION_SYSTEM  ESCHERICHIA COLI; EXPRESSION_SYSTEM_TAXID  562; EXPRESSION_SYSTEM_VECTOR_TYPE  PLASMID; EXPRESSION_SYSTEM_PLASMID  PET11D: 1.82:-1.00
SSVLQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEELTIAGMTFTTFDLGRVWKNYLPAI
NGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEERLREMFGLYGQTTGKGSVSLKEL
NARPLEVFMCSVLKRQGYGEGFRWMAQYID/SSVLQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDPTLHPTSEE
LTIAGMTFTTFDLGVWKNYLPAINGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEE
RLREMFGLYGQTTGKGSVSLKELNARPLEVFMCSVLKRQGYGEGFRWMAQYID/.x/.$x/wwwwwwwwwwwwww
wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww
wwwwwwwww/wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww
wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww*

第四步：用产生的序列文件产生一个alignment.ali脚本文件：把上面的序列复制到脚本中

>P1;2FMX
structureX:2FMX.pdb: 13 :A:4872 :B:undefined:undefined:-1.00:-1.00
---------SSVLQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEELTIAGMTFTTFDLG------RVWKNYLPAI
NGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEERLREMFGLYGQTTGKGSVSLKEL
NARPLEVFMCSVLKRQGYGEGFRWMAQYID/---------SSVLQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKD-------PTLHPTSEE
LTIAGMTFTTFDLG-------VWKNYLPAINGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEE
RLREMFGLYGQTTGKGSVSLKELNARPLEVFMCSVLKRQGYGEGFRWMAQYID/.x/.$x/wwwwwwwwwwwwww
wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww
wwwwwwwww/wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww
wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww*

>P1;2FMX_fill
sequence:2FMX_fill:    :A  :    :B ::: 0.00: 0.00
HHHHHHSGFSSVLQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEELTIAGMTFTTFDLGGHIQARRVWKNYLPAI
NGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEERLREMFGLYGQTTGKGSVSLKEL
NARPLEVFMCSVLKRQGYGEGFRWMAQYID/HHHHHHSGFSSVLQFLGLYKKTGKLVFLGLDNAGKTTLLHMLKDDRLGQHVPTLHPTSEE
LTIAGMTFTTFDLGGHIQARRVWKNYLPAINGIVFLVDCADHERLLESKEELDSLMTDETIANVPILILGNKIDRPEAISEE
RLREMFGLYGQTTGKGSVSLKELNARPLEVFMCSVLKRQGYGEGFRWMAQYID/.x/.$x/wwwwwwwwwwwwww
wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww
wwwwwwwww/wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww
wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww*
第五步建立一个select.py文件，在这个文件中我们只优化补全的那部分残基。

from modeller import *
from modeller.automodel import * # Load the automodel class

log.verbose()
env = environ()

# directories for input atom files
env.io.atom_files_directories = ['.', '../atom_file']
env.io.water = True
env.io.hetatm = True

class MyModel(loopmodel):
def select_atoms(self):
      return selection(self.residue_range('1:A', '9:A'),
                     self.residue_range('68:A', '73:A'),
                     self.residue_range('196:B', '204:B'),
                     self.residue_range('240:B', '246:B'),
                     self.residue_range('270:B', '276:B'))


a = MyModel(env,
         alnfile = 'ali-ligand.ali',
         knowns = '2FMX',
         sequence = '2FMX_fill')

a.starting_model = 1
a.ending_model = 2

a.make()
注释:这里因为有两条链所以编号时B链的是在A链的基础上编的
第六步运行 mod9v7 select.py 结果提取pdb文件即可上面说的只是一些基本的东西，要想深入的学好Modeller 还要花点时间看看在线的tutorial
http://salilab.org/modeller/tutorial/

[ Last edited by lei0736 on 2009-11-24 at 22:03 ]

回复此楼

» 收录本帖的淘帖专辑推荐

» 本帖已获得的红花（最新10朵）

liuwei4410

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

似水年华

1楼2009-10-12 10:01:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiaoling0523

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2162.3
帖子: 83
在线: 37小时
虫号: 536980
注册: 2008-03-31

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

非常感谢您的资料，正在学习，太有帮助了

赞一下(1人)

回复此楼

2楼2009-10-12 15:02:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sijifengsd

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 273.7
散金: 50
红花: 2
帖子: 94
在线: 12.1小时
虫号: 794988
注册: 2009-06-15
性别: GG
专业: 物理化学

呵呵

希望如此，共同交流

赞一下

回复此楼

似水年华

3楼2009-10-12 16:03:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

mophyworld

金虫 (正式写手)

应助: 3 (幼儿园)
金币: 1282.1
散金: 10
红花: 4
帖子: 459
在线: 88.2小时
虫号: 631894
注册: 2008-10-20
性别: MM
专业: 理论和计算化学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

看上去很专业啊。先收藏了

赞一下(1人)

回复此楼

Ineedsomeinspirationandsomeconfidence

4楼2009-10-14 15:55:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

gogen

金虫 (小有名气)

应助: 4 (幼儿园)
金币: 1335.9
红花: 2
帖子: 160
在线: 66小时
虫号: 589241
注册: 2008-08-29
性别: GG
专业: 环境微生物学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

基础差啊，看到我头晕死了。。LZ这类软件有没有最简单的教学版本啊。。想先入门。。那个官方版的介绍基本看不懂

赞一下(1人)

回复此楼

5楼2010-03-07 22:33:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tyl681108

谢谢分享

6楼2010-03-08 17:22:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tjegg

铁杆木虫 (著名写手)

应助: 26 (小学生)
金币: 6869.3
红花: 16
帖子: 2615
在线: 448小时
虫号: 338983
注册: 2007-04-05
性别: GG
专业: 抗肿瘤药物药理

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

不错，谢谢搂住，非常感谢您的资料，太有帮助了

赞一下(1人)

回复此楼

除了你的亲人，没有人应该对你好，对你好的人，一定要珍惜。

7楼2010-03-08 19:09:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

patent

至尊木虫 (正式写手)

应助: 3 (幼儿园)
金币: 13980
散金: 20
帖子: 458
在线: 845.7小时
虫号: 57321
注册: 2005-01-29
性别: GG
专业: 合成药物化学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

第四步：用产生的序列文件产生一个alignment.ali脚本文件：把上面的序列复制到脚本中

请问这一步具体怎样操作呀，怎样用产生的序列文件产生一个alignment.ali脚本文件

赞一下(1人)

回复此楼

8楼2010-03-09 20:52:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sijifengsd

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 273.7
散金: 50
红花: 2
帖子: 94
在线: 12.1小时
虫号: 794988
注册: 2009-06-15
性别: GG
专业: 物理化学

★
lei0736(金币+1):谢谢 2010-04-10 10:14

引用回帖:

Originally posted by patent at 2010-03-09 20:52:22:
第四步：用产生的序列文件产生一个alignment.ali脚本文件：把上面的序列复制到脚本中

请问这一步具体怎样操作呀，怎样用产生的序列文件产生一个alignment.ali脚本文件

这个很简单，就是先建一个alignment.ali文件，然后把序列复制进去就行了

赞一下(1人)

回复此楼

似水年华

9楼2010-04-10 08:36:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

奋斗1s

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 1258.1
散金: 109
红花: 1
帖子: 415
在线: 106.7小时
虫号: 980571
注册: 2010-03-24
专业: 有机合成

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by sijifengsd at 2009-10-12 10:01:04:
最近专注于分子动力学的模拟研究。发现好多的蛋白质晶体结构或多或少的有一部分残基丢失了，要做动力学模拟，就需要把丢失的残基给补上。看文献发现比较公认的、又是免费的软件当属Modeller了。所以本人花了两个周 ...

你好，我现在在学习modeller，但是有很多地方不懂，我想问你一下：env.io.atom_files_directory = ['.', '../atom_files']中，① ['.', '../atom_files']中应该写什么？路径？还是什么都不用写？要写的话你能给我一个模板吗？②atom_files应该怎样准备？
麻烦你了，谢谢

赞一下(1人)

回复此楼

10楼2011-04-05 10:00:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 sijifengsd 的主题更新

返回列表