<div class="gmail_quote">On Thu, May 14, 2009 at 1:32 AM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im"><br>
<br>
Manik Mayur wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Users,<br>
<br>
I constructed a .pdb file for TIP4P water model which I tried to convert to .gro using pdb2gmx using:<br>
<br>
$ pdb2gmx -f Sim1.pdb -o Sim1.gro -water tip4p<br>
<br>
</blockquote>
<br></div>
There is no need to use pdb2gmx.<br>
<br>
If you just need a TIP4P solvent box, one exists in the /share/top subdirectory of your Gromacs installation (tip4p.gro).<br>
</blockquote><br>Thanks for the reply. I needed to use pdb2gmx as it automatically checks for errors in the configuration file. Moreover apart from water, my system has silicon walls enclosing water. <br><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
If you want to use your configuration, use editconf to convert it to .gro and #include "tip4p.itp" in your topology (which you can write by hand for a simple water system), or leave it as .pdb; Gromacs can use many different file formats for structures.</blockquote>
<div class="im"><br>yes, editconf will convert it to .gro but I was wondering why pdb2gmx is behaving in such unexpected manner?<br><br></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im">
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
with the choice of forcefield as 5 (OPLS-AA/L all-atom force field). After conversion I found that the position of HW3(MW), the dummy atom is not the same relative to OW, what I kept in the .pdb file. Even during the conversion pdb2gmx gave warnings as "Warning: Long Bond (1773-1776 =....".<br>
</blockquote>
<br></div>
What values does it find? (... does not help)<br>
</blockquote><br>e.g.<br>Warning: Long Bond (8637-8640 = 8.73449 nm)<br> .. and so on for many more pairs (all ~8 nm)<br><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Interestingly, if one runs pdb2gmx on the existing tip4p.gro that comes pre-installed, the program reports "Short bond" warnings.<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
My .pdb file goes like this:<br>
............<br>
ATOM 1773 OW SOL 108 45.944 29.409 58.762 1.00 0.00<br>
ATOM 1774 HW1 SOL 108 46.530 29.409 59.519 1.00 0.00<br>
ATOM 1775 HW2 SOL 108 46.530 29.409 58.005 1.00 0.00<br>
ATOM 1776 HW3 SOL 108 46.094 29.409 58.762 1.00 0.00<br>
............<br>
note that the OW - HW3 is 0.15 angs as it should be.<br>
<br>
So why is pdb2gmx recalculating the co-ordinates of HW3? How is the OW-HW3 bond is getting stretched? Also, what should be the correct procedure for the conversion?<br>
<br>
</blockquote>
<br></div>
The residue name "SOL" is interpreted as "HOH," or a three-point water molecule. OPLS calls tip4p HO4 (per the .rtp file). This could be a potential source of error.<br>
</blockquote><br>with "-water tip4p" in pdb2gmx, SOL is interpreted as HO4, the default one however is HOH.<br><br>Thanks,<br><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
-Justin<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im">
Thanks,<br>
<br>
Manik Mayur<br>
Graduate student<br>
Microfluidics Lab<br>
Dept. of Mechanical Engg.<br>
IIT Kharagpur<br>
INDIA<br>
<br>
<br></div>
------------------------------------------------------------------------<br>
<br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">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't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can'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>
</blockquote>
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<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>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">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't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can'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>
</blockquote></div><br>