Dear professors,
Thanks very much for your reply. For convenience, I built a system including two amino acids and a MG2+ and then performed a short MD simulation. The protein and MG2+ were destined to the receptor and ligand, respectively. Then, I calculated MMPBSA using the 'MMPBSA.py' program and got the following errors:
[tony.tony py_mmpbsa]$ MMPBSA.py -O -i mmpbsa.in -o MMPBSA_result.dat -sp ../com_solvated.prmtop -cp ../com.prmtop -rp ../rec.prmtop -lp ../lig.prmtop -y ../prod1.mdcrd
ptraj found! Using /home/tony/program/amber11/exe/ptraj
sander found! Using /home/tony/program/amber11/exe/sander (serial only!)
Warning: igb=2 should be used with mbondi2 pbradii set. Yours are modified Bondi radii (mbondi)
Preparing trajectories with ptraj...
checkCoordinates(): Could not predict number of frames for AMBER trajectory file: _MMPBSA_complex.mdcrd
If this is not a compressed file then there is a problem
checkCoordinates(): Could not predict number of frames for AMBER trajectory file: _MMPBSA_complex.mdcrd
If this is not a compressed file then there is a problem
checkCoordinates(): Could not predict number of frames for AMBER trajectory file: _MMPBSA_complex.mdcrd
If this is not a compressed file then there is a problem
readAmberRestart(): topology/coordinates file are inconsist with
NATOMS = -2 (1)
Error! Ptraj failed. Check coordinate and topology files for the complex.
NOTE: All files have been retained for debugging purposes. Type MMPBSA.py --clean to erase these files.
I checked the file named "_MMPBSA_ptraj6.out" and found following informations:
\-/ Residue labels:
MG2
PTRAJ: Processing input from file _MMPBSA_ligandinpcrd.in
PTRAJ: trajin _MMPBSA_ligand.mdcrd 1 1 1
Checking coordinates: _MMPBSA_ligand.mdcrd
Could not process trajectory _MMPBSA_ligand.mdcrd
I also run the mm_pbsa.pl using the input file similar to the MMPBSA tutorial but the script interrupted, the snapshot_lig.all file showed as flowing:
MM
GB
PB
MS
PB_SURFTEN 0.0072
PB_SURFOFF 0.00
GB_SURFTEN 0.0072
GB_SURFOFF 0.00
1
BOND = 0.0000 ANGLE = 0.0000 DIHED = 0.0000
VDWAALS = 0.0000 EEL = 0.0000 EGB = -465.1086
1-4 VDW = 0.0000 1-4 EEL = 0.0000 RESTRAINT = 0.0000
corrected reaction field energy: -733.423395
surface area = 83.647
ECAVITY = 83.647
EDISPER = 0.0000
2
BOND = 0.0000 ANGLE = 0.0000 DIHED = 0.0000
VDWAALS = 0.0000 EEL = 0.0000 EGB = -465.1086
1-4 VDW = 0.0000 1-4 EEL = 0.0000 RESTRAINT = 0.0000
corrected reaction field energy: -733.423227
surface area = 83.647
ECAVITY = 83.647
EDISPER = 0.0000
and the snapshot_statistics.out file showed as flowing:
# COMPLEX RECEPTOR LIGAND
# ----------------------- ----------------------- -----------------------
# MEAN STD MEAN STD MEAN STD
# ======================= ======================= ======================
ELE -137.41 163.23 0.00 0.00 -76.08 56.90
VDW 50986.49 82978.20 0.00 0.00 50985.07 82979.48
INT 280795.09 191390.71 0.00 0.00 280795.09 191390.71
GAS 331644.18 249493.62 0.00 0.00 331704.08 249428.77
PBSUR 6.77 2.50 0.60 0.00 6.64 2.38
PBCAL -652.88 232.67 -437.21 0.00 -299.30 127.24
PBSOL -646.11 230.94 -436.61 0.00 -292.66 125.52
PBELE -790.29 77.59 -437.21 0.00 -375.38 77.39
PBTOT 330998.07 249337.32 -436.61 0.00 331411.42 249344.52
GBSUR 6.77 2.50 0.60 0.00 6.64 2.38
GB -734.84 266.96 -465.11 0.00 -345.33 150.53
GBSOL -728.07 265.16 -464.51 0.00 -338.69 148.65
GBELE -872.25 113.65 -465.11 0.00 -421.41 104.22
GBTOT 330916.11 249304.55 -464.51 0.00 331365.39 249317.52
# DELTA
# -----------------------
# MEAN STD
# =======================
ELE -61.33 118.71
VDW 1.42 3.26
INT 0.00 0.00
GAS -59.90 115.91
PBSUR -0.47 0.14
PBCAL 83.64 115.53
PBSOL 83.17 115.47
PBELE 22.31 8.87
PBTOT 23.26 9.42
GBSUR -0.47 0.14
GB 75.60 124.91
GBSOL 75.13 124.86
GBELE 14.27 10.33
GBTOT 15.23 13.23
thank you very much!
------------------ ÔʼÓʼþ ------------------
·¢¼þÈË: "Jason Swails"<jason.swails.gmail.com>;
·¢ËÍʱ¼ä: 2011Äê1ÔÂ28ÈÕ(ÐÇÆÚÎå) ÍíÉÏ9:55
ÊÕ¼þÈË: "AMBER Mailing List"<amber.ambermd.org>;
Ö÷Ìâ: Re: [AMBER] MM-PBSA error about protein + Mg2+
http://archive.ambermd.org/201101/0449.html
On Fri, Jan 28, 2011 at 2:47 AM, ФÕýÌÎ <xzt41.126.com> wrote:
> Dear amber professors,
> I am doing the free energy calculation between the protein and Mg2+
> using MM-PBSA with Amber10(both mm_pbsa.pl and MMPBSA.py),
> but many trials failed, when i check the temporary files, the following
> results were found in the "snapshot_lig.all.out":
> MM
> GB
> PB
> MS
> 1
> BOND = 0.0000 ANGLE = 0.0000 DIHED =
> 0.0000
> VDWAALS = 0.0000 EEL = 0.0000 EGB =
> -465.1086
> 1-4 VDW = 0.0000 1-4 EEL = 0.0000 RESTRAINT =
> 0.0000
> corrected reaction field energy: -733.413162
> surface area = 83.647
> ECAVITY = 83.647
> EDISPER = 0.0000
>
>
>
>
> I am confused with this, then i calculated the MM-PBSA energy of protein
> and the Na+, the same results was produced. Little information about this
> can be found by google, so i hope you can help me and give me some advices.
> thank you and best wishes to you!!
>
>
> Xiao Zhengtao ,
>
> Department of bioinformatics in Northwest Sci-Tech
> University of Agriculture and Forestry
>
>
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
Jason M. Swails
Quantum Theory Project,
University of Florida
Ph.D. Graduate Student
352-392-4032
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Jan 30 2011 - 01:30:03 PST