On 9/01/2012 4:53 PM, John Ladasky wrote:
Hello everyone,

I've used the rot+trans option in GROMACS trjconv to superimpose groups of atoms within a single molecular dynamics simulation. I am now interested in modeling a protein, and a rather thoroughly scrambled circular permutation of that same protein. I want to construct a superimposition which optimizes the alignment of selected residues within the two structures. But unlike in a standard trjconv superimposition, I will be trying to align atoms in one structure which correspond to atoms with different numbers in the other structure. So the idea of groups as numbered atoms, as conventionally defined in a GROMACS .ndx file, would not seem to apply here.

Is this a practical task within GROMACS, or should I be looking elsewhere for the tool to accomplish the job?


Doesn't trjconv -fit rot+trans prompt for a group for each structure for the superposition? If so, you should be able to construct an index file with groups like

[ normal ]

3 45 223

[ scrambled ]

451 97 345

and fit those.

Mark
-- 
gmx-users mailing list    [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [email protected].
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to