tooluniverse-binder-discovery
Discover novel small-molecule binders for protein targets using structure-based and ligand-based screening. Covers druggability assessment, known-ligand mining (ChEMBL, BindingDB), similarity expansion, ADMET filtering, and synthesis feasibility. Use for hit identification, virtual screening, target
By mims-harvard · 379 installs
npx skills add mims-harvard/tooluniverse --skill tooluniverse-binder-discovery
Source repository · Upstream listing
Small Molecule Binder Discovery Strategy
Systematic discovery of novel small molecule binders using 60+ ToolUniverse tools across druggability assessment, known ligand mining, similarity expansion, ADMET filtering, and synthesis feasibility.
LOOK UP DON'T GUESS Always retrieve actual data from tools before drawing conclusions. Do not assume druggability, binding sites, or compound properties based on target class alone.
KEY PRINCIPLES :
1. Report first approach Create report file FIRST, then populate progressively
2. Target validation FIRST Confirm druggability before compound searching
3. Multi strategy approach Combine structure based and ligand based methods
4. ADMET aware filtering Eliminate poor compounds early
5. Evidence grading Grade candidates by supporting evidence
6. Actionable output Provide prioritized candidates with rationale
7. English first queries Always use English terms in tool calls. Respond in the user's language
Binding Site Reasoning (Start Here)
Before any tool call, reason about the target's structural biology:
Is the binding site a well defined pocket (small molecule accessible) or a flat protein protein interface (needs peptide/macrocycle)? This determines your screening strategy.
Enzymes with active sites (proteases, kinases, ATPases): deep, well defined pockets. Classic small molecule territory. Prioritize co crystal structure search and known inhibitor scaffold analysis.
GPCRs and ion channels : transmembrane pockets. Structure often available; start with GPCRdb and GtoPdb for known pharmacology.
Nuclear receptors : deep hydrophobic pockets. Excellent small molecule tractability; ligand based methods are well powered.
Protein protein interfaces : flat, large contact surface. Small molecules rarely compete effectively unless there is a "hot spot" cavity. Check whether any allosteric pockets exist before committing to small molecule strategy. Warn the user if no pocket is found.
Intrinsically disordered regions : essentially no small molecule approach. Redirect to peptide or degrader strategies.
Scaffolding / adaptor proteins : assess co crystal structures for unexpected pockets before declaring undruggable.
Use this reasoning to select phases and warn the user about challenges before executing a full workflow.
Critical Workflow Requirements
1. Report First Approach (MANDATORY)
DO NOT show search process or tool outputs to the user. Instead:
1. Create the report file FIRST Before any data collection:
File name: [TARGET] binder discovery report.md
Initialize with all section headers from the template (see REPORT TEMPLATE.md)
Add placeholder text: [Researching...] in each section
2. Progressively update the report As you gather data, update each section immediately.
3. Output separate data files :
[TARGET] candidate compounds.csv Prioritized compounds with SMILES, scores
[TARGET] bibliography.json Literature references (optional)
2. Citation Requirements (MANDATORY)
Every piece of information MUST include its source:
Example: Source: ChEMBL via ChEMBL get target activities (CHEMBL203)
Workflow Overview
Phases in order:
Phase 0 : Tool verification (check parameter names with get tool info )
Phase 1 : Target validation — resolve IDs, assess druggability, identify binding sites, predict structure if needed
Phase 2 : Known ligand mining — ChEMBL, BindingDB, GtoPdb, PubChem BioAssay, chemical probes; SAR analysis
Phase 3 : Structure analysis — PDB co crystals, EMDB (membrane targets), binding pocket characterization
Phase 3.5 : Docking validation — dock reference inhibitor to validate pocket geometry
Phase 4 : Compound expansion — similarity/substructure search (seeds: 3 5 diverse actives) + de novo generation
Phase 5 : ADMET filtering — physicochemical, bioavailability, toxicity, CYP, structural alerts
Phase 6 : Candidate docking and prioritization — score and rank top 20
Phase 6.5 : Literature evidence — PubMed, EuropePMC, OpenAlex
Phase 7 : Report synthesis and delivery
Phase 0: Tool Verification
CRITICAL : Verify tool parameters before calling unfamiliar tools.
Common parameter corrections (verify with get tool info if uncertain):
OpenTargets : ensemblId (camelCase); ADMETAI : smiles must be a list
NvidiaNIM alphafold2 (requires NVIDIA API KEY env var; free key at build.nvidia.com) : sequence not seq ; NvidiaNIM genmol (requires NVIDIA API KEY env var; free key at build.nvidia.com) : SMILES must contain [ {min max}]
NvidiaNIM boltz2 (requires NVIDIA API KEY env var; free key at build.nvidia.com) : polymers=[{"molecule type": "protein", "sequence": "..."}]
Phase 1: Target Validation
1.1 Identifier Resolution
Resolve all IDs upfront and store for downstream queries:
1.2 Druggability Assessment
Use multi source triangulation:
OpenTargets get target tractability by ensemblID(ensemblId) tractability bucket
DGIdb get gene druggability(genes=[gene symbol]) druggability categories
OpenTargets get target classes by ensemblID(ensemblId) target class
For GPCRs: GPCRdb get protein + GPCRdb get ligands + GPCRdb get structures
For antibody landscape: TheraSAbDab search by target(target=target name)
Decision Point : If no tractability data and binding site reasoning suggests PPI or disordered region, explicitly warn the user before proceeding.
1.3 Binding Site Analysis
ChEMBL search binding sites(target chembl id)
get binding affinity by pdb id(pdb id) for co crystallized ligands
InterPro get protein domains(accession) for domain architecture
1.4 Structure Prediction (NVIDIA NIM)
Requires NVIDIA API KEY . Two options:
AlphaFold2 : NvidiaNIM alphafold2(sequence, algorithm="mmseqs2") high accuracy, 5 15 min
ESMFold : ESMFold predict structure(sequence) fast (~30s), max 1024 AA
pLDDT guidance: =90 very high confidence, 70 90 confident, <70 use with caution. Low pLDDT in the putative binding region undermines docking reliability.
Phase 2: Known Ligand Mining
Priority order for bioactivity data:
1. ChEMBL get target activities curated, SAR ready
2. BindingDB get ligands by uniprot direct Ki/Kd with literature links
3. GtoPdb search ligands pharmacology focus (GPCRs, channels)
4. PubChem search assays by target gene HTS screens, novel scaffolds
5. OpenTargets get chemical probes by target ensemblID validated probes
Key steps:
1. Filter to IC50/Ki/Kd < 10 uM; retrieve molecule details for top actives
2. Identify chemical probes and approved drugs
3. Analyze SAR: common scaffolds, key modifications
4. Check off target selectivity: BindingDB get targets by compound
Phase 3: Structure Analysis
Tools:
PDB search similar structures(query=uniprot, type="sequence") find PDB entries
get protein metadata by pdb id(pdb id) resolution, method
get binding affinity by pdb id(pdb id) co crystal ligand affinities
get ligand smiles by chem comp id(chem comp id) ligand SMILES from PDB
EMDB search structures(query) cryo EM structures (prefer for GPCRs, ion channels)
alphafold get prediction(qualifier) AlphaFold DB fallback
Phase 3.5: Docking Validation (NVIDIA NIM)
If PDB + SDF available: use get diffdock info(protein=PDB, ligand=SDF, num poses=10) .
If only sequence + SMILES: use NvidiaNIM boltz2(polymers=[...], ligands=[...]) .
Dock a known reference inhibitor first to validate the binding pocket geometry before running candidates.
Phase 4: Compound Expansion
4.1 4.3 Search Based Expansion
Use 3 5 diverse actives as seeds, similarity threshold 70 85%:
ChEMBL search similar molecules(molecule=SMILES, similarity=70)
PubChem search compounds by similarity(smiles, threshold=0.7)
ChEMBL search substructure(smiles=core scaffold)
STITCH get chemical protein interactions(identifier=gene, species=9606)
4.4 De Novo Generation (NVIDIA NIM)
GenMol scaffold hopping with masked regions:
MolMIM controlled analog generation:
Phase 5: ADMET Filtering
Apply sequentially (all tools accept smiles=[list] ):
1. Physicochemical : ADMETAI predict physicochemical properties Lipinski violations <= 1, QED 0.3, MW 200 600
2. Bioavailability : ADMETAI predict bioavailability oral bioavailability 0.3
3. Toxicity : ADMETAI predict toxicity AMES < 0.5, hERG < 0.5, DILI < 0.5
4. CYP : ADMETAI predict CYP interactions flag CYP3A4 inhibitors
5. Alerts : ChEMBL search compound structural alerts no PAINS
Include a filter funnel summary in the report showing pass/fail counts at each stage.
Phase 6: Candidate Docking & Prioritization
Composite score: docking confidence (40%) + ADMET score (30%) + similarity to known active (20%) + novelty (10%, not in ChEMBL + novel scaffold bonus).
Evidence tiers for candidates:
T1 (3 stars): Experimental IC50/Ki < 100 nM
T2 (2 stars): Docking within 5% of reference OR IC50 100 1000 nM
T3 (1 star): 80% similarity to T1 compound
T4 (0 stars): 70 80% similarity, scaffold match only
T5 (no stars): Generated molecule, ADMET passed, no docking
Deliver top 20 candidates with: Rank, ID, SMILES, docking score, ADMET score, overall score, source, evidence tier.
Phase 6.5: Literature Evidence
PubMed search articles(query="[TARGET] inhibitor SAR") peer reviewed
EuropePMC search articles(query, source="PPR") preprints (not peer reviewed)
openalex search works(query) citation analysis
Fallback Chains
Programmatic Access (Beyond Tools)
When ToolUniverse tools return limited compound sets, access chemical databases directly:
See tooluniverse data wrangling skill for format cookbook and pagination patterns.
NVIDIA NIM Runtime Notes
AlphaFold2: 5 15 min (async, max ~2000 AA). ESMFold: ~30 sec (max 1024 AA). DiffDock: ~1 2 min/ligand. Boltz2: ~2 5 min. GenMol/MolMIM: ~1 3 min.
Always check: import os; nvidia available = bool(os.environ.get("NVIDIA API KEY"))
For large expansions ( 500 compounds): batch in chunks of 100, prioritize top candidates for docking.
Reference Files
[WORKFLOW DETAILS.md](./WORKFLOW DETAILS.md) Phase by phase procedures, code patterns, screening protocols
[TOOLS REFERENCE.md](./TOOLS REFERENCE.md) Complete tool reference with parameters and fallback chains
[REPORT TEMPLATE.md](./REPORT TEMPLATE.md) Report file template and evidence grading system
[EXAMPLES.md](./EXAMPLES.md) End to end workflow examples (EGFR, novel target, lead optimization)
[CHECKLIST.md](./CHECKLIST.md) Pre delivery verification checklist