计算化学公社

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

[GROMACS] gmx2pdb报错提示的残基编号对应的氨基酸与pdb中对不上

[复制链接 Copy URL]

262

帖子

0

威望

635

eV
积分
897

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 KazusaT 于 2025-9-19 01:12 编辑

各位老师好,最近在进行含有非标准氨基酸蛋白的模拟,使用sobtop工具构建了rtp条目后随便搞了一个pdb用来测试rtp条目和hdb文件编写是否正确,使用的命令如下:
gmx pdb2gmx -f .\testrtp.pdb -o test.gro -ff amber19sb -ignh
之后报错:
Fatal error:
Residue 3 named LMA of a molecule in the input file was mapped
to an entry in the topology database, but the atom C used in
that entry is not found in the input file. Perhaps your atom
and/or residue naming needs to be fixed.


LMA为编辑过的赖氨酸,即我使用的非标准残基,然而我用于测试的文件是GLY-LMA-GLY,residue3对应的根本不是LMA,请问这种情况怎么解决?

测试pdb作为附件上传,rtp中LMA相关条目为:
[ LMA ]
  [ atoms ]
  ; Name        type          charge    chg grp
     N1          N         -0.3976005430      1
     H2          H          0.2652278046      2
     C3          CT        -0.2221055101      3
     H4          H1         0.1297590841      4
     C5          CT         0.1628087828      5
     C6          C          0.6515984075      6
     H7          HC         0.0003253387      7
     H8          HC         0.0003253387      8
     C9          CT        -0.0512850643      9
     O10         O         -0.5624298351     10
     H11         HC         0.0143695187     11
     H12         HC         0.0143695187     12
     C13         CT        -0.1434419752     13
     H14         HC         0.0351074812     14
     H15         HC         0.0351074812     15
     C16         CT         0.4698012210     16
     H17         H1        -0.0462972567     17
     H18         H1        -0.0462972567     18
     N19         N         -0.6239073109     19
     H20         H          0.3238520352     20
     C21         CT         0.0789216254     21
     H22         HC         0.0180339240     22
     C23         CT        -0.1503586640     23
     H24         HC         0.0417268001     24
     H25         HC         0.0417268001     25
     H26         HC         0.0417268001     26
     C27         CT        -0.1620758962     27
     H28         HC         0.0443269901     28
     H29         HC         0.0443269901     29
     H30         HC         0.0443269901     30
     C31         C          0.5074289788     31
     O32         O         -0.5593985989     32

后面为bond、angle等

testrtp.pdb

3.88 KB, 下载次数 Times of downloads: 2

aminoacids.hdb

20.17 KB, 下载次数 Times of downloads: 0

aminoacids.rtp

98.12 KB, 下载次数 Times of downloads: 0

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 15:33 , Processed in 0.250194 second(s), 23 queries , Gzip On.

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