Hi,
Short answer is unfortunately no. The 'watershell' command (or really
any other command in cpptraj) doesn't have a concept of "inside" a
protein. It works by counting waters that are within certain distance
cutoffs from each atom. I suppose if your protein were roughly
spherical you could pick 1 atom or residue near the center of the
protein, choose a short first cutoff, then a second cutoff that's
roughly the radius of your protein, but that seems pretty problematic.
Would be interesting in the future to implement a command with concept
of "inside" a protein, but I don't see myself getting to it anytime
soon unfortunately...
-Dan
On Wed, May 6, 2020 at 9:18 AM 石娜 <shina14.mails.jlu.edu.cn> wrote:
>
>
> Dear amber developers:
>
> In the process of analyse the trajectory, I need to calculate the number of
> water molecules inside the protein in each snapshot. Can watershell achieve
> this goal? If not, can you recommend a better method. I am looking forward
> to your reply! Thanks very much!
> Yours sincerely!
>
> Na Shi
> _______________________________________________
> 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 Wed May 20 2020 - 09:30:04 PDT
