<br><br>
<div class="gmail_quote">On Tue, Oct 18, 2011 at 2:50 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote"><br><br>Steven Neumann wrote:<br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<div class="im"><br><br>On Tue, Oct 18, 2011 at 1:36 PM, lloyd riggs <<a href="mailto:lloyd.riggs@gmx.ch" target="_blank">lloyd.riggs@gmx.ch</a> <mailto:<a href="mailto:lloyd.riggs@gmx.ch" target="_blank">lloyd.riggs@gmx.ch</a>>> wrote:<br>
<br> Message: 1<br> Date: Tue, 18 Oct 2011 06:41:55 -0400<br></div> From: "Justin A. Lemkul" <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>>
<div class="im"><br> Subject: Re: [gmx-users] Interaction energy<br> To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><u></u>><br>
Message-ID: <<a href="mailto:4E9D57F3.10905@vt.edu" target="_blank">4E9D57F3.10905@vt.edu</a> <mailto:<a href="mailto:4E9D57F3.10905@vt.edu" target="_blank">4E9D57F3.10905@vt.edu</a>><u></u>>
<div class="im"><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<br> <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>><br>
</div>
<div>
<div></div>
<div class="h5"> > <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">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<br>
ligands<br> > in the system). I would like to see what is the<br> contribution of<br> > electrostatic and vdW interactions between ligand and<br> 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<br> suggest<br> > how to do this?<br> ><br> ><br> ><br> > For such information, you have to specify these groups as<br>
energygrps<br> > in the .mdp file. You can rerun the trajectory using mdrun<br> -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<br>
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<br>
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> >=============================<u></u>===========<br><br> Dear Steven,<br><br> I did the following, but like Justin mentioned the files get large.<br>
Basically, I made a very complex ndx file from the beggining. I<br> have the whole proteins and peptides involved, then a series of<br> loops, and protein minus loop, and then individual amino acids and<br> protein minus the amino acids. Basically every time you make the<br>
ndx for a particular amino acid/chain/loop you also have to generate<br> the ndx for the rest of your protein minus what you were looking at.<br> Otherwise you run out of groups to assign, and will start getting<br>
"residue or atom in multiple groups" when it is referenced more than<br> one time.<br><br> This might not be necessary in your case, but lets you divie up<br> energy contributions easier, although multiple re-runs using varied<br>
indexes might be necessary, ie for the sum total, then varied parts,<br> 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/<u></u>freephone</a><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/<u></u>mailman/listinfo/gmx-users</a><br> Please search the archive at<br> <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>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></div> <mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>>.
<div class="im"><br> Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br><br><br> Thank you Stephan! I have faced the problem you mentioned:<br>
Fatal error: Atom 1 in multiple Energy Mon. groups (1 and 3)<br> So as I understood what you have written I should specify in my INDEX.ndx:<br> [r_91] - residue 91<br>1 2 3.... 7<br>[Protein-r91] - Protein without residue 91<br>
8 9 10 11... 382<br>[r_92] - resideu 92<br>8 9 10 11...17<br>[Protein-r92] - Protein without residue 92<br>1 2..7 18 19 20...38<br> etc. etc.<br> And the same with ligands but substracting from my [System] group?<br> <br>
</div></blockquote><br>This is not a correct approach, per the fatal error. If you want to monitor energies involving different residues, the only requirement is that those residues be in their own distinct groups, i.e.:<br>
<br>energygrps = r_91 r_92 r_93 ... Ligand<br><br>The energy groups cannot be redundant. What Stephan was suggesting is that you cannot, for instance, specify:<br><br>energygrps = Protein r_91 r_92 ...<br><br>because "Protein" encompasses residues 91, 92, etc. You can, however, specify energygrps as in my first example, such that the whole "Protein" group is not one of the energygrps. If you want energies involving those residues, you have to do as Stephan suggested and create an "incomplete" protein group that excludes the residues that you want to specifically analyze.<br>
<br>-Justin<br></blockquote>
<div> </div>
<div>Thank you both!!!!</div>
<div> </div>
<div>Steven</div>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote"><br>-- <br>==============================<u></u>==========<br><br>Justin A. Lemkul<br>Ph.D. Candidate<br>ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>Department of Biochemistry<br>Virginia Tech<br>Blacksburg, VA<br>jalemkul[at]<a href="http://vt.edu/" target="_blank">vt.edu</a> | <a href="tel:%28540%29%20231-9080" target="_blank" value="+15402319080">(540) 231-9080</a><br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a><br><br>==============================<u></u>==========
<div>
<div></div>
<div class="h5"><br>-- <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/<u></u>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/<u></u>Support/Mailing_Lists/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/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br></div></div></blockquote></div><br>