Poojari, Chetan wrote:
Hi Justin,
I used the below g_select command to get the list of phosphorus atoms (from DPPC lipid) within 0.5 nm of protein.
g_select -f ordered.xtc -n initial.ndx -on select.ndx -select "name P and within 0.5
of group Protein" -s md_500ns.tpr -b 100000 -dt 1000
My output is an index file which contains different phosphorus atom number at
each frame within 0.5nm of protein and this number changes for every frame.
Please can I know if its possible to use this index file to convert each frame
into trajectory file and do analysis for lipids within 0.5 of protein.
Yes, but not in an automatic way. Each index group corresponds to the atoms
that satisfy the given criterion in a given frame. You can parse these atoms
from the trajectory by writing individual snapshots by selecting each of the
indices for each of the frames. Then you might be able to use trjcat to stitch
a trajectory back together, but if there are different numbers of atoms in
different frames, the trajectory will be invalid.
-Justin
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
--
gmx-users mailing list gmx-users@gromacs.org
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 gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
