deeprank2 icon indicating copy to clipboard operation
deeprank2 copied to clipboard

Fix hydrogens unkown to the forcefield

Open Max1461 opened this issue 2 years ago • 3 comments

@LilySnow @gcroci2 @DaniBodor I run into a situation that when I use queries.process() to make and save my graph(s) the following type of warning messages are printed for several hydrogens atoms:

is unknown to the forcefield; charge is set to 0.0
is unknown to the forcefield; vanderwaals_parameters set to (0.0, 0.0, 0.0, 0.0)

The test situation in which this occurs is when I use the following pdb file(I have changed the extension otherwise I could not upload it to git) 1AKJ.txt

The hydrogen atoms are then added to this pdb file using:

reduce -FLIP -Quiet $pdbFL > $pdbFL.reduce_flip

Resulting in the following pdb file with hydrogens present: 1AKJ.txt

From this pdb file with added hydrogens a micro-environment is selected and written to a new pdb file, in this case: 1AKJ_1_ILE.txt

This pdb structure is added to queries = QueryCollection()

And the graph is then generated and saved with

queries.process()

During the process step the warning message is printed, which i have added in the following text file: warning_message.txt

This occurs because the hydrogen naming conventions in the pdb that has the hydrogens added through reduce do not match the hydrogen naming conventions present in the forcefield topology file used by DeepRank core.

I have already discussed this issue today during a meeting with @LilySnow and have made this issue after this. Hopefully I have provided enough information through the files, and if not I would gladly help out further.

Max1461 avatar Mar 08 '23 20:03 Max1461

Many thanks @Max1461, we'll look soon further into this :)

gcroci2 avatar Mar 09 '23 12:03 gcroci2

I believe that this will be dealt with in #393

DaniBodor avatar Mar 20 '23 09:03 DaniBodor

This issue is stale because it has been open for 30 days with no activity.

github-actions[bot] avatar Apr 24 '23 03:04 github-actions[bot]