<font size="4"><font face="verdana,sans-serif">Hi Eric,<br>
<br>
I tried the dockprep.py file you pointed to me a few emails ago. It
seems like its exactly what I am looking for. The only problem is I get
an error at the moment; I also have two other questions about it.</font></font><br><br><font size="4"><font face="verdana,sans-serif">ep.py file you pointed to me a few emails ago. It
seems like its exactly what I am looking for. The only problem is I get
an error at the moment; I also have two other questions about it.<br><br>Here is the command I used:<br><br>./chimera --nogui pdb:2h4f $DESKTOP/dockprep.py<br><br>The program then returns a bunch of text on what its doing. When adding the charges to the ligands with the AM1-BCC method is when the error pops up. Here is the (slightly abridged) error message, plus a few lines back:<br>
<br><font size="4"><b><span style="font-family:courier new,monospace">Running ANTECHAMBER for residue NAD</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">Running ANTECHAMBER command: /home/abhay/Desktop/rcg-research/personal-work/CHIMERA/bin/amber11/bin/antechamber -ek qm_theory='AM1', grms_tol=0.02, tight_p_conv=0, scfconv=1.d-10, -i /tmp/tmplP08co/ante.in.mol2 -fi mol2 -o /tmp/tmplP08co/ante.out.mol2 -fo mol2 -c bcc -nc -1 -j 5 -s 2</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace"></span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">(NAD) Total number of electrons: 344; net charge: -1</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace"></span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">(NAD) Running: /home/abhay/Desktop/rcg-research/personal-work/CHIMERA/bin/amber11/bin/sqm -O -i <a href="http://sqm.in">sqm.in</a> -o sqm.out</span><br style="font-family:courier new,monospace">
<br style="font-family:courier new,monospace;color:rgb(0,0,0)"><span style="font-family:courier new,monospace;color:rgb(0,0,0)">(NAD) Error: cannot run "/home/abhay/Desktop/rcg-research/personal-work/CHIMERA/bin/amber11/bin/sqm -O -i <a href="http://sqm.in">sqm.in</a> -o sqm.out" of bcc() in charge.c properly, exit</span><br style="font-family:courier new,monospace">
<br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">(NAD) Running: /home/abhay/Desktop/rcg-research/personal-work/CHIMERA/bin/amber11/bin/bondtype -j part -i ANTECHAMBER_BOND_TYPE.AC0 -o <a href="http://ANTECHAMBER_BOND_TYPE.AC">ANTECHAMBER_BOND_TYPE.AC</a> -f ac</span><span style="font-family:courier new,monospace"></span><br style="font-family:courier new,monospace">
<br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">(NAD) Running: /home/abhay/Desktop/rcg-research/personal-work/CHIMERA/bin/amber11/bin/atomtype -i ANTECHAMBER_AC.AC0 -o <a href="http://ANTECHAMBER_AC.AC">ANTECHAMBER_AC.AC</a> -p gaff</span><br style="font-family:courier new,monospace">
<br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">Error while processing /home/abhay/Desktop/dockprep.py:</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">Failure running ANTECHAMBER for residue NAD</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">Check reply log for details</span></b></font><br><br>I am not sure what the problem is. $AMBERHOME and $DOCK_HOME are not defined, but that only seemed to cause problems before. Perhaps they need to be defined now?<br>
I would give you the reply log, but I don't know how to find it. I tried using:<br><br></font></font><font size="4"><font face="verdana,sans-serif">$ ./chimera --start Reply_Log --nogui pdb:2h4f $DESKTOP/dockprep.py<br>
<br>I also tried reply_log, Reply_log, reply_Log, ReplyLog, replyLog, Replylog, and transposing --start and --nogui. None of those worked either. The program immediately returns:<br><br><b style="font-family:courier new,monospace">Starting extension 'Reply_Log' failed</b></font></font><font size="4"><font face="verdana,sans-serif"> (or whatever the pertinent permutation of "reply_log" happens to be.)</font></font><br>
<br><font size="4"><font face="verdana,sans-serif">Assuming this can be fixed easily, I have one more question: this appears to only prepare the receptor for docking. What about prepping the ligand? Does that need its own new script?<br>
<br>Thank you so much for your help,<br><br>-Brad<br></font></font><br><div class="gmail_quote">On Wed, Feb 29, 2012 at 7:08 PM, Eric Pettersen <span dir="ltr"><<a href="mailto:pett@cgl.ucsf.edu">pett@cgl.ucsf.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">This previous posting to the dock-fans list seems to be more directly pertinent:<div>
<br></div><div><a href="http://mailman.docking.org/pipermail/dock-fans/2007-May/001043.html" target="_blank">http://mailman.docking.org/pipermail/dock-fans/2007-May/001043.html</a></div><div><br></div><div>--Eric</div><div>
<br><div><span style="border-collapse:separate;color:rgb(0,0,0);font-family:Helvetica;font-size:16px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div style="word-wrap:break-word">
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px"><font style="font:normal normal normal 16px/normal Helvetica" size="5" face="Helvetica"> Eric Pettersen</font></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">
<font style="font:normal normal normal 16px/normal Helvetica" size="5" face="Helvetica"> UCSF Computer Graphics Lab</font></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">
<font style="font:normal normal normal 16px/normal Helvetica" size="5" face="Helvetica"> <a href="http://www.cgl.ucsf.edu" target="_blank">http://www.cgl.ucsf.edu</a></font></div><div><font style="font:normal normal normal 16px/normal Helvetica" size="5" face="Helvetica"><br>
</font></div></div></span></div><div><div><div></div><div class="h5"><div>On Feb 29, 2012, at 2:47 PM, Trent E. Balius wrote:</div><br></div></div><blockquote type="cite"><div><div></div><div class="h5"><font color="black" face="Arial, Helvetica, sans-serif"><font color="black" face="Arial, Helvetica, sans-serif"> <div>
Hi Brad, <br> <br> I am cc'ing this to the chimera list. <br> <br> I think that you should write a python script to run DOCK prep on the command line. See the following:<br> <br> <a href="http://www.cgl.ucsf.edu/pipermail/chimera-dev/2011/000784.html" target="_blank">http://www.cgl.ucsf.edu/pipermail/chimera-dev/2011/000784.html</a><br>
<br> For other example scripts see the following:<br> <br> <a href="http://www.cgl.ucsf.edu/pipermail/chimera-users/2009-February/003539.html" target="_blank">http://www.cgl.ucsf.edu/pipermail/chimera-users/2009-February/003539.html</a><br>
</div> <div> <br> <font color="black" face="Arial, Helvetica, sans-serif">You should not have to specify the AMBERHOME path. This is done internally to chimera.</font><br> <br> <font color="black" face="Arial, Helvetica, sans-serif">as an aside, It is unnecessary to run antichamber with sudo.</font><br>
<br> I hope this helps,<br> <br> </div> <div style="clear:both">Trent E. Balius<br> Graduate Student, Rizzo Group,<br> Department of Applied Mathematics and Statistics,<br> Stony Brook University.<br> Office: Math Tower 3-129, Phone: <a href="tel:%28631%29%20632-8519" value="+16316328519" target="_blank">(631) 632-8519</a><br>
<font><font face="Arial, Helvetica, sans-serif">URL: <a href="http://www.ams.sunysb.edu/%7Etbalius/" target="_blank">http://www.ams.sunysb.edu/~tbalius/</a></font></font><br> </div> <div> <br> </div> <div> <br> </div> <div style="font-family:arial,helvetica;font-size:10pt;color:black">
-----Original Message-----<br> From: Brad Ridder <<a href="mailto:clustro@gmail.com" target="_blank">clustro@gmail.com</a>><br> To: dock-fans <<a href="mailto:dock-fans@docking.org" target="_blank">dock-fans@docking.org</a>><br>
Sent: Wed, Feb 29, 2012 4:38 pm<br> Subject: [Dock-fans] Help with antechamber<br> <br> <div> <font size="4"><font face="verdana,sans-serif">Hi Dockfans,<br> <br> My name is Brad; I am a grad student in the dept. of chemical engineering at Purdue University.<br>
<br> We are experimenting with DOCK for various drug discovery work. At the moment, I am trying to do execute all of the "dock prep" commands through the Terminal, since it is not practical to do large volumes of Dock Prep by hand through the CHIMERA GUI.<br>
<br> Corresponding to what is in the DOCK tutorial, I am on "Receptor and Ligand Preparation." I got reduce to work, in order to add hydrogens.<br> <br> Now, I am trying to get antechamber to work, in order to add the charges, but it is not working. Here is the command I am running:<br>
<br> <b><span style="font-family:courier new,monospace">sudo ./antechamber -i $DESKTOP/3GLR.pdb -fi pdb -o $DESKTOP/3GLRcharged.mol2 -fo mol2</span></b><br> <br> Which returns the error:<br> <b><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">Warning: $AMBERHOME and $DOCK_HOME enviornment strings are not set, use "mopac.sh" in the work directory</span><br style="font-family:courier new,monospace"><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">Cannot open CONNECT.TPL , exit</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">The atom number exceeds the MAXATOM, reallocate memory</span></b><br>
<br> I looked around in the dock6 directory, and found a folder called "parameters/antechamber", which contains CONNECT.TPL. I figured this is what should be the $AMBERHOME directory.<br> For DOCK_HOME, I figured it should be the dock6 folder.<br>
<br> Here is what I have in my .bash_profile file in my home folder.<br> <br> <b style="font-family:courier new,monospace">export AMBERHOME='/home/.../DOCK6.5/dock6/parameters/antechamber'<br> export DOCK_HOME='/home/.../DOCK6.5/dock6'</b><br>
<br> I also tried using the "amber11" folder that comes with CHIMERA as AMBERHOME, but to no effect.<br> <br> I don't actually have AMBER installed on my computer, but since CHIMERA can add charges without actually requiring the purchase of AMBER, I do not think that is the problem.<br>
<br> I did not know what to do with "mopac.sh." I ran the script, but it says it needs input files.<br> <br> What am I doing wrong that is preventing antechamber from executing correctly?<br clear="all"></font></font><br>
-- <br> Bradley James Ridder<br> Chakrabarti Group<br> Graduate Student<br> School of Chemical Engineering<br> Purdue University<br> <br> <br> </div> <div style="margin:0px;font-family:Tahoma,Verdana,Arial,Sans-Serif;font-size:12px;color:#000;background-color:#fff">
<pre style="font-size:9pt"><tt>_______________________________________________
Dock-fans mailing list
<a>Dock-fans@docking.org</a>
<a>http://mailman.docking.org/mailman/listinfo/dock-fans</a>
</tt></pre> </div> </div> </font></font></div></div>_______________________________________________<br>Chimera-users mailing list<br><a href="mailto:Chimera-users@cgl.ucsf.edu" target="_blank">Chimera-users@cgl.ucsf.edu</a><br>
<a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" target="_blank">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br></blockquote></div><br></div><div><span style="border-collapse:separate;color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div style="word-wrap:break-word">
<br></div></span> </div><br></div></blockquote></div><br><br clear="all"><br>-- <br>Bradley James Ridder<br>Chakrabarti Group<br>Graduate Student<br>School of Chemical Engineering<br>Purdue University<br><br><br>