Charmm-note

charmm_executable < charmm_script.inp > charmm_output.out

charmm –i 1.inp –o 1.out

RTF(residue topology file):包含了成键和电荷信息。Charmm包含了核酸,磷脂、蛋白和糖类的标准RTF文件

PARM(Parameter File):决定了能量:键长键角二面角力常数和范德华参数。Charmm包含了核酸,磷脂、蛋白、糖类和水的标准PARM文件

CRD(coordinates):体系的标准笛卡尔坐标

PSF(protein structure file):hold了成键和非键列表,键长键角二面角。对计算能量至关重要

标题必须有:

* This would be a short title
* 

注释:

! this is a comment on a line by itself
 
结束:
Stop
 
 
 
蛋白,核酸,磷脂等的力场参数:
Charmm的力场参数要自行下载,放到输入文件所在文件夹下面供其调用https://www.charmm.org/charmm/showcase/news/new-protein-force-field/
 
小分子的力场参数:
1、通过amber的力场参数转化……………………blabla很难受!!但是,因为目前的任务是想把amber动力学之后的结构,用charmm来做计算,所以最好的方法是做力场转化。如果不转化的话,可能charmm的方法生成的小分子的原子类型之类的对不上,导致出错。(比如说加氢,amber是用的reduce,但是charmm我用的是gv保存的mol2)
2、charmm-gui(http://www.charmm-gui.org/?doc=input/ligandrm),提供了小分子的名字,可以自动生成其rtf文件和prm文件。PS:所用方法是CGenFF1.0版本。
3、直接用CgenFF(charmm小分子制作的标准方法,相当于antechamber)(https://cgenff.paramchem.org/initguess/#20181120_3/rg1.str
 
(1)使用方法是先用gv保存分子的mol2文件,然后将mol2文件中的所有Ar替换为ar,在Molecule Name一行改为小分子的名字,然后导入分子,就可以直接生成。
(2)默认使用最新版本,目前的版本是2.2.0(使用的CgenFF版本是4.0),可以选择使用1.0版本(使用的CgenFF版本是3.0.1),输出结果与charmm-gui结果是一样的。
 

如:全部使用1.0版本:

带有小分子的蛋白的动力学:

1、下载最新版本的蛋白质力场

2、下载与生成小分子相匹配的CGenFF力场,如:使用的CgenFF版本是2.2.0,则需要使用CgenFF的4.0版本。

3、下载水和离子的力场toppar_water_ions_namd.str

4、得到小分子力场my_ligands.str

5、读取参数文件和坐标

! read topology and parameter files

(1)With recent CHARMM force field releases, you can just read in the topology and parameter files as follows:

read rtf card name @TOPPAR/top_all36_prot.rtf

read param card flex name @TOPPAR/par_all36_prot.prm

read rtf card append name @TOPPAR/top_all36_cgenff.rtf

read para card flex append name @TOPPAR/par_all36_cgenff.prm

stream toppar_water_ions.str

(2)(optional) read in your own stream file:

stream my_ligands.str

! Read PSF and coordinates from file

(1)read psf card name all.psf

      read coor card name all.crd

PS:可通过这个例子中的方法生成蛋白整体的psf和crd(https://www.charmmtutorial.org/index.php/Full_example)

原文地址:https://www.cnblogs.com/jszd/p/11178806.html