<P>
&nbsp; <BR>
Thank you very much Justin<BR>
<BR>
On Mon, 02 Jun 2008 Justin A.Lemkul wrote :<BR>
&gt;Quoting minnale &lt;minnale_gnos@rediffmail.com&gt;:<BR>
&gt;<BR>
&gt; &gt;<BR>
&gt; &gt; Thanks for your prompt reply to Erik and Florian<BR>
&gt; &gt; i have done according to Erik's reply but there is no change in the RMSD<BR>
&gt; &gt; values<BR>
&gt; &gt; iam mentioning what i have done<BR>
&gt; &gt;<BR>
&gt; &gt; 1.make_ndx -f protein.gro -o calpha_re20_90.ndx<BR>
&gt; &gt;&nbsp; selected 3 &amp; r 20-90 ( 3 = calpha)<BR>
&gt; &gt;<BR>
&gt; &gt; this ndx gave input for g_rms<BR>
&gt; &gt; g_rms -f 7nsprot.xtc -n calpha_re20_90.ndx -pbc -s minwat -o rms_calpha.xvg<BR>
&gt; &gt;&nbsp;  selected calpha_res20-90 for least square fit<BR>
&gt; &gt;&nbsp;  selected calpha for RMSD<BR>
&gt; &gt; but there is no difference in RMSD values of protein, selected<BR>
&gt; &gt; protein&amp;res20-90(index for segment) and calpha&amp;res20-90(suggested one)<BR>
&gt; &gt;<BR>
&gt; &gt; have i made any mistake in mentioning options in commands<BR>
&gt;<BR>
&gt;Yes, as Erik already mentioned, you must use the calpha_res20-90 group for both<BR>
&gt;the fitting and the calculation.&nbsp; What you're showing is essentially the same<BR>
&gt;thing you already showed in your previous message - you are calculating the<BR>
&gt;RMSD of the entire 'calpha' group for the protein.<BR>
&gt;<BR>
&gt;-Justin<BR>
&gt;<BR>
&gt; &gt;<BR>
&gt; &gt; Thanks in advance.<BR>
&gt; &gt;<BR>
&gt; &gt; On Mon, 02 Jun 2008 Erik Marklund wrote :<BR>
&gt; &gt; &gt;minnale skrev:<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;&nbsp; Hi all,<BR>
&gt; &gt; &gt;&gt;I have run the protein in water for 7ns and plotted RMSD for whole protein<BR>
&gt; &gt; with 250 residues and segmented 70 residues out of 250 residues by using<BR>
&gt; &gt; index file, result in two cases values are nearly same and showing almost<BR>
&gt; &gt; same graph. I think it should be different in two cases.<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;the commands I have given for whole protein and segmented protein are<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;g_rms -f 7ns_pro.xtc -pbc -s min_pro.tpr -o rms_whol7nsprot<BR>
&gt; &gt; &gt;&gt;&nbsp;  selected c-alpha for least square fit<BR>
&gt; &gt; &gt;&gt;&nbsp;  selceted c-alpha for RMSD calculation<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;g_rms -f 7ns_pro.xtc -pbc -s min_pro.tpr -n 20_90res.ndx -o rms_segm7nsprot<BR>
&gt; &gt; &gt;&gt;&nbsp;  selected 20_90res for least square fit<BR>
&gt; &gt; &gt;&gt;&nbsp;  selected c-alpha for RMSD calculation<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;Any suggestions are appreciated<BR>
&gt; &gt; &gt;&gt;Thanks in advance<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;Swiss<BR>
&gt; &gt;<BR>
&gt;&lt;http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-home.htm/1050715198@Middle5/2133997_2126411/2134144/1?PARTNER=3&amp;OAS_QUERY=null&gt;<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;------------------------------------------------------------------------<BR>
&gt; &gt; &gt;&gt;<BR>
&gt; &gt; &gt;&gt;_______________________________________________<BR>
&gt; &gt; &gt;&gt;gmx-users mailing list&nbsp; &nbsp; gmx-users@gromacs.org<BR>
&gt; &gt; &gt;&gt;http://www.gromacs.org/mailman/listinfo/gmx-users<BR>
&gt; &gt; &gt;&gt;Please search the archive at http://www.gromacs.org/search before posting!<BR>
&gt; &gt; &gt;&gt;Please don't post (un)subscribe requests to the list. Use the www interface<BR>
&gt; &gt; or send it to gmx-users-request@gromacs.org.<BR>
&gt; &gt; &gt;&gt;Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>
&gt; &gt; &gt;Hi,<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;You have in fact calculated RMSD for all calphas in both cases, but only<BR>
&gt; &gt; used a subset of them for fitting in the second case. You need to create an<BR>
&gt; &gt; index group for the calphas in your segment and use that for your calculation<BR>
&gt; &gt; (and fitting).<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;-- -----------------------------------------------<BR>
&gt; &gt; &gt;Erik Marklund, PhD student<BR>
&gt; &gt; &gt;Laboratory of Molecular Biophysics,<BR>
&gt; &gt; &gt;Dept. of Cell and Molecular Biology, Uppsala University.<BR>
&gt; &gt; &gt;Husargatan 3, Box 596,&nbsp; &nbsp; 75124 Uppsala, Sweden<BR>
&gt; &gt; &gt;phone:&nbsp; &nbsp; +46 18 471 4537&nbsp; &nbsp; &nbsp; &nbsp; fax: +46 18 511 755<BR>
&gt; &gt; &gt;erikm@xray.bmc.uu.se&nbsp; &nbsp; http://xray.bmc.uu.se/molbiophys<BR>
&gt; &gt; &gt;<BR>
&gt; &gt;<BR>
&gt;<BR>
&gt;<BR>
&gt;<BR>
&gt;========================================<BR>
&gt;<BR>
&gt;Justin A. Lemkul<BR>
&gt;Graduate Research Assistant<BR>
&gt;Department of Biochemistry<BR>
&gt;Virginia Tech<BR>
&gt;Blacksburg, VA<BR>
&gt;jalemkul@vt.edu | (540) 231-9080<BR>
&gt;http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/<BR>
&gt;<BR>
&gt;========================================<BR>

</P>
<br><br>
<Table border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style='font-family:Verdana;font-size:11px;line-height:15px;'><TR><td><img src ='http://imadworks.rediff.com/cgi-bin/AdWorks/adimage.cgi/2129697_2122129/creative_2129852.gif'  alt='Amity'  border=0></td></TR></Table>