Hi all,
I’ve run into 2 issues related to the check_structure
tool.
-
keep_h
, listed here is not actually supported with version biobb version 3.7.6 (my check_structure version is 3.9.11). I get this warning when I try to use it:
/home/pbarletta/anaconda3/envs/locuaz/lib/python3.7/site-packages/biobb_common/generic/biobb_object.py:71: UserWarning: Warning: keep_h is not a recognized property. The most similar property is: ph
error_property, close_property))
And indeed, when I check the docstrings for StrCheckAddHydrogens, it’s not there.
- I can’t understand the selection syntax for the
add_hydrogen
command.
Say I want to add hydrogens to the 3rd residue of the B chain, which is a glutamate.
If I use this properties dictionary:
prop = {
'charges': False,
'mode': 'list',
'list': "GLU3GLU"
}
I get an error:
h_rules[ion_res_list[res]]
KeyError: False
And if I use this one:
prop = {
'charges': False,
'mode': 'list',
'list': "GLU3GLU"
}
it runs, but it modifies every residue. I tried some other variations, sometimes I’d get this:
#ERROR: no mutations available for *:GLUB3
I couldn’t get it working on the terminal either (calling check_structure
directly).
I’d appreciate any help on these issues.
Cheers,