Global Matching Sphere Optimization

From DISI
Revision as of 06:58, 30 April 2023 by Iamkaant (talk | contribs) (Created page with "== Goal == To optimize your matching sphere (MS) setups getting more enrichment with fewer spheres. == Description == The program performs optimization of matching spheres using genetic algorithm. It selects spheres from two sets: * heavy atoms of xtal-lig * spheres prepared by SPHGEN program At each generation, N matching sphere sets are created, containing a maximum of M spheres each. Then retrospective docking is done for each set, and sets are ranked by the enrichme...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Goal

To optimize your matching sphere (MS) setups getting more enrichment with fewer spheres.

Description

The program performs optimization of matching spheres using genetic algorithm. It selects spheres from two sets:

  • heavy atoms of xtal-lig
  • spheres prepared by SPHGEN program

At each generation, N matching sphere sets are created, containing a maximum of M spheres each. Then retrospective docking is done for each set, and sets are ranked by the enrichment, RMSD of the docked pose to the experimental one. After that, a quarter of sets "survive" and produce a new generation by direct transfer, mutations, and crossover. This process is repeated until enrichment, RMSD and minimum number of spheres do not change substantially in 10 generations. The program consists of two main modules: a Py

  • thon script (juggler.py) that performs MS generation, optimization, and ranking.
  • a Bas

h script, that watches created directory structure, runs docking and processes docking results

Setup

So far, the program is running on Wynton. LMK if you are interested in launching it on Gimel or other clusters. Prepare dockfiles directory with any tools of your liking (blastermaster, dockopt etc). You will also need rec.pdb, rec.crg.pdb, xtal-lig.pdb, ligands.names, decoys.names and a "sdi" directory with the paths to ligand .tgz files. To get RMSD of xtal-lig docked poses to the experimental pose, your xtal-lig.pdb must have correct bond orders and atom valences. You can edit it in Schrodinger and save as xtal-lig.pdb. Prepare juggler_config.yml file. Put it into an empty directory. Enter a screen environment so your run is not interrupted if you disconnect your SSH session. Then do: source /wynton/group/bks/soft/python_envs/python3.8.5.sh python ../../juggler-v0.7.py You can detach from the screen (Ctrl-A d). Open a new screen. In the same directory launch a docking daemon sh ../../rundockd-wynton-taskid.sh