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

Error to process neutral Cter of ASN #386

Open
micmar15sr opened this issue Jan 24, 2024 · 2 comments
Open

Error to process neutral Cter of ASN #386

micmar15sr opened this issue Jan 24, 2024 · 2 comments

Comments

@micmar15sr
Copy link

When i process a pdb file obtained from a protein modelled with Charmm36m topology extracted from a trajectory by the gmx trjconv software I obtain the following error:

ValueError: Biomolecular structure is incomplete:  Found gap in biomolecule structure for atom ATOM   2654  HT2 ASN   166      72.292  57.589  76.447  0.0000 0.0000

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/home/marcello/anaconda3/envs/propka/lib/python3.11/site-packages/pdb2pqr/main.py", line 792, in main_driver
    results = non_trivial(
              ^^^^^^^^^^^^
  File "/home/marcello/anaconda3/envs/propka/lib/python3.11/site-packages/pdb2pqr/main.py", line 631, in non_trivial
    raise ValueError(err)
ValueError: Unable to debump biomolecule. Biomolecular structure is incomplete:  Found gap in biomolecule structure for atom ATOM   2654  HT2 ASN   166      72.292  57.589  76.447  0.0000 0.0000

The above exception was the direct cause of the following exception:

Traceback (most recent call last):
  File "/home/marcello/anaconda3/envs/propka/bin/pdb2pqr30", line 10, in <module>
    sys.exit(main())
             ^^^^^^
  File "/home/marcello/anaconda3/envs/propka/lib/python3.11/site-packages/pdb2pqr/main.py", line 829, in main
    if main_driver(args) == 1:
       ^^^^^^^^^^^^^^^^^
  File "/home/marcello/anaconda3/envs/propka/lib/python3.11/site-packages/pdb2pqr/main.py", line 802, in main_driver
    raise RuntimeError from err
RuntimeError

It seems that renaming the "HT2" atom to "H ", make the program work (properly??). I saw older issue that are related to the *.xml files, maybe there are some changes needed but i'm not able to undertand if is just related to add another <altname> or create another group.

@sobolevnrm
Copy link
Member

Thank you for bringing this up. I need to fix #367 before I can address this problem.

@micmar15sr
Copy link
Author

Ok, thank you, does modifying the PDB as I did make it work properly? Is there a guide to understand and maybe try to modify the xml file on my own properly and then suggest a commit? Could you suggest some clue on the file to modify?

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