I noticed that this max_distance filter is not being applied to all relevant protein - ligand atom distances.
|
max_distance = 0.8 * (unit_cell[0] / 2) |
The intent of this filter is to select protein atoms that are less than half the box size away from 1. the selected ligand atom(s) and 2. the other protein atoms. Currently, it e.g. doesn't check the distance between P1-L1 (in the case where no protein atoms had been provided by the user) which could potentially lead to pbc wrapping issues.