Re: Fwd: AMBER: How to calculate the mmpbsa for explicit MD ?

From: Thomas Steinbrecher <steinbrt.scripps.edu>
Date: Mon, 12 Jun 2006 11:50:30 -0700 (PDT)

Hi Jitrayut

> Thank you for your kind answer but I am still confused about MM-PBSA script
> that you say strips waters and counterions from the trajectorie. Please
> explain in more specific details or how I can do with this script.

When I did MM-PBSA calculations, I always did two steps one after the
other:

For the processing of your trajectory set GC=1 and AS,DC,MM,GB,PB,MS=0 in
your input script and insert the details of your trajectory into the
.MAKECRD section (see amber 9 manual p273 for the meaning of the
parameters in MAKECRD) as well as the filename of your crd-file at
.TRAJECTORY

run mm_pbsa.pl [your input file] > mm_pbsa.log

This creates a lot of coordinate files in your directory.
Create a second input script in which you set GC=0 and e.g. MM=1, PB=1

run mm_pbsa.pl [your second input file]

and it should process your snapshots.

You have to set COMPT,RECPT and LIGPT to the filenames of your parameter
files without water molecules.

Regards,

Thomas
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jun 14 2006 - 06:07:11 PDT
Custom Search