[AMBER] error while running mmpbsa.py in amber12

From: Alessandro Contini <alessandro.contini.unimi.it>
Date: Wed, 30 Jan 2013 12:47:41 +0100

Dear community,
I've recently installed amber12 (composer XE 2013.1.117) and I'm
experiencing an error while running mmpbsa.py on a system where I've
previously run the very same calculation with amber11.

here follows the input:

GB/PBSA analysis
&general
   startframe=1, endframe=1000, interval=20, verbose=2, keep_files=0,
   strip_mask=:WAT,MG2,Na+,GTP,GDP,:1-438
/
&gb
  igb=2, saltcon=0.150,
/
&pb
  istrng=0.15, fillratio=4.0
/


I've launched mmpbsa.py trough the following command:

/usr/local/amber12/bin/MMPBSA.py -O -i mmpbsa_py.in -o MMPBSA_PB-GB.txt
-sp solvated_complex.top -cp complex.top -rp rec.top -lp lig.top -y
traj.mdcrd.bz2 >& mmpbsa_py.out

and this is the output:

Reading command-line arguments and input files...
Loading and checking parameter files for compatibility...
mmpbsa_py_energy found! Using /usr/local/amber12/bin/mmpbsa_py_energy
cpptraj found! Using /usr/local/amber12/bin/cpptraj
Preparing trajectories for simulation...
ValueError: invalid literal for int() with base 10: 'an'
Exiting. All files have been retained.

I did a quick google search and I understood that it probably is a
python error, although I'm not experienced with python programming.

I also found something on amber mailing list,
http://archive.ambermd.org/201106/0285.html and related messages, but my
calculations are run on standard topologies generated by tleap found in
previous amber versions. Hereafter follows the head of each top

solvated_complex.top
%VERSION VERSION_STAMP = V0001.000 DATE = 06/21/12
14:17:25
%FLAG
TITLE
%
FORMAT(20a4)
                                                                                
%FLAG
POINTERS
%FORMAT(10I8)
  254878 21 240786 14374 29476 19866 58636 45348 0
0
  452886 77696 14374 19866 45348 86 194 85 55
1
       0 0 0 0 0 0 0 1 119
0
       0

complex.top
%VERSION VERSION_STAMP = V0001.000 DATE = 03/04/11
10:09:56
%FLAG
TITLE
%
FORMAT(20a4)
                                                                                
%FLAG
POINTERS
%
FORMAT(10I8)
   13454 15 6623 6989 15021 9497 29072 21848 0
0
   74154 865 6989 9497 21848 67 156 67 44
0
       0 0 0 0 0 0 0 0 119
0
       0

rec.top
%VERSION VERSION_STAMP = V0001.000 DATE = 03/02/11
11:56:55
%FLAG
TITLE
%
FORMAT(20a4)
                                                                                
%FLAG
POINTERS
%
FORMAT(10I8)
    6742 14 3325 3495 7550 4746 14555 10953 0
0
   37194 438 3495 4746 10953 43 99 50 29
0
       0 0 0 0 0 0 0 0 26
0
       0
lig.top
%VERSION VERSION_STAMP = V0001.000 DATE = 03/04/11
10:09:56
%FLAG
TITLE
%
FORMAT(20a4)
                                                                                
%FLAG
POINTERS
%
FORMAT(10I8)
    6712 15 3298 3494 7471 4751 14517 10895 0
0
   36960 427 3494 4751 10895 67 156 67 44
0
       0 0 0 0 0 0 0 0 119
0
       0


All test run (serial) from the amber test suite run successfully.

Any suggestion will be greatly appreciated.

Best regards,

Alessandro Contini

 

-- 
Alessandro Contini, PhD
Dipartimento di Scienze Farmaceutiche
Sezione di Chimica Generale e Organica "A. Marchesini"
Via Venezian, 21 20133 Milano
tel. +390250314480
e-mail alessandro.contini.unimi.it
skype alessandrocontini
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jan 30 2013 - 04:00:02 PST
Custom Search