Zuzana Benkova wrote:
Dear GROMACS users,
I have a system of 64 chains and want to calculate averaged values of
end-to-end distance of all chains as well as the averaged radius of
gyration from 1ns trajectory. In the case of end-to-end distance I have
prepared index file containing one group of the first atoms in ascending
order and the second one containing the terminal atoms in ascending
order. I have used g_dist but it looks like it combined all pairs
between both groups. I have looked at the mailing list but did not find
the answer.
g_dist should only be measuring the COM distance between the two selected
groups, so I don't know what you mean by "combined all pairs" - perhaps you can
provide an example (index groups, relevant output, what you were expecting, and
why you think it's wrong)? If your groups contain all starting residues in one,
then all ending residues in another, then yes, you're going to get weird
results. If you have one chain that is residues 1-100, you should have a group
for r_1 and a group for r_100; select those for g_dist analysis. If another
chain is residue 101-200, r_101 and r_200 are what you want, etc.
Is it possible to do such analysis using directly GROMACS tools or I
need to prepare some scripts with loops running over all chains. Thank
you for response in advance.
It certainly should be possible to iteratively use g_dist to measure distances
between whatever groups you select. Fairly straightforward, actually, so if
you're not getting what you expect, provide the information I asked for above
and hopefully someone should be able to figure out what's going wrong.
-Justin
Zuzana Benkova
--
========================================
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 [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php