[Biopython] Problem with ResidueDepth class

Carlos Henrique carlosmr12 at gmail.com
Mon Nov 9 22:23:12 UTC 2015


Hi Lenna,
Thank you so much for your quick answer!

And that was exactly the problem here. I extracted only the *msms*,
*pdb_to_xyrz* and *pdb_to_xyrzn* files to my* /usr/local/bin*.

I just moved the *atmtypenumbers* and changed the path in *pdb_to_xyrn* and
worked like a charm.
:)
Thanks once more.

Cheers,

Carlos

2015-11-09 20:03 GMT-02:00 Lenna Peterson <arklenna at gmail.com>:

> Hi Carlos,
>
> I am able to get ResidueDepth and surface from 1yrk (although I notice
> neither function is covered by the test suite, perhaps due to the soft
> dependencies).
>
> The error means that pdb_to_xyzr is not able to find an entry for a water
> oxygen in the file "atmtypenumbers".
>
> In pdb_to_xyzr, the path to atmtypenumbers is specified relatively, so my
> suspicion is that you're not running it from the directory where
> atmtypenumbers is located and it's not being found at all. If that is the
> case, I would recommend editing pdb_to_xyzr to specify an absolute path
> (e.g. change the line 'numfile = "./atmtypenumbers"' to something like
> 'numfile = /home/carlos/...').
>
> Please let me know if that fixes your problem.
>
> Cheers,
>
> Lenna
>
> On Mon, Nov 9, 2015 at 3:58 PM, Carlos Henrique <carlosmr12 at gmail.com>
> wrote:
>
>> Hello everyone,
>> I have been struggling with the calculation of Residue depth using the
>> class ResidueDepth as mentioned in [1
>> <http://biopython.org/wiki/The_Biopython_Structural_Bioinformatics_FAQ#How_do_I_calculate_residue_depth.3F>
>> ].
>> I already set everything up and installed the msms from [2
>> <http://mgltools.scripps.edu/downloads#msms>].
>> However, I cannot make it run as it should. At first I thought it could
>> be a problem with my .pdb files, but I tested with 3 different files
>> downloaded directly from the protein data bank (1YRK, 1B4F, 1U46) and I
>> keep getting the same error message which is something like this:
>>
>> A bunch of this:
>> pdb_to_xyzr: error, file 1YRK.pdb line 1687 residue   37 atom pattern HOH
>> O was not found in ./atmtypenumbers
>> pdb_to_xyzr: error, file 1YRK.pdb line 1688 residue   38 atom pattern HOH
>> O was not found in ./atmtypenumbers
>> pdb_to_xyzr: error, file 1YRK.pdb line 1689 residue   39 atom pattern HOH
>> O was not found in ./atmtypenumbers
>> pdb_to_xyzr: error, file 1YRK.pdb line 1690 residue   40 atom pattern HOH
>> O was not found in ./atmtypenumbers
>> pdb_to_xyzr: error, file 1YRK.pdb line 1691 residue   41 atom pattern HOH
>> O was not found in ./atmtypenumbers
>> pdb_to_xyzr: error, file 1YRK.pdb line 1692 residue   42 atom pattern HOH
>> O was not found in ./atmtypenumbers
>>
>> and the traceback error message right after:
>> Traceback (most recent call last):
>>   File "generate_measurements.py", line 82, in <module>
>>     rd = ResidueDepth(struct[0], '1YRK.pdb')
>>   File
>> "/home/carlos/.virtualenvs/pdb/local/lib/python2.7/site-packages/Bio/PDB/ResidueDepth.py",
>> line 144, in __init__
>>     surface=get_surface(pdb_file)
>>   File
>> "/home/carlos/.virtualenvs/pdb/local/lib/python2.7/site-packages/Bio/PDB/ResidueDepth.py",
>> line 90, in get_surface
>>     "Failed to generate surface file using command:\n%s" % make_surface
>> AssertionError: Failed to generate surface file using command:
>> msms -probe_radius 1.5 -if /tmp/tmpwodW16 -of /tmp/tmpCbMTwN >
>> /tmp/tmprkvaQv
>>
>> I have been trying to solve this since yesterday in vain. Does anyone
>> have any clue or maybe I can get the residue depth using something else?
>>
>> Cheers,
>>
>> [1]
>> http://biopython.org/wiki/The_Biopython_Structural_Bioinformatics_FAQ#How_do_I_calculate_residue_depth.3F
>> [2] http://mgltools.scripps.edu/downloads#msms
>>
>> --
>> Carlos Rodrigues
>>
>>
>> _______________________________________________
>> Biopython mailing list  -  Biopython at mailman.open-bio.org
>> http://mailman.open-bio.org/mailman/listinfo/biopython
>>
>
>


-- 
Carlos Rodrigues
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.open-bio.org/pipermail/biopython/attachments/20151109/ab953294/attachment-0001.html>


More information about the Biopython mailing list