Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

ligpargen running problem #2

Closed
zhubd11 opened this issue Dec 29, 2021 · 1 comment
Closed

ligpargen running problem #2

zhubd11 opened this issue Dec 29, 2021 · 1 comment

Comments

@zhubd11
Copy link

zhubd11 commented Dec 29, 2021

I tried with BOSS 4.9 and ligpargen compiled with conda python3.7, but when I run it, I get this error, below are the error output:

ligpargen -i part2.pdb -n phenol -p phenol -r MOL -c 0 -o 0 -cgen CM1A
[16:59:22] Explicit valence for atom # 25 O, 5, is greater than permitted
Traceback (most recent call last):
File "/public1/home/scb1100/miniconda3/envs/test/bin/ligpargen", line 33, in
sys.exit(load_entry_point('LigPargen', 'console_scripts', 'ligpargen')())
File "/public1/home/scb1100/ligpargen/ligpargen-main/ligpargen/ligpargen.py", line 494, in main
molname=args.molname, workdir= args.path, debug= args.debug)
File "", line 17, in init
File "/public1/home/scb1100/ligpargen/ligpargen-main/ligpargen/ligpargen.py", line 174, in post_init
moleculeRDkit, newIndexToOriginalIndex, atomsNameOriginal, residueNameOriginal = utilities.generateRDkitMolecule(self.ifile, self.smile, self.workdir, molnameA, self.debug)
File "/public1/home/scb1100/ligpargen/ligpargen-main/ligpargen/tools/utilities.py", line 165, in generateRDkitMolecule
atomsIndexLstWithRightOrder = buildProperOrderAtomIndexLst(molecule, molname, workdir)
File "/public1/home/scb1100/ligpargen/ligpargen-main/ligpargen/tools/utilities.py", line 513, in buildProperOrderAtomIndexLst
heavyAtomsIndexLst = [atom.GetIdx() for atom in molecule.GetAtoms() if atom.GetSymbol() != 'H']
AttributeError: 'NoneType' object has no attribute 'GetAtoms'

Any ideas?

@Isra3l
Copy link
Owner

Isra3l commented Feb 1, 2022

Hey zhubd11!

By default, BOSS checks the molecule integrity as a sanity check of the inputs. Probably your input PDB (part2.pdb) has something wrong (atoms too close, wrong charge, ...). You should take a look in pymol/chimera/VMD/... and see where is the problem.

Good luck!

Israel

@Isra3l Isra3l closed this as completed Feb 1, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants