<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>
<div>I did this with a small molecule and the ions were in the solvent, but associated with the ligand, and conversly if the site has a Mg or something etc...it wouldn&#39;t be restrained in the normal posres.itp unless you made one for them.&nbsp; It is in the end a matter of view, but I am assuming the change in the binding sites energy would reflect also the ions or the ions being attached to the ligand would not represent a completly solvated ligand (ie energy missing).&nbsp; The positional restraint I stated meant you can just make a two ion -posres_ion and include it in an .mdp however it might screw your system up if your protein is not positionally restrained, or calphas etc...thus I dont know based on your experiment.&nbsp; If the ions move into solution it probably doesnt matter (a quick visual of the simulation).&nbsp; But somone else might have better ideas to help out.</div>

<div>&nbsp;</div>

<div>Stephan</div>

<div>&nbsp;
<div name="quote" style="margin:10px 5px 5px 10px; padding: 10px 0 10px 10px; border-left:2px solid #C3D9E5; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div style="margin:0 0 10px 0;"><b>Gesendet:</b>&nbsp;Mittwoch, 31. Juli 2013 um 12:46 Uhr<br/>
<b>Von:</b>&nbsp;&quot;Steven Neumann&quot; &lt;s.neumann08@gmail.com&gt;<br/>
<b>An:</b>&nbsp;&quot;Discussion list for GROMACS users&quot; &lt;gmx-users@gromacs.org&gt;<br/>
<b>Betreff:</b>&nbsp;Re: [gmx-users] Umbrella Sampling</div>

<div name="quoted-content">They do not dissociate...Are you sure? My mdp specifies only ligand as a<br/>
pull_group1. I think it would change having ions in this group included.<br/>
<br/>
<br/>
On Wed, Jul 31, 2013 at 11:01 AM, lloyd riggs &lt;lloyd.riggs@gmx.ch&gt; wrote:<br/>
<br/>
&gt;<br/>
&gt; &gt;will get the PMF profile for my<br/>
&gt; &gt;ligand binding or ligand and two ions binding?<br/>
&gt;<br/>
&gt; It would be the ligand and two ions unless the ions also at some point<br/>
&gt; discossiate from the ligand once in solvent. Could add positional restraint<br/>
&gt; for them, but dont know how that effects the calculation?<br/>
&gt; *Gesendet:* Mittwoch, 31. Juli 2013 um 09:29 Uhr<br/>
&gt; *Von:* &quot;Steven Neumann&quot; &lt;s.neumann08@gmail.com&gt;<br/>
&gt; *An:* &quot;Discussion list for GROMACS users&quot; &lt;gmx-users@gromacs.org&gt;<br/>
&gt; *Betreff:* [gmx-users] Umbrella Sampling<br/>
&gt; Dear Gmx Users,<br/>
&gt;<br/>
&gt; I run SMD to extract the windows for US calculations. The system involves<br/>
&gt; negatively charged ligand and protein. I generated the protein-ligand<br/>
&gt; complex within self assembly MD simulations.<br/>
&gt;<br/>
&gt; I pulled my molecule away and two ions were also detached from the protein<br/>
&gt; surface being attached to my ligand.<br/>
&gt;<br/>
&gt; My question: if I run my US caluclation and combine windows by WHAM (I<br/>
&gt; specified in my umbrella.mdp my ligand as a pull_group1 and same protein<br/>
&gt; residues I pulled it from as pull_group0) will get the PMF profile for my<br/>
&gt; ligand binding or ligand and two ions binding?<br/>
&gt;<br/>
&gt; Steven<br/>
&gt; --<br/>
&gt; gmx-users mailing list gmx-users@gromacs.org<br/>
&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br/>
&gt; * Please search the archive at<br/>
&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br/>
&gt; * Please don&#39;t post (un)subscribe requests to the list. Use the<br/>
&gt; www interface or send it to gmx-users-request@gromacs.org.<br/>
&gt; * 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/>
&gt;<br/>
&gt; --<br/>
&gt; gmx-users mailing list gmx-users@gromacs.org<br/>
&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br/>
&gt; * Please search the archive at<br/>
&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br/>
&gt; * Please don&#39;t post (un)subscribe requests to the list. Use the<br/>
&gt; www interface or send it to gmx-users-request@gromacs.org.<br/>
&gt; * 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/>
&gt;<br/>
--<br/>
gmx-users mailing list gmx-users@gromacs.org<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 gmx-users-request@gromacs.org.<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></div>
</div>
</div>
</div></div></body></html>