计算化学公社

标题: 求助：如何将L型氨基酸转为D氨基酸进行模拟 [打印本页]

作者
Author: 小熊也有梦想 时间: 2023-7-16 12:05
标题: 求助：如何将L型氨基酸转为D氨基酸进行模拟
老师们好，我目前已经有一个8-L型氨基酸聚合物的模拟体系在NAMD上跑模拟，现在想做一个镜像体系，请问要如何获得对应的8-D型氨基酸聚合物的PDB、PSF、拓扑文件和参数文件呢，有这方面的教程吗？望各位老师不吝赐教！谢谢。

作者
Author: ionexchangeC 时间: 2023-7-16 22:34
所有原子的X坐标和Y坐标不变，Z坐标乘以-1

作者
Author: 喵星大佬 时间: 2023-7-17 00:49
pdb文件直接用gromacs转一下就行，gmx editconf -scal 1 1 -1

拓扑和参数动都不用动

作者
Author: 小熊也有梦想 时间: 2023-7-30 08:56
谢谢各位老师，目前问题已经解决了。后来在一篇文献中看到说，在Charmm力场中D氨基酸和L氨基酸的CMAP需要倒转，charmm提供了倒转后的力场文件。

因此我用老师们的方法改变坐标后，用新的力场文件重新构建了体系。希望给其他有同样问题的小伙伴提供一些参考。

作者
Author: coisiniiiiiii 时间: 2025-12-24 15:49

小熊也有梦想发表于 2023-7-30 08:56
谢谢各位老师，目前问题已经解决了。后来在一篇文献中看到说，在Charmm力场中D氨基酸和L氨基酸的CMAP需要倒 ...

您好，方便分享一下文献嘛？

作者
Author: student0618 时间: 2025-12-24 18:32
在CHARMM36m力场中，D-氨基酸的残基名是DXXX。可以手动编辑PDB文件中的残基名直接使用。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)