计算化学公社

 找回密码 Forget password
 注册 Register
Views: 20|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] mmPBSA计算氨基酸残基贡献能太大正常吗?

[复制链接 Copy URL]

5

帖子

0

威望

35

eV
积分
40

Level 2 能力者

请教:我在通过gmxPBSA计算蛋白与蛋白的氨基酸残基贡献能太大正常吗?
这是计算的蛋白与蛋白的pb结合能:

Delta (Complex - Receptor - Ligand):
Energy Component       Average     SD(Prop.)         SD   SEM(Prop.)        SEM
-------------------------------------------------------------------------------
螖BOND                     0.00         18.96       0.00         1.89       0.00
螖ANGLE                    0.00         23.59       0.00         2.35       0.00
螖DIHED                    0.00         12.00       0.00         1.19       0.00
螖UB                      -0.00          3.94       0.00         0.39       0.00
螖IMP                      0.00          4.58       0.00         0.46       0.00
螖CMAP                    -0.00          5.44       0.00         0.54       0.00
螖VDWAALS               -285.55         12.42       9.85         1.24       0.98
螖EEL                    142.08         57.00      48.51         5.67       4.83
螖1-4 VDW                  0.00         13.07       0.00         1.30       0.00
螖1-4 EEL                 -0.00         27.37       0.00         2.72       0.00
螖EPB                     -0.59         29.34      44.47         2.92       4.43
螖ENPOLAR                -34.84          0.69       0.64         0.07       0.06
螖EDISPER                  0.00          0.00       0.00         0.00       0.00

螖GGAS                  -143.47         58.91      43.86         5.86       4.36
螖GSOLV                  -35.44         29.35      44.63         2.92       4.44

螖TOTAL                 -178.91         65.81      16.12         6.55       1.60
这是计算的氨基酸残基贡献能:
        Avg.        Std. Dev.        Std. Err. of Mean
R:A:ARG:205        -246.40303        14.97752        4.51589
L:B:ARG:1262        -245.95876        5.67422        1.71084
R:A:ARG:863        -245.45585        12.04191        3.63077
R:A:ARG:904        -238.4797        7.84836        2.36637
L:B:ARG:1137        -237.73359        8.08867        2.43883
R:A:ARG:869        -236.88081        15.29091        4.61038
R:A:ARG:864        -236.293        16.65653        5.02213
R:A:ARG:192        -233.16172        13.47289        4.06223
R:A:ARG:191        -233.1054        13.67005        4.12168
L:B:ARG:1134        -225.61067        6.64485        2.0035
R:A:ARG:546        -223.92876        22.24628        6.7075
L:B:ASP:1275        -96.02517        21.38084        6.44657
R:A:ASP:528        -91.22533        22.29914        6.72344
R:A:ASP:193        -87.77585        12.7511        3.8446
R:A:ASP:868        -87.56817        26.10986        7.87242
L:B:GLU:1115        -76.40435        23.61997        7.12169
L:B:GLU:1124        -75.1394        21.51473        6.48694
L:B:TYR:1287        -63.23057        21.39722        6.45151
R:A:ASN:940        -62.18934        2.40581        0.72538
L:B:ASN:1177        -58.71483        4.09924        1.23597
R:A:ASN:871        -55.64511        4.02612        1.21392
R:A:ASN:982        -51.09245        4.29388        1.29465
R:A:ASN:194        -50.61462        4.74084        1.42942
R:A:GLN:872        -37.12881        3.28985        0.99193
L:B:GLN:1059        -34.23394        4.05239        1.22184
R:A:GLN:937        -33.73124        3.78774        1.14205
L:B:GLN:1117        -32.81811        4.47229        1.34845
L:B:GLN:1060        -32.32123        5.28173        1.5925
L:B:GLN:1267        -28.67152        3.93577        1.18668
L:B:LYS:1274        -27.85976        21.51852        6.48808
R:A:LYS:867        -26.71306        22.98879        6.93138
L:B:LYS:1062        -25.4495        10.20154        3.07588
L:B:LYS:1090        -25.18064        17.89742        5.39628
L:B:LYS:1239        -24.98858        14.34851        4.32624
L:B:LYS:1272        -23.31091        8.26301        2.49139
R:A:LYS:198        -18.80256        8.83462        2.66374
R:A:HIS:201        -9.02112        4.29438        1.29481
L:B:LYS:1181        -8.24005        13.81447        4.16522
R:A:HIS:260        -6.03911        3.65348        1.10157
L:B:THR:1171        -2.86182        2.444        0.73689
L:B:MET:1056        -1.90861        13.38203        4.03483
这是mmpbsa_pb.in文件参数,计算的是平衡的后10ns的结合能:
&general
sys_name="Prot-Prot-CHARMM",
startframe=4500,
endframe=99999,
interval=5,
PBRadii=7,
verbose=2,
/

&pb
istrng=0.15,
fillratio=4.0,
radiopt=0,
inp=1,
use_sav=0,
/

&decomp
idecomp=2,
dec_verbose=3,
print_res="within 5",
csv_format=1,
/



202608181928542523..png (36.28 KB, 下载次数 Times of downloads: 0)

202608181928542523..png

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-8-19 00:43 , Processed in 0.494310 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list