Binding affinity within 1 kcal/mol for novel scaffolds across protein families
Decades of effort and an experimental noise floor close to the target (about 0.9 kcal/mol assay reproducibility in public data) make scaffold-agnostic 1 kcal/mol accuracy research-grade and harder than the neighbouring drug-discovery problems (co-folding generalisation 47, residence time 61), though well below the defining open problems of mathematics. Application value 76 is the highest in this area because affinity ranking is the central economic lever of structure-based discovery and would serve antibody design (74) and generative chemistry alike; it sits above near-memory LLM decode (66) and below LWE (80). Leverage 56 since virtual screening, co-folding and generative chemistry all wait on a trustworthy scoring function, comparable to magic-state cultivation (55).
The problem
Predicting, within about a five-fold error, how tightly a brand-new kind of drug molecule sticks to its target.
- 1Candidate molecule and target
- 2Predict pose, then grip strength
- 3Compare with lab value
- 4Error near five-fold is the goal
- 5Needs new chemistry, blind
Why 1 kcal/mol. Binding free energy . At 298 K, kcal/mol, so an error of 1 kcal/mol changes the predicted by . Experimental reproducibility between laboratories is about 0.9 kcal/mol pairwise RMSE in public data, so a target of 1 to 1.2 kcal/mol RMSE is close to the limit.
Built from
earlier problems → this oneThe problems this one is built from: progress or tools it needs (builds on), special cases that come first (generalizes), and weaker results on the way (implies). On the right, what builds on it.
- Start with
- Converged reference interaction energies for 100-atom noncovalent complexes: reconciling CCSD(T) and diffusion Monte CarloComputational Chemistry · D40
- One classical force field accurate for folded proteins, disordered proteins and weak protein–protein associationComputational Chemistry · D51.5
- Predicting fold-switching proteins and both of their folds from sequenceComputational Medicine · D58
- Step 2
- A transferable exchange–correlation functional free of delocalization and self-interaction errorComputational Chemistry · D72
- Boltzmann-weighted conformational ensembles of multi-state proteins from sequenceComputational Medicine · D64
- Step 3
Path to understanding
foundations → problemThe research entries this problem rests on, in the order to read them. Each links to its explainer.
- Start with
- Reverse-mode automatic differentiationAlgorithm
- Gravitational N-body simulationAlgorithm
- Step 2
- The transformer architectureAlgorithm
- Molecular dynamics for biomoleculesAlgorithm
- Step 3
- Step 4
- Molecular docking and binding free-energy calculationAlgorithm · host
- Then
- Binding affinity within 1 kcal/mol for novel scaffolds across protein familiesThe problem
State of the art
- Physics-based relative FEP. On the curated FEP+ benchmark, FEP+ achieves pairwise RMSE of about 1.25 kcal/mol, against 0.91 kcal/mol for experimental reproducibility in public data (Ross et al., Communications Chemistry 2023). It performs well for congeneric series on well-behaved targets, much less well for charge-changing perturbations, scaffold hops, flexible loops, membrane proteins, and nucleic-acid targets (about 1.39 kcal/mol RMSE for DNA/RNA binders in a vendor-sponsored preprint). Open implementations (OpenFE) are somewhat less accurate than FEP+ in a 59-dataset comparison by OMSF, with atom-mapping failures producing outliers.
- CASP16 (2024). In the first large blind affinity test in the CASP series, 30 groups predicted 140 affinities on two targets: the best Kendall's was 0.42, below the noise ceiling of about 0.73; providing experimental structures did not help, which the assessors take as evidence that scoring functions are the limiting factor.
- Learned affinity heads. Boltz-2 (preprint, June 2025) reports average Pearson 0.66 on a four-target FEP+ subset and approaches OpenFE on 876 measurements at more than 1,000 times lower cost, outperforming all CASP16 participants out of the box; on eight blinded internal industrial assays it exceeded Pearson 0.55 on only three, and the authors state that performance "varies strongly between assays". These are developer-reported numbers without independent prospective replication in sources reviewed.
- Hit discovery. The same preprint reports an enrichment factor of 18.4 at the top 0.5% in a retrospective screen (MF-PCBA); this measures ranking against decoys, not accurate ΔG on scaffolds.
- Not shown. No method meets the benchmark above prospectively across families without target-specific information.
Realistic applications
Accurate scaffold-agnostic affinity prediction would let medicinal chemists rank virtual compounds before synthesis, cut the number of design-make-test cycles in lead optimisation, make large virtual screens reliable, and allow design of selectivity (predicting off-target binding). It would also be the scoring function needed for generative chemistry.
What counts as solved
A prospective, pre-registered, blinded evaluation (a D3R/CASP-style or industry-sponsored challenge with unpublished measured affinities) in which one method, with no target-specific data, meets the RMSE, unsigned-error and Kendall's thresholds above on at least 10 protein families and at least 200 compounds with novel scaffolds, with the analysis published in a refereed journal and the method accessible to independent groups.
How it gets verified
ObservationSettled by new measurements or observations.
Needs unpublished measured affinities: one method with no target-specific data must meet the RMSE, unsigned-error and Kendall tau thresholds on at least 10 protein families and 200 compounds with novel scaffolds, in a prospective, pre-registered, blinded challenge (D3R/CASP-style or industry-sponsored).
Claims under verification
+ add a claimNo claimed results on the Verification board for this problem yet.
Where the work is
Gains need better energy functions and sampling, much more compute-efficient free-energy protocols, and prospective blinded data for scaffold-diverse compounds.
Bottlenecks
ranked by votes- #1EmpiricalNoise floor and lack of prospective blinded data
Assay reproducibility is about 0.9 kcal/mol pairwise RMSE; public sets are leaky and few prospective blinded affinity tests exist (CASP16 had two targets).
Community grades
Sign in to add your own grade.
Connected problems
Progress
- Oct 1, 2025CASP16 affinity categorysource
Best Kendall's tau 0.42 over 140 affinities; experimental structures did not help.
- Jun 18, 2025Boltz-2 (preprint)source
First open co-folding model with an affinity head reporting approach to FEP-level correlation on selected sets; large variance between assays.
- Oct 1, 2023Maximal accuracy of rigorous FEPsource
Ross et al. quantify experimental reproducibility (about 0.91 kcal/mol pairwise RMSE) against FEP+ performance (about 1.25).
- Feb 18, 2015FEP+ prospective protocolsource
Wang et al. report mean unsigned errors near 1 kcal/mol for relative potencies on eight targets.