bioexcel / biobb_chemistry

Biobb_chemistry is the Biobb module collection to perform chemistry over molecular dynamics simulations.
https://mmb.irbbarcelona.org/biobb/
Apache License 2.0
4 stars 1 forks source link

adding hydrogens ignores the pH option #6

Closed amjjbonvin closed 3 months ago

amjjbonvin commented 9 months ago

Describe the bug

When running the ligand parametrisation example and part adding the hydrogens, obabel reports that it will ignore the defined pH because both the -p and -p options are given. Extract from a log file:

*** Open Babel Warning  in DoTransformations
  Both -p and -h options are set. All implicit hydrogens (-h) will be added without considering pH.

To Reproduce

Run the ligand parametrisation notebook, and check the log file generated by the Add Hydrogens step

Expected behavior

the -h and -p options should be exclusive, allowing to control the pH to define the protonation state.

Desktop (please complete the following information):

amjjbonvin commented 9 months ago

Tested manually with obabel that when using only the -p options hydrogens are added considering the defined pH. With both -h and -p all hydrogens are added by default, irrespective of the defined pH

gbayarri commented 9 months ago

Fixed in PR https://github.com/bioexcel/biobb_chemistry/pull/7, it will be available in the next release.

Thanks!

github-actions[bot] commented 4 months ago

This issue is stale because it has been open for 30 days with no activity. If there is no more activity in the following 14 days it will be automatically closed.

github-actions[bot] commented 3 months ago

This issue was closed because it has been inactive for 14 days since being marked as stale.

amjjbonvin commented 3 months ago

This strange to close an issue as completed when nothing has happened… it is not because it was stale that it is completed.

gbayarri commented 3 months ago

Hi, @amjjbonvin , it's an automatic action. It's reopened until next release.

Best.

amjjbonvin commented 3 months ago

Hi Genis

Thanks. Then you could close it as “not planned” or something else, but not as “completed” (might be easier to extract stats of what has been done) :-)

Cheers Alexandre

gbayarri commented 3 months ago

Fixed in the new 4.2.0 version