前段时间我想做小分子化合物与蛋白共价结合的蛋白复合物的动力学,发现amber可以做,但我一直在学gromacs,而gromacs关于此方面的东西比较少,最近我突然想到可以把共价结合的小分子化合物当做非标准氨基酸处理,就看了一些网上的一些知识,还是一头雾水。举个例子,如果我想用amber99SB力场在gromacs下跑,那这个非标准氨基酸我应该怎么得到它的力场参数?
我在gromacs官网中看到有关于添加新力场参数的做法,http://www.gromacs.org/Documentation/How-tos/Adding_a_Residue_to_a_Force_Field,但这个对我来说太笼统,让我关注的就只有这样一句话,You may need to use an external topology generation tool and adapt the results to the .rtp format. 这句话的意思是不是,举个例子,我可用GAFF力场写出这个非标准氨基酸的top文件,然后根据你想要得到的amber99SB力场参数形式修改前面得到的原子类型等参数,最后按照amber99SB力场的参数格式paste到力场文件中就可以了?