This is a tough issue and there are some unresolved questions here. <div><br></div><div>I recommend seeing the work of Kastenholz and Hunenberger on ion solvation in JCP a couple years ago (two companion papers) as these explain the key issues.</div>
<div><br>Thanks.<br><br><div class="gmail_quote">On Wed, Aug 17, 2011 at 3:24 AM, Marcelino Arciniega Castro <span dir="ltr"><<a href="mailto:macastro@mytum.de">macastro@mytum.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi there,<br>
<br>
Could someone please share some ideas about the next situation?:<br>
<br>
I am trying to estimate the relative binding affinity of two small molecules to a protein. In the setup process, ions are added in the solvent to neutralize the system. Then as everybody suggest, the ligand partial charges should be decouple to zero. But my specific case, the small molecules are charged, that means that after turning off the ligand partial charges, the system will have a net charge (because the ions in the solvent). Turning off the ions charges does not solve the problem, because of the protein charges.<br>
<br>
The concrete question are:<br>
Should I continue to atom mutation step using PME in a non-neutral system?<br>
It would make sense to turn off the partial charges of the some atoms in the small molecule (just those that are going to be mutated)?<br>
<br>
Thanks in advance<br>
Marcelino<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><br><br clear="all"><div><br></div>-- <br>David Mobley<br><a href="mailto:dmobley@gmail.com" target="_blank">dmobley@gmail.com</a><br>504-383-3662<br><br><br>
</div>