Dear All,
I am using AMBER10... I am using MM PBSA for calculating Binding free energy of a ligand...
i ran 1 ns production MD..
i used the following commands:
mm_pbsa.pl extract_coords.mmpbsa > extract_coords.log and
mm_pbsa.pl binding_energy.mmpbsa > binding_energy.log
After i run mm_pbsa.pl binding_energy.mmpbsa > binding_energy.log i get the following error
/opt/amber10/exe/sander -O
-i pbs a_lig.in -o pbsa_lig.8.out -c
/home/bharatam/amberresults1/vikamb_q1results/snap shot_lig.crd.8 -p
/home/bharatam/amberresults1/vikamb_q1/mol.prmtop not successful
mol is the name of my ligand...
How can i overcome this error?????
Please help...
the input files(extract_coords.mmpbsa and binding_energy.mmpbsa) are given below..
extract_coords.mmpbsa:
.GENERAL
PREFIX snapshot
PATH
/home/bharatam/amberresults1/vikamb_q1results
#
COMPLEX 1
RECEPTOR 1
LIGAND 1
#
COMPT
/home/bharatam/amberresults1/vikamb_q1/Awr.prmtop
RECPT
/home/bharatam/amberresults1/vikamb_q1/A.prmtop
LIGPT
/home/bharatam/amberresults1/vikamb_q1/wr.prmtop
#
GC 1
AS 0
DC 0
#
MM 0
GB 0
PB 0
MS 0
#
NM 0
.MAKECRD
BOX YES
NTOTAL 31303
NSTART 1
NSTOP 100
NFREQ 1
#
NUMBER_LIG_GROUPS 1
LSTART 3716
LSTOP 3757
NUMBER_REC_GROUPS 1
RSTART 1
RSTOP 3715
.TRAJECTORY
TRAJECTORY /home/bharatam/amberresults1/vikamb_q1/prod1.mdcrd
TRAJECTORY /home/bharatam/amberresults1/vikamb_q1/prod2.mdcrd
binding_energy.mmpbsa
.GENERAL
PREFIX snapshot
PATH
/home/bharatam/amberresults1/vikamb_q1results
#
COMPLEX 1
RECEPTOR 1
LIGAND 1
#
COMPT
/home/bharatam/amberresults1/vikamb_q1/Amol.prmtop
RECPT
/home/bharatam/amberresults1/vikamb_q1/A.prmtop
LIGPT
/home/bharatam/amberresults1/vikamb_q1/mol.prmtop
#
GC 0
AS 0
DC 0
#
MM 1
GB 1
PB 1
MS 1
#
NM 0
.PB
PROC 2
REFE 0
INDI 1.0
EXDI 80.0
SCALE 2
LINIT 1000
PRBRAD 1.4
ISTRNG 0.0
RADIOPT 0
NPOPT 1
CAVITY_SURFTEN 0.0072
CAVITY_OFFSET 0.00
#
SURFTEN 0.0072
SURFOFF 0.00
.MM
DIELC 1.0
.GB
IGB 2
GBSA 1
SALTCON 0.00
EXTDIEL 80.0
INTDIEL 1.0
#
SURFTEN 0.0072
SURFOFF 0.00
.MS
PROBE 0.0
.PROGRAMS
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed May 20 2009 - 17:19:06 PDT
