Autogrow4: tutorial failure

Computer-Aided Drug Design at the Durrant Lab Forums AutoGrow4 Autogrow4: tutorial failure

Viewing 7 reply threads
  • Author
    • #16740

      HI! We are running autogrow4 on Linux cluster with python 3.7.6 and mgltools. When running the tutorial with the correct paths, we keep getting an error with the input files.No changes have been made to the file or file.

      python /home/witte.t/autogrow4/ -j /home/witte.t/autogrow4/tutorial_variables.json

      File "/home/witte.t/autogrow4/", line 657, in <module>
      vars, printout = load_in_commandline_parameters(INPUTS)
      File "/home/witte.t/autogrow4/autogrow/", line 964, in load_in_commandline_parameters
      File "/home/witte.t/autogrow4/autogrow/", line 438, in check_for_required_inputs
      "Receptor file can not be found. File must be a .PDB file."
      NotImplementedError: Receptor file can not be found. File must be a .PDB file.

      JSON file
      "filename_of_receptor": "./home/witte.t/autogrow4/tutorial/PARP/4r6eA_PARP1_prepared.pdb",
      "center_x": -70.76,
      "center_y": 21.82,
      "center_z": 28.33,
      "size_x": 25.0,
      "size_y": 16.0,
      "size_z": 25.0,
      "source_compound_file": "./home/witte.t/autogrow4/source_compounds/naphthalene_smiles.smi",
      "root_output_folder": "./home/witte.t/autogrow4/de_novo_design/autogrow4_results/",
      "number_of_mutants_first_generation": 10,
      "number_of_crossovers_first_generation": 10,
      "number_of_mutants": 3,
      "number_of_crossovers": 3,
      "number_elitism_advance_from_previous_gen": 3,
      "top_mols_to_seed_next_generation": 5,
      "diversity_mols_to_seed_first_generation": 3,
      "diversity_seed_depreciation_per_gen": 0,
      "num_generations": 1,
      "mgltools_directory": "./home/witte.t/autogrow4/mgltools_x86_64Linux2_1.5.6",
      "number_of_processors": -1,
      "selector_choice": "Rank_Selector",
      "max_variants_per_compound": 1,
      "filter_source_compounds": false,
      "use_docked_source_compounds": false,
      "LipinskiLenientFilter": true,
      "docking_exhaustiveness": 1,
      "start_a_new_run": true,
      "generate_plot": true,
      "rescore_lig_efficiency": true,
      "scoring_choice": "VINA",
      "dock_choice": "QuickVina2Docking",
      "gypsum_timeout_limit": 1,
      "debug_mode": false,
      "rxn_library": "all_rxns"

    • #16747

      hello. What are the permissions of the receptor.pdb? I.e.
      ls -l /home/witte.t/autogrow4/tutorial/PARP/4r6eA_PARP1_prepared.pdb

    • #16780

      Hi! Thank you! Changing the access to the folder did help! However it came up with another error in plotting.

      Finished generation 1


      Overall Scoring Average for all Compounds

      Traceback (most recent call last):

      File "/home/witte.t/autogrow4/", line 675, in <module>


      File "/home/witte.t/autogrow4/autogrow/", line 140, in main_execute


      File "/home/witte.t/autogrow4/autogrow/plotting/", line 442, in generate_figures

      dict_of_averages = print_data_table(infolder, folder_list) 

      File "/home/witte.t/autogrow4/autogrow/plotting/", line 388, in print_data_table

      average_affinity_dict = get_average_score_per_gen(infolder, folder_list) 

      File "/home/witte.t/autogrow4/autogrow/plotting/", line 104, in get_average_score_per_gen

      gen_affinity_average = gen_affinity_sum / num_lines_counter 

      ZeroDivisionError: float division by zero

    • #16835

      Still looks like a permission error to me. It’s either the ranked file(s) that was not created or it’s empty. Would you mind checking?
      ls -lR | grep "_ranked.smi"

    • #16842

      The smi files only allow the user to write which could be leading to an issue. However to hault this from happening on all codes. Where in the code would you suggest i put the "chmod u=rwx,g=rwx,o=rwx $PATH/generation1_ranked.smi" to change the permissions errors.

      I also noticed that the code is failing to "Convert Ligand to PDBQT format Begun" on all of the files.

      Failed to convert 0 times: /home/witte.t/autogrow4/de_novo_design/autogrow4_results/Run_4/generation_1/PDBs/naphthalene_60__1.pdb

      Thank you for your help!

    • #16942

      I wouldn’t suggest doing this, to begin with. Yet if no other solution is available, I would assume it makes sense to change the perms before the files are accessed. For example, in:
      You could add
      import os
      after line 8
      Then, after line 86 (after the loop for rank_file in ranked_file:)
      os.chmod(rank_file , 0777)
      Disclaimer: I’m not the author of autogrow project. This all looks like a poor man’s solution. Ideally would be to deal with the OS setup, so that the files would be accessible with no tweaks.
      Let’s try and call someone else to help you.

    • #17086
      Jacob Durrant

      Hi all. I noticed in your JSON file that "filename_of_receptor" is set to "./home/witte.t/autogrow4/tutorial/PARP/4r6eA_PARP1_prepared.pdb". Did you mean to include the initial "."? I suspect that’s the problem. Take care.

    • #17088
      Jacob Durrant

      Re. the second error, can you confirm that your ""*_ranked.smi" files exist, and that they are not empty? It may be that rerunning AutoGrow with the corrected JSON file with correct this second issue as well.

Viewing 7 reply threads
  • The topic ‘Autogrow4: tutorial failure’ is closed to new replies.