<p>Hi,</p>
<p>Do mind that calcium binding may involve (quantum) effects that are ill captured by classical force fields. If the binding site plays a central role in your research question, this may be problematic.</p>
<p>Cheers,</p>
<p>Tsjerk</p>
<p><blockquote type="cite">On Oct 16, 2010 2:34 PM, &quot;Justin A. Lemkul&quot; &lt;<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>&gt; wrote:<br><br><p><font color="#500050">

leila karami wrote:
&gt;
&gt; Dear mohsen ramezanpour/ /
&gt;
&gt; yes. calcium is a typical ligand in my pdb ...</font></p>
There is no need to invoke PRODRG for any of this.  Calcium ions are present in nearly all the force fields in Gromacs, and require no special intervention - pdb2gmx can handle them.  PRODRG is only useful for small molecule topologies, (which are usually not so accurate) so passing it an entire protein is not going to work (as you can see).  It will also fail with an error if you pass it a single Ca ion.<br>

<br>
-Justin<br>
<br>
-- <br>
========================================<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> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br><font color="#888888">
-- <br></font><p><font color="#500050">gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a>
...</font></p></blockquote></p>