AMBER Archive (2005)

Subject: Re: AMBER: Can I get each corresponding sanpshots' delta delta energy term from mm_pbsa

From: JunJun Liu (ljjlp03_at_gmail.com)
Date: Mon Sep 19 2005 - 21:51:32 CDT


Hi Zhang:

you can extract energy terms for each snapshots from prefix_com.all.out(or
rec.all.out, lig.all.out) into a seperated file. For example:

1). The content of prefix_com.all.out is:
=====prefix_com.all.out========
MM
PB
MS
1
  BOND = 207.2842 ANGLE = 886.4916 DIHED =
2458.2773
  VDWAALS = -2609.7641 EEL = -25189.6876 HBOND =
0.0000
  1-4 VDW = 1073.6154 1-4 EEL = 12345.0770 RESTRAINT =
0.0000
corrected reaction field energy: -4923.742000
surface area = 14783.808
2
  BOND = 205.2730 ANGLE = 882.7424 DIHED =
2452.5204
  VDWAALS = -2625.4418 EEL = -25127.6069 HBOND =
0.0000
  1-4 VDW = 1072.0818 1-4 EEL = 12354.1791 RESTRAINT =
0.0000
corrected reaction field energy: -5061.096000
surface area = 14620.204
=============================

2). extracted seperated files are as follows:
======prefix_com.1.out=======
MM
PB
MS
1
  BOND = 207.2842 ANGLE = 886.4916 DIHED =
2458.2773
  VDWAALS = -2609.7641 EEL = -25189.6876 HBOND =
0.0000
  1-4 VDW = 1073.6154 1-4 EEL = 12345.0770 RESTRAINT =
0.0000
corrected reaction field energy: -4923.742000
surface area = 14783.808
==============================

=====prefix_com.2.out========
MM
PB
MS
1
  BOND = 205.2730 ANGLE = 882.7424 DIHED =
2452.5204
  VDWAALS = -2625.4418 EEL = -25127.6069 HBOND =
0.0000
  1-4 VDW = 1072.0818 1-4 EEL = 12354.1791 RESTRAINT =
0.0000
corrected reaction field energy: -5061.096000
surface area = 14620.204
=============================

3). doing similar operation for prefix_rec.all.out and prefix_lig.all.out

4). write input file for mm_pbsa_statistics.pl
===prefix_statistics.1.in===
prefix_com.1.out
prefix_rec.1.out
prefix_lig.1.out
============================

5). get statistical results for each snapshot
mm_pbsa_statistics.pl 1 0 prefix_statistics.1.in prefix_statistics.1.out
mm_pbsa_statistics.pl 1 0 prefix_statistics.2.in prefix_statistics.2.out
....
...

You can write a simple shell script for this job.

Best Regards!

Liu

On Tue, 13 Sep 2005 06:43:13 +0800, Suxin Zheng <sxzheng_at_u.washington.edu>
wrote:

> Hi,
> Can I output each snapshots' energy term from mm_pbsa, such as PBTOTAL,
> vdw,elec term of delta delta G = com_x - rec_x - lig_x, where x
> represents the snaphosts' number.
> Thanks
> suxin
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu

-- 
JunJun Liu

College of Chemistry Central China Normal University WuHan 430079 P.R. China ----------------------------------------------------------------------- The AMBER Mail Reflector To post, send mail to amber_at_scripps.edu To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu