计算化学公社

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

[Amber] pdb转化top和crd文件出错

[复制链接 Copy URL]

41

帖子

0

威望

181

eV
积分
222

Level 3 能力者

想问一下各位老师,这种情况应该怎么解决呢
我的pdb文件增加了二硫键,第一个氨基酸便是半胱氨酸,显示必须是atom,回去查看了自己的pdb文件,不管把二硫键的位置放在pdb文件的开头还是末尾,都显示此错误


207d85c8f77208f66871674d37099cc.png (44.55 KB, 下载次数 Times of downloads: 19)

错误信息

错误信息

d53c097bcdad4759f179b04b41c89e9.png (106.61 KB, 下载次数 Times of downloads: 17)

pdb文件

pdb文件

90

帖子

0

威望

1686

eV
积分
1776

Level 5 (御坂)

2#
发表于 Post on 2023-10-20 11:00:46 | 只看该作者 Only view this author
本帖最后由 Serious 于 2023-10-20 11:01 编辑

1)Error信息写得听清楚的,不需要加”SSBOND“;要不改残基名,要不就在tleap加bond(反正也会自动检测);

41

帖子

0

威望

181

eV
积分
222

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-10-20 14:44:35 | 只看该作者 Only view this author
因为我在分析的时候需要分析二硫键,所以我把CYS已经改为CYX,同时tleap文件也已经修改了,但是就是显示此错误

310

帖子

0

威望

1503

eV
积分
1813

Level 5 (御坂)

4#
发表于 Post on 2023-10-20 16:18:31 | 只看该作者 Only view this author
molly85 发表于 2023-10-20 14:44
因为我在分析的时候需要分析二硫键,所以我把CYS已经改为CYX,同时tleap文件也已经修改了,但是就是显示此 ...

错误提示指出说tleap中bond命令用法不对,所以你可以把tleap的输入文件或者命令发出来看看。

41

帖子

0

威望

181

eV
积分
222

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-10-21 13:50:05 | 只看该作者 Only view this author
本帖最后由 molly85 于 2023-10-21 13:52 编辑
Frozen-Penguin 发表于 2023-10-20 16:18
错误提示指出说tleap中bond命令用法不对,所以你可以把tleap的输入文件或者命令发出来看看。

老师,这是我的tleap文件。

Screenshot from 2023-10-21 13-43-38.png (74.1 KB, 下载次数 Times of downloads: 17)

Screenshot from 2023-10-21 13-43-38.png

310

帖子

0

威望

1503

eV
积分
1813

Level 5 (御坂)

6#
发表于 Post on 2023-10-21 20:43:35 | 只看该作者 Only view this author
molly85 发表于 2023-10-21 13:50
老师,这是我的tleap文件。

大多数模拟软件会重新设置残基编号以方便后续计算,不会沿用PDB中的残基编号。
在使用bond定义键时,中间的数字是残基编号,要用重新设定之后的。
可以读取pdb后直接保存一个pdb,这样就可以看到重新设定后的残基编号。

41

帖子

0

威望

181

eV
积分
222

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-10-21 21:23:57 | 只看该作者 Only view this author
Frozen-Penguin 发表于 2023-10-21 20:43
大多数模拟软件会重新设置残基编号以方便后续计算,不会沿用PDB中的残基编号。
在使用bond定义键时,中 ...

老师,用什么命令读取原始pdb文件之后会改变残基序号呀,之前没遇到过,然后新的pdb文件生成之后需要改动嘛?就是半胱氨酸。

310

帖子

0

威望

1503

eV
积分
1813

Level 5 (御坂)

8#
发表于 Post on 2023-10-22 10:00:03 | 只看该作者 Only view this author
molly85 发表于 2023-10-21 21:23
老师,用什么命令读取原始pdb文件之后会改变残基序号呀,之前没遇到过,然后新的pdb文件生成之后需要改动 ...

直接读取然后保存pdb就行了,残基号会自动重新设定好
  1. source leaprc.protein.ff19SB
  2. mol = loadpdb input.pdb
  3. savepdb mol output.pdb
  4. quit
复制代码

41

帖子

0

威望

181

eV
积分
222

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-10-22 12:20:21 | 只看该作者 Only view this author
Frozen-Penguin 发表于 2023-10-22 10:00
直接读取然后保存pdb就行了,残基号会自动重新设定好

多谢老师,问题在您的指导下已经解决了

3

帖子

0

威望

79

eV
积分
82

Level 2 能力者

10#
发表于 Post on 2023-12-8 14:55:50 | 只看该作者 Only view this author
请问这是在什么软件输入的命令?

41

帖子

0

威望

181

eV
积分
222

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2023-12-18 16:23:46 | 只看该作者 Only view this author
liangxf 发表于 2023-12-8 14:55
请问这是在什么软件输入的命令?

amber

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 08:51 , Processed in 0.542697 second(s), 23 queries , Gzip On.

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