<br><br>
<div class="gmail_quote">On Tue, Oct 18, 2011 at 1:36 PM, lloyd riggs <span dir="ltr"><<a href="mailto:lloyd.riggs@gmx.ch">lloyd.riggs@gmx.ch</a>></span> wrote:<br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">Message: 1<br>Date: Tue, 18 Oct 2011 06:41:55 -0400<br>From: "Justin A. Lemkul" <<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>><br>
Subject: Re: [gmx-users] Interaction energy<br>To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>><br>Message-ID: <<a href="mailto:4E9D57F3.10905@vt.edu">4E9D57F3.10905@vt.edu</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br><br><br><br>Steven Neumann wrote:<br>><br>><br>> On Mon, Oct 17, 2011 at 6:02 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a><br>
> <mailto:<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>>> wrote:<br>><br>><br>><br>> Steven Neumann wrote:<br>><br>> Dear Gmx Users,<br>> I would like to calculate the interaction energy (LJ and<br>
> electrostatic) between each residue and my ligands (10 ligands<br>> in the system). I would like to see what is the contribution of<br>> electrostatic and vdW interactions between ligand and each of my<br>
> residue. I thought to use g_energy and specify each of my<br>> residues in index files but it is not possible. Will you suggest<br>> how to do this?<br>><br>><br>><br>> For such information, you have to specify these groups as energygrps<br>
> in the .mdp file. You can rerun the trajectory using mdrun -rerun<br>> and a new .tpr file specifying these groups, but depending on the<br>> output frequency, the result may not be as accurate as you'd like.<br>
><br>> -Justin<br>><br>><br>> Thank you Justin. Now I have two groups sepcified in my mdp file:<br>><br>> energygrps = Protein LIG<br>><br>> How can I specify each residue of my protein separately and each ligand?<br>
> In my md.gro file I have residues:<br>><br>> 91GLY 92TYR ..... 161LIG 162LIG...<br>><br>><br>> Will it be correct like this<br>><br>> energygrps = 91GLY 92 TYR ... 161LIG 162LIG...<br>><br>
<br>No. The names must correspond to valid groups in an index file.<br><br>> If yes, will this simulation take longer? Thank you<br><br>Perhaps, but certainly your energy file will be considerably larger.<br><br>-Justin<br>
<br>><br>> Steven<br>><br>><br>><br>><br>><br><br>>--<br>>========================================<br><br>Dear Steven,<br><br>I did the following, but like Justin mentioned the files get large. Basically, I made a very complex ndx file from the beggining. I have the whole proteins and peptides involved, then a series of loops, and protein minus loop, and then individual amino acids and protein minus the amino acids. Basically every time you make the ndx for a particular amino acid/chain/loop you also have to generate the ndx for the rest of your protein minus what you were looking at. Otherwise you run out of groups to assign, and will start getting "residue or atom in multiple groups" when it is referenced more than one time.<br>
<br>This might not be necessary in your case, but lets you divie up energy contributions easier, although multiple re-runs using varied indexes might be necessary, ie for the sum total, then varied parts, etc...<br><br>At least that is a suggestion from my experience..<br>
<br>Freundlichen Grüsse<br><br>Stephan Watkins<br>--<br>NEU: FreePhone - 0ct/min Handyspartarif mit Geld-zurück-Garantie!<br>Jetzt informieren: <a href="http://www.gmx.net/de/go/freephone" target="_blank">http://www.gmx.net/de/go/freephone</a><br>
<font color="#888888">--<br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">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">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>
<div><br> </div>
<div>Thank you Stephan! I have faced the problem you mentioned:</div>
<div> </div>
<div>Fatal error: Atom 1 in multiple Energy Mon. groups (1 and 3)</div>
<div> </div>
<div>So as I understood what you have written I should specify in my INDEX.ndx:</div>
<div> </div>
<div>[r_91] - residue 91</div>
<div>1 2 3.... 7</div>
<div>[Protein-r91] - Protein without residue 91</div>
<div>8 9 10 11... 382</div>
<div>[r_92] - resideu 92</div>
<div>8 9 10 11...17</div>
<div>[Protein-r92] - Protein without residue 92</div>
<div>1 2..7 18 19 20...38</div>
<div> </div>
<div>etc. etc. </div>
<div> </div>
<div>And the same with ligands but substracting from my [System] group?</div>
<div> </div>
<div>Then</div>
<div>[Ligand] - all Ligands atoms</div>
<div>[Protein] - all Protein atoms</div>
<div>[Water] - all Water atoms.</div>
<div> </div>
<div>Please let me whether I am correct. Did such simulation take much longer? I run 100ns simulations for 2 days on the cluster with specified energy groups for Protein and Ligand only. My system consists of app. 4000 atoms only (water+protein+ligands).</div>
<div>Thank you,</div>
<div> </div>
<div>Steven</div>
<div> </div>