求助：修改蛋白质pdb金属离子文件后金属离子位置发生变化

查看数: 3625 | 评论数: 6 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2022-3-16 10:43

使用GROMOS96 54a7力场对含有两个镁离子的蛋白质进行分子动力学模拟时，因为力场的aminoacids.rtp文件中镁离子的名称是MG2+，从而将蛋白质pdb文件中的镁离子名称改成MG2+，但发现镁离子跑到蛋白质外面，并且用该文件 ...

回复 Reply

wangdafeixh 发表于 Post on 2022-3-21 15:28:34

Frozen-Penguin 发表于 2022-3-21 13:08
虽然坐标的位置对了，但是可能残基名称还没有对齐，A（链编号）要和上面对齐，MG2+左边和上面的残基名称 ...

已经解决问题，太感谢您了！

Frozen-Penguin 发表于 Post on 2022-3-21 13:08:57

wangdafeixh 发表于 2022-3-21 09:59
谢谢老师的回答，我按照老师的方法将pdb文件的镁离子调整到蛋白质内了，但是用修改好的pdb文件生成gro文 ...

虽然坐标的位置对了，但是可能残基名称还没有对齐，A（链编号）要和上面对齐，MG2+左边和上面的残基名称对齐，右边和A连在一起，中间没有空格。

具体的pdb格式要求：https://www.cgl.ucsf.edu/chimera ... rials/pdbintro.html

Frozen-Penguin 发表于 Post on 2022-3-18 00:54:26

pdb文件是按列读取的，而不是看空格，加了“2+”以后后面的坐标往后移了2格，所以读取的时候都歪了
比如pdb里面的32.892在gro里面变成了3.28和9.2
把后面的内容移回来就可以了

wangdafeixh 发表于 Post on 2022-3-17 09:44:20

sobereva 发表于 2022-3-16 16:36
把盒子显示出来，然后用VMD观看轨迹，看是否是越过了盒子导致的

谢谢老师回答，可是这是在进行模拟之前，准备蛋白质结构时就发生的情况，应该怎么导出轨迹呢？

sobereva 发表于 Post on 2022-3-16 16:36:05

把盒子显示出来，然后用VMD观看轨迹，看是否是越过了盒子导致的

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register