Hi!<br>It is actually a phosphorylated serine!! So I am wondering if GMX has any special building block for handling such special cases?<br>Jayant<br><br><br><br><div class="gmail_quote">On Wed, Apr 9, 2008 at 12:26 AM, Tsjerk Wassenaar &lt;<a href="mailto:tsjerkw@gmail.com">tsjerkw@gmail.com</a>&gt; wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Darn, forget that post, I wasn&#39;t reading clearly. Too early here. It&#39;s<br>
part of a phosphorylated residue, you could&#39;ve at least mentioned the<br>
residue to go with it. There&#39;s some stuff on the contributions pages<br>
of the gromacs site about phosphorylated residues. Have a look there.<br>
<br>
Sorry for the previous babbling... The reference to the wiki still goes though<br>
<br>
Cheers,<br>
<font color="#888888"><br>
Tsjerk<br>
</font><div><div></div><div class="Wj3C7c"><br>
<br>
<br>
On Wed, Apr 9, 2008 at 9:24 AM, Tsjerk Wassenaar &lt;<a href="mailto:tsjerkw@gmail.com">tsjerkw@gmail.com</a>&gt; wrote:<br>
&gt; Hi Jayant,<br>
&gt;<br>
&gt; &nbsp;PO3 seems to me the anion of metaphosphoric acid, not much to do with<br>
&gt; &nbsp;phosphoric acid. The latter could be sort of extracted from phosphoric<br>
&gt; &nbsp;acid containing residues, although that doesn&#39;t give guarantees for<br>
&gt; &nbsp;good behaviour. Also, I wouldn&#39;t count on the PO4 parameters thus<br>
&gt; &nbsp;extracted to have much to do with PO3. Very likely that the atomic<br>
&gt; &nbsp;charges are quite different, but maybe LJ parameters too.<br>
&gt;<br>
&gt; &nbsp;Have a look at<br>
&gt;<br>
&gt; &nbsp;<a href="http://wiki.gromacs.org/index.php/Parameterization" target="_blank">http://wiki.gromacs.org/index.php/Parameterization</a><br>
&gt;<br>
&gt; &nbsp;Cheers,<br>
&gt;<br>
&gt; &nbsp;Tsjerk<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; &nbsp;On Wed, Apr 9, 2008 at 3:44 AM, jayant james &lt;<a href="mailto:jayant.james@gmail.com">jayant.james@gmail.com</a>&gt; wrote:<br>
&gt; &nbsp;&gt; Hi !<br>
&gt; &nbsp;&gt; I am attempting to simulate a protein that is phosphorylated below is a part<br>
&gt; &nbsp;&gt; of the PDF file<br>
&gt; &nbsp;&gt; HETATM &nbsp;460 &nbsp;P &nbsp; PO3 A 200 &nbsp; &nbsp; &nbsp;14.995 &nbsp; 1.523 &nbsp; 4.011 &nbsp;1.00 &nbsp;0.89<br>
&gt; &nbsp;&gt; P<br>
&gt; &nbsp;&gt; HETATM &nbsp;461 &nbsp;O1 &nbsp;PO3 A 200 &nbsp; &nbsp; &nbsp;14.287 &nbsp; 1.882 &nbsp; 5.272 &nbsp;1.00 &nbsp;1.16<br>
&gt; &nbsp;&gt; O<br>
&gt; &nbsp;&gt; &nbsp;HETATM &nbsp;462 &nbsp;O2 &nbsp;PO3 A 200 &nbsp; &nbsp; &nbsp;15.841 &nbsp; 0.179 &nbsp; 4.048 &nbsp;1.00 &nbsp;0.77<br>
&gt; &nbsp;&gt; O<br>
&gt; &nbsp;&gt; HETATM &nbsp;463 &nbsp;O3 &nbsp;PO3 A 200 &nbsp; &nbsp; &nbsp;14.134 &nbsp; 1.500 &nbsp; 2.694 &nbsp;1.00 &nbsp;0.97<br>
&gt; &nbsp;&gt; O<br>
&gt; &nbsp;&gt; HETATM &nbsp;464 &nbsp;P &nbsp; PO3 A 201 &nbsp; &nbsp; &nbsp;13.207 &nbsp;-1.138 &nbsp;-1.342 &nbsp;1.00 &nbsp;1.31<br>
&gt; &nbsp;&gt; P<br>
&gt; &nbsp;&gt; &nbsp;HETATM &nbsp;465 &nbsp;O1 &nbsp;PO3 A 201 &nbsp; &nbsp; &nbsp;12.997 &nbsp;-2.551 &nbsp;-0.914 &nbsp;1.00 &nbsp;1.43<br>
&gt; &nbsp;&gt; O<br>
&gt; &nbsp;&gt; HETATM &nbsp;466 &nbsp;O2 &nbsp;PO3 A 201 &nbsp; &nbsp; &nbsp;11.918 &nbsp;-0.211 &nbsp;-1.324 &nbsp;1.00 &nbsp;1.47<br>
&gt; &nbsp;&gt; O<br>
&gt; &nbsp;&gt; HETATM &nbsp;467 &nbsp;O3 &nbsp;PO3 A 201 &nbsp; &nbsp; &nbsp;13.927 &nbsp;-0.914 &nbsp;-2.727 &nbsp;1.00 &nbsp;1.41<br>
&gt; &nbsp;&gt; O<br>
&gt; &nbsp;&gt;<br>
&gt; &nbsp;&gt; but upon executing pdb2gmx of the PDB file I get an error<br>
&gt; &nbsp;&gt; -----------------------------------------------------------------<br>
&gt; &nbsp;&gt; Fatal error:<br>
&gt; &nbsp;&gt; Residue &#39;PO3&#39; not found in residue topology database<br>
&gt; &nbsp;&gt; ------------------------------------------------------------------<br>
&gt; &nbsp;&gt;<br>
&gt; &nbsp;&gt; I do know that DNA can be simulated in GROMOS (option 0 of pdb2gmx). So I am<br>
&gt; &nbsp;&gt; wondering how to overcome the above error.<br>
&gt; &nbsp;&gt;<br>
&gt; &nbsp;&gt; Awaiting suggestions<br>
&gt; &nbsp;&gt; Thanks<br>
&gt; &nbsp;&gt; Jayant<br>
&gt; &nbsp;&gt;<br>
&gt; &nbsp;&gt;<br>
&gt; &nbsp;&gt; --<br>
&gt; &nbsp;&gt; Jayasundar Jayant James<br>
&gt; &nbsp;&gt;<br>
&gt; &nbsp;&gt; <a href="http://www.chick.com/reading/tracts/0096/0096_01.asp" target="_blank">www.chick.com/reading/tracts/0096/0096_01.asp</a>)<br>
&gt; &nbsp;&gt; Residence -24935864, cell-9841042164<br>
&gt; &nbsp;&gt; _______________________________________________<br>
&gt; &nbsp;&gt; &nbsp;gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt; &nbsp;&gt; &nbsp;<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt; &nbsp;&gt; &nbsp;Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
&gt; &nbsp;&gt; &nbsp;Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt; &nbsp;&gt; &nbsp;www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt; &nbsp;&gt; &nbsp;Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
&gt; &nbsp;&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; &nbsp;--<br>
&gt; &nbsp;Tsjerk A. Wassenaar, Ph.D.<br>
&gt; &nbsp;Junior UD (post-doc)<br>
&gt; &nbsp;Biomolecular NMR, Bijvoet Center<br>
&gt; &nbsp;Utrecht University<br>
&gt; &nbsp;Padualaan 8<br>
&gt; &nbsp;3584 CH Utrecht<br>
&gt; &nbsp;The Netherlands<br>
&gt; &nbsp;P: +31-30-2539931<br>
&gt; &nbsp;F: +31-30-2537623<br>
&gt;<br>
<br>
<br>
<br>
--<br>
Tsjerk A. Wassenaar, Ph.D.<br>
Junior UD (post-doc)<br>
Biomolecular NMR, Bijvoet Center<br>
Utrecht University<br>
Padualaan 8<br>
3584 CH Utrecht<br>
The Netherlands<br>
P: +31-30-2539931<br>
F: +31-30-2537623<br>
_______________________________________________<br>
gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Jayasundar Jayant James<br><br><a href="http://www.chick.com/reading/tracts/0096/0096_01.asp">www.chick.com/reading/tracts/0096/0096_01.asp</a>) <br>Residence -24935864, cell-9841042164