I have a set of homologous bacterial proteins, where experimental evidence exists that a drug ligand binds to some of them, but not to others. The working hypothesis is that this is dependent on the protein structure.
Now I am looking for a suitable descriptor of the respective binding pockets that might inform how to mutate either of them to reproduce the binding behaviour of the other one.
Note that for the entire set of proteins only AlphaFold predictions of structure are available.
Is that a sensible endeavour to begin with? What would be a suitable descriptor, or better yet a tool that implements a comparison like that?