Hi Ray,
Thank you for your reply  
I used Amber18/AmberTools18. 
The final system parameters are: 
 NSTEP = 49990000   TIME(PS) =  111980.000  TEMP(K) =   309.04  PRESS =     0.0
 Etot   =   -264550.3893  EKtot   =     68926.3047  EPtot      =   -333476.6940
 BOND   =      2659.8247  ANGLE   =      7371.0874  DIHED      =     11181.7492
 1-4 NB =      3195.8521  1-4 EEL =     32048.3876  VDWAALS    =     36214.6452
 EELEC  =   -426148.2401  EHBOND  =         0.0000  RESTRAINT  =         0.0000
 EKCMT  =         0.0000  VIRIAL  =         0.0000  VOLUME     =   1091009.1860
                                                    Density    =         1.0239
And here is the input file used for the production: 
prod 
 &cntrl
   imin=0,
   irest=1,
   ntx=5,
   nstlim=50000000,
   dt=0.002,
   ntc=2,
   ntf=2, 
   tol=0.00001,
   cut=10.0, 
   ntb=2,
   ntp=1,
   taup=2.0,
   ntpr=10000, 
   ntwx=10000, 
   ntwr = 10000, 
   ntwe = 10000, 
   iwrap=1,
   ntt=3, 
   barostat=2, 
   gamma_ln=2.0,
   temp0=310.15, 
   tempi=310.15, 
   ig=-1,
/
Best regards, 
Harutyun 
-------- Forwarded message -------
From: "Ray Luo" <rluo.uci.edu>
To: "AMBER Mailing List" <amber.ambermd.org>
Sent: October 1, 2020 2:57 AM
Subject: Re: [AMBER] controversy in PBSA and GBSA results
Harutyun,
Which version of Amber are you using? Also, it's better for us to see
your detailed energy output ...
All the best,
Ray
--
Ray Luo, Ph.D.
Professor of Structural Biology/Biochemistry/Biophysics,
Chemical and Materials Physics, Chemical and Biomolecular Engineering,
Biomedical Engineering, and Materials Science and Engineering
Department of Molecular Biology and Biochemistry
University of California, Irvine, CA 92697-3900
On Wed, Sep 30, 2020 at 4:44 AM <h_sahakyan.mb.sci.am> wrote:
> Dear All,
> 
> I used MM-PBSA/GBSA methods to calculate the binding free energy between two subunits of a protein.
> The binding energies were -101 and -94 kcal/mol for GBSA and PBSA respectively. However, in the
> presents of a ligand (between two subunits), binding energies were -124/-67 (GBSA/PBSA). As you can
> see, according to GBSA the ligand makes the interaction stronger and vice versa according to the
> PBSA calculation. What can be the reason for this controversy?
> 
> I used MMPBSA.py with following input:
> 
> pbsa_gbsa
> &general
> startframe=1, interval=50, endframe=5000,
> verbose=2, keep_files=0, netcdf=0,
> strip_mask = ":WAT,Cl-,Na+,K+,",
> ligand_mask = ":1-437,867,869",
> receptor_mask = ":438-865,866,868,870",
> /
> &gb
> saltcon=0.150,
> /
> &pb
> inp=1 # I tried inp=2 it gives more weird results, generally with positive energy values
> istrng=0.150
> /
> Best wishes,
> Harutyun
> Senior Assistant at Laboratory of
> Computational Modeling of Biological Processes,
> Institute of Molecular Biology, NAS RA
> 7 Hasratyan St, Yerevan, Armenia
> Tel: +374 93 323990
> E-mail: h_sahakyan.mb.sci.am (h_sahakyan.mb.sci.am)
> _______________________________________________
> 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
-- 
Senior Assistant at Laboratory of
Computational Modeling of Biological Processes,
Institute of Molecular Biology, NAS RA
7 Hasratyan St, Yerevan, Armenia
Tel: +374 93 323990
E-mail: h_sahakyan.mb.sci.am
October 1, 2020 2:57 AM, "Ray Luo" <rluo.uci.edu> wrote:
> Harutyun,
> 
> Which version of Amber are you using? Also, it's better for us to see
> your detailed energy output ...
> 
> All the best,
> Ray
> --
> Ray Luo, Ph.D.
> Professor of Structural Biology/Biochemistry/Biophysics,
> Chemical and Materials Physics, Chemical and Biomolecular Engineering,
> Biomedical Engineering, and Materials Science and Engineering
> Department of Molecular Biology and Biochemistry
> University of California, Irvine, CA 92697-3900
> 
> On Wed, Sep 30, 2020 at 4:44 AM <h_sahakyan.mb.sci.am> wrote:
> 
>> Dear All,
>> 
>> I used MM-PBSA/GBSA methods to calculate the binding free energy between two subunits of a protein.
>> The binding energies were -101 and -94 kcal/mol for GBSA and PBSA respectively. However, in the
>> presents of a ligand (between two subunits), binding energies were -124/-67 (GBSA/PBSA). As you can
>> see, according to GBSA the ligand makes the interaction stronger and vice versa according to the
>> PBSA calculation. What can be the reason for this controversy?
>> 
>> I used MMPBSA.py with following input:
>> 
>> pbsa_gbsa
>> &general
>> startframe=1, interval=50, endframe=5000,
>> verbose=2, keep_files=0, netcdf=0,
>> strip_mask = ":WAT,Cl-,Na+,K+,",
>> ligand_mask = ":1-437,867,869",
>> receptor_mask = ":438-865,866,868,870",
>> /
>> &gb
>> saltcon=0.150,
>> /
>> &pb
>> inp=1 # I tried inp=2 it gives more weird results, generally with positive energy values
>> istrng=0.150
>> /
>> Best wishes,
>> Harutyun
>> Senior Assistant at Laboratory of
>> Computational Modeling of Biological Processes,
>> Institute of Molecular Biology, NAS RA
>> 7 Hasratyan St, Yerevan, Armenia
>> Tel: +374 93 323990
>> E-mail: h_sahakyan.mb.sci.am (h_sahakyan.mb.sci.am)
>> _______________________________________________
>> 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 Wed Sep 30 2020 - 21:30:02 PDT