Please does anyone know how is calculated the B factor in gromacs? What is the formula that gives the B factor with the average coordinates with g_rmsf?<br><br>Carla<br><br><div class="gmail_quote">On Fri, Jun 11, 2010 at 10:53 AM, Erik Marklund <span dir="ltr"><<a href="mailto:erikm@xray.bmc.uu.se">erikm@xray.bmc.uu.se</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi,<br>
<br>
yes it can be done with g_covar, but such an average structure is often of little physical significace. Imagine for instance the acerage structure of a rotating methyle group (it's a bunch of atoms on a line).<br>
<br>
Erik<br>
<br>
Carla Jamous skrev:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im">
Hi Tsjerk,<br>
<br>
thank you for your answer.<br>
Actually, for the initial structure, I took the values of the B factor, and calculated the mean square displacement per atom. This is what I meant by saying RMSF of initial structure. Anyway, thanks for the explanation.<br>
<br>
But I have another question: I need to take 3 structures and make an average structure of these 3. Is there a way to do it with gromacs?<br>
<br>
Cheers,<br>
Carla<br>
<br>
<br></div><div class="im">
On Thu, Jun 10, 2010 at 12:15 PM, Tsjerk Wassenaar <<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a> <mailto:<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a>>> wrote:<br>
<br>
Hi Carla,<br>
<br>
On Thu, Jun 10, 2010 at 12:03 PM, Carla Jamous<br></div><div><div></div><div class="h5">
<<a href="mailto:carlajamous@gmail.com" target="_blank">carlajamous@gmail.com</a> <mailto:<a href="mailto:carlajamous@gmail.com" target="_blank">carlajamous@gmail.com</a>>> wrote:<br>
> Hi Everyone,<br>
><br>
> please I have a question concerning g_rmsf.<br>
> I need to compare the RMSF from my initial structure to the RMSF<br>
of my<br>
> average structure.<br>
<br>
Single structures (initial c.q. average) do not have an RMSF.<br>
<br>
> When I did<br>
><br>
> g_rmsf -s .tpr -f .xtc -n .ndx -b xx -e yy -o rmsf<br>
><br>
> or<br>
><br>
> g_rmsf -s .tpr -f .xtc -n .ndx -b xx -e yy -res -o rmsf<br>
<br>
This is not what you did. Please copy/paste command lines.<br>
<br>
> I got the same result, when choosing C-alpha for root mean square<br>
> calculation.<br>
<br>
Sure, when selecting C-alphas, averaging the RMSF per residue<br>
(sum_over_calphas_in_residue/number_of_calphas_in_residue) will<br>
evidently be identical to calculating the RMSF on an atom basis for<br>
each Calpha.<br>
<br>
> So please can anyone explain how can I get the average per<br>
residue over a<br>
> period of time?<br>
<br>
Select 'protein' (and use the -res flag).<br>
<br>
Cheers,<br>
<br>
Tsjerk<br>
<br>
<br>
--<br>
Tsjerk A. Wassenaar, Ph.D.<br>
<br>
post-doctoral researcher<br>
Molecular Dynamics Group<br>
Groningen Institute for Biomolecular Research and Biotechnology<br>
University of Groningen<br>
The Netherlands<br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div class="im"><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<div class="im"><br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br>
</div></blockquote>
<br>
<br>
-- <br>
-----------------------------------------------<br>
Erik Marklund, PhD student<br>
Dept. of Cell and Molecular Biology, Uppsala University.<br>
Husargatan 3, Box 596, 75124 Uppsala, Sweden<br>
phone: +46 18 471 4537 fax: +46 18 511 755<br>
<a href="mailto:erikm@xray.bmc.uu.se" target="_blank">erikm@xray.bmc.uu.se</a> <a href="http://folding.bmc.uu.se/" target="_blank">http://folding.bmc.uu.se/</a><br><font color="#888888">
<br>
-- <br></font><div><div></div><div class="h5">
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br>