[AMBER] MMPBSA ATP-Mg

From: Chiara Pallara <chiara.pallara.bsc.es>
Date: Wed, 13 Mar 2013 18:54:38 +0100

Hi,
/I have a protein/Mg-ATP structure and I'm trying to calculate //the
binding free energy of this protein-cofactor system.
//My input is:

&general
// strip_mask=:WAT:Cl-:Cs+:IB:K+:Li+:Na+:Rb+,
    startframe=1, endframe=5, interval=1,
    verbose=1,
/
&gb
   igb=5, saltcon=0.100,
/

The calculation run until the end but the results are really senseless.
The .pdb files for Ligand, Rec and Complex are like this
MODEL 1
ATOM 1 HH31 ACE 1 -nan -nan -nan 0.00 0.00
ATOM 2 CH3 ACE 1 -nan -nan -nan 0.00 0.00
ATOM 3 HH32 ACE 1 -nan -nan -nan 0.00 0.00
ATOM 4 HH33 ACE 1 -nan -nan -nan 0.00 0.00
ATOM 5 C ACE 1 -nan -nan -nan 0.00 0.00
ATOM 6 O ACE 1 -nan -nan -nan 0.00 0.00
ATOM 7 N GLU 2 -nan -nan -nan 0.00 0.00
ATOM 8 H GLU 2 -nan -nan -nan 0.00 0.00
ATOM 9 CA GLU 2 -nan -nan -nan 0.00 0.00
ATOM 10 HA GLU 2 -nan -nan -nan 0.00 0.00
ATOM 11 CB GLU 2 -nan -nan -nan 0.00 0.00
ATOM 12 HB2 GLU 2 -nan -nan -nan 0.00 0.00

And in the "FINAL_RESULTS_MMPBSA.dat", I found...

GENERALIZED BORN:

Complex:
Energy Component Average Std. Dev. Std. Err.
of Mean
-------------------------------------------------------------------------------
EGB -66427.3241 0.0000 0.0000

G gas 0.0000 0.0000 0.0000
G solv -66427.3241 0.0000 0.0000

TOTAL -66427.3241 0.0000 0.0000


Receptor:
Energy Component Average Std. Dev. Std. Err.
of Mean
-------------------------------------------------------------------------------
EGB -64081.8425 0.0010 0.0004

G gas 0.0000 0.0000 0.0000
G solv -64081.8425 0.0010 0.0004

TOTAL -64081.8425 0.0010 0.0004


Ligand:
Energy Component Average Std. Dev. Std. Err.
of Mean
-------------------------------------------------------------------------------
EGB -2345.4816 0.0000 0.0000

G gas 0.0000 0.0000 0.0000
G solv -2345.4816 0.0000 0.0000

TOTAL -2345.4816 0.0000 0.0000


Differences (Complex - Receptor - Ligand):
Energy Component Average Std. Dev. Std. Err.
of Mean
-------------------------------------------------------------------------------
VDWAALS 0.0000 0.0000 0.0000
EEL 0.0000 0.0000 0.0000
EGB 0.0000 0.0000 0.0000
ESURF 0.0000 0.0000 0.0000

DELTA G gas 0.0000 0.0000 0.0000
DELTA G solv 0.0000 0.0000 0.0000

DELTA TOTAL 0.0000 0.0000 0.0000


-------------------------------------------------------------------------------
-------------------------------------------------------------------------------

/
Can anyone tell me what i'm doing wrong? or some suggestion?
Thanks in advance,
/Chiara

/


WARNING / LEGAL TEXT: This message is intended only for the use of the
individual or entity to which it is addressed and may contain
information which is privileged, confidential, proprietary, or exempt
from disclosure under applicable law. If you are not the intended
recipient or the person responsible for delivering the message to the
intended recipient, you are strictly prohibited from disclosing,
distributing, copying, or in any way using this message. If you have
received this communication in error, please notify the sender and
destroy and delete any copies you may have received.

http://www.bsc.es/disclaimer
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Mar 13 2013 - 11:00:03 PDT
Custom Search