[gmx-users] bad box error

amna khan amnakhan369 at gmail.com
Tue Jan 29 22:04:11 CET 2013


yes i have following this tutorial , the only difference i made is that i
didn't change the charges for ligand file.

it seems it is necessary to do this ?
am i rite ?


On Wed, Jan 30, 2013 at 1:51 AM, Justin Lemkul <jalemkul at vt.edu> wrote:

>
>
> On 1/29/13 3:46 PM, amna khan wrote:
>
>> i am using the best conformation of protein ligand file obtained after
>> docking results
>>
>> i made the .gro file  by
>>   pdb2gmx -f pro.pdb -o proceed.gro -water spc -missing -ignh
>>
>> then the ligand codinates from ligand,gro was obtained from prodrug sever,
>>
>
> Beware the quality of PRODRG topologies - they are never sufficiently
> accurate without manual reparameterization.  Using the coordinates from
> PRODRG may or may not cause problems.  Make sure you don't allow PRODRG to
> run EM on the structure, or else the docked pose will be altered.
>
>
>  then i pasted ligand coordinates in this file
>> the first three lines are
>>
>> Gallium Rubidium Oxygen Manganese Argon Carbon Silicon
>>   2059
>>    705SER      N    1  -3.373  -1.131   1.295
>>    705SER     H1    2  -3.442  -1.166   1.232
>>    705SER     H2    3  -3.339  -1.043   1.261
>>    705SER     H3    4  -3.298  -1.196   1.303
>>
>> last three lines are
>>      1LIG  OAR     33 0.425 -0.911   1.194
>>      1LIG  CAC     34   0.593  -1.073   1.167
>>      1LIG  CAT     35   0.500  -1.191   1.198
>>     9.21266   2.49507   2.95990
>>
>>
>>
>>
>> and the output for command is
>> wc -l proceed.gro
>> 2100 proceed.gro
>>
>>
> This shows your problem.  You didn't account for the fact that you pasted
> in a ligand by properly incrementing the atom count on the second line of
> the file. Thus, when editconf reaches atom 2059, it expects to find box
> vectors on the next line, but instead finds another atom and gives you the
> "bad box" message. Your file actually has 2097 atoms (2100 - 3).
>
> Note that there is a protein-ligand tutorial available if these concepts
> are unfamiliar to you.
>
> http://www.bevanlab.biochem.**vt.edu/Pages/Personal/justin/**
> gmx-tutorials/complex/index.**html<http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/index.html>
>
> -Justin
>
>
> --
> ==============================**==========
>
> Justin A. Lemkul, Ph.D.
> Research Scientist
> Department of Biochemistry
> Virginia Tech
> Blacksburg, VA
> jalemkul[at]vt.edu | (540) 231-9080
> http://www.bevanlab.biochem.**vt.edu/Pages/Personal/justin<http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin>
>
> ==============================**==========
> --
> gmx-users mailing list    gmx-users at gromacs.org
> http://lists.gromacs.org/**mailman/listinfo/gmx-users<http://lists.gromacs.org/mailman/listinfo/gmx-users>
> * Please search the archive at http://www.gromacs.org/**
> Support/Mailing_Lists/Search<http://www.gromacs.org/Support/Mailing_Lists/Search>before posting!
> * Please don't post (un)subscribe requests to the list. Use the www
> interface or send it to gmx-users-request at gromacs.org.
> * Can't post? Read http://www.gromacs.org/**Support/Mailing_Lists<http://www.gromacs.org/Support/Mailing_Lists>
>



More information about the gromacs.org_gmx-users mailing list