转载请联系,且保留链接

翻译自:MOLINFO

目标

此命令主要是提供使用信息在已存在的系统中

流程

plumed提供了三种识别不同的残基模式:MOLTYPE=protein,MOLTYPE=rna,MOLTYPE=dna
MOLINFO进行蛋白或者核酸pdb文件操作可以使用@特征符:

1
2
@"definition"-chainresiduenum
@"definition"-residuenum

例子如下:

1
2
@psi-1 will select the atoms defining the psi torsion of residue 1
@psi-C1 will define the same torsion for residue 1 of chain C.

下面列出现在可用的 “definition”
对于蛋白残基

1
2
3
4
@phi-#
@psi-#
@omega-#
@chi1-#

选择合适的残基,定义每个残基#的二面角

对于DNA或者RNA,下面的组是可用的:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
# quadruplets for backbone dihedral angles
@alpha-#
@beta-#
@gamma-#
@delta-#
@epsilon-#
@zeta-#
# quadruplets for sugar dihedral angles
@v0-#
@v1-#
@v2-#
@v3-#
@v4-#
# quadruplet corresponding to the chi torsional angle
@chi-#
# backbone, sugar, and base heavy atoms
@back-#
@sugar-#
@base-#
# ordered triplets of atoms on the 6-membered ring of nucleobases
# namely:
# C2/C4/C6 for pyrimidines
# C2/C6/C4 for purines
@lcs-#

当然也可以使用单原子@atom-chainresiduenum or @atom-residuenum.

例子

此命令提供蛋白哪些原子为骨架:

1
2
MOLINFO STRUCTURE=reference.pdb
ALPHARMSD RESIDUES=all TYPE=DRMSD LESS_THAN={RATIONAL R_0=0.08 NN=8 MM=12} LABEL=a

下面的例子为打印GC Watson-Crick 配对的氢键练习

1
2
3
4
5
MOLINFO STRUCTURE=reference.pdb
hb1: DISTANCE ATOMS=@N2-1,@O2-14
hb2: DISTANCE ATOMS=@N1-1,@N3-14
hb3: DISTANCE ATOMS=@O6-1,@N4-14
PRINT ARG=hb1,hb2,hb3

下面的例子为计算转角:

1
2
3
4
MOLINFO MOLTYPE=protein STRUCTURE=myprotein.pdb
t1: TORSION ATOMS=@phi-3
t2: TORSION ATOMS=@psi-4
PRINT ARG=t1,t2 FILE=colvar STRIDE=10

简单说明:@phi-3告诉plumed你需要计算第三个蛋白残基的phi角,类似的@psi-4告诉plumed你需要计算第四个蛋白残基的psi角