I checked the docked structure and the structure obtained after adding the ligand coordinates to the processed file obtained after using pdb2gmx command. It's very surprising that in the docked structure ligand is at the correct place but the processed gromacs file the ligand lies far apart from the protein. Any clue what could be the reason for this ??<br>
<br><div class="gmail_quote">On Fri, Dec 2, 2011 at 11:00 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="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><br>
<br>
bharat gupta wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Actually I check the file newbox.gro in VMD and I found that the phosphate ion instead of being docked to my protein lies somewhere far away from the protein. The docked complex was taken from autodock's docking result. So what could have wrong . I guess this could be the reason for the box being too large .<br>
<br>
</blockquote>
<br></div>
Sounds like you chose Autodock's undocked complex where the protein and ligand are far apart, but that's more of an Autodock question rather than a Gromacs one.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
On Fri, Dec 2, 2011 at 10:54 AM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br>
<br>
<br>
bharat gupta wrote:<br>
<br>
Hi,<br>
<br>
I am trying the simulation of a docked complex of my protein .<br>
While solvating the box using the following command :-<br>
genbox -cp newbox.gro -cs spc216.gro -p topol.top -o solv.gro<br>
<br>
<br>
The processing does not stop and continues to run . Here's the<br>
output that I got while solvating the box , which is still<br>
continuing :-<br>
<br>
<br>
Reading solute configuration<br>
GRoups of Organic Molecules in ACtion for Science<br>
Containing 3620 atoms in 231 residues<br>
Initialising van der waals distances...<br>
<br>
WARNING: masses and atomic (Van der Waals) radii will be determined<br>
based on residue and atom names. These numbers can deviate<br>
from the correct mass and radius of the atom type.<br>
<br>
Reading solvent configuration<br></div>
"216H2O,WATJP01,SPC216,SPC-__<u></u>MODEL,300K,BOX(M)=1.86206NM,__<u></u>WFVG,MAR.<div class="im"><br>
1984"<br>
solvent configuration contains 648 atoms in 216 residues<br>
<br>
Initialising van der waals distances...<br>
Will generate new solvent configuration of 46x46x46 boxes<br>
Generating configuration<br>
Sorting configuration<br>
Found 1 molecule type:<br>
SOL ( 3 atoms): 21024576 residues<br>
Calculating Overlap...<br>
box_margin = 0.315<br>
Removed 1047843 atoms that were outside the box<br>
<br>
What have gone wrong for solvation to take this much time on 12<br>
processors CPU.<br>
<br>
<br>
genbox is not parallelized; you can only use 1 CPU.<br>
<br>
You've added 21 million water molecules (yikes!) to the box and<br>
genbox is trying to remove over 1 million atoms - I'd say your box<br>
is much too large for a solute of 3620 atoms, by several orders of<br>
magnitude. You're probably running out of memory to do this operation.<br>
<br>
-Justin<br>
<br></div>
-- ==============================<u></u>__==========<div class="im"><br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br></div>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>__<a href="http://vt.edu/Pages/Personal/justin" target="_blank"><u></u>vt.edu/Pages/Personal/justin</a><br>
<<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a>><br>
<br>
==============================<u></u>__==========<div class="im"><br>
<br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<a href="http://lists.gromacs.org/__mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/__<u></u>mailman/listinfo/gmx-users</a><div class="im"><br>
<<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a>><br>
Please search the archive at<br></div>
<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a><div class="im"><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a>> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>>.<br>
Can't post? Read <a href="http://www.gromacs.org/__Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists</a><div class="im"><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a>><br>
<br>
<br>
<br>
<br>
-- <br>
Bharat<br>
Ph.D. Candidate<br>
Room No. : 7202A, 2nd Floor<br>
Biomolecular Engineering Laboratory<br>
Division of Chemical Engineering and Polymer Science<br>
Pusan National University<br>
Busan -609735<br>
South Korea<br>
Lab phone no. - +82-51-510-3680, +82-51-583-8343<br>
Mobile no. - 010-5818-3680<br></div>
E-mail : <a href="mailto:monu46010@yahoo.com" target="_blank">monu46010@yahoo.com</a> <mailto:<a href="mailto:monu46010@yahoo.com" target="_blank">monu46010@yahoo.com</a>><br>
<br>
</blockquote><div class="HOEnZb"><div class="h5">
<br>
-- <br>
==============================<u></u>==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<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.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>
==============================<u></u>==========<br>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/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/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Bharat<br>Ph.D. Candidate<br>Room No. : 7202A, 2nd Floor<br>Biomolecular Engineering Laboratory<br>Division of Chemical Engineering and Polymer Science<br>
Pusan National University<br>Busan -609735<br>South Korea<br>Lab phone no. - +82-51-510-3680, +82-51-583-8343<div>Mobile no. - 010-5818-3680<br>E-mail : <a href="mailto:monu46010@yahoo.com" target="_blank">monu46010@yahoo.com</a></div>
<br>