|
|
请教:我在通过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,
/
|
|