Hello! I have been working with PyRosetta for a while to model protein-ligand interactions. I was getting hydrogen-bond information from other software, but wanted to use PyRosetta's Hbond set for my pose. Rosetta did not find all the expected hydrogen bonds, so I wanted to check what criteria was used by Rosetta to find the HBs, and also see whether I could tweak these criteria in any way (things like maximum distance, or maximum angles, etc.).
I have looked up and down the documentations, even scripts, and found no way of easily doing this. The closest I got was to find an "HBondOptions" class, which seems to be private, so I have no way of checking and/or modifying it, to the best of my understanding.
So this is my question: Is there a way to check/modify or otherwise relax the geometric criteria that PyRosetta/Rosetta uses to find/determine hydrogen bonds, or am I stuck with whatever Rosetta decides?
Thanks in advance for any help.