amber中非标准氨基酸残基的参数生成

在amber教程中给出了相关的教程，具体流程都是根据这个教程来制作的。

http://ambermd.org/tutorials/basic/tutorial5/

但是他的电荷使用的是bcc电荷，但是并没有介绍如何把电荷转换为resp。

流程：

0.在pdb数据库中搜索残基名称，选择 ligand id 进行下载对应的cif文件；

1.antechamber -fi ccif -i CRO.cif -bk CRO -fo gcrt -o cro.gjf -ch "cro.chk" -gm "%mem=2048MB" -gn "%nproc=24" -nc 0 #直接将下载的cif转为高斯输入文件，进行高斯计算；

2.antechamber -fi ccif -i CRO.cif -bk CRO -fo ac -o cro_step1.ac -at amber #和教程上一致，目的为了生成参照的ac文件，-at 必须制定，否则原子类型默认gaff；

3.antechamber -fi gout -i cro.log -rn CRO -fo ac -o cro_step2.ac -c resp -at amber #赋予resp电荷到新的ac文件中；

4.文本操作将cro_step2.ac中的坐标以及电荷信息列覆盖到cro_step1.ac中，同时修改主链N原子的原子类型从NT到N，保存为cro_step3.ac；

5.制作mc文件，格式如下图，参照cro_step3.ac文件制作；