On Sat, 5 Jul 2008 10:40:21 -0400
"rams rams" <[EMAIL PROTECTED]> wrote:
Dear users,
Is it possible to evaluate the rotational diffusion of proteins using
gromacs tools ??
No directly. However you can use g_rotacf to generate the autocorrelation
function of vectors (option -d). By defining vectors representing your
molecule/protein you can access the rotational correlation time of your
representative vector. You can imagine different way to get a statistically
significant value. One would be to define many vectors between backbone
atoms and average your results. Another would be to again define many
vectors but this time between the center of mass of the protein and each
Ca atoms and average ...
You can also hack the g_rms code to extract the rotational matrix during
the overlay of your protein to a reference structure and apply it to
a unit vector from whose trajectory you can again use g_rotacf to get
the autocorrelation function of that vector ...
An important point in the comparison of your result to experimental
values is the way the rotational correlation time is extracted
experimentally. They select different mode of relaxation (1 or 2) and thus
you have to use the corresponding Legendre polynomial when calculating the
autocorrelation function. From NMR relaxation l=2.
XAvier.
-----------------------------------------------------
XAvier Periole - PhD
Molecular Dynamics Group
- NMR and Computation -
University of Groningen
The Netherlands
-----------------------------------------------------
_______________________________________________
gmx-users mailing list [email protected]
http://www.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
