Difference between revisions of "De novo Developer Progress"

Latest revision as of 11:53, 4 February 2019

This is the Rizzo lab wiki page for coordinating bugs and progress on the de novo project.

Valgrind clean version of the code on cluster that Rizzo lab should be using:

Lauren:

/gpfs/projects/rizzo/zzz.programs/dock6.9_release
This version includes all changes of the merge.

Path to Generic Fragment Library:

/gpfs/projects/rizzo/leprentis/gen-frags-12

Path to Frequency Anchors:

/gpfs/projects/rizzo/leprentis/zinc1_ancs_freq

Current Coding Progress:

Working on these currently:

1. Lauren: Check MGS+(-50)TAN before and after fingerprinting fix for 663 systems
2. Lauren: Implement Roulette fragment picking into graph and random as an option
3. Lauren: Implement Adjacency Matrix into fraglib/dn (initialize matrix and utilize matrix for graph and random fragment picking)
4. Lauren: Testing with sb2012 default values for molecules passed on to the next layer and root (looking for timing and efficiency).
5. Lauren&John: Rework VS protocol to integrate de novo protocol more smoothly

Need to be fixed:

1. score_molecules and internal_energy problem (for simple_build)
2. HMS needs to fixed when no heavy atoms matching
3. Get rid of some scoring functions and only use descriptor score

Not working on these right now:

1. Lauren: Adjacency matrix vs tors env
2. Lauren: Addition of "3mer" combination fragment check (post tors check)
3. Lauren: Min and Max for charge to replace absolute value of charge.(Broke everything)
4. Lauren: Capping groups for post growth process (halogens and methyls)
5. Lauren: Fix Frag_String output into chimera for Refinement situations (current space can remove the spaces in the mol2 file - temp fix)
6. Stephen: Change scaling factor to a function of decay (currently a straight line to lowest score cutoff)
7. Lauren&John: SMILEs and ZINC script (for dn and ga)
8. Lauren: incorporate tan pruning as final step (post growth) as user option (replace make_unique script)

Completed:

1. Lauren: determine if random seed is reset for each aps
2. Lauren: Create testset for each dn function
3. Lauren: Test simple build function with merged de novo
4. Lauren&Stephen: clean make_unique script for release
5. Lauren: merge GA into dock/dn
6. Dwight & Lauren: MPI wrapper for 192 processors (8 nodes) for testsets on rizzo cluster
7. Lauren: Create short testsets for denovo frag gen, focused fragment generic for DOCK6.9 release
8. Dwight+Lauren: merge parameter files of de novo with DOCK
9. Lauren: add dn_defn file for separate defn with Hydrogens
10. Lauren: Implement csingleton fix for orienting fragments with less than 3 heavy atoms
11. Lauren: Test bfochtman fix for rotatable bonds within an user defined anchor
12. Lauren: Test csingleton fix for orienting fragments with Du
13. Lauren: Test MGS focused fragment library results with dn paper
14. Stephen: editting script to calculate SMILE string of de novo molecules in OpenBable

List of features that we definitely want for the 6.9 release:

List of features/ideas for future releases:

• Using different references for different layers of dn growth
• Stereo centers / volume overlap pruning
• Capping group functions (H, CH3, Halogen)
• Incorporate GA at the end of each layer
• Overhaul the simple-build function
• Monte carlo algorithm that checks bond frequency
• Scaling max root / layer size with layer
• Select torenv before selecting fragment. Will need to overhaul fraggraph, will keep us from needing to assemble mols that will not pass torenv.
• Add fragname string to restart and dump files, already done for final and fraglib files.
• Add ZINC name to torenv table
• Unusual behavior during library generation when frequency cutoff == 0
• Print out how many molecules cannot be capped. (Difference between ensemble size and dump.)
• building from anchor 0 -> building from scf.98
• Possible torenv check for dump molecules after capping before printing.
• keep tables of what fragments (and torsion types) are already included in a growing molecule (i.e.e the name string has this info) and only accept a new fragment (or torsion type) within certain ranges and probabilities. In other words use knowledge of chemical makeup probabilities to keep from over including or under including certain fragment and bond types (essentially use datamining to help us only build molecules within certain boundaries)

List of SB2012 systems that we will use for tests:

For now, let's use 5-15 rotatable bonds inclusive; total = 709 systems ("drug-like" size molecules). De novo paper only used 663 systems that removed 46 systems where the cognate ligand did not fall with a +/-2 formal charge. (5through15 = 709, 5through15_ch2 = 663)

{5RB = 107; 6RB = 96; 7RB = 103; 8RB = 75; 9RB = 66; 10RB = 75; 11RB = 57; 12RB = 41; 13RB = 38; 14RB = 26; 15RB = 25}