← projects

Action Identifiability in Protein-Design Agents, and Evaluating Biological Agents

Saanvi S. Subramanian, Gabriela Lobinska

AITHYRA

MoML 2026 (MIT) AgenticLS @ NeurIPS 2026 WiML @ NeurIPS 2026 AI for Science · Trust-AI-Eval · GEM-Bio · Verify-Agents @ NeurIPS 2026

A closed-loop molecular agent can only choose between actions its scorer keeps apart. This report asks whether a standard scorer does, builds a five-category battery for auditing any homology-free claim, and turns the same lens on the proxy, benchmark, and verifier an agent leans on.

A scorer that cannot resolve competing actions turns every reported mechanism into a guess. Across all three parts, a number read in a learned representation stops meaning what it claims the moment an optimizer selects on it.

Abstract

Molecular design agents choose interventions by scoring predicted outcomes, which requires the scorer to distinguish the available actions. This joint report tests that requirement with 78 COSMIC mutational processes as actions and distance in ESM-2 space as the score, then turns the same discipline on the proxy, benchmark, and verifier a design agent delegates to.

Median pairwise distance falls from 0.815 in COSMIC spectrum space to 0.0068 between ESM-2 outcome regions, a 120x compression, and a composition-matched control reproduces the destinations. A learned policy identifies the true action in 7.55% of held-out cases, split-conformal prediction reaches 92.2% coverage only by returning 66.8 of 78 actions, and an independent Pfam profile HMM finds no PF17041 hits at the curated threshold. A five-category audit battery separates a reproducible embedding effect from the mechanistic claim attached to it, and the agent's parts fail the same way: a zero-shot score ranks below counting, a benchmark correlation stops predicting selection, and a trust gate certifies a scorer that memorized the answers.

120×
Compression of the action space, 0.815 in COSMIC space to 0.0068 between ESM-2 destinations
7.55%
Top-1 accuracy of a policy naming the action, against 1.28% chance and 92.45% still wrong
0.991
Gate AUROC of a scorer that memorized the outcomes, against 0.663 for the truthful scorer
A useful eval answersIn this report
What it measuresMovement in a learned representation, a benchmark correlation, a gate's in-sample AUROC
What it leaves outBiological function, selection value at the top of a ranking, a supplied scorer's implementation
How it can be gamedA composition-matched direction, a memorizing scorer, a single-residue metric exploit at 5 s.d.
When it should change a decisionOnly after matched controls and an independent check survive; otherwise widen the set or abstain

Part I · Action identifiability

A scorer can only choose between actions it keeps apart

MoML 2026 (MIT)AgenticLS @ NeurIPS 2026

A closed-loop molecular agent proposes an intervention, scores the resulting molecule, and uses the score to choose what to do next. This searches large spaces efficiently, but only if the scorer resolves the actions available to the policy. If several interventions produce indistinguishable scored outcomes, the agent cannot infer which intervention caused a result, so any mechanism it reports afterwards is unidentified.

Action identifiability The ability to distinguish the outcomes of competing actions using the agent's own scoring model.

Identifiability and biological admissibility are different properties. An action can represent a perfectly valid biological process even when the scorer discards exactly the information that distinguishes it from other processes. The two come apart, and only one is usually checked.

To test the requirement I need an action space whose ground-truth distinctions are known independently of any learned model. The 78 COSMIC single-base-substitution signatures give exactly that: each is a distribution over 96 trinucleotide substitution types, catalogued from real mutational processes. An action is a signature, executed on a coding sequence to produce a mutant, which is translated, embedded with ESM-2, and scored by cosine distance to a target Pfam-domain centroid. That gives both halves of the comparison: how far apart two actions are as processes, and how far apart their outcomes are after scoring.

The scorer compresses the action space

Across all 3,003 pairs of COSMIC signatures, the median distance between native 96-channel spectra is 0.815; between the same actions' ESM-2 outcome regions it is 0.0068. The scored action space is roughly two orders of magnitude tighter than the space the actions actually live in. 1,780 of 3,003 pairs are more than 0.5 apart in COSMIC space yet fall within 0.011 of each other after scoring. SBS1 and SBS7a, clock-like ageing and UV exposure, are 0.952 apart as spectra and 0.0055 apart as outcomes.

A composition-matched control pins down why. Scrambling trinucleotide context and mutation placement while preserving the realized substitution composition reproduces each action's ESM-2 destination to within cosine distance 0.02, and the same pattern appears at both ESM-2 8M and 650M. The agent has 78 nominal actions and a far coarser effective action space underneath.

COSMIC spectrum spaceESM-2 outcome space
median pairwise distance 0.815
Interactive. The scorer compresses the action space. Scrub from the 78 mutational actions as they sit in COSMIC spectrum space, spread apart, to where their outcomes land after ESM-2 scoring, one tight region. Median pairwise distance falls from 0.815 to 0.0068, a 120x compression, and SBS1 and SBS7a, far apart as processes, land almost on top of each other. Point positions are schematic; the collapse magnitude is the measured 120x.
Left, a steering cross-matrix showing each mutational process moves held-out proteins toward its own outcome region, with SBS1, SBS7a and SBS2 overlapping and SBS17a farther out. Right, the same pattern in pairwise destination distance.
Figure 1. Distinct processes reach the same place. Left, enacting each process moves held-out proteins toward its own outcome region, but SBS1 (clock-like), SBS7a (UV) and SBS2 (APOBEC) overlap while SBS17a sits farther out. Right, the same in pairwise destination distance. SBS2 and SBS1 are near-orthogonal in COSMIC space at cosine 0.02 and still merge after scoring, so the representation collapses a distinction the catalogue makes.
Finding 1

Biological grounding leaves scorer-level identifiability open: the effective action space is far coarser than the nominal one.

Matched controls remove the gain

On the PF17041 target the structured action looks like it works, reducing ESM-2 distance by 0.0736 ± 0.0030 across five seeds against 0.0218 ± 0.0051 for edit-count-matched random mutation. That is the result a design paper would report and stop. Then composition matching erases it: reapplying the same amino-acid substitution multiset at random positions produces 0.0764, against 0.0736 for the structured action. The paired difference is −0.0027, 95% CI [−0.0063, 0.0012]. Trinucleotide context and mutation placement, the things that define a signature, add no steering power beyond the substitution composition ESM-2 already preserves.

Scorer feedback adds nothing either. An unselected-drift control scores every proposal and accepts it regardless, improving target distance by 0.0737 against 0.0736 under Metropolis selection. Whatever the loop gets, it gets from the proposal distribution rather than the feedback. More compute also leaves it where it is: at 100 ESM-2 scorer calls, five-action one-step lookahead reaches mean target distance 0.177, worse than a fixed policy at 0.165 and random selection at 0.173.

Two panels. Left, paired seed-level improvements for structured, composition-matched, and random arms, with structured and composition-matched overlapping and both above random. Right, mean final distance to target against ESM-2 scorer-call budget for fixed, random, and adaptive policies, adaptive above the others at every budget.
Figure 2. Matched controls and matched compute. (a) Each point is one paired seed; bars are means. The structured arm matches composition-matched mutation, while both beat edit-count-matched random. (b) Under equal scorer-call budgets, five-action adaptive lookahead stays farther from the target than fixed or random selection across 60 trajectories per budget.
Finding 2

Every matched control removes a result that would otherwise be reported, and the scorer's resolution, not the agent's effort, sets what is recoverable.

The policy recovers the action in a small fraction of cases

A supervised policy trained on 35,403 realized Pfam-domain transitions to predict which of the 78 signatures produced a given transition reaches 7.55% top-1 accuracy on a protein-grouped split, against 6.12% for a nearest-neighbor baseline and 1.28% for uniform chance, an edge of +0.93 points over a train-only frequency lookup (95% CI [+0.71, +1.27]). So the policy carries real information and still predicts the wrong action in 92.45% of held-out cases. Holding out complete processes drops top-1 to 0.8%, below the 1.3% baseline, and holding out SBS17a yields 0% recovery.

Uncertainty-aware interfaces return broad sets rather than a specific action. Adaptive prediction sets reach 92.2% coverage but return 66.8 of 78 actions on average, with a median set size of 78, the entire catalogue; a direct ESM-2 neighborhood reaches 82.6% coverage at mean set size 58.2. The uncertainty quantification works as intended, reporting that the evidence spreads across many actions.

Two panels. Left, mean returned set size against true-action coverage for ESM-2 neighborhood and adaptive prediction sets, both rising steeply. Right, conditional accuracy against the fraction of cases acted on under validation-frozen abstention thresholds.
Figure 3. Uncertainty-aware interfaces return broad sets. (a) Direct ESM-2 neighborhoods and adaptive prediction sets raise coverage by returning much of the 78-action catalogue; the dashed line marks 90%. (b) Validation-frozen abstention improves conditional accuracy by acting on fewer cases.

The representation still carries a signal. A probe that reads an ESM-2 displacement and names which of five processes produced it reaches 66% five-way accuracy against 20% chance, and 89% on SBS17a. Its errors follow substitution composition: SBS1 is confused with SBS2 even though their COSMIC spectra are near-orthogonal at cosine 0.02. The representation decodes a real signal and still merges processes the catalogue separates. Decoding differs from identifiability.

Confusion matrix for a probe decoding which of five mutational processes produced an ESM-2 displacement, showing 66 percent five-way accuracy with off-diagonal confusions following substitution composition.
Figure 4. Decoding differs from identifiability. A probe naming which of five processes produced a displacement reaches 66% against 20% chance, so the representation carries signal. The off-diagonal is the point: the confusions follow substitution composition rather than the catalogue's own structure.
Finding 3

The policy carries a real signal yet recovers only a small fraction of specific actions, and set-valued prediction buys coverage by returning most of the catalogue.

The objective moves proteins; an independent check tests acquisition

The optimization moves proteins measurably closer to the PF17041 centroid. Checked against a model that shares nothing with the objective, an independent Pfam profile HMM evaluated with HMMER 3 against PF17041.11, zero evaluated sequences pass the curated gathering threshold, and under a permissive threshold the structured-action sequences produce no match at all (best E = 0.34 elsewhere).

The objective is also exploitable. An agent making unconstrained single-residue edits drives cosine distance to the centroid from 0.26 to 0.009 over 200 edits, landing 5.0 standard deviations below the natural-protein reference mean of 0.194 ± 0.037 across 25 measured proteins, while ESM-2 naturalness falls alongside it from −0.41 to −0.50. The plausibility signal detects the departure, but because both the objective and that self-check live in ESM-2, a check in the same representation counts as independent evidence only from outside it.

Finding 4

Under the ESM-2 objective, proteins move toward the centroid while acquiring no PF17041-like properties, and an independent check outside that representation is what tells the two apart.

Part II · The audit battery, five failure modes

Five failure modes of a homology-free claim

WiML @ NeurIPS 2026

Deciding whether a mutation changes protein function drives variant interpretation and design, so a method that answers without a database match must carry its own evidence. Sequence-homology searches provide both a prediction and an external biological reference, but 30 to 50% of genes lack annotated homologs. Protein language models offer a database-free alternative by treating movement in embedding space as evidence of functional change; removing the database match also removes the reference that would expose a misleading score. The battery below tests whether a homology-free claim survives an adversarial audit of its data, generator, generalization, attribution, and metric, using PRISM's strongest reported result, linking SBS17a to the PF17041 domain family across two structurally distinct protein clans, as the subject.

Failure modeObserved resultWhat it changes
Data integrity1.8× hit inflationDuplicate jobs and self-pairs alter discovery counts.
GenerationSpurious UV signal; 6× temperature shiftDecoder collapse and alignment create the apparent biology.
GeneralizationAUC 0.94 in-distribution; 0.8% top-1 held outLeave-signatures-out accuracy falls below the 1.3% frequency baseline.
Metric attack5.0 s.d. closer than natural proteinsOptimizing centroid distance reduces protein naturalness.
AttributionSBS4 recovered first in 82% of trialsCorrelated COSMIC signatures recover interchangeably.

Length correction reduces the SBS17a to PF17041 enrichment from 60.7x to 18.7x while preserving its trajectory rank, #1 of 78. A composition-matched direction has cosine 0.86 ± 0.05 with the SBS17a direction and reproduces 72% of its displacement toward the PF17041 centroid. Because the direction and the centroid occupy the same ESM-2 space, independent validation needs a model outside it. The evidence supports a reproducible, composition-associated property of ESM-2 geometry rather than SBS17a as a biological driver of PF17041 function.

Finding 5

The five audit categories separate a reproducible embedding effect from the mechanistic claim attached to it, and the battery applies wherever a learned metric replaces an external reference.

Part III · Evaluating biological agents

The proxy, benchmark, and verifier an agent delegates to

AI for Science @ NeurIPS 2026Trust-AI-Eval @ NeurIPS 2026GEM-Bio @ NeurIPS 2026Verify-Agents @ NeurIPS 2026

A closed-loop agent measures only part of what it decides on. It ranks thousands of candidates with a cheap proxy, trusts that proxy because a benchmark reported a correlation, and delegates ranking to a supplied scorer only if a trust gate clears it. Each is a measurement read off a static corpus, and each fails the same way action identifiability does.

The proxy ranks below counting. On 504 measured PET-hydrolase variants across 25 studies, a zero-shot ESM-2 650M score ranks activity at a within-study Spearman of 0.01, below a hydropathy feature and well below a plain mutation count at 0.25. The expensive model is the weakest of four scorers on the property a campaign would use it for.

Within-study rank correlation for four scorers on PET-hydrolytic activity and melting temperature; the ESM-2 zero-shot score sits near zero on activity, below a mutation count.
Figure 5. A score that ranks below counting. Within-study rank correlation with measured outcomes for four scorers, on PET-hydrolytic activity (504 variants, 25 studies) and melting temperature (376 variants, 15 studies). Within study, so it is not confounded by between-study offsets.

Benchmark correlation stops predicting where campaigns act. Across 70 benchmark units in four corpora, reported correlation explains half the variance in top-10% selection value and 16% at top-1%; three of 41 ProteinGym assays give worse-than-random top-1% selections despite positive reported correlations. The same correlation means different things across families: at matched correlation, SSMuLA libraries deliver +0.265 more top-1% value than ProteinGym (p = 0.0001), so a leaderboard scalar is comparable only within a family mixture.

Reported correlation against selection value at two cutoffs; the relationship weakens at the sharper cutoff and three assays fall below zero at top-1%.
Figure 6. Correlation decays where campaigns act. Each point is one benchmark unit. The relationship weakens sharply at the tighter cutoff, and three ProteinGym assays fall below zero at top-1% despite positive reported correlations.

A trust gate certifies memorization. Across 4,920 trials, a scorer that memorizes the observed outcomes clears the gate at AUROC 0.991 and passes every trial, against 0.663 for the truthful scorer worth +1.082. Forward validation closes that surface, and a pool-side attacker that touches only unmeasured candidates still returns identical gate statistics by construction and costs 1.24 good candidates per batch. A probe catches 99.2% of naive attacks yet recovers about 10% of the lost value with a random fallback, because catching an attack and replacing it are different problems.

Three panels. Left, gate scores rise with leakage. Center, only the truthful scorer improves selection. Right, forward validation preserves the truthful scorer while cutting leaked scorers' pass rate.
Figure 7. The gate rewards the wrong signal. Left, the gate's score rises with leakage. Center, only the truthful scorer helps. Right, forward validation holds the truthful scorer's value while pushing the leaked scorers to chance, the defense that works on leakage, with the pool-side attack covered in the text.
Finding 6

The proxy, benchmark, and verifier each fail where the agent uses them, because a static number is not the quantity an optimizer pushes on.

Takeaways

What the joint report adds up to

Distinct actions become indistinguishable
Biological grounding leaves scorer-level identifiability open, and ESM-2 maps many distinct processes to nearly the same place.
Resolution sets the ceiling, not agency
Once the scorer discards a distinction, adaptive selection, lookahead, and uncertainty quantification work only with what remains.
Optimization needs independent validation
Movement under an objective comes apart from acquiring the property, and a check sharing the objective's representation leaves the two tangled.
One battery, three subjects
The five-category audit that separates effect from mechanism in a homology-free claim also exposes the proxy, benchmark, and verifier an agent leans on.
An eval earns the decision
A result should move a model or a campaign only after matched controls and an independent check survive; until then the honest output is a wider uncertainty set or an abstention.

A three-step evaluation

  1. Measure whether the scorer distinguishes the actions available to the policy.Compare distances in the action's own parameterization against distances between scored outcomes. Cheap, and diagnostic.
  2. Compare optimization against matched controls.Match proposal distribution, edit count, sequence length, and scorer-call budget. Each matched control here removed a result that would otherwise have been reported.
  3. Validate the molecular claim with a measurement independent of the objective.If the objective and the check share a representation, independence is what makes the check count as evidence.

Limitations

The identifiability conclusions apply to the action-scorer pair studied here; other actions or scorers may preserve more mechanistic information, and the recipe treats this as a measurable question rather than an assumption. The steering and decoding experiments cover a subset of signatures while the destination analysis covers all 78, so catalogue-wide overlap is established while per-signature steering specificity remains open. ESM-2 8M and 650M share training data and tokenization, so the scale comparison is a within-family check rather than an independent-model replication, and the 0.011 collapse bar is inherited from an earlier four-signature experiment. The evaluating-biological-agents results are each scoped to their own data: PET-hydrolytic activity for the proxy, four corpora with the family offset estimated from two of them for the benchmark, and an offline counterfactual evaluation on pre-selected pools for the verifier, so absolute performance there is comparative rather than a population estimate. Methods that steer protein representations are dual-use; the action space is restricted to characterized processes and the work audits failure modes rather than claiming successful molecular design.

Cite

@inproceedings{subramanian2026action, title = {Action Identifiability in Protein-Design Agents}, author = {Subramanian, Saanvi S. and Lobinska, Gabriela}, booktitle = {Workshop on Machine Learning for Molecules (MoML)}, year = {2026} } @inproceedings{subramanian2026noblast, title = {No BLAST, New Problems: Auditing Homology-Free Protein-Function Prediction}, author = {Subramanian, Saanvi S. and Lobinska, Gabriela}, booktitle = {Women in Machine Learning (WiML) @ NeurIPS}, year = {2026} } @inproceedings{subramanian2026benchmark, title = {Does Benchmark Correlation Predict Selection Quality for Protein Design?}, author = {Subramanian, Saanvi S.}, booktitle = {NeurIPS Workshops (Trust-AI-Eval; GEM-Bio)}, year = {2026}, url = {https://openreview.net/forum?id=KxSFpyRXzp} } @inproceedings{subramanian2026verifier, title = {When the Verifier Is Wrong: Red-Teaming Agent Verification Where Ground Truth Exists}, author = {Subramanian, Saanvi S.}, booktitle = {NeurIPS Verify-Agents Workshop}, year = {2026}, url = {https://openreview.net/forum?id=RyzGAFDM8S} }

← all projects