<br><div>This is exactly what I thought.</div><div>I found out that saving the file with .mol ending in PyMol does not give the unwanted bonds.</div><div><br></div><div>Thank you all.</div><div><br class="webkit-block-placeholder">
</div><div><br><br><div class="gmail_quote">On Thu, Nov 12, 2009 at 8:49 AM, Ran Friedman <span dir="ltr">&lt;<a href="mailto:r.friedman@bioc.uzh.ch">r.friedman@bioc.uzh.ch</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
With VMD it&#39;s even simpler: use &quot;dynamic bonds&quot;.<br>
<font color="#888888"><br>
Ran.<br>
</font><div><div></div><div class="h5"><br>
Nicolas Sapay wrote:<br>
&gt;<br>
&gt;<br>
&gt; Dallas B. Warren a écrit :<br>
&gt;&gt;<br>
&gt;&gt; Coordinate files like pdb and gro aren’t used by GROMACS to provide<br>
&gt;&gt; any bonding information. That is what the topology files are for. So<br>
&gt;&gt; their “presence” in your pdb isn’t an issue.<br>
&gt;&gt;<br>
&gt;&gt; Actually, what is probably happening is that PyMol is guessing the<br>
&gt;&gt; bonds presence, based on the distance between atoms, and displaying<br>
&gt;&gt; it (which is what VMD does as well). So the bonds aren’t actually<br>
&gt;&gt; there at all in the pdb file.<br>
&gt;&gt;<br>
&gt; VMD can read/write CHARMM/NAMD topology files (namely psf files). If<br>
&gt; your problem is *just* a visualization artefact, you can load your<br>
&gt; structure in VMD, write a psf file and delete the unwanted bonds (this<br>
&gt; does not require a specific forcefield). You can also combine 2 psf<br>
&gt; files (1 for your protein and 1 for your ligand), that require some<br>
&gt; basic knowledge of TCL though. After that, you just have to load the<br>
&gt; topology abd the coordinates:<br>
&gt;<br>
&gt;    vmd -psf topology.psf -pdb coordinate.pdb<br>
&gt;<br>
&gt; You won&#39;t see any weird bonds.<br>
&gt;<br>
&gt; Nicolas<br>
&gt;&gt;<br>
&gt;&gt; Catch ya,<br>
&gt;&gt;<br>
&gt;&gt; Dr. Dallas Warren<br>
&gt;&gt; Drug Delivery, Disposition and Dynamics<br>
&gt;&gt; Monash Institute of Pharmaceutical Sciences, Monash University<br>
&gt;&gt; 381 Royal Parade, Parkville VIC 3010<br>
&gt;&gt; <a href="mailto:dallas.warren@pharm.monash.edu.au">dallas.warren@pharm.monash.edu.au</a><br>
&gt;&gt; +61 3 9903 9167<br>
&gt;&gt; ---------------------------------<br>
&gt;&gt; When the only tool you own is a hammer, every problem begins to<br>
&gt;&gt; resemble a nail.<br>
&gt;&gt;<br>
&gt;&gt; *From:* <a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a><br>
&gt;&gt; [mailto:<a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a>] *On Behalf Of *Gunnar<br>
&gt;&gt; Widtfeldt Reginsson<br>
&gt;&gt; *Sent:* Thursday, 12 November 2009 10:13 AM<br>
&gt;&gt; *To:* <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;&gt; *Subject:* [gmx-users] Docking with PyMol and using Gromacs<br>
&gt;&gt;<br>
&gt;&gt; Hi.<br>
&gt;&gt;<br>
&gt;&gt; I am a new user of Gromacs.<br>
&gt;&gt;<br>
&gt;&gt; My question is both PyMol and Gromacs related.<br>
&gt;&gt;<br>
&gt;&gt; I tried the PyMol users mailing list but couldn&#39;t find anything.<br>
&gt;&gt;<br>
&gt;&gt; I have a small organic molecule that I am inserting into DNA in pymol.<br>
&gt;&gt;<br>
&gt;&gt; I have the DNA as one pdb file and the organic molecule as another<br>
&gt;&gt; pdb file. I open the DNA file in pymol and then load the organic<br>
&gt;&gt; molecule. After docking the organic molecule I write &quot;save name.pdb&quot;<br>
&gt;&gt;<br>
&gt;&gt; When viewing the name.pdb file in pymol there are some bonds between<br>
&gt;&gt; the organic molecule and the DNA that I don&#39;t want. Somehow pymol<br>
&gt;&gt; creates them and I don&#39;t see those bonds in the name.pdb file when I<br>
&gt;&gt; open it in a text reader.<br>
&gt;&gt;<br>
&gt;&gt; I then create a .gro file with pdb2gmx of the DNA.pdb and a .gro file<br>
&gt;&gt; of the organic molecule with topolbuilder 1.2 , and unite those .gro<br>
&gt;&gt; files and convert into a .pdb with editconf<br>
&gt;&gt;<br>
&gt;&gt; The newly created pdb file still has those unwanted bonds.<br>
&gt;&gt;<br>
&gt;&gt; My question is:<br>
&gt;&gt;<br>
&gt;&gt; Can I ignore those bonds?<br>
&gt;&gt;<br>
&gt;&gt; If not, how can I prevent pymol making those bonds?<br>
&gt;&gt;<br>
&gt;&gt; Thanks.<br>
&gt;&gt;<br>
<br>
</div></div><div><div></div><div class="h5">--<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/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></div>