I'm having some problems using ligand_identification. I would like to restrict the search to a specific density peak, even if it is not the highest unmodeled peak or the highest Fo-Fc peak in the map. I tried using options: search_center="67.5 18.3 11.8" or ligand_near_res=S2063, Will the search be restricted to that region, or if a particuar ligand doesn't fit that blob, will it search through the rest of the map? If it is restricted, in what radius? Is this radius affected by the "search_dist" or "local_search" parameters? and local_search = True is default? Is there a threshold level for density level, below which building a ligand in a blob will not be attempted? I've tried with both search_center= and ligand_near_res=, and something gets built far away from that site. But there were errors, so I may have something wrong: phenix.ligand_identification mtz_in=sqr2803or13_031.mtz input_labels="2FOFCWT PH2FOFCWT" \ model=sqr2803or13_031.pdb ligand_near_res=S2063 nproc=2 After preparing the ligand library, then: Running LigandFit process 1... Number of atoms in ligand suc.pdb is 23 Running job sequence 1, ligand 2, in /tb/sb/usr20c/berry/ref/sqrdep/sqr2803dep2/phenix2/Temp_1... Evaluating all ligands in ligand-lib now...and placing fittedligand ### in resolve_ligand_###.pdb Number of atoms in ligand 2pe.pdb is 28 Running job sequence 0, ligand 1, in /tb/sb/usr20c/berry/ref/sqrdep/sqr2803dep2/phenix2/Temp_0... Process Process-2: Traceback (most recent call last): File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/multiprocessing/process.py", line 258, in _bootstrap self.run() File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/multiprocessing/process.py", line 114, in run self._target(*self._args, **self._kwargs) File "/sw/lnx/phenix-1.12-2829/build/../modules/phenix/phenix/command_line/ligand_identification.py", line 1382, in RunLigandFit shutil.rmtree(ligandfit_dir) File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/shutil.py", line 247, in rmtree rmtree(fullname, ignore_errors, onerror) File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/shutil.py", line 256, in rmtree onerror(os.rmdir, path, sys.exc_info()) File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/shutil.py", line 254, in rmtree os.rmdir(path) OSError: [Errno 39] Directory not empty: '/tb/sb/usr20c/berry/ref/sqrdep/sqr2803dep2/phenix2/Temp_1/LigandFit_run_1_/TEMP0' Process Process-1: Traceback (most recent call last): File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/multiprocessing/process.py", line 258, in _bootstrap self.run() File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/multiprocessing/process.py", line 114, in run self._target(*self._args, **self._kwargs) File "/sw/lnx/phenix-1.12-2829/build/../modules/phenix/phenix/command_line/ligand_identification.py", line 1382, in RunLigandFit shutil.rmtree(ligandfit_dir) File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/shutil.py", line 247, in rmtree rmtree(fullname, ignore_errors, onerror) File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/shutil.py", line 256, in rmtree onerror(os.rmdir, path, sys.exc_info()) File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/shutil.py", line 254, in rmtree os.rmdir(path) OSError: [Errno 39] Directory not empty: '/tb/sb/usr20c/berry/ref/sqrdep/sqr2803dep2/phenix2/Temp_0/LigandFit_run_1_/TEMP0' Evaluating LigandFit results ... The run continues, but it does not test any more ligands after those first two but goes on to evaluate the results. With nproc = 1, only 1 ligand gets tested. In all cases the first ligand in the library (2PE.pdb) is evaluated as the best. It is placed in density, but density that has already been built out with (and looks more like) a string of water molecules. And this is far from the selected residue or coordinates specified. The directory that raised the error when attempting to be deleted does eventually get removed: after the run there is no TEMP_N in the parent directory. Any suggestions would be welcome. Ed P.S. - running with .eff file: ['--show_defaults'] ligand_identification { mtz_in = sqr2803or13_031.mtz mtz_type = *F diffmap model = sqr2803or13_031.pdb ncpu = 1 n_indiv_tries_min = 30 n_indiv_tries_max = 300 n_group_search = 4 search_dist = 10 local_search = True search_center = "67.5 18.3 11.8" # ligand_near_res = S2063 verbose = False debug = False use_ligandfit = True search_mode = *default LigandFit temp_dir = Auto dry_run = False # number_of_ligands = 1 cc_min = 0.75 open_in_coot = False non_bonded = True keep_all_files = False # cif_def_file_list = real_space_target_weight = 10 # job_title = None ligandfit { } } gives: [['2pe.pdb', 'suc.pdb', . . . 'upl.pdb']] /tb/sb/usr20c/berry/ref/sqrdep/sqr2803dep2/phenix2 ******************************************************************************* Sorry, the protein model file None does not seem to exist? ******************************************************************************* Running LigandFit process 0... Process Process-1: Traceback (most recent call last): File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/multiprocessing/process.py", line 258, in _bootstrap self.run() File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/multiprocessing/process.py", line 114, in run self._target(*self._args, **self._kwargs) File "/sw/lnx/phenix-1.12-2829/build/../modules/phenix/phenix/command_line/ligand_identification.py", line 1122, in RunLigandFit shutil.copyfile(mtz_in,data_local) File "/sw/lnx/phenix-1.12-2829/base/lib/python2.7/shutil.py", line 82, in copyfile with open(src, 'rb') as fsrc: IOError: [Errno 2] No such file or directory: 'None' Evaluating LigandFit results ... Lig_seq Placed/total cc_all cc cc_adj score Code HBscore Cannot find overall_ligand_scores.log0. This could mean that none of that (sub)set of ligand fitted well. None of the ligand fit the difference desity well enough. Please try the following -- 1) if you input a custom library, try to use the default library (no extra keywords needed), or 2) if you used the default library already, ususlly this means that the density is too small (> 6 atome or more is needed.) Exiting ...... No good ligand found.