Fix ligand selection when cofactors share resname UNK#119
Fix ligand selection when cofactors share resname UNK#119siddhantkhatod wants to merge 4 commits into
Conversation
In OpenFE hybrid topologies, both ligands and cofactors can be assigned resname UNK. gather_rms_data() now auto-detects the actual ligand by checking for hybrid tempfactors (0.25, 0.5, 0.75) when the default 'resname UNK' selection matches multiple residues. Fixes OpenFreeEnergy#118
In OpenFE hybrid topologies, both ligands and cofactors can be assigned resname UNK. gather_rms_data() now auto-detects the actual ligand by checking for hybrid tempfactors (0.25, 0.5, 0.75) when the default 'resname UNK' selection matches multiple residues. Fixes OpenFreeEnergy#118
|
pre-commit.ci autofix |
for more information, see https://pre-commit.ci
|
Hi @siddhantkhatod, thank you for this contribution, the tempfactor-based detection is a clever way to find the ligand without relying on residue names. |
Description
Fixes #118.
In OpenFE hybrid topologies, both the ligand and cofactors can be assigned the residue name UNK. The gather_rms_data function currently selects all atoms with resname UNK, which incorrectly includes cofactors in the ligand RMSD calculation.
Changes
Tests