[AMBER] MMPBSA: Per Residue vs Pairwise decomposition Calculation(idecomp= 2 o idecomp=4)

From: Rodriguez, Yoel <yoel.rodriguez.mssm.edu>
Date: Tue, 27 May 2014 17:22:06 +0000

Hi Amber Users,

We have been working on some MMPB(GB)SA calculations using the new implementation of AmberTools14 and we are getting strange results.

We have used two options for decomposing the total energy into components: idecomp=2 (per residue) or idecomp=4 (pairwise). They should not affect the total energy in the calculation. However, the total energies differ depending whether we use idecomp=2 or idecomp=4. Using GB, the results are the same for both idecomp options, but using PB the result differ.

Please, find below an example of this calculation on exactly the same system (only using idecomp=2 or 4).


FINAL RESULTS - Per Residue Decomposition (idecomp=2)


   NSTEP ENERGY RMS GMAX NAME NUMBER

      1 -2.0693E+03 1.8571E+01 1.2681E+02 C 1456


 BOND = 348.5332 ANGLE = 949.0414 DIHED = 1280.5748

 VDWAALS = -899.3273 EEL = -7849.1436 EPB = -1458.7532

 1-4 VDW = 433.1028 1-4 EEL = 4767.9374 RESTRAINT = 0.0000

 ECAVITY = 954.5187 EDISPER = -595.7736

FINAL RESULTS - Pairwise Decomposition (idecomp=4)


   NSTEP ENERGY RMS GMAX NAME NUMBER

      1 -6.2420E+02 1.8571E+01 1.2681E+02 C 1456


 BOND = 348.5332 ANGLE = 949.0414 DIHED = 1280.5748

 VDWAALS = -899.3273 EEL = -7849.1436 EPB = -13.6645

 1-4 VDW = 433.1028 1-4 EEL = 4767.9374 RESTRAINT = 0.0000

 ECAVITY = 954.5187 EDISPER = -595.7736

As it can be seen above, all the values are the same with the exception of the EPB (in red).

Could you please help us to clarify this issue?

Thanks very much!

Best,

Yoel
Mount Sinai School of Medicine, NYC


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber

Received on Tue May 27 2014 - 10:30:02 PDT
Custom Search