AMBER Archive (2006)Subject: Re: AMBER: GBSA EGB = NaN
From: William Wei (william.wei_at_utoronto.ca)
Date: Mon Aug 21 2006 - 11:04:36 CDT
Dear Ross,
I tried to run the mdin file generated by mm-pbsa, and did not find any
abnormal output. That give the same EGB=NaN. I generted the structure from
my .crd and .top files, I checked using sybyl, I did not see any improper
contact or bond. The VDWAALS and BOND items are okay. (I will attach the
output from the seperate run).
For the IGB, I used 4. The amberversion I am using is amber 7. Is this
version too old? I do not know if this version have option IGB=5 or not.
Thank you again for your help,
William
###Complex sander input file:
File generated by mm_pbsa.pl. Using MM GB
&cntrl
ntf = 1, ntb = 0, dielc = 1,
idecomp= 0,
igb = 4, saltcon= 0.00,
offset = 0, extdiel= 80.0,
gbsa = 0, surften= 1.0,
cut = 99999.0, nsnb = 99999,
scnb = 2.0, scee = 1.2,
imin = 1, maxcyc = 1, ncyc = 0,
&end
### sander output file:
-------------------------------------------------------
Amber 7 SANDER Scripps/UCSF 2002
-------------------------------------------------------
| Mon Aug 21 11:41:16 2006
[-O]verwriting output
File Assignments:
| MDIN: ./sander_com.in
| MDOUT: ./sander_com_test.out
|INPCRD: ../sn_test_20/ODC_OMP_com.crd
| PARM: ../ODCase_Com.top
|RESTRT: restrt
| REFC: refc
| MDVEL: mdvel
| MDEN: mden
| MDCRD: mdcrd
|MDINFO: mdinfo
|INPDIP: inpdip
|RSTDIP: rstdip
Here is the input file:
File generated by mm_pbsa.pl. Using MM GB
&cntrl
ntf = 1, ntb = 0, dielc = 1,
idecomp= 0,
igb = 4, saltcon= 0.00,
offset = 0, extdiel= 80.0,
gbsa = 0, surften= 1.0,
cut = 99999.0, nsnb = 99999,
scnb = 2.0, scee = 1.2,
imin = 1, maxcyc = 1, ncyc = 0,
&end
----------------------------------------------------------------------------
----
1. RESOURCE USE:
----------------------------------------------------------------------------
----
| Flags: SGIFFT MEM_ALLOC MPI RLE ROWAT HAS_10_12 HAS_FTN_ERFC
| New format PARM file being parsed.
| Version = 1.000 Date = 08/01/06 Time = 15:42:24
NATOM = 6430 NTYPES = 18 NBONH = 3216 MBONA = 3262
NTHETH = 7318 MTHETA = 4412 NPHIH = 13512 MPHIA = 8232
NHPARM = 0 NPARM = 0 NNB = 35364 NRES = 420
NBONA = 3262 NTHETA = 4412 NPHIA = 8232 NUMBND = 44
NUMANG = 94 NPTRA = 47 NATYP = 30 NPHB = 0
IFBOX = 0 NMXRS = 34 IFCAP = 0 NEXTRA = 0
| Memory Use Allocated
| Real 257515
| Hollerith 39002
| Integer 384888
| Max Pairs 1
| Max Rstack 96450
| Max Istack 32150
| Duplicated 0 dihedrals
| Duplicated 0 dihedrals
----------------------------------------------------------------------------
----
2. CONTROL DATA FOR THE RUN
----------------------------------------------------------------------------
----
General flags:
imin = 1, nmropt = 0
Nature and format of input:
ntx = 1, irest = 0, ntrx = 1
Nature and format of output:
ntxo = 1, ntpr = 50, ntrx = 1, ntwr =
500
iwrap = 0, ntwx = 0, ntwv = 0, ntwe =
0
ioutfm = 0, ntwprt = 0, idecomp = 0, rbornstat=
0
Potential function:
ntf = 1, ntb = 0, igb = 4, nsnb =
99999
ipol = 0, gbsa = 0
dielc = 1.00000, cut =**********, intdiel = 1.00000
saltcon = 0.00000, offset = 0.00000, gbalpha= 1.00000
gbbeta = 0.00000, gbgamma = 0.00000, surften = 1.00000
scnb = 2.00000, scee = 1.20000
Frozen or restrained atoms:
ibelly = 0, ntr = 0
Energy minimization:
maxcyc = 1, ncyc = 0, ntmin = 1
dx0 = 0.01000, dxm = 0.50000, drms = 0.00010
| PLEVEL = 1: runmd parallelization, no EKCMR
----------------------------------------------------------------------------
----
3. ATOMIC COORDINATES AND VELOCITIES
----------------------------------------------------------------------------
----
begin time read from input coords = 0.000 ps
| Reorder water Setup:
| ----- Reordering water will NOT be done
| Does not work with minimization
|
Number of triangulated 3-point waters found: 0
| Atom division among processors:
| 0 6430
| Running AMBER/MPI version on 1 nodes
----------------------------------------------------------------------------
----
4. RESULTS
----------------------------------------------------------------------------
----
NSTEP ENERGY RMS GMAX NAME NUMBER
1 NaN NaN 0.0000E+00 N 1
BOND = 1198.5268 ANGLE = 3183.2750 DIHED =
3858.7614
VDWAALS = -3497.1212 EEL = -31925.1169 EGB =
NaN
1-4 VDW = 1484.4497 1-4 EEL = 21392.8129 RESTRAINT =
0.0000
MAXIMUM NUMBER OF F EVALUATION EXCEEDED
FINAL RESULTS
NSTEP ENERGY RMS GMAX NAME NUMBER
1 NaN NaN 0.0000E+00 N 1
BOND = 1198.5268 ANGLE = 3183.2750 DIHED =
3858.7614
VDWAALS = -3497.1212 EEL = -31925.1169 EGB =
NaN
1-4 VDW = 1484.4497 1-4 EEL = 21392.8129 RESTRAINT =
0.0000
----------------------------------------------------------------------------
----
5. TIMINGS
----------------------------------------------------------------------------
----
| Read coords time 0.04 ( 0.07% of Total)
| Calc gb radii 9.40 (19.04% of Gen B)
| Calc gb diag 18.47 (37.43% of Gen B)
| Calc gb off-diag 21.47 (43.52% of Gen B)
| Gen Born time 49.34 (100.0% of Nonbo)
| Nonbond force 49.34 (99.86% of Force)
| Angle energy 0.02 ( 0.05% of Force)
| Dihedral energy 0.04 ( 0.09% of Force)
| Force time 49.41 (100.0% of Runmd)
| Runmd Time 49.41 (98.35% of Total)
| Other 0.79 ( 1.58% of Total)
| Total time 50.24 (100.0% of ALL )
| Highest rstack allocated: 96450
| Highest istack allocated: 12860
| Setup wallclock 1 seconds
| Nonsetup wallclock 49 seconds
##############
----- Original Message -----
From: "Ross Walker" <ross_at_rosswalker.co.uk>
To: <amber_at_scripps.edu>
Sent: Friday, August 18, 2006 6:11 PM
Subject: RE: AMBER: GBSA EGB = NaN
> Dear William,
>
> I was more interested in seeing the mdin file greated by mmpbsa. Can you
run
> the calculation manually outside of the mm-pbsa script. What happens in
this
> case?
>
> There is nothing that looks strange in your output files so it is unlikely
> that you have a bad structure. Do the tests cases all pass? Have you tried
> setting up a simple GB calculation of your system independently of mmpbsa
> and verified that this give the same problem?
>
> I would also try recreating the prmtop file and then also taking a look at
> the section labelled:
>
> %FLAG RADII
> %FORMAT(5E16.8)
>
> Does anything look weird here. Most of these numbers should be around 0.1
to
> 2.0. I would pay close attention to anything that is not.
>
> Also which version of GB are you using? I know there can be some issues
with
> IGB=1 or 2 although I forget specifically what the issues are - perhaps
> someone else on the list can contribute here. If you are using IGB=1 or 2
> try using IGB=5 and see if this fixes the problem.
>
> All the best
> Ross
>
> /\
> \/
> |\oss Walker
>
> | HPC Consultant and Staff Scientist |
> | San Diego Supercomputer Center |
> | Tel: +1 858 822 0854 | EMail:- ross_at_rosswalker.co.uk |
> | http://www.rosswalker.co.uk | PGP Key available on request |
>
> Note: Electronic Mail is not secure, has no guarantee of delivery, may not
> be read every day, and should not be used for urgent or sensitive issues.
>
> > -----Original Message-----
> > From: owner-amber_at_scripps.edu
> > [mailto:owner-amber_at_scripps.edu] On Behalf Of William Wei
> > Sent: Friday, August 18, 2006 08:12
> > To: amber_at_scripps.edu
> > Subject: Re: AMBER: GBSA EGB = NaN
> >
> > Dear Ross,
> >
> > At the beginning of my calculation, it was running sander with the
> > snapshots. I did not see anything unusual. My system is not
> > big, as about
> > 6430atoms in my complex. I am using SGI Onyx 3800 with IRIX
> > 6.5 system.
> > Here is my input file and the output of my complex.
> >
> > Thank you very much,
> > William
> > GBSA input:
> > @GENERAL
> > #
> > PREFIX ODC_OMP
> > PATH ../sn_test_20/
> > #
> > COMPLEX 1
> > RECEPTOR 1
> > LIGAND 1
> >
> > #
> > #
> >
> > COMPT ../ODCase_Com.top
> >
> > RECPT ../ODCase_Rec.top
> >
> > LIGPT ../ODCase_Lig.top
> >
> > #
> >
> > GC 0
> >
> > AS 0
> >
> > DC 0
> >
> > #
> >
> > MM 1
> >
> > GB 1
> >
> > PB 0
> >
> > MS 1
> >
> > #
> >
> > NM 0
> >
> > #
> >
> > ##############################################################
> > ##############
> > ####
> >
> > @DELPHI
> >
> > #
> >
> > #
> >
> > FOCUS 0
> >
> > INDI 1.0
> >
> > EXDI 80.0
> >
> > PERFIL 80.0
> >
> > SCALE 2
> >
> > LINIT 1000
> >
> > BNDCON 4
> >
> > #CHARGE ../my_amber94_delphi.crg
> >
> > #SIZE ../my_parse_delphi.siz
> >
> > CHARGE ../amber94_delphi_OMP.crg
> >
> > SIZE ../parseres_OMP.siz
> >
> > #
> >
> > SURFTEN 0.00542
> >
> > SURFOFF 0.92
> >
> > #
> >
> > ##############################################################
> > ##############
> > ####
> >
> > @GB
> >
> > #
> >
> > # GB parameters
> >
> > #
> >
> > IGB 4
> >
> > SALTCON 0.00
> >
> > EXTDIEL 80.0
> >
> > #
> >
> > SURFTEN 0.0072
> >
> > SURFOFF 0.00
> >
> > #
> > ##############################################################
> > ##############
> > ####
> >
> > @MS
> > #
> > # Molsurf parameters
> > #
> > # PROBE - Radius of the probe sphere used to calculate the SAS.
> >
> > # RADII - Name of the radii file.
> > #
> > PROBE 1.4
> >
> > #RADII ../atmtypenumbers
> >
> > RADII ../atmtypenumbers_OMP
> > #
> > ##############################################################
> > ##############
> > #####
> >
> > @PROGRAMS
> > #
> > # Program executables and paths
> > #
> > #DELPHI /home/gohlke/src/delphi.98/exe.R10000/delphi
> > DELPHI /aspirin/william/Software/delphi/IRIX/delphi
> > #
> > ##############################################################
> > ##############
> > ####
> >
> >
> >
> > Output of my complex:
> > MM
> > GB
> > MS
> > 1
> > BOND = 1198.5268 ANGLE = 3183.2750 DIHED =
> > 3858.7614
> > VDWAALS = -3497.1212 EEL = -31925.1169 EGB =
> > NaN
> > 1-4 VDW = 1484.4497 1-4 EEL = 21392.8129 RESTRAINT =
> > 0.0000
> > surface area = 15877.709
> > 2
> > BOND = 1249.4284 ANGLE = 3112.2372 DIHED =
> > 3852.2178
> > VDWAALS = -3529.7154 EEL = -31904.5266 EGB =
> > NaN
> > 1-4 VDW = 1476.8527 1-4 EEL = 21318.8379 RESTRAINT =
> > 0.0000
> > surface area = 15876.309
> > 3
> > BOND = 1206.7192 ANGLE = 3145.3704 DIHED =
> > 3880.0414
> > VDWAALS = -3512.8353 EEL = -31834.1365 EGB =
> > NaN
> > 1-4 VDW = 1508.3817 1-4 EEL = 21321.6349 RESTRAINT =
> > 0.0000
> > surface area = 15878.668
> > 4
> > BOND = 1236.1746 ANGLE = 3158.5354 DIHED =
> > 3869.2900
> > VDWAALS = -3485.6651 EEL = -31855.4143 EGB =
> > NaN
> > 1-4 VDW = 1461.7513 1-4 EEL = 21376.0072 RESTRAINT =
> > 0.0000
> > surface area = 15798.397
> > 5
> > BOND = 1247.7303 ANGLE = 3273.7533 DIHED =
> > 3876.9796
> > VDWAALS = -3534.9529 EEL = -31871.1197 EGB =
> > NaN
> > 1-4 VDW = 1458.7460 1-4 EEL = 21307.7514 RESTRAINT =
> > 0.0000
> > surface area = 15895.754
> > 6
> > BOND = 1241.4778 ANGLE = 3222.0305 DIHED =
> > 3925.6273
> > VDWAALS = -3488.7917 EEL = -31868.4341 EGB =
> > NaN
> > 1-4 VDW = 1480.5582 1-4 EEL = 21310.2085 RESTRAINT =
> > 0.0000
> > surface area = 15913.447
> > 7
> > BOND = 1246.3852 ANGLE = 3210.2564 DIHED =
> > 3940.9817
> > VDWAALS = -3559.4844 EEL = -31873.9004 EGB =
> > NaN
> > 1-4 VDW = 1471.1604 1-4 EEL = 21342.0148 RESTRAINT =
> > 0.0000
> > surface area = 15736.376
> > 8
> > BOND = 1239.2338 ANGLE = 3248.2204 DIHED =
> > 3892.2518
> > VDWAALS = -3497.6510 EEL = -31832.1607 EGB =
> > NaN
> > 1-4 VDW = 1449.2936 1-4 EEL = 21270.9055 RESTRAINT =
> > 0.0000
> > surface area = 15745.638
> >
> >
> > ----- Original Message -----
> > From: "Ross Walker" <ross_at_rosswalker.co.uk>
> > To: <amber_at_scripps.edu>
> > Sent: Thursday, August 17, 2006 10:57 PM
> > Subject: RE: AMBER: GBSA EGB = NaN
> >
> >
> > > Dear William,
> > >
> > > Can you post the input files you used both for
> > equilibration and the GB
> > part
> > > of the run. Also include the first few steps of the output
> > file so we can
> > > look more carefully. What happens if you run the
> > calculation from the
> > > beginning with GB turned on? E.g. during the initial
> > minimisation - do you
> > > get a problem here? Also do all the GB test cases work? And
> > how big is
> > your
> > > molecule? What is you hardware (compilers, math library
> > etc). I have come
> > > across bugs in IBM's MASSV library on Blue Gene systems
> > where you would
> > get
> > > NaN for GB simulations on more than about 20,000 atoms
> > (vexp on vector
> > > lengths of > 10^8 gives NaN's) but this is a very specific
> > case. There
> > could
> > > be other systems that have bugs in their vector math
> > libraries though. But
> > > first of all before we jump to conclusions it would be good
> > to see all
> > your
> > > input and output files.
> > >
> > > All the best
> > > Ross
> > >
> > > /\
> > > \/
> > > |\oss Walker
> > >
> > > | HPC Consultant and Staff Scientist |
> > > | San Diego Supercomputer Center |
> > > | Tel: +1 858 822 0854 | EMail:- ross_at_rosswalker.co.uk |
> > > | http://www.rosswalker.co.uk | PGP Key available on request |
> > >
> > > Note: Electronic Mail is not secure, has no guarantee of
> > delivery, may not
> > > be read every day, and should not be used for urgent or
> > sensitive issues.
> > >
> > > > -----Original Message-----
> > > > From: owner-amber_at_scripps.edu
> > > > [mailto:owner-amber_at_scripps.edu] On Behalf Of William Lianhu Wei
> > > > Sent: Thursday, August 17, 2006 19:31
> > > > To: amber_at_scripps.edu
> > > > Subject: AMBER: GBSA EGB = NaN
> > > >
> > > > Dear all,
> > > >
> > > > I finished molecular dynamics simulation of my system about
> > > > 3ns. From the
> > > > energy profile, the system reached equilibrium at about 1.5
> > > > ns. I was trying
> > > > to calculate the binding free energy. The trajectories I
> > > > used was after
> > > > 1.5ns. No matter which trajectory I used, I always got
> > > > EGB=NaN (on Complex
> > > > and receptor), and the program stopped at statistical
> > > > calculation. The VDW
> > > > contribution is okay. Does anyone meet this before, or give
> > > > me any clue for
> > > > that. In 2004 archive, there is somebody asked this, but
> > > > there was no reply
> > > > for this. Thanks in advance.
> > > >
> > > > William
> > > >
> > > > --
> > > > --------------------------------------------------------------
> > > > William Wei Tel: +1-416-946-7551
> > > > Faculty of Pharmacy Email: william.wei_at_utoronto.ca
> > > > University of Toronto weilianhu_at_hotmail.com
> > > > 19 Russell St.
> > > > Toronto, Ontario
> > > > M5S 2S2,Canada
> > > > --------------------------------------------------------------
> > > > ---------
> > > > The AMBER Mail Reflector
> > > > To post, send mail to amber_at_scripps.edu
> > > > To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
> > > >
> > >
> > >
> > >
> > --------------------------------------------------------------
> > ---------
> > > The AMBER Mail Reflector
> > > To post, send mail to amber_at_scripps.edu
> > > To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
> > >
> >
> > --------------------------------------------------------------
> > ---------
> > The AMBER Mail Reflector
> > To post, send mail to amber_at_scripps.edu
> > To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
> >
>
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|