<br>
<div class="gmail_quote">On Tue, Mar 6, 2012 at 11:44 AM, Mark Abraham <span dir="ltr">&lt;<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>&gt;</span> wrote:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">
<div text="#000000" bgcolor="#FFFFFF">
<div class="im">On 6/03/2012 10:40 PM, Steven Neumann wrote: 
<blockquote type="cite"><br><br>
<div class="gmail_quote">On Tue, Mar 6, 2012 at 11:35 AM, Mark Abraham <span dir="ltr">&lt;<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>&gt;</span> wrote:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">
<div>
<div>On 6/03/2012 8:21 PM, Steven Neumann wrote:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">Dear Gmx Users,<br>I am running my protein with ligand. at the begining of my simulation the protein is far away from the ligand. Then in about 5-10 ns my ligand binds to the protein surface remaining there for rest 90-95 ns. Would you adjust coupling groups in thermostat as:<br>
a) Protein_LIg and Water_and_ions<br>b) Protein Water_ions_LIG<br>c) Protein LIG Water_ions<br><br></blockquote><br></div></div>I&#39;d consider a single coupling group before any of those. See final part of manual 3.4.8 and ref therein. c) is unsound because the ligand is likely too small. The lack of T-group atom contiguity in a) might be an issue, I&#39;m not sure. b) is probably what people do, but you know this from your literature searching about similar work, right? :-)</blockquote>

<div> </div>
<div>Well...I used a) approach as Ligand becaome a part of the protein quite fast.</div></div></blockquote><br></div>Sure, but if the path that took it there might be unsound, then you can&#39;t refute the argument that the fact that that path occurred might be a numerical artefact. 
<div class="im"><br><br>
<blockquote type="cite">
<div class="gmail_quote">
<div>The b) just was my own tought as at the begining of the simulation is in water only, however just for 500 ps as far as I see out of 100ns. How would you tackle this?</div></div></blockquote><br></div>Didn&#39;t I answer this already? :-)</div>
</blockquote>
<div> </div>
<div>My ligand will stack to everything and I put it around my protein in different positions and run separate simulations. Then the final coordinates will be chosen for binding free energy calculations (umbrella sampling) so think it is better approach.</div>

<div> </div>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">
<div text="#000000" bgcolor="#FFFFFF"><span class="HOEnZb"><font color="#888888"><br><br>Mark<br></font></span></div><br>--<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&#39;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&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</blockquote></div><br>