计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1049|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[Amber] 求助:怎么搭建有2条链的蛋白?

[复制链接 Copy URL]

113

帖子

0

威望

377

eV
积分
490

Level 3 能力者

老师好,

我想用amber软件跑一个蛋白和配体和下游蛋白的模拟。建模的时候我分别把受体、下游蛋白、配体按顺序粘在了PDB里,中间加TER。然后用tleap生成prmtop、inpcrd。之后存出复合物的PDB。

但是我看到tleap把受体和下游蛋白当成一个蛋白了,中间加了一段loop。如下图:

老师,请问怎么让tleap把受体和下游蛋白当成2个蛋白处理?

谢谢老师。

202312241747328292..png (29.21 KB, 下载次数 Times of downloads: 31)

202312241747328292..png

33

帖子

0

威望

697

eV
积分
730

Level 4 (黑子)

2#
发表于 Post on 2023-12-25 08:54:05 | 只看该作者 Only view this author
tleap中直接导出复合物试试

113

帖子

0

威望

377

eV
积分
490

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-12-25 14:24:05 | 只看该作者 Only view this author
你好,您是说用tleap的combine指令吗

1

帖子

0

威望

160

eV
积分
161

Level 3 能力者

4#
发表于 Post on 2023-12-25 17:41:28 | 只看该作者 Only view this author
这应该是个pymol的显示问题。要确定有没有出问题的话,可以看tleap的输出,比如你的两个蛋白都是G开始A结束,那么可以看到NGLY和CALA都是2。那这样应该没问题。还可以使用vmd打开top和crd,看看虚线处是否有成键。

113

帖子

0

威望

377

eV
积分
490

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-12-25 18:41:46 | 只看该作者 Only view this author
本帖最后由 灵芝5 于 2023-12-25 18:45 编辑

好的,谢谢老师解答。

另外,老师想请问您一下,如果想对这2个蛋白的氨基酸进行限制,是不是在in文件里连续写残基号就可以了。
比如,第一个蛋白残基号是1-395,第2个蛋白残基号是396-543,只需要在in文件里写restraintmask = '(:1-543)’就可以了。

241

帖子

0

威望

1077

eV
积分
1318

Level 4 (黑子)

6#
发表于 Post on 2023-12-27 10:16:12 | 只看该作者 Only view this author
灵芝5 发表于 2023-12-25 18:41
好的,谢谢老师解答。

另外,老师想请问您一下,如果想对这2个蛋白的氨基酸进行限制,是不是在in文件里 ...

113

帖子

0

威望

377

eV
积分
490

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-12-27 16:31:14 | 只看该作者 Only view this author
好的,谢谢老师解答

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 02:58 , Processed in 0.205652 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list