<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div style="RIGHT: auto"><SPAN style="RIGHT: auto">Dear Justin</SPAN></div>
<div style="RIGHT: auto"><SPAN style="RIGHT: auto"></SPAN> </div>
<div style="RIGHT: auto"><SPAN style="RIGHT: auto">Thank you very muchh</SPAN></div>
<div style="RIGHT: auto"><SPAN style="RIGHT: auto"></SPAN> </div>
<div style="RIGHT: auto"><SPAN style="RIGHT: auto">Sogol<VAR id=yui-ie-cursor></VAR></SPAN></div>
<div><BR></div>
<DIV style="FONT-FAMILY: times new roman, new york, times, serif; FONT-SIZE: 12pt">
<DIV style="FONT-FAMILY: times new roman, new york, times, serif; FONT-SIZE: 12pt">
<DIV dir=ltr><FONT size=2 face=Arial>
<DIV style="BORDER-BOTTOM: #ccc 1px solid; BORDER-LEFT: #ccc 1px solid; PADDING-BOTTOM: 0px; LINE-HEIGHT: 0; MARGIN: 5px 0px; PADDING-LEFT: 0px; PADDING-RIGHT: 0px; HEIGHT: 0px; FONT-SIZE: 0px; BORDER-TOP: #ccc 1px solid; BORDER-RIGHT: #ccc 1px solid; PADDING-TOP: 0px" class=hr contentEditable=false readonly="true"></DIV><B><SPAN style="FONT-WEIGHT: bold">From:</SPAN></B> Justin A. Lemkul <jalemkul@vt.edu><BR><B><SPAN style="FONT-WEIGHT: bold">To:</SPAN></B> Kowsar Bagherzadeh <kw_bagherzadeh@yahoo.com>; Discussion list for GROMACS users <gmx-users@gromacs.org> <BR><B><SPAN style="FONT-WEIGHT: bold">Sent:</SPAN></B> Thursday, May 24, 2012 9:59 AM<BR><B><SPAN style="FONT-WEIGHT: bold">Subject:</SPAN></B> Re: [gmx-users] g_rms -bm<BR></FONT></DIV><BR><BR><BR>On 5/24/12 7:24 AM, Kowsar Bagherzadeh wrote:<BR>> <BR>> <BR>> Dear Users,<BR>> I am trying to analyze a ligand-protein simulation results. I read in the manual<BR>>
that using g_rms command with –bm option produces a matrix of average bond angle<BR>> deviations. And only bonds between atoms in the comparison groups are<BR>> considered. Does it mean that it is for the bonds and their angles that are<BR>> already in existence? (Not the ones that may be formed throughout simulation, I<BR>> mean the ligand may for example interact with residues through H-bonds) .I have<BR><BR>In this context, a "bond" means an actual chemical bond. A hydrogen bond is a nonbonded interaction.<BR><BR>> made a group in my index file named Active site (including only the active site<BR>> residues), and I have a LIG group as well. If I choose these two groups for<BR>> g_rms with this command:<BR>> /g_rms –s *.tpr –f *.trr –o rmsd.xvg –bm bond.xpm –n *.ndx**/<BR>> Does it show me how the ligand affects the active site residues bond angles?<BR><BR>Potentially.<BR><BR>> And one more question, how
can I study the ligand orientation in the active site?<BR><BR>That depends on how you define orientation - internal metrics like dihedrals or angles between planes of groups in the ligand, relative measurements like its position with respect to protein residues, etc. All analysis tools are listed in the manual, Chapter 8 and Appendix D. It's quite a lot to read, but you'll be able to identify all the various things you can analyze and how the information might be connected across different analysis routines.<BR><BR>-Justin<BR><BR>-- ========================================<BR><BR>Justin A. Lemkul, Ph.D.<BR>Research Scientist<BR>Department of Biochemistry<BR>Virginia Tech<BR>Blacksburg, VA<BR>jalemkul[at]<A href="http://vt.edu/" target=_blank>vt.edu</A> | (540) 231-9080<BR>http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<BR><BR>========================================<BR><BR><BR></DIV></DIV></div></body></html>