I searched the user list in order to get some idea about calculating cosine content of first PCA .. but I am not able to do so ??<div><br></div><div>Can u direct me to a good thread or discussion..</div><div><br></div><div>
<br><div><div class="gmail_quote">On Tue, Feb 1, 2011 at 2:04 AM, bharat gupta <span dir="ltr"><<a href="mailto:bharat.85.monu@gmail.com">bharat.85.monu@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Thanks for the reply ...<div><br></div><div>I am following this thread - <a href="http://oldwww.gromacs.org/pipermail/gmx-users/2004-May/010438.html" target="_blank">http://oldwww.gromacs.org/pipermail/gmx-users/2004-May/010438.html</a></div>
<div><a href="http://oldwww.gromacs.org/pipermail/gmx-users/2004-May/010438.html" target="_blank"></a>Also I read the following paper <a href="http://pubs.acs.org/doi/abs/10.1021/jp983120q" target="_blank">http://pubs.acs.org/doi/abs/10.1021/jp983120q</a> (Internal dynamics of GFP) in which they have carried out such projections of the eigen vector in order to find the regions showing minimum and maximum amplitudes..</div>
<div><br></div><div>What I am thinking is shall I take those c-alpha atoms that are well stable and equilibrated after referring to RMSD graph generated after simulation ...</div><div><br></div><div>I am stating here what I have done so far and what I have to analyze :-</div>
<div><br></div><div>I exchanged the loops regions of GFP with loops of longer length taken from pdb .. I simulated both the crystal str. and modeled strs. to compare how the loop incorporation leads to the change in the structure of GFP .... Since the RMSD of the GFP crystal str and that of model are in comparable range .. I thought of looking at the motions of loops in the structure that are giving such fluctuations in RMSD of model and for that I decided to generate eigen vectors and project them in 3d (like the one given in paper mentioned above) to see how that particular loop behaves .. I don't know whether what I have hypothesized is correct or not .. since I am not able to find any other insilico method apart from MDS for such a study .. Also if my prediction comes out to be good enough I have to carry out such insertion in the wet lab and has to check for fluorescence ..any help in this regard will be appreciated... </div>
<div><div><div></div><div class="h5"><br><br><div class="gmail_quote">On Tue, Feb 1, 2011 at 1:28 AM, Tsjerk Wassenaar <span dir="ltr"><<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Bharat,<br>
<div><br>
On Tue, Feb 1, 2011 at 8:33 AM, bharat gupta <<a href="mailto:bharat.85.monu@gmail.com" target="_blank">bharat.85.monu@gmail.com</a>> wrote:<br>
> Hi,<br>
> After searching though gmxuserlist I found the relevant thread for ED... I<br>
> followed the following steps :-<br>
<br>
</div>What thread are you referring to?<br>
<div><br>
> 1) g_covar - to generate a covariance matrix and diagonalize it (for c-alpha<br>
> atoms only)<br>
> 2) g_anaeig - to generate eigen vectors<br>
<br>
</div>g_covar calculates the eigenvectors. It's what you end up with through<br>
diagonalization.<br>
<div><br>
> 3)g_rmsf - for calculating RMSD of first 8 eigen vectors , I got an error at<br>
> this step that -<br>
<br>
</div>g_rmsf is not for calculating RMSD, and hasn't got much to do with<br>
eigenvector analysis. That's what g_anaeig is for (like with the<br>
option -rmsf). Using g_anaeig will avoid the error you observe.<br>
<br>
As a side note, 3 ns is rather short for this sort of thing. You have<br>
to check the cosine content of the first principal components to see<br>
if you've reached equilibrium already.<br>
<br>
Cheers,<br>
<br>
Tsjerk<br>
<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>
* Zernike Institute for Advanced Materials<br>
University of Groningen<br>
The Netherlands<br>
<font color="#888888">--<br>
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/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before 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>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</font></blockquote></div><br><br clear="all"><br></div></div><div class="im">-- <br>Bharat<br>Ph.D. Candidate<br>Room No. : 7202A, 2nd Floor<br>Biomolecular Engineering Laboratory<br>Division of Chemical Engineering and Polymer Science<br>
Pusan National University<br>
Busan -609735<br>South Korea<br>Lab phone no. - +82-51-510-3680, +82-51-583-8343<div>Mobile no. - 010-5818-3680<br>E-mail : <a href="mailto:monu46010@yahoo.com" target="_blank">monu46010@yahoo.com</a></div><br>
</div></div>
</blockquote></div><br><br clear="all"><br>-- <br>Bharat<br>Ph.D. Candidate<br>Room No. : 7202A, 2nd Floor<br>Biomolecular Engineering Laboratory<br>Division of Chemical Engineering and Polymer Science<br>Pusan National University<br>
Busan -609735<br>South Korea<br>Lab phone no. - +82-51-510-3680, +82-51-583-8343<div>Mobile no. - 010-5818-3680<br>E-mail : <a href="mailto:monu46010@yahoo.com" target="_blank">monu46010@yahoo.com</a></div><br>
</div></div>