masif_seed icon indicating copy to clipboard operation
masif_seed copied to clipboard

select specific sites for masif_seed_search

Open slieped opened this issue 1 year ago • 1 comments

I was trying to find an option to select a specific site (not residueID) to run the masif_seach docking algorithm. I was checking the code at masif_seed_search_nn.py and realized that the conditional block at line 103 is coded to return the specific topN instead of a specific site.

if 'selected_site_ixs' in params:
    site_ixs = params['selected_site_ixs']
    site_vixs = [vix for ix,vix in enumerate(target_vertices) if ix in site_ixs]

Is this the expected outcome? Should the variable params["selected_site_ixs"] in params_peptide.py allow selecting a specific site? If thats the case, a simple edit in the code as follows should solve it!

site_ixs = [ix for ix,vix in enumerate(params['selected_site_ixs'])]
site_vixs = [vix for vix in params['selected_site_ixs'] if vix in target_vertices]

This will return the desired site, if and only if its part of the TOP_N vertices

Thanks for this amazing tool! Victor

slieped avatar Jul 12 '23 09:07 slieped