Yes, I am not using the same surface tension for both because I am following the recommendations of the manual for GB5 and PB4 respectively (page 699). Should I used different flavours of MMPBSA/MMGBSA? These settings seem to work for all other systems, it is just this one that has given me a peculiar result with PB specifically.
________________________________
From: Ray Luo <rluo.uci.edu>
Sent: 06 December 2016 17:03:46
To: AMBER Mailing List
Subject: Re: [AMBER] MMGBSA vs MMPBSA
Pietro,
For the NP part, you are not using the same coefficient ... 5
cal/mol-A^2 for GB and 7.3 cal/mol-A^2 for PB.
The delta G EPB and EGB are 519 kcal/mol and 505 kcal/mol,
respectively. They are different but apparently close.
All the best,
Ray
--
Ray Luo, Ph.D.
Professor
Biochemistry, Molecular Biophysics, Chemical Physics,
Chemical and Biomedical Engineering
University of California, Irvine, CA 92697-3900
On Mon, Dec 5, 2016 at 8:15 PM, Hai Nguyen <nhai.qn.gmail.com> wrote:
> hi
>
> I don't see anything wrong either. You can try with igb=8 too (recommend to
> use mbondi2 for GBSA in this case).
>
> Hai
>
> On Mon, Dec 5, 2016 at 11:10 PM, Aronica, Pietro <
> pietro.aronica07.imperial.ac.uk> wrote:
>
>> Hello,
>>
>> I have recently run a simulation of a simple protein receptor + peptide
>> ligand system and tried to calculate the MMGBSA and MMPBSA difference in
>> energy between the bound and unbound state.
>>
>> In my experience, the absolute value for the PB tends to be a few kcal/mol
>> lower than the GB one, but in this system this was not the case. While the
>> GB was -67.7195 ± 5.4461, the PB was -54.7742 ± 11.3891. This is not an
>> isolated result and the same pattern repeats across different runs. The
>> discrepancy and the much higher error of the PB give me pause and I have
>> been looking at the system to determine if something was wrong or if I made
>> a mistake.
>>
>> As far as I can tell I have used all the input files correctly, or at
>> least in the same way I usually do. What can I be on the lookout for? What
>> should I check in order to determine the cause of this bizarre result?
>> Should I even be worried?
>>
>> I also print the entire GB and PB output files, in case these are of any
>> help.
>>
>> If there are suggestions or comments, please do let me know.
>>
>> Cheers
>>
>> Pietro Aronica
>>
>>
>> GB:
>>
>> | Run on Tue Dec 6 00:03:58 2016
>> |
>> |Input file:
>> |--------------------------------------------------------------
>> |&general
>> |/
>> |&gb
>> | igb=5, surften=0.0050
>> |/
>> |&decomp
>> | idecomp=2, dec_verbose=1
>> |/
>> |--------------------------------------------------------------
>> |MMPBSA.py Version=14.0
>> |Solvated complex topology file: /home/pietroa/Charge_
>> derivation/GLU/MMPBSA/FF99SB_KEK_mbondi.prmtop
>> |Complex topology file: /home/pietroa/Charge_
>> derivation/GLU/MMPBSA/FF99SB_KEK_com_mbondi2.prmtop
>> |Receptor topology file: /home/pietroa/Charge_
>> derivation/GLU/MMPBSA/FF99SB_KEK_rec_mbondi2.prmtop
>> |Ligand topology file: /home/pietroa/Charge_
>> derivation/GLU/MMPBSA/FF99SB_KEK_lig_mbondi2.prmtop
>> |Initial mdcrd(s): ../Total_run.mdcrd
>> |
>> |Receptor mask: ":1-517"
>> |Ligand mask: ":518-520"
>> |
>> |Calculations performed using 100 complex frames.
>> |
>> |Generalized Born ESURF calculated using 'LCPO' surface areas
>> |
>> |All units are reported in kcal/mole.
>> ------------------------------------------------------------
>> -------------------
>> ------------------------------------------------------------
>> -------------------
>>
>> GENERALIZED BORN:
>>
>> Complex:
>> Energy Component Average Std. Dev. Std. Err. of
>> Mean
>> ------------------------------------------------------------
>> -------------------
>> VDWAALS -4461.2579 33.9886
>> 3.3989
>> EEL -39522.6040 144.1058
>> 14.4106
>> EGB -6004.7227 119.4563
>> 11.9456
>> ESURF 99.6832 2.0092
>> 0.2009
>>
>> G gas -43983.8620 147.5266
>> 14.7527
>> G solv -5905.0395 118.7491
>> 11.8749
>>
>> TOTAL -49888.9015 69.0225
>> 6.9023
>>
>>
>> Receptor:
>> Energy Component Average Std. Dev. Std. Err. of
>> Mean
>> ------------------------------------------------------------
>> -------------------
>> VDWAALS -4403.5819 33.5620
>> 3.3562
>> EEL -38747.6831 144.3997
>> 14.4400
>> EGB -6221.3471 120.6745
>> 12.0674
>> ESURF 103.4116 1.9760
>> 0.1976
>>
>> G gas -43151.2651 147.2434
>> 14.7243
>> G solv -6117.9355 120.0168
>> 12.0017
>>
>> TOTAL -49269.2005 68.8915
>> 6.8892
>>
>>
>> Ligand:
>> Energy Component Average Std. Dev. Std. Err. of
>> Mean
>> ------------------------------------------------------------
>> -------------------
>> VDWAALS -6.5483 1.2008
>> 0.1201
>> EEL -258.8440 13.4295
>> 1.3430
>> EGB -288.8132 12.4559
>> 1.2456
>> ESURF 2.2239 0.0509
>> 0.0051
>>
>> G gas -265.3922 12.9580
>> 1.2958
>> G solv -286.5893 12.4173
>> 1.2417
>>
>> TOTAL -551.9815 3.7165
>> 0.3717
>>
>>
>> Differences (Complex - Receptor - Ligand):
>> Energy Component Average Std. Dev. Std. Err. of
>> Mean
>> ------------------------------------------------------------
>> -------------------
>> VDWAALS -51.1277 4.2314
>> 0.4231
>> EEL -516.0770 17.5053
>> 1.7505
>> EGB 505.4375 15.2402
>> 1.5240
>> ESURF -5.9523 0.1398
>> 0.0140
>>
>> DELTA G gas -567.2047 16.5366
>> 1.6537
>> DELTA G solv 499.4852 15.2016
>> 1.5202
>>
>> DELTA TOTAL -67.7195 5.4461
>> 0.5446
>>
>>
>>
>> PB:
>>
>> | Run on Tue Dec 6 08:33:43 2016
>> |
>> |Input file:
>> |--------------------------------------------------------------
>> |&general
>> |/
>> |&pb
>> | inp=1, radiopt=0, cavity_surften=0.0072
>> |/
>> |&decomp
>> | idecomp=2, dec_verbose=1
>> |/
>> |--------------------------------------------------------------
>> |MMPBSA.py Version=14.0
>> |Solvated complex topology file: /home/pietroa/Charge_
>> derivation/GLU/MMPBSA/FF99SB_KEK_mbondi.prmtop
>> |Complex topology file: /home/pietroa/Charge_
>> derivation/GLU/MMPBSA/FF99SB_KEK_com_mbondi.prmtop
>> |Receptor topology file: /home/pietroa/Charge_
>> derivation/GLU/MMPBSA/FF99SB_KEK_rec_mbondi.prmtop
>> |Ligand topology file: /home/pietroa/Charge_
>> derivation/GLU/MMPBSA/FF99SB_KEK_lig_mbondi.prmtop
>> |Initial mdcrd(s): ../Total_run.mdcrd
>> |
>> |Receptor mask: ":1-517"
>> |Ligand mask: ":518-520"
>> |
>> |Calculations performed using 100 complex frames.
>> |Poisson Boltzmann calculations performed using internal PBSA solver in
>> sander.
>> |
>> |All units are reported in kcal/mole.
>> ------------------------------------------------------------
>> -------------------
>> ------------------------------------------------------------
>> -------------------
>>
>> POISSON BOLTZMANN:
>>
>> Complex:
>> Energy Component Average Std. Dev. Std. Err. of
>> Mean
>> ------------------------------------------------------------
>> -------------------
>> VDWAALS -4461.2579 33.9886
>> 3.3989
>> EEL -39522.6040 144.1058
>> 14.4106
>> EPB -5966.1358 119.2844
>> 11.9284
>> ENPOLAR 156.8866 1.7105
>> 0.1710
>>
>> G gas -43983.8620 147.5266
>> 14.7527
>> G solv -5809.2492 118.3545
>> 11.8354
>>
>> TOTAL -49793.1112 75.5242
>> 7.5524
>>
>>
>> Receptor:
>> Energy Component Average Std. Dev. Std. Err. of
>> Mean
>> ------------------------------------------------------------
>> -------------------
>> VDWAALS -4403.5819 33.5620
>> 3.3562
>> EEL -38747.6831 144.3997
>> 14.4400
>> EPB -6196.9759 123.0006
>> 12.3001
>> ENPOLAR 158.8439 1.7092
>> 0.1709
>>
>> G gas -43151.2651 147.2434
>> 14.7243
>> G solv -6038.1320 122.1447
>> 12.2145
>> TOTAL -49189.3971 72.9118
>> 7.2912
>>
>>
>> Ligand:
>> Energy Component Average Std. Dev. Std. Err. of
>> Mean
>> ------------------------------------------------------------
>> -------------------
>> VDWAALS -6.5483 1.2008
>> 0.1201
>> EEL -258.8440 13.4295
>> 1.3430
>> EPB -288.9018 12.0239
>> 1.2024
>> ENPOLAR 5.3541 0.0626
>> 0.0063
>>
>> G gas -265.3922 12.9580
>> 1.2958
>> G solv -283.5477 11.9841
>> 1.1984
>>
>> TOTAL -548.9399 3.7137
>> 0.3714
>>
>>
>> Differences (Complex - Receptor - Ligand):
>> Energy Component Average Std. Dev. Std. Err. of
>> Mean
>> ------------------------------------------------------------
>> -------------------
>> VDWAALS -51.1277 4.2314
>> 0.4231
>> EEL -516.0770 17.5053
>> 1.7505
>> EPB 519.7418 15.9509
>> 1.5951
>> ENPOLAR -7.3113 0.1295
>> 0.0129
>> EDISPER 0.0000 0.0000
>> 0.0000
>>
>> DELTA G gas -567.2047 16.5366
>> 1.6537
>> DELTA G solv 512.4305 15.9443
>> 1.5944
>>
>> DELTA TOTAL -54.7742 11.3891
>> 1.1389
>>
>> _______________________________________________
>> AMBER mailing list
>> AMBER.ambermd.org
>> http://lists.ambermd.org/mailman/listinfo/amber
>>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Dec 06 2016 - 18:30:02 PST