STRC h01 Phase 8g-v2 — APBS Off-Target Rescore 2026-04-26
Replaces failed Phase 8g Vina off-target panel + inconclusive Phase 8h-lite #5 Coulomb proxy. Motivated by Hydroxamic Acid Divalent Cation Liability and the heavy-queue (2a) item from STRC h01 DR3 distillation log entry 2026-04-26.
Method
For each of 5 receptors (STRC E1659A K1141 reference + 4 cochlear off-targets), with crystallographic metals retained:
pdb2pqr --ff=PARSE(preserves HETATM Zn/Ca/Mg/Fe with PARSE ion charges).- APBS
mg-auto npbe, fine grid 65³, fine box 28 Å on lead-pose centroid, coarse box ≥ protein-extent + 40 Å, εp=2.0 / εs=78.54 / 0.15 M NaCl / T=310 K. - Sample φ at each hydroxamate-O atom of the Phase 8g best-pose (
v5.2__aq3__adamantyl__CONHOH__-ClMODEL 1). - Sample φ in 6×6×6 Å pocket-local box around lead centroid.
- Record metal proximity to anion-O.
Script: ~/STRC/hypotheses/h01-pharmacochaperone/scripts/phase8g_v2/apbs_offtarget_rescore.py. Output: artifacts/phase8g_v2_apbs_offtarget.json. Wall: ~5 min, 3-way parallel.
Result
Mean φ at the lead’s hydroxamate-O atoms (high → low):
| Receptor | anion-O ⟨φ⟩ (kT/e) | metals | min metal–anion dist |
|---|---|---|---|
| TRPM4 (8RD9) | +9.36 | 0 | — |
| TMEM16A (7ZK3) | +4.89 | 6 Ca²⁺ | 17.5 Å |
| STRC E1659A K1141 | +2.51 | 0 | — |
| KCNQ4 (7BYM) | +0.91 | 0 | — |
| Cx50 / Cx36 (8QOJ) | −3.02 | 0 | — |
Margin (STRC − max-off-target) = −6.85 kT/e (STRC is 6.85 kT/e LESS attractive than TRPM4 at the lead’s anion site). Gate verdict: FAIL.
Two surprises that rewrite the story
(1) Metals are not the failure mode
DR3 §I.6 predicted Zn²⁺/Ca²⁺ chelation as the hydroxamate liability. But TMEM16A’s 6 crystallographic Ca²⁺ sit 17.5 to 62.7 Å from the lead’s anion-O — far outside coordination range (~2.4 Å). The +4.89 kT/e at TMEM16A is residue-driven, not metal-driven. TRPM4 has zero crystallographic metals and the highest φ.
Implication: bioisosteric escape via tetrazole or CONHSO₂Me (hypothesised in Hydroxamic Acid Divalent Cation Liability) may not save selectivity — any anionic warhead will be pulled into a positive pocket.
(2) STRC reference is sub-optimal in the docked pose
Phase 5k MD-ensemble-averaged K1141 pocket centre: +5.99 ± 1.37 kT/e. Phase 8g-v2 lead anion-O on STRC: +2.51 kT/e. The lead’s hydroxamate-O is sitting at less than half the maximum-φ point of its own pocket. Either (a) the Phase 8e Vina pose for STRC is sub-optimally placed for electrostatic complementarity, or (b) the K1141 maximum-φ region is not where the lead’s anion lands geometrically.
If (a): re-rank Phase 8e poses by APBS φ at anion-O instead of Vina ΔG. Expect rescue. If (b): the lead is mechanistically under-optimised for the K1141 pocket; v5.3 design needs anion-position-aware geometry, not just Coulomb-aware composition.
Caveat — Phase 8g pose quality
The off-target poses are Vina blind-dock results from Phase 8g, which itself was flagged as Vina-failure. Vina without site definition often parks ligands in the most-electrostatically-favourable surface site, not the biologically-relevant binding pocket. The +9.36 kT/e on TRPM4 may reflect a non-functional surface basin, not the actual TRPM4 ligand-binding site.
This means the FAIL verdict is conditional: “the lead’s anion is electrostatically attracted to alternative cochlear-channel surfaces with stronger Coulomb pull than its target K1141 pocket as currently posed.” That is bad — but the strength of “bad” depends on whether the off-target pockets we found are real binding sites or surface artefacts.
What this changes
Mech 4 holds. The K1141-pocket binding mechanism is intact; this proof says nothing against it. What it says is that other proteins also have positive pockets attractive to this lead.
Deliv 4 → 3. Selectivity is no longer “wet-lab to determine” — it is “computationally suspect” until corrected. A pediatric ear drop whose lead has stronger Coulomb attraction to TRPM4 than to its target cannot claim selectivity for paper §3.
Tier A held, but on warning. If the next two compute jobs (Phase 8g-v3 site-defined off-target re-dock, Phase 8e-v2 STRC pose-anion-O re-rank) do not reverse this, tier drops.
Ranking delta
- tier: A (held, on warning)
- mech: 4 → 4 (held)
- deliv: 4 → 3 (selectivity gap quantified — first FAIL verdict, not just inconclusive)
- misha_fit: 4 → 4 (held — Misha’s specific need for the lead is unchanged)
- next_step: heavy queue reordered:
- (NEW) Phase 8g-v3 — site-defined off-target re-dock. Use UniProt/PDB-annotated functional pockets (TRPM4 cytoplasmic Ca²⁺/CaM site, TMEM16A inner-pore lipid-binding site, KCNQ4 retigabine pocket, Cx50 hemichannel pore). Re-dock lead at each defined site, then re-run Phase 8g-v2 APBS on those poses. Expected: drop the Vina-blind-dock surface artefact contribution.
- (NEW) Phase 8e-v2 — APBS re-rank of v5.2 STRC poses. Score all 14 ADMET-clean v5.2 poses on STRC by anion-O φ at the K1141 pocket; see whether any pose has anion-O above +5 kT/e. If yes → that lead becomes new mech-anchor; if no → v5.3 design must add anion-position-aware geometry constraint.
- (HELD) Phase 5l fragment ladder — substrate for τRAMD; AF3 jobs scaffolded at
~/STRC/models/af3_jobs_2026-04-26_phase5l/, awaiting Egor’s AlphaFold Server upload. - (HELD) Phase 5m τRAMD — skeleton at
scripts/phase5m_taramd_skeleton.py; depends on Phase 5l winner + ramd-openmm install. - (DEMOTED) v5.3 design (tetrazole/quinoxaline) — defer until Phase 8g-v3 confirms whether anion-warhead choice can rescue selectivity, or whether the issue is geometric (different question).
Connections
[part-of]h01 hub[supersedes]STRC h01 Phase 8g v5.2 Off-Target Panel 2026-04-25[supersedes]STRC h01 Phase 8h-lite Light Computational Evidence Package 2026-04-26 (selectivity proxy #5 only — other 4 proofs unaffected)[motivated-by]Hydroxamic Acid Divalent Cation Liability[falsifies-prediction]Hydroxamic Acid Divalent Cation Liability — metal-chelation predicted as primary mechanism; observed mechanism is residue-driven Coulomb on metal-free pockets[informs]STRC Pharmacochaperone Virtual Screen E1659A