<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
&gt; Hi!<br>
&gt; I&#39;m working with kinase with ATP in active site. ATP topology was<br>
&gt; created by <a href="http://swissparam.ch" target="_blank">swissparam.ch</a> &lt;<a href="http://swissparam.ch" target="_blank">http://swissparam.ch</a>&gt; server, as well as I&#39;m<br>
&gt; working with charmm27 in gromacs 4.5.3. Using superposition of the<br>
&gt; protein models I&#39;ve found  coordinates for ions which have to stabilize<br>
&gt; ATP (Mg2+). For this purpose we used MgCl in the real experiments. I<br>
&gt; know, that Zn atoms in protein I can bind changing protonation type of<br>
&gt; aminoacid CYS and HIS which formed 4 bonds with Zn. And I need your help<br>
&gt; in this question! Of course I shouldn&#39;t create &quot;bond&quot; and &quot;angel&quot;<br>
&gt; between etheric O and magnesium. May be I have to PULL this Mg2+ atoms?<br>
&gt; and if it is so - which PULLING method I need (distance or else)?<br>
&gt;<br>
<br>
I don&#39;t see how pulling is relevant here.  You could use it to enforce a<br>
restraint, but it&#39;s not necessary.  Distance restraints or type 6 bonds are<br>
probably more useful and straightforward.  You do certainly have the<br>
complication of whether or not the force field&#39;s charges for the amino acids are<br>
correct, which most likely they aren&#39;t, as most QM/MM studies would indicate.<br>
<br>
-Justin<br>
<br>
&gt; Thank you very much<br></blockquote></div><br>Thank you for a reply )))<br><br>First of all I&#39;ve understood that Pulling is a bad idea. I tried to apply several manuals, including this one <a href="http://www.gromacs.org/Documentation/How-tos/Distance_Restraints">http://www.gromacs.org/Documentation/How-tos/Distance_Restraints</a>. But what about typing of 6 bonds? Bonds between ATP and Mg, Mg aminoacids? <br>
This part of your message is unclear &quot;whether or not the force field&#39;s charges for the amino acids are
correct, which most likely they aren&#39;t&quot; what aa did you mean, those I described with Zn binding or supposed Mg binding?<br>Sorry for I repeat the question<br>