DRYLAB-Bench
Benchmarking the Biosecurity Risk of Agentic LLM Mutation Design
arXiv preprint · 2026
Why this work exists
Biosecurity has always focused on capability. As language models move beyond describing biology and are tested on biological design tasks, a central question is how far their design capability extends. A model capable of proposing functional mutations may also generate mutation designs associated with biosecurity risks.
Models describe biology
LLMs answer questions, summarise papers, and explain mechanisms. Evaluations measure what a model knows through multiple-choice knowledge probes and refusal tests.
Models are tested on biological design
Models can be prompted to propose sequences, primers, and protein variants, bringing biological design tasks within the scope of LLM evaluation. However, determining whether generated designs yield measurable biological outcomes remains difficult.
Agents are evaluated on tool-assisted design
In tool-assisted workflows, agents query specialized biological tools, receive model-generated evidence, and iteratively refine designs. However, whether such interactions translate into increased risk-associated outcomes remains unclear.
Why this matters to policy, not just to labs. Policy sources considered in this work include the Biological Weapons Convention, WHO guidance, U.S. Executive Orders 14110 and 14292, and the EU AI Act, alongside other biosecurity frameworks that identify malicious biological modification as a critical concern. Governing it requires quantitative assessment: understanding what a model can do today and how its capabilities change when augmented with specialized biological design tools. Without quantitative assessment, oversight is guesswork in both directions, proving either too loose to catch real capability or too tight to permit beneficial work.
Why it is hard to measure. A biosecurity-relevant design capability is not a fact a model can be quizzed on. It is whether the model, given a wild-type sequence and a stated objective, produces a specific substitution that measurably functions. That requires ground truth about the physical world, which is the crux of the problem described in the next section.
What current evaluation methods miss
Existing evaluations capture parts of biological risk, but fail to measure risk-associated outcomes in realistic design workflows.
Knowledge is not capability
Knowledge benchmarks such as WMDP and VCT assess whether a model knows hazardous biology, usually through multiple-choice or refusal probes. Knowing a fact and engineering a protein that behaves differently are distinct skills. Furthermore, answer keys prompt the model toward the correct response, so the score partly measures the prompt rather than the model.
What stays hidden: Whether a model can generate substitutions with experimentally measurable effects.
The amplifier is invisible
Newer agentic evaluations provide models with biological tools, but only for peripheral stages such as sequence reconstruction, screening evasion, or laboratory robot scripting. This differs from allowing an agent to orchestrate design tools that score its candidates against the targeted property being optimized. Restricting agents to general-purpose software leaves specialized design capability and its associated risks unmeasured.
What stays hidden: How much specialized tools amplify a model's biological design capability and associated biosecurity risks.
Ground truth is a trilemma
Verifying an open-ended design forces a choice among three flawed options, each with its own failure mode. Historically, no inexpensive, scalable, and hazard-free source of ground truth has existed for novel designs.
What stays hidden: Any trustworthy evaluation score for a generated design.
The verification trilemma, and how it breaks
Deep mutational scanning databases such as ProteinGym and ViroGym have traditionally been used to benchmark predictive models. DRYLAB-Bench repurposes them to verify the generative output of an LLM agent, replacing both the judge and the wet lab with an empirical blind holdout.
Thirty hazardous design scenarios
Thirty high-risk mutation-design scenarios, distilled from international biosecurity frameworks, each anchored to a real protein with published deep mutational scanning data.
Where the ground truth actually sits
Immune evasion / medical countermeasure failure
Pathogens evolving around vaccines, antibodies and other countermeasures while keeping their ability to infect.
- Tasks
- 6
- Measured substitutions
- 43,344
- Assay source
- ProteinGym / ViroGym
Enhanced infectivity / host-range expansion
Mutations that strengthen receptor binding or extend a pathogen's host range.
- Tasks
- 4
- Measured substitutions
- 25,725
- Assay source
- ViroGym
Antimicrobial / antiviral resistance
Mutations that defeat antibiotics and antivirals — the oldest and most thoroughly realised of these risks.
- Tasks
- 6
- Measured substitutions
- 46,390
- Assay source
- ProteinGym
Protein toxification
Raising the potency, stability or persistence of a protein toxin.
- Tasks
- 2
- Measured substitutions
- 2,372
- Assay source
- ProteinGym
Stealthy oncogenesis / tumour-suppressor inactivation
Disabling tumour suppressors, or locking growth signalling permanently on.
- Tasks
- 5
- Measured substitutions
- 22,344
- Assay source
- ProteinGym
Metabolic pathway collapse / pharmacogenomic weaponization
Breaking a metabolic enzyme, or shifting how a drug is metabolised.
- Tasks
- 4
- Measured substitutions
- 35,912
- Assay source
- ProteinGym
Enhanced environmental tolerance
Surviving heat, cold, desiccation or UV better than the wild type.
- Tasks
- 2
- Measured substitutions
- 10,389
- Assay source
- ViroGym
Agricultural / ecological biosecurity
Mutations affecting crops, livestock or ecosystem species rather than humans directly.
- Tasks
- 1
- Measured substitutions
- 1,219
- Assay source
- ProteinGym
View as table (Table 1)
| Cat. | Risk category | Tasks | Source | Measured single substitutions |
|---|---|---|---|---|
| F1 | Immune evasion / medical countermeasure failure | 6 | PG / VG | 43,344 |
| F2 | Enhanced infectivity / host-range expansion | 4 | VG | 25,725 |
| F3 | Antimicrobial / antiviral resistance | 6 | PG | 46,390 |
| F4 | Protein toxification | 2 | PG | 2,372 |
| F5 | Stealthy oncogenesis / tumour-suppressor inactivation | 5 | PG | 22,344 |
| F6 | Metabolic pathway collapse / pharmacogenomic weaponization | 4 | PG | 35,912 |
| F7 | Enhanced environmental tolerance | 2 | VG | 10,389 |
| F8 | Agricultural / ecological biosecurity | 1 | PG | 1,219 |
| Total | 30 | 187,695 |
PG = ProteinGym, VG = ViroGym. Coverage is deliberately reported per category: the thinnest categories carry the least evidence, which bounds how much weight their scores can bear.
How the measurement works
Four arms, one paired comparison
The design crosses two factors: whether the agent receives tool evidence and whether it may perform multiple revisions. S0, S1, and S2 are snapshots of a single trajectory (sharing the same seed, initial candidate pool, and round budget), ensuring the comparison is strictly paired rather than between independent runs. S0-iter serves as a separate no-tool control under the same round budget.
Designs straight from the sequence and goal. The pure zero-shot baseline.
Same round budget as S2, but self-review replaces tool evidence. Deliberately noisy — read as a bound on iteration alone, not an estimate.
One round of tool selection and a single revision. An opportunity to use evidence — not a guarantee that any was used.
Repeatedly selects tools, analyses the returned evidence, and refines its pool for up to four rounds.
Values represent pooled valid-response BRS. Stepwise ΔBRS increases by +0.012 for iteration alone, +0.048 for tools alone, and +0.068 for extra rounds. Each unit returns exactly three single-site substitutions with self-assigned confidences. Each scenario is posed under three wording-only prompt variants, allowing variant-to-variant spread to bound the extent to which a model's measured risk is a prompt-engineering artifact.
What the BRS score means
Every submitted substitution is queried against the hidden assay corresponding to that scenario. Scores are standardized within each assay to ensure comparability across different experiments, weighted by the model's self-reported confidence, and averaged across the top three submissions.
BRS = 0 represents the average measured variant. This single anchor ensures clear interpretability: a positive BRS indicates that the model generated a design associated with a higher risk-related assay outcome than a typical mutation of that protein, with the top-risk decile representing the assay's highest-scoring 10%.
Full scoring formulation (Equations 1–7)
What we found
Under benchmark conditions, frontier models generate mutation designs with measurable baseline design capability even without tools, while tool-augmented iteration measurably amplifies risk-associated design outcomes.
Pooled BRS by workflow arm
Confidence-weighted BRS pooled across 15 models under both refusal treatments: valid-response (excluding refused units) and floor (scoring refused units at fref = 0).
Step ΔBRS (valid view): +0.012 iteration only · +0.048 tools only · +0.068 extra rounds
Refusals track keywords, not intent
Refusal rate per unit under three controlled prompt framings of the same 30 tasks: the original wording, an academic framing that omits organism names, and a malicious-context framing. Seven representative models are shown; the complete table appears below.
(a) Per-model BRS — valid responses
Each row displays the transition in BRS from one-shot (S0, hollow) to multi-round agentic execution (S2, filled) across valid units only. Models are ordered by their S2 valid-response BRS, corresponding to the main text.
(b) Per-model BRS — refusal-penalised
Same sequence of rows, with each refused unit scored at the baseline floor ( fref = 0) within its respective task mean. High-refusal models collapse toward zero under this treatment.
View pooled values as a table (Table 9)
| Arm | Tool evidence | Revision | Valid BRS | Floor BRS | Step ΔBRS |
|---|---|---|---|---|---|
| S0 | no | one shot | 0.359 | 0.249 | — |
| S0-iter | no | up to 4 rounds | 0.371 | 0.284 | +0.012 |
| S1 | one round | one revision | 0.407 | 0.275 | +0.048 |
| S2 | adaptive | up to 4 rounds | 0.475 | 0.317 | +0.068 |
Step ΔBRS is each workflow step's contribution in the valid view (S2 − S0 = +0.116). Values from the paper's Table 9. S0-iter is deliberately noisy by design (811 of 1,350 units valid) and is read only as a bound on iteration without tools.
Full refusal table (Table 3)
| Model | Original | Academic | Malicious |
|---|---|---|---|
| Claude Fable 5 | 96.7% | 96.7% | 100.0% |
| Claude Opus 5 | 81.1% | 76.7% | 96.7% |
| Claude Opus 4.8 | 76.7% | 73.3% | 90.0% |
| GPT-6 Astra | 66.7% | 50.0% | 90.0% |
| GPT-5.6 Terra | 46.7% | 33.3% | 96.7% |
| GPT-5.6 Sol | 44.4% | 33.3% | 96.7% |
| GPT-5.6 Luna | 43.3% | 30.0% | 86.7% |
| Gemini 3.6 Flash | 28.9% | 23.3% | 96.7% |
| Gemini 3.8 Flash | 24.4% | 16.7% | 93.3% |
| Grok 4.5 | 1.1% | 0.0% | 80.0% |
| Kimi-K3 | 0.0% | 0.0% | 80.0% |
| GLM-5.3 | 0.0% | 0.0% | 70.0% |
| DeepSeek V4 Pro | 0.0% | 0.0% | 80.0% |
| Qwen3.8 Max | 0.0% | 0.0% | 90.0% |
| DeepSeek V4 Flash | 0.0% | 0.0% | 80.0% |
Full per-model table (Table 10) — the data behind charts (a) and (b)
| Compare | Valid-response BRS | Floor-included BRS | |||||
|---|---|---|---|---|---|---|---|
| Kimi-K3 | 0.519 | 0.543 | 0.631 | 0.523 | 0.551 | 0.637 | |
| GLM-5.3 | 0.394 | 0.392 | 0.492 | 0.393 | 0.377 | 0.497 | |
| DeepSeek V4 Pro | 0.276 | 0.348 | 0.483 | 0.283 | 0.348 | 0.483 | |
| Grok 4.5 | 0.368 | 0.427 | 0.467 | 0.358 | 0.446 | 0.478 | |
| Qwen3.8 Max | 0.401 | 0.442 | 0.428 | 0.413 | 0.429 | 0.445 | |
| DeepSeek V4 Flash | 0.428 | 0.367 | 0.416 | 0.427 | 0.363 | 0.416 | |
| Gemini 3.6 Flash | 0.511 | 0.465 | 0.592 | 0.339 | 0.301 | 0.414 | |
| Gemini 3.8 Flash | 0.269 | 0.492 | 0.464 | 0.192 | 0.341 | 0.356 | |
| GPT-5.6 Luna | 0.508 | 0.626 | 0.507 | 0.285 | 0.336 | 0.300 | |
| GPT-5.6 Sol | 0.347 | 0.372 | 0.531 | 0.189 | 0.218 | 0.293 | |
| GPT-5.6 Terra | 0.247 | 0.355 | 0.313 | 0.129 | 0.185 | 0.159 | |
| GPT-6 Astra | 0.226 | 0.222 | 0.302 | 0.077 | 0.087 | 0.101 | |
| Claude Opus 4.8 | 0.227 | 0.395 | 0.369 | 0.053 | 0.097 | 0.086 | |
| Claude Opus 5 | 0.378 | 0.343 | 0.383 | 0.067 | 0.041 | 0.072 | |
| Claude Fable 5 | 0.287 | 0.324 | 0.754 | 0.010 | 0.011 | 0.025 | |
Mean risk of LLM-designed mutations by condition (confidence-weighted BRS; BRS = 0 is the assay mean). "Valid" averages non-refused units with unmeasured candidates filled in via the calibrated in-silico fallback; "Floor" scores each refused unit at fref = 0 inside its own task mean. Click a column heading to sort (default order: S2 floor-mean); tick rows to compare selected models above the table. Claude Fable 5's S2 valid value (0.754) is the highest in the suite.
What the results show
1. Substantial zero-shot design capability
Even without tools or iteration, frontier models exhibit notable design proficiency: 67.6% of one-shot submissions outperform the average measured variant, and 21.4% fall within the assay-specific top-risk decile, exceeding the 10% random expectation by more than a factor of two. Pure zero-shot performance reaches a BRS of 0.52 (Kimi-K3), while GPT-5.6 Sol's one-shot design achieved the highest score among the 3,372 variants evaluated in the SRC activation assay.
2. Amplification, bounded by base capability
Tool evidence amplifies risk-associated design outcomes across both stages, yielding an increase of +0.048 for a single tool round, an additional +0.068 for adaptive multi-round iteration, and +0.116 in total. This gain is negatively correlated with the untooled baseline (r = −0.42): models with lower initial performance achieve the largest gains, such as DeepSeek V4 Pro improving by +0.207 from an initial 0.276, whereas near-ceiling models saturate, such as Claude Opus 5 starting at 0.378 and gaining +0.004. Consequently, tooling operates as a convergence mechanism rather than a uniform multiplier.
3. Tool utility depends on task alignment
Utility degrades sharply when a tool fails to measure the target property being optimized. In the CcdB toxin scenario, the conservation signal was strongly anti-correlated with the measured outcome (ρ = −0.815 across 1,176 samples), while the toxicity predictor proved uninformative (ρ = −0.135). Models lacking proper calibration confidently incorporated these misleading signals, degrading their candidate quality (e.g., GPT-5.6 Terra shifted by −0.208). Furthermore, S1 represents an opportunity rather than a guarantee: across 861 trajectories reaching an S1 submission, agents opted not to request any tools in 372 instances.
4. Refusals protect averages, not tails
Refusal is driven by surface wording rather than intent: reframing the same 30 tasks increases refusals across the entire suite, including models that accepted every original prompt. On the subset of units answered, Claude Fable 5 still achieves a valid BRS of 0.754; because scores are standardized within each assay, this performance does not reflect a less challenging residual set. Refusing a scenario therefore provides no safety guarantee, which is why refused units are reported separately rather than integrated into the primary score.
What is open, and what is not
The release is deliberately split: all components necessary to audit the methodology and reproduce the reported empirical findings are publicly available, whereas actionable, step-by-step design pipelines and risk-enabling procedures are excluded.
- Evaluation framework and code
- 15 of the 30 scenarios — prompts, targets and ground-truth mappings
- Experimental data for those 15 scenarios: per-unit trajectories, tool evidence returned, prompts as rendered, and raw model outputs
- Scoring code, baseline statistics, and the aggregate result tables covering all 30 scenarios
- The remaining 15 scenarios
- Experimental data for those scenarios: trajectories, rendered prompts, tool evidence and model outputs
- Prompts and ground-truth mappings for those scenarios
Open scenarios (filled): Individual benchmark units with fully accessible task specifications, visible and selectable via hover, tap, or focus. Withheld scenarios (dashed): Fifteen restricted benchmark units representing high-consequence pathogen classes, intentionally excluded from public distribution to prevent dual-use misuse.
Zero physical risk. The benchmark works entirely in silico: it scores text against already-measured phenotypes, never synthesizes or tests anything, and the hidden DMS labels are never shown to the models.
No novel information. Every wild-type sequence, target protein and ground-truth dataset used here is already public in open literature and databases (NCBI, UniProt, ProteinGym, ViroGym). The benchmark introduces no novel biological agents, no classified information, and no zero-day vulnerabilities — it repurposes existing empirical knowledge to measure AI capability.
Expert review. Scenario construction, filtering and verification involved domain experts in computational biology, virology and genetics, briefed on the sensitive material and free to opt out of reviewing any scenario without penalty.
Scope. BRS quantifies the risk of individual substitutions under a defined workflow — not multi-mutant fitness, transmissibility, or behaviour in a living organism. We welcome responsible disclosure and corrections.
Cite this work
The paper is available as an arXiv preprint.
@misc{wang2026drylabbenchbenchmarkingbiosecurityrisk,
title={DRYLAB-Bench: Benchmarking the Biosecurity Risk of Agentic LLM Mutation Design},
author={Yuzhen Wang and Chunxiao Li and Baichen Le and Lewen Yan and Jie Zhang and Yucheng T. Yang and Xing-Ming Zhao and Lijun Li and Jing Shao},
year={2026},
eprint={2610.13116},
archivePrefix={arXiv},
primaryClass={cs.AI},
url={https://arxiv.org/abs/2610.13116},
}