The hit identification and drug discovery strategy consist in high-throughput docking for the identification of binder to the TKB domain of CBL-B that (1) bind to the same pocket and complete with the co-crystallized compound, (2) represent novel chemical templates and (3) bind with a KD below 30 uM.
Large library of commercially available compounds (Enamine) will be downloaded and prepared using the LigPrep preparation workflow from Maestro (Schrödinger, Inc). Briefly, the main tautomers from each compound will be generated and only compounds with reasonable physico-chemical properties will be considered (e.g., molecular weight, HBAcc/HBDon, number of rings, log P/D, PSA, ligand flexibility, as well as number of undefined chiral centers. We will also use our proprietary filters to eliminate reactive and chemically unstable compounds, compounds with undesirable functional groups and Pan-Assay INterfeering Structures (PAINS). The docking simulations will be performed on the PDB structure 8GCY.
The virtual screening stage will consist in a several-step molecular docking workflow using the popular Glide (Schrödinger, Inc) docking tool. First, all the compounds will be docked using the fast HTVS scoring function. Then, the top 500K compounds will be docked using the SP function. The top 50K compounds will then be docking using the more accurate XP function. For each considered PDB structure, the top 20K poses from SP function and 5K poses from XP scoring functions will be finally merged before the refinement stage.
The extracted poses will be subsequently re-evaluated using the independent SeeSAR (BioSolveIT) tool using the HYDE scoring function with key structural waters added. This structure-based software will be mainly used to discard putative false positives from virtual screening (e.g., poses with either geometry or energy warnings).
We will explore application of a MM/ML approach that approximate high level ab initio binding energies with a fast machine learning approach to rescore and prioritize the final list of virtual hits for ordering.
Finally, a thorough visual inspection by experienced computational and medicinal chemists over all kept poses with the aforementioned recommendations from a variety of methods (Glide SP & XP, HYDE scores and MM/ML score that make into account structural waters) will be conducted to select the most promising compounds to be experimentally tested.