Mutation-aware docking

Find the compounds that
shift toward the mutant.

Pick a clinically relevant mutation, pick your compounds, and Liganx docks them against wild-type and the mutant in parallel — then ranks the shifts so you can prioritise what to test next. A fast first pass, not a final verdict. No PyMOL, no FoldX setup, no AutoDock wrangling.

No install Vina · GNINA CNN · Boltz-2 co-folding — all live Mutation-aware virtual screening Free for academic use

Built for proprietary work. No account required, no third-party trackers or cookies, and nothing leaves North America. Your compounds are sent to a GPU worker only for the run itself, then discarded there — and we'll delete any job on request. How we handle your data →

Sample
EGFR · selectivity matrix
example
CompoundWTT790MC797S
Osimertinib-8.8-9.3-5.2
Gefitinib-7.2-5.1-4.9
Erlotinib-7.6-5.4-5.3
Compound X-6.1-7.8-7.1
Liganx flags Compound X — its docking score shifts 1.7 kcal/mol toward the T790M mutant — as a hypothesis worth testing.
Vina scoring noise is roughly ±1 kcal/mol at default exhaustiveness — Δs above ~1 kcal/mol are interpretable, smaller deltas live near the noise floor.
Benchmarked, not claimed

Every result is a fresh GPU run with the cache off — and the whole benchmark is public.

We publish the full pose-accuracy benchmark, the per-pose confidence-badge validation, and the per-target scores — so you can audit the numbers instead of taking a slogan on faith.

GNINA 69% top-1 < 2 Å on Astex-42, automated · confidence badge 12/13 right when green · see Validation →

Built on tools the community already trusts

AutoDock VinaRDKitFoldXMol*ProLIFRunPod
New · Flagship · Beta

Resistance Radar — flag which mutations are most likely to weaken your compound

Rank a target's known mutations by how likely each is to weaken your compound — in one click. No other self-serve tool packages this.

Bring your own molecule. In one pass, Resistance Radar docks it across a target's entire known resistance panel and returns a ranked liability map — per-variant Δ-vs-wild-type docking plus a forecast score for how likely each mutation is to confer resistance (calibration not yet validated). It's the resistance question asked about your compound, not a generic drug: where will it hold, and where will it fail?

  • One pass, whole panel. Submit once; every clinically-relevant variant is co-docked against wild-type under identical conditions.
  • A ranked forecast, not just a docking Δ. The same 2-signal model as the Atlas (docking Δ + ESM-2 fitness), 5-fold cross-validated ROC-AUC 0.81 (95% CI 0.62–0.96, n=25). AUC measures ranking; calibration as a probability is not yet validated.
  • A shareable liability map. Every scan gets its own page and link — hand a collaborator the resistance forecast, no re-running.

Early-access beta. Panel coverage is growing target by target; request access from the Studio.

KRAS · your compound
liability map
  • Q61H
    likely
  • G13D
    borderline
  • G12C
    borderline
  • G12V
    unlikely
Illustrative preview — run your own compound for live, scored results.
How it works

From mutation to selectivity matrix in three clicks.

1

Pick a target

Choose from clinically actionable kinases or upload your own PDB. Pocket boxes are pre-defined.

2

Pick mutations

Click EGFR T790M, KRAS G12C, BRAF V600E — anything from the curated library, or type your own.

3

Read the matrix

We dock every compound against WT and each mutant. The Δ-score highlights how each compound's binding shifts between variants.

What's new

Recently shipped.

Just shipped

Model confidence on every Boltz-2 prediction

Every Boltz-2 pose now carries the model's own confidence (ligand pLDDT), binned green / amber / red and checked on a 30-target holdout of post-2024 structures Boltz-2 never trained on. There, green poses were right 12 of 13 times; low-confidence poses are flagged, so you don't chase a bad hit. A small-sample trust signal, not a calibrated probability.

See it in Studio
Just shipped

Resistance Atlas — rank the mutations most likely to break a drug

For every FDA-approved targeted cancer drug, the Atlas ranks which mutations are most likely to emerge as clinical resistance. Combines docking Δ + ESM-2 protein-language-model fitness in a 2-signal logistic model, fit on 25 published clinical-resistance events (16 resistance / 9 non-resistance): ROC-AUC 0.90 in-sample, 0.81 5-fold cross-validated (95% CI 0.62–0.96 — small-n, indicative). 15 drugs covered today. Novel (gene, position, mutant) lookups now run real ESM-2 on our GPU pod on demand.

Open the Atlas
Pro beta

Test your own (drug, mutation) data against our model

Upload up to 10 (gene, position, wt, mutant, drug) rows as CSV. We score each through the same 2-signal model the Atlas uses — real ESM-2 inference for novel mutations, instant cache hits for the 49-event ESM-2 cache — and return a forecast score (uncalibrated), verdict, and AUC if you provide ground truth. Free tier: 10 rows / day.

Calibrate your data
Just shipped

Pre-computed FDA-drug screenings — no setup, no GPU wait

We pre-ran 30 oncology kinase inhibitors against every resistance mutation in our catalog — KRAS G12C/G12D/Q61H, EGFR T790M/L858R/C797S, BCR-ABL T315I/E255K. Hit a public URL and see ranked selectivity hits in 1 second. Click any compound to see its 3D pose.

Browse pre-computed screenings
What you get

Everything an early-discovery med-chemist actually wants.

Selectivity matrix

N compounds × M mutants in one view, cells colored by Δ-score so resistance and selectivity gain pop out instantly. The whole product on one screen.

MM/GBSA affinity ranking

Single-snapshot rescore of the docked pose (ff14SB / OpenFF Sage / OBC2, ε=1, no salt). Ranks compounds within a target — not a Kd. Short MD is a pose-stability badge, never averaged in. Median Spearman ρ 0.71 across 8 congeneric series (184 ligands; range 0.46–0.83), rescoring reference poses. Available by request.

Pre-computed FDA-drug screenings

30 FDA-launched oncology kinase inhibitors pre-docked against every catalog resistance mutation — public landing pages, no login. Open any card on /library and see the ranked selectivity hits in one second.

Mutation-aware virtual screening

Drop up to 1000 compounds against a (target, mutation) pair; we dock each against WT and the mutant in parallel and return a hit list ranked by selectivity index. Promote the top hits to a full job in one click — no re-dock.

Resistance Atlas + test your own data

A public per-drug atlas ranking the mutations most likely to break each FDA-approved targeted drug, fit on 25 clinical-resistance events (16 resistance / 9 non-resistance; out-of-fold ROC-AUC 0.81, 95% CI 0.62–0.96 — small-n). Upload your own (drug, mutation) CSV to score against the same ESM-2-backed model.

Validated poses on GPU Vina

Free jobs run on QuickVina2-GPU, with every pose checked by PoseBusters, Vinardo re-score, and RDKit strain analysis — most free tools give you no validation at all. GNINA CNN rescoring and Boltz-2 ML co-folding run side-by-side on the same job on Pro.

Honest by default

Every Δ near the ±1 kcal/mol Vina noise floor gets a within-noise badge; mutations outside the pocket are flagged, not scored; every pose comes with a plain-English readout and an inline ADMET panel (hERG, DILI, CYP, BBB). We publish our method limits on purpose.

MM/GBSA

Physics-based ranking, one rung above docking.

Docking is fast triage; MM/GBSA is the next rung of confidence. We take the docked pose and do a single-snapshot rescore with implicit-solvent molecular mechanics (Amber ff14SB / OpenFF Sage 2.2 / OBC2, εin = 1, no salt), then rank your compounds by the binding-energy estimate. A separate short MD run checks the pose actually holds — reported as a stability badge, never averaged into the score.

Validated across 8 congeneric series (184 ligands): median Spearman ρ 0.71, Pearson r 0.67, ranging from HIF2A ρ 0.46 to TYK2 ρ 0.83 — rescoring reference poses. We lead with rank correlation because that is the honest measure of a rescoring signal, and we publish every target, including the weak ones.

TYK2 · MM/GBSA vs experimentSpearman ρ 0.8313 ligands · reference poses
−8−9−10−11−12Experimental ΔG (kcal/mol)−42−44−46−48Predicted MM/GBSA

Each dot is one ligand; the dashed line is the least-squares fit. More-negative = stronger binding. Curated reference poses.

Case study (TYK2, n = 12): our template-anchored docking — which uses the co-crystal ligand as a reference but docks the query itself rather than borrowing a crystal pose — followed by MM/GBSA rescoring reaches ρ 0.72, versus ρ 0.83 when the same ligands are rescored from reference poses. One series so far; other targets still run the de-novo docking path.

MM/GBSA magnitudes run several-fold larger than experimental ΔG and omit configurational entropy, so read this as a way to order compounds within a target, not an absolute affinity. For the close calls, FEP is available by request.

MM/GBSA scoring is available by request.

Where we sit

The missing middle.

Free serversLiganxSchrödinger Maestro
Compound resistance forecast across a variant panelpartial
Runs in the browser, no install
Mutation-aware WT-vs-mutant matrixpartial
Public resistance-mutation atlas
Pre-computed FDA-drug screenings
Bulk virtual screening, selectivity-rankedpartial
Inline ADMET (hERG / DILI / CYP / BBB)partial
Multiple scoring engines side-by-sidepartial
Ensemble / flexible-receptor docking
Published, reproducible validation report

Reflects publicly known features as of May 2026. Free-server and Schrödinger capabilities vary by version, license tier, and module.

Method limitations we publish on purpose
See full validation report →
  • Eleven literature-anchored controls, public verdict. ABL T315I, EGFR T790M, BRAF V600E, KIT D816V, BTK C481S, KRAS G12C, EGFR C797S, EGFR L858R — five of eleven PASS at above-noise magnitude in the published direction. Five NOISE results sit in documented method-limit territory (covalent acrylamides, active-conformation selectivity, conformational activation). One FAIL (EGFR L858R + Gefitinib) is explained candidly with the structural reason — rigid-receptor docking can't capture L858R's conformational activation. The full per-case verdict and the open-source script that re-derives it are public.
  • Vina noise floor. Vina/QuickVina2 scoring has roughly ±1 kcal/mol noise at default exhaustiveness. We surface a "within-noise" badge for any Δ inside that band so a reader doesn't over-interpret 0.3 kcal/mol shifts.
  • Mutant-receptor build path. The hosted app builds the mutant with PDBFixer's residue substitution, which applies the new identity but does not energy-minimise the structure; FoldX BuildModel is used only on environments where an academic licence is present. Drastic side-chain changes (e.g. small→large) can introduce clash signal in the Δ that isn't pure binding affinity. Submitting the same mutation in both modes and comparing flags this when it matters.
  • FoldX academic licensing. FoldX is free for academic use under its own EULA but requires a commercial licence for industry workflows. Liganx ships the FoldX call path; users running commercial work should verify their licence with the FoldX team at the Centre for Genomic Regulation directly.
  • Single-conformation rigid-receptor docking. We dock against one PDB conformation per (target, mutation). Mutations far from the binding pocket — typical for activation-loop, allosteric, or distant-domain residues — get an "outside pocket" badge instead of a possibly-misleading Δ, and the matrix shows them as not-scored rather than zero.

Stop hand-rolling mutation-aware docking.

One UI. Real Vina under the hood. Selectivity matrix in minutes, not days.