Towards Autonomous Mechanistic Reasoning in Virtual Cells
Organizations: Korea Advanced Institute of Science and Technology (KAIST) · Valence Labs · Recursion · University College London
Abstract
Large language models (LLMs) have recently gained significant attention as a promising approach to accelerate scientific discovery. However, their application in open-ended scientific domains such as biology remains limited, primarily due to the lack of factually grounded and actionable explanations. To address this, we introduce a structured explanation formalism for virtual cells that represents biological reasoning as mechanistic action graphs, enabling systematic verification and falsification. Building upon this, we propose VCR-Agent, a multi-agent framework that integrates biologically grounded knowledge retrieval with a verifier-based filtering approach to generate and validate mechanistic reasoning autonomously. Using this framework, we release VC-TRACES dataset, which consists of verified mechanistic explanations derived from the Tahoe-100M atlas. Empirically, we demonstrate that training with these explanations improves factual precision and provides a more effective supervision signal for downstream gene expression prediction. These results underscore the importance of reliable mechanistic reasoning for virtual cells, achieved through the synergy of multi-agent and rigorous verification.
Figures & tables
| Format | Verifier | ||||
| Model | Validity | Verifiability | Drug-target interaction | Differential expression | |
| Open | DeepSeek-R1-8B | 0.003 0.004 | 1.000 0.081 | 0.000 0.000 | 0.000 0.000 |
| Qwen3-30B | 0.860 0.016 | 0.965 0.007 | 0.528 0.004 | 0.272 0.030 | |
| Llama3.3-70B | 0.841 0.046 | 0.887 0.017 | 0.322 0.034 | 0.090 0.010 | |
| Closed | Claude-Sonnet-4 | 1.000 0.000 | 0.867 0.008 | 0.657 0.005 | 0.504 0.064 |
| VCR-Agent (Ours) | 1.000 0.000 | 0.945 0.011 | 0.725 0.017 | 0.528 0.060 | |
Appendix figures & tables9 assets
Supplementary material from the paper’s appendix.
Appendix
| Category | Sub-category | Action name | Arguments |
| system initialization | set_context | set_context ({cell_type, genotype, disease, prior_perturbation, extras}) | |
| metabolic | converts_substrate | converts_substrate (id, enzyme, substrate, product, {via, confidence}) | |
| regulation | activity modulation | modulates_molecule_activity | modulates_molecule_activity (id, target, direction, {via, confidence}) |
| modulates_pathway_activity | modulates_pathway_activity (id, pathway, direction, {via, confidence}) | ||
| protein regulation | modulates_complex | modulates_complex (id, members, complex, direction, {stoichiometry, via, confidence}) | |
| post_translational_modification | post_translational_modification (id, protein, mod_type, site, direction, {via, confidence}) | ||
| Task | Model | C32 | HOP62 | HepG2/C3A | Hs 766T | PANC-1 | Average | Union |
| Differential expression | Random | 0.269 | 0.254 | 0.232 | 0.239 | 0.301 | 0.259 | 0.286 |
| Mean | 0.405 | 0.419 | 0.345 | 0.465 | 0.330 | 0.393 | 0.401 | |
| K-neighbor | 0.459 | 0.419 | 0.343 | 0.453 | 0.355 | 0.406 | 0.403 | |
| STATE (ST) | 0.273 | 0.280 | 0.199 | 0.219 | 0.313 | 0.257 | 0.251 | |
| Zero-shot | 0.363 | 0.423 | 0.316 | 0.332 | 0.399 | 0.367 | 0.338 | |
| SFT | 0.344 | 0.353 | 0.163 | 0.363 | 0.236 | 0.292 | 0.291 |
| Format | Verifier | |||
|---|---|---|---|---|
| Retrieval source | Validity | Verifiability | Drug-target interaction | Differential expression |
| PubMed | 1.000 0.000 | 0.300 0.099 | 0.432 0.059 | 0.435 0.050 |
| StarkPrimeKG | 1.000 0.000 | 0.296 0.089 | 0.382 0.031 | 0.449 0.054 |
| Harmonizome | 1.000 0.000 | 0.297 0.092 | 0.399 0.013 | 0.457 0.051 |
| Wikipedia | 1.000 0.000 | 0.292 0.091 | 0.469 0.028 | 0.485 0.063 |
| All (Ours) | 1.000 0.000 | 0.601 0.009 | 0.766 0.008 | 0.487 0.054 |
| Format | Verifier | |||
|---|---|---|---|---|
| Backbone | Validity | Verifiability | Drug-target interaction | Differential expression |
| GPT-4.1 | 0.970 0.002 | 0.299 0.105 | 0.399 0.063 | 0.527 0.076 |
| Gemini-2.5-Flash | 0.060 0.016 | 0.963 0.078 | 0.338 0.036 | 0.207 0.036 |
| Claude (Ours) | 1.000 0.000 | 0.601 0.009 | 0.766 0.008 | 0.487 0.054 |
| Format | Verifier | |||
|---|---|---|---|---|
| Step | Validity | Verifiability | Drug-target interaction | Differential expression |
| One-step | 0.999 0.002 | 0.534 0.067 | 0.329 0.026 | 0.267 0.051 |
| Two-step (Ours) | 1.000 0.000 | 0.601 0.009 | 0.766 0.008 | 0.487 0.054 |
| Format | Verifier | |||
|---|---|---|---|---|
| Model | Validity | Verifiability | Drug-target interaction | Differential expression |
| DeepSeek | 0.046 0.051 | 0.711 0.269 | 0.185 0.099 | 0.246 0.145 |
| Qwen3 | 0.937 0.055 | 0.340 0.173 | 0.081 0.056 | 0.556 0.152 |
| Llama3.3 | 0.789 0.105 | 0.468 0.227 | 0.079 0.098 | 0.231 0.146 |
| Two-step (Ours) | 1.000 0.000 | 0.601 0.009 | 0.766 0.008 | 0.487 0.054 |
| Scientific accuracy | Logical consistency | Mechanistic clarity | |
|---|---|---|---|
| Without verifiers | 0.641 | 0.720 | 0.725 |
| +DE | 0.642 | 0.721 | 0.725 |
| +DTI | 0.643 | 0.721 | 0.725 |
| +DE&DTI | 0.645 | 0.724 | 0.725 |
| Scientific accuracy | Logical consistency | Mechanistic clarity | |
|---|---|---|---|
| Human–LLM correlation ( ) | 0.72 | 0.69 | 0.76 |
| Average LLM-judge score (full) | 0.65 | 0.72 | 0.73 |
Scientific Reasoning in Language Models
Scientific reasoning in language models studies how AI systems use evidence, domain knowledge, and tools to solve problems and support discovery.
See all 144 papers