Hi Tsjerk,<br>Thanks for your response, <br>I am finding problem in using original coordinates with the PRODRG topology because of mismatch in number of Hydrogen atoms in original coordinate file and PRODRG topology.<br>Can you tell me how can I try 2nd option you suggested i.e. Fit the structure obtained from PRODRG to original structure?<br>
<br>Thanks &amp; Regards,<br>Vivek<br><br><div class="gmail_quote">On 27 January 2010 10:35, Tsjerk Wassenaar <span dir="ltr">&lt;<a href="mailto:tsjerkw@gmail.com">tsjerkw@gmail.com</a>&gt;</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;">
Hi Vivek,<br>
<br>
Either<br>
<br>
1. Use the original ligand coordinates with the PRODRG topology<br>
<br>
or<br>
<br>
2. Fit the structure obtained from PRODRG to the original structure<br>
<br>
Cheers,<br>
<font color="#888888"><br>
Tsjerk<br>
</font><div><div></div><div class="h5"><br>
On Wed, Jan 27, 2010 at 5:58 AM, vivek sharma<br>
&lt;<a href="mailto:viveksharma.iitb@gmail.com">viveksharma.iitb@gmail.com</a>&gt; wrote:<br>
&gt; Hi Dallas,<br>
&gt; I am trying to run MD simulation over a docked complex (protein+ligand), to<br>
&gt; confirm their dynamic stability in water media.<br>
&gt; For the same I am using PRODRG server to generate topologies for ligand<br>
&gt; molecule as gromacs can generate topology for 20 standard residues. As<br>
&gt; mentioned in tutorial for drug-enzyme complex, I am editing the .top and<br>
&gt; .gro files to include the PRODRG generated files (DRGGMX.ITP in .top and<br>
&gt; DRGAPH.GRO in .gro file).<br>
&gt; I observe that their are changes in co-ordinate of ligand after processing<br>
&gt; them through PRODRG server. So these new co-ordinates for ligand are placing<br>
&gt; ligand away from the protein while the ligand molecule was in protein pocket<br>
&gt; in original docked complex.<br>
&gt;<br>
&gt; I hope it gives what I am trying to do, and where I am getting stuck.<br>
&gt;<br>
&gt; I am looking for some suggestions and more insight of the problem to solve<br>
&gt; it.<br>
&gt; Earlier I have done same procedure successfully for a different docked<br>
&gt; complex.<br>
&gt;<br>
&gt; Regards,<br>
&gt; Vivek<br>
&gt;<br>
&gt; 2010/1/27 Dallas B. Warren &lt;<a href="mailto:Dallas.Warren@pharm.monash.edu.au">Dallas.Warren@pharm.monash.edu.au</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; So, what EXACTLY are you doing?<br>
&gt;&gt;<br>
&gt;&gt;<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 resemble a<br>
&gt;&gt; nail.<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; From: <a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a> [mailto:<a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a>]<br>
&gt;&gt; On Behalf Of vivek sharma<br>
&gt;&gt; Sent: Monday, 25 January 2010 7:38 PM<br>
&gt;&gt; To: Discussion list for GROMACS users<br>
&gt;&gt; Subject: Re: [gmx-users] Ligand coming out while trying Drug-enzyme<br>
&gt;&gt; tutorial<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; HI Tsjerk,<br>
&gt;&gt;<br>
&gt;&gt; Thanks for your reply. But, I can&#39;t see if it is going suddenly or<br>
&gt;&gt; gradually.<br>
&gt;&gt; What i can see is the ligand is away from the molecule after editing the<br>
&gt;&gt; gro file with PRODRG output.<br>
&gt;&gt;<br>
&gt;&gt; It seems liek PRODRG has modified the co-ordinates that places ligand away<br>
&gt;&gt; from the protein.<br>
&gt;&gt;<br>
&gt;&gt; ~Vivek<br>
&gt;&gt;<br>
&gt;&gt; 2010/1/25 Tsjerk Wassenaar &lt;<a href="mailto:tsjerkw@gmail.com">tsjerkw@gmail.com</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; Hi Vivek,<br>
&gt;&gt;<br>
&gt;&gt; &gt; Now when I am processing the modified .gro file to generate box, the<br>
&gt;&gt; &gt; ligand<br>
&gt;&gt; &gt; and cofactor are going away from the protein molecule and I am not able<br>
&gt;&gt; &gt; to<br>
&gt;&gt; &gt; analyze the complex.<br>
&gt;&gt;<br>
&gt;&gt; Gradually going away, or suddenly jumping?<br>
&gt;&gt;<br>
&gt;&gt; In the latter case, read up on periodic boundary conditions.<br>
&gt;&gt;<br>
&gt;&gt; Tsjerk<br>
&gt;&gt;<br>
&gt;&gt; --<br>
&gt;&gt; Tsjerk A. Wassenaar, Ph.D.<br>
&gt;&gt;<br>
&gt;&gt; Computational Chemist<br>
&gt;&gt; Medicinal Chemist<br>
&gt;&gt; Neuropharmacologist<br>
&gt;&gt; --<br>
&gt;&gt;<br>
&gt;&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt;&gt; Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
&gt;&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt;&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt;&gt; 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;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; --<br>
&gt;&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt;&gt; Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
&gt;&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt;&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt;&gt; 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;<br>
&gt;<br>
&gt; --<br>
&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt; Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt; 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;<br>
<br>
<br>
<br>
--<br>
Tsjerk A. Wassenaar, Ph.D.<br>
<br>
Computational Chemist<br>
Medicinal Chemist<br>
Neuropharmacologist<br>
--<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>