Re: [AMBER] ptraj usage for clustering of protein ensembles manual or tutorial

From: Andrew Voronkov <drugdesign.yandex.ru>
Date: Fri, 04 Dec 2009 22:09:58 +0300

I actually want to divide trajectory of 5000 snapshots from 5 nanoseconds into 10-20 representative structures as it is done in relaxed complex scheme, but in my case i dont have a ligand, just want to simulate the protein flexibility by usage of most common structures with rigid docking. So what procedure is best for such a purpose and is there anywhere more detailed description of these procedures?

Sincerely yours,
Andrey

04.12.09, 10:01, "Thomas Cheatham III" <tec3.utah.edu>:

>
> > trajin fzd2_w_md5.mdcrd
> > cluster out avg-1 representative none average none all none averagelinkage clusters 1 mass rms
> > "(:1-240) & .CA,C,N,O"
> > go
>
> The "none" keyword means you are choosing not to dump any coordinate
> output for the "representative" snapshot (that closest to the centroid),
> the "average" over each cluster, and dumping of "all" the frames in each
> cluster...
>
> Also, I do not not understand why you would want to cluster into a single
> (cluster 1) cluster as this will simply give back your original
> trajectory?
>
>
>
> cluster out c-out representative pdb average pdb averagelinkage \
> clusters 10 mass rms :1-240.CA,C,N,O
>
> will create
>
> c-out.avg.n <- average structures for each cluster as PDB
> c-out.rep.n <- representative " "
> c-out.n <- trajectories in mdcrd format for each cluster
>
>
> -- tec3
>
>
> _______________________________________________
> 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 Fri Dec 04 2009 - 11:30:02 PST
Custom Search