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><div class="gmail_quote">On Fri, Dec 2, 2011 at 10:54 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="HOEnZb"><div class="h5"><br>
<br>
bharat gupta wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
<br>
I am trying the simulation of a docked complex of my protein . 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 output that I got while solvating the box , which is still 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>
"216H2O,WATJP01,SPC216,SPC-<u></u>MODEL,300K,BOX(M)=1.86206NM,<u></u>WFVG,MAR. 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 processors CPU.<br>
</blockquote>
<br></div></div>
genbox is not parallelized; you can only use 1 CPU.<br>
<br>
You've added 21 million water molecules (yikes!) to the box and genbox is trying to remove over 1 million atoms - I'd say your box is much too large for a solute of 3620 atoms, by several orders of magnitude. You're probably running out of memory to do this operation.<br>
<br>
-Justin<br>
<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>==========<div class="HOEnZb"><div class="h5"><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>