Hi John,
Try something like this:
-------------------------------------------
solvent byres :WAT
closest 100 : 1-81 oxygen noimage
trajout near100.pdb pdb
-------------------------------------------
May be it helps.
Good luck!
Igor
On Wed, Dec 8, 2010 at 7:44 AM, Beale, John <John.Beale.stlcop.edu> wrote:
>
> I am trying to determine the 100 closest water molecules to my protein.
> This is my ptraj input file:
>
> #!/bin/csh
> ptraj file1.prmtop <<EOF
> trajin file2.nc 1 1
> center :1-81 mass origin
> image origin center familiar
> solvent byres :WAT
> closest 100 : 1-81 first
> average near100.pdb pdb
> EOF
>
> When I run the script, the pdb file shows the water molecules clustered in
> a sphere in the nucleotide binding site of the protein and not distributed
> over the entire protein. The starting x-ray structure shows that water
> molecules should be distributed around the protein and not clustered in the
> binding site.
>
> Is there something wrong with my input file (above)?
>
> Thanks!
>
> John
>
>
> _______________________________________________
> 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 Mon Dec 13 2010 - 08:00:04 PST
