code: def extract_bs(self, cutoff, ligcentroid, resis):
          
          return [obres.GetIdx() for obres in resis if self.res_belongs_to_bs(obres, cutoff, ligcentroid)]
docs: Return list of ids from residues belonging to the binding site