[AMBER] ptraj closest water

From: George Tzotzos <gtzotzos.me.com>
Date: Mon, 18 Jul 2011 18:00:59 +0200

Hi everybody,

The sequence of actions below is self explanatory:

1. strip all but closest water molecules to a given protein residue (:124) and generate a pdb file (new.pdb)

2. process a 2ns trajectory file by keeping only the closest water molecules and generate a new trajectory file (new.mdcrd)

3. create new topology files from new.pdb using tleap

The above scheme works with the following hiccup.

In step 1., the script

trajin prod_2ns.mdcrd
center :1-124 mass origin
image origin center familiar
solvent byres :WAT
closest 4 :124.N first no image
strip :Na+
trajout test.pdb pdb 1 200 1

produces 200 individual pdb files.

My questions:

Is there a way to produce only one new.pdb file? Furthermore, if the working directory contains other trajectories (e.g. prod_4ns, prod_6ns, etc.), the script scans them all - although they are not specified in the first line of the script. Any suggestions why so?

Is there a way to merge the script of step 1 with that of step 2, below?

trajin prod_2ns.mdcrd
trajout prod_closest_WAT_2.mdcrd netcdf
center :1-142 mass origin
image origin center familiar
solvent byres :WAT
closest 4 :124.N first no image
strip :Na+





_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jul 18 2011 - 09:30:02 PDT
Custom Search