Peer Revenue Publish or profit

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Quantitative Biology Quantitative Methods

q-bio.QM Watch Quantitative Methods (sign in to use)

Wednesday 16 September

Quantitative Methods · Cross-listed from cs.LG

Sharpness-Aware Minimization (SAM) Improves Classification Accuracy of Bacterial Raman Spectral Data Enabling Portable Diagnostics

Antimicrobial resistance is expected to claim 10 million lives per year by 2050, and resource-limited regions are most affected. Raman spectroscopy is a novel pathogen diagnostic approach promising rapid and portable antibiotic resistance testing within a few hours, compared to days when using gold standard methods. However, current algorithms for Raman spectra analysis 1) are unable to generalize well on limited datasets across diverse patient populations and 2) require increased complexity due to the necessity of non-trivial pre-processing steps, such as feature extraction, which are essential to mitigate the low-quality nature of Raman spectral data. In this work, we address these limitations using Sharpness-Aware Minimization (SAM) to enhance model generalization across a diverse array of hyperparameters in clinical bacterial isolate classification tasks. We demonstrate that SAM achieves accuracy improvements of up to 10.5% on a single split, and an increase in average accuracy of 2.7% across all splits in spectral classification tasks over the traditional optimizer, Adam. These results display the capability of SAM to advance the clinical application of AI-powered Raman spectroscopy tools.
Read the full abstract

All authors Kaitlin Zareno, Jarett Dewbury, Siamak K. Sorooshyari, Hossein Mobahi, Loza F. Tadesse

Quantitative Methods

Automatic denoising and differentiation based on Savitzky-Golay filtering and Homogeneous Differentiators for attractor reconstruction via differential embedding

Differential embedding methods aim to reconstruct attractors of dynamical systems from noisy measured time series, but require accurate estimates of signal derivatives. We introduce SHADED (Savitzky-Golay and Homogeneous-differentiator based Automatic DEnoising and Differentiation), a novel methodology for denoising and estimation of derivatives up to an arbitrary order, which enables attractor reconstruction via differential embedding from noisy time series data. Homogeneous Differentiators (HD) guarantee finite-time derivative estimates in the presence of noise, while subsequent Savitzky-Golay (SG) filtering attenuates chattering. Crucially, SHADED extracts all parameters required for application of both HD and SG automatically from the data, without requiring manual tuning that may lead to inaccurate reconstruction, and can also incorporate prior knowledge, if available, thereby yielding a flexible tool for data-driven numerical differentiation of noisy signals. The obtained differential embeddings can reveal features of the underlying dynamics that are useful, e.g., for system identification, pattern recognition and discrimination between dynamic regimes; the latter application is particularly important in biomedical settings, to help distinguish between different physiological and pathological states. We demonstrate the efficacy of SHADED by testing it on computational neuroscience models, LTspice-simulated chaotic electronic circuits, and photoplethysmography and arterial blood pressure experimental recordings: across all these case studies, SHADED produces accurate derivative estimates and accurate attractor reconstructions via differential embedding (whenever a ground truth is available) or geometrically coherent and reproducible reconstructions consistent with the expected dynamics (in the absence of a ground truth), without the need for manual parameter tuning.
Read the full abstract

All authors Uros Sutulovic, Daniele Proverbio, Rami Katz, Giulia Giordano

Quantitative Methods · Cross-listed from quant-ph

Quantum Probability Current Guided Reduction of Coupling Control Degrees of Freedom for Excitation Transport

Time-dependent coherent control can enhance excitation transport in open quantum networks, but independently controlling every inter-site coupling creates a control space of high dimension and leads to difficult optimization problems. We introduce an edge-ranking strategy based on the control-induced change in the gradient component of the time-integrated quantum probability current, which is obtained via a graph Hodge decomposition. When our strategy is applied to the seven-site Fenna-Matthews-Olson (FMO) model, the six-edge set retains $99.83\%$ of the enhancement achieved by full control, and the four-edge set retains $97.60\%$ while reducing the pulse fluence—used here as a proxy for control effort—by $41.55\%$ relative to full control. Dephasing scans and comparisons with random edge sets and random networks provide numerical support for the relevance and potential broader utility of the ranking. These results show that edge selection guided by the quantum probability current can substantially reduce the control space while preserving high transport performance with lower control effort.
Read the full abstract

Tuesday 15 September

Quantitative Methods · Cross-listed from cs.AI

BioPhys-Bridge: A Benchmark for Interdisciplinary Scientific Reasoning in Physics-Grounded Biological Research

Language models face unique challenges in analyzing interdisciplinary scientific research literature. In biophysics research, faithful answers require grounding observed data in source evidence, interpreting it through a quantitative physics model, and linking it to a biological mechanism. To address this challenge, we introduce BioPhys-Bridge, a novel benchmark dataset for evidence-grounded scientific reasoning over biophysical literature. Each case contains evidence blocks, stable evidence IDs, quantitative values, units, equations, assumptions, mechanisms, and next decisions as grounding targets for question answering (QA) and retrieval-augmented generation (RAG). The initial release contains 500 cases, 1,517 agent-facing tasks, and covers six biological domains and nine physical model families, including three sparse families reserved for future expansion. We enforce strict quality gates for all cases in schema, evidence-integrity, quantitative-grounding, source-license, duplicate, unit-normalization, with domain expert review and annotation for 81 cases. Preliminary evaluations show that DeepSeek-V4-Flash obtain the highest evidence-ID $F_1$ score (0.360), followed by Qwen3.7-Max (0.316) and GPT-4o-mini (0.294). BioPhys-Bridge is an interdisciplinary benchmark for evaluating attribution, faithfulness, hallucination reduction, and biological experiment design with complex, multi-step scientific reasoning. Future works will increase the size and complexity of the dataset and perform comprehensive evaluations. Code and data are available in the GitHub repository and on Hugging Face.
Read the full abstract

Quantitative Methods

FlowLOT: Linearized Optimal Transport for Flow Cytometry Analysis

Multiparameter flow cytometry generates high-dimensional, unordered single-cell mea- surement data for disease diagnosis and monitoring, yet analysis often remains dependent on manual gating, limiting scalability and reproducibility. Existing machine-learning ap- proaches can reduce annotation burden but frequently require large training cohorts and offer limited interpretability. To address these challenges, we introduce FlowLOT , an optimal-transport-based framework that models the single-cell measurement data of a pa- tient sample as an empirical cellular distribution and maps it directly into a fixed-length feature vector. Within a single transparent architecture, FlowLOT unifies high-dimensional classification, interpretable visualization, and continuous quantitative inference. In few-shot regimes, using as few as 16 patients per class on FlowCAP-II and 8 patients per class on BLAST110, it accurately distinguishes healthy from acute myeloid leukemia (AML) sam- ples, reaching 94.3% and 98.0% balanced accuracy, respectively. The underlying embedding exposes marker-level variation driving disease-associated population shifts and enables quantitative measurable residual disease (MRD) estimation, achieving a Pearson correlation of 0.82 on held-out samples and 0.79 under cross-dataset transfer. Furthermore, at the clinically relevant 0.1% threshold for leukemia-associated immunophenotype (LAIP) residual disease, FlowLOT detects positivity with 72% sensitivity at 100% specificity. By replacing subjective manual gating and black-box deep learning with a distribution-aware framework, FlowLOT offers a sample-efficient, scalable, and interpretable solution for high- dimensional cytometry under realistic clinical and experimental constraints.
Read the full abstract

All authors Naqib Sad Pathan, Mohammad Shifat-E-Rabbi, Kristofor E. Pas, Ivan Medri, Bartek Rajwa, Gustavo K. Rohde

Quantitative Methods · Cross-listed from stat.ML

Bracketing Uncertainty in Clustering Under the Manifold Hypothesis

The manifold hypothesis suggests a natural criterion for clustering: partition data according to the manifold component from which each point is drawn. Whether two components are separable depends on a geometric tradeoff: the ambient separation between components versus the largest gap in sampling. In practice, this tradeoff is rarely assessed explicitly, leading standard methods to over-commit to a single clustering assignment even when the data do not support a unique answer. We formalize this tradeoff by combining intrinsic manifold geometry (volume growth and reach) with sample-level quantities (fill distance and density), yielding a threshold phenomenon for mutual-$k$-nearest-neighbor graphs: when the offset-to-fill ratio exceeds a conservative upper threshold, component separation is preserved; below a lower threshold, components fuse. The gap between these thresholds defines a geometric uncertainty zone in which the number of clusters is not identifiable from the data. Nevertheless, conventional approaches still seek one: sweeping parameters (an engineering approach) or fitting a generative mixture model (a model-based approach). Rather than forcing a single estimate of the number of clusters, we propose Manifold-Based Clustering (MBC), which returns an explicit bracket interval to quantify the underlying data uncertainty. This bracket acts as an empirically calibrated diagnostic: it narrows when a single resolution is supported, widens when multiple resolutions coexist, and collapses to one when no separated structure is detectable. Empirically, we find that many real datasets lie within the uncertainty zone rather than admitting one clear answer. Our results suggest that ambiguity in cluster number is often intrinsic, and should be quantified rather than resolved.
Read the full abstract

Quantitative Methods

GIA: Germline-Informed Aging with AlphaGenome Finds Genetically Regulated CpGs

Epigenetic clocks estimate age and aging-related phenotypes from DNA methylation at selected CpG sites, but the extent to which these inputs are influenced by germline genetic variation is unclear. Because methylation at many CpGs is genetically regulated, some between-person variation in clock estimates may reflect inherited genetic differences rather than aging-related change alone. Here we developed GIA (Germline-Informed Aging), a framework that maps CpGs selected from 13 published epigenetic clocks to blood methylation quantitative trait loci (meQTLs) and scores associated genetic variants with AlphaGenome. We show that clock CpGs were enriched for blood meQTLs relative to matched unused Illumina 450k probes (62.7% versus 39.6%; OR 2.57), across multiple clock families, suggesting that age-informative methylation sites are heavily influenced by germline genetic variation. Ranking by predicted chromatin effect isolated rs10190186, a cis-acting variant at FHL2 predicted to increase blood chromatin accessibility (ATAC +1.00; DNase +1.64) and FHL2 RNA (+0.30). This locus illustrates how inherited variation may shape methylation features repeatedly used by epigenetic clocks, motivating direct tests of whether such variants shift baseline clock estimates or longitudinal aging trajectories.
Read the full abstract

Quantitative Methods · Cross-listed from stat.ME

Information Set Emulation: Causal Certificates for AI Derived EHR Features

AI and large language models can recover clinically meaningful features from electronic health records (EHRs), but predictive usefulness does not establish admissibility for causal inference. We introduce information set emulation: an AI typed lift attaches source evidence, clinical and recording times, decision-time availability, representation version, proposed causal roles, and unresolved ambiguity to extracted features under a locked target trial. Causal certificates record auditable evidence for those roles. Features with unresolved downstream roles are routed to compatible reporting or separate analyses. Typed evidence defines an observational fiber of causal worlds consistent with the observed law. The locked scalar estimand maps this fiber to a compatible image whose squared Chebyshev radius equals the residual minimax mean squared error when the image is nonempty and compact. This classical identity provides a target-specific measure of information ambiguity. The contribution is its integration with a joint EHR observation map and an auditable certificate architecture. Under explicit exchangeability, positivity, and nuisance-consistency conditions, we give identification and cross-fitted augmented inverse probability weighted estimation, distinguishing empirical and population targets. An EHR compression-drift identity separates the roles of frame presence, treatment assignment, and outcome observation. Artificial simulations and a common-law finite-world example illustrate estimation failures and information-radius reduction. Synthetic Phase 0 notes demonstrate audit diagnostics; a separate role-specific analysis spread illustrates routing and is not an exact fiber radius. All experiments are synthetic. The framework specifies when reconstructed information can support a point claim and when compatible reporting is required.
Read the full abstract

Quantitative Methods · Cross-listed from stat.ME

When AI Generates Covariates: Causal Typing and Estimand Drift in Sequential Experiments

AI-generated covariates from notes, conversations, images, and wearable streams can change the causal question when their roles are left unspecified. A generated feature may represent a treatment version, pre-action state, history, design variable, mediator, outcome proxy, observation process, or intercurrent event; these roles are not interchangeable. We formulate a causal type discipline for sequential experiments: a versioned representation map, a causal role classifier, a claim-status filter, and an estimand lock. The lock fixes a standardized proximal effect before generated covariates enter the analysis. Under audit correctness and standard identification assumptions, admissible role assignments preserve this estimand. We apply the established conditional-covariance characterization of compression bias to substitution of generated representations for design-relevant states. A standardized decomposition separates compression, conditional-law, and standardization drift. Further results cover mediator adjustment, post-action leakage, marker-intervention conflation, outcome-guided discovery, and state-measurement error. Cluster-level orthogonal estimators distinguish empirical and superpopulation targets under repeated sessions and missing outcomes. Simulations show that refinement helps when it retains design-relevant information, whereas design erasure, leakage, and same-data marker selection can produce bias or undercoverage. The framework places causal semantics and claim status before confirmatory inference with generated representations.
Read the full abstract

Quantitative Methods

Decoding Extrahepatic Targeting of Lipid Nanoparticles with Interpretable Machine Learning

Lipid nanoparticles (LNPs) have transformed RNA medicine, yet their clinical utility remains constrained by predominant hepatic accumulation after systemic administration. Redirecting LNPs to extrahepatic tissues requires understanding of how lipid chemistry and formulation composition jointly govern in vivo biodistribution. Here, we develop an interpretable machine learning framework to predict hepatic versus extrahepatic LNP accumulation and identify molecular design rules for extrahepatic RNA delivery. A literature-derived dataset of 476 intravenous LNP formulations was curated from 81 studies, integrating formulation composition, lipid chemical structures, and IVIS-based biodistribution profiles. Standardized SMILES representations of ionizable lipids, helper lipids, sterols, PEGylated or polymer-conjugated lipids, additional lipids, and polymer repeat units were converted into RDKit Expert descriptors and combined with formulation-level variables to generate an 808-dimensional feature representation. Logistic regression, random forest, and XGBoost achieved ROC-AUC values of 0.839, 0.866, and 0.874, respectively. SHAP-based interpretation and consensus feature ranking revealed that ionizable-lipid descriptors dominate biodistribution prediction, while formulation composition, particularly ionizable lipid, sterol, and PEGylated/polymer-conjugated lipid fractions, contributes substantially. The top 20 consensus features retained nearly all predictive information in tree-based models. The most informative features implicated electrotopological surface properties, charge- and hydrophobicity-weighted surface areas, molecular topology, and amide/alkyl structural motifs as drivers of extrahepatic accumulation. This study establishes an interpretable, data-driven strategy for decoding LNP biodistribution and provides actionable design principles for engineering LNPs beyond the liver.
Read the full abstract

All authors Asal Mehradfar, Mohammad Shahab Sepehri, Owen Antholine, Varun Shankar, Glen S. Kwon, Salman Avestimehr, Morteza Rasoulianboroujeni

Quantitative Methods · Cross-listed from cs.SI

Graphlets as structural fingerprints of complex networks

Complex networks are often compared using selected graph-theoretical measures that capture a selected set of properties with effects ranging from local to global, such as degree, clustering or betweenness centrality. Here we introduce a structural fingerprinting framework based on graphlets: small rooted subgraphs whose distributions provide a systematic description of local-to-mesoscale topology. Across synthetic networks generated from several random graph models, graphlet fingerprints capture parameter-dependent structural differences, outperform standard graph-theoretical measures, and identify even subtle local patterns driving discrimination. We then apply the framework to empirical resting-state functional connectomes, documenting that while graphlets show superior sensitivity also to controlled topological perturbations of brain connectivity, specifically in schizophrenia-control classification they perform only comparably to classical graph-theoretical features. This is in line with the notion that schizophrenia-related alterations are dominated by spatially localized connectivity changes rather than general topological reorganization. Altogether, the generative modeling, targeted perturbations and real-world neuroimaging classification challenge position graphlets as flexible structural fingerprints of complex networks, while carefully outlining their strength and weaknesses compared to more classical graph theoretical features.
Read the full abstract

All authors Anna Pidnebesna, David Hartman, Aneta Pokorna, Daniel Trlifaj, Jaroslav Hlinka

Quantitative Methods · Cross-listed from eess.SY

Development of a 4D Cerebral Microvascular Imaging Platform for Mouse Stroke Model

In ischemic stroke, changes in cerebral hemodynamics during both the ischemic and reperfusion phases strongly influence stroke outcomes. However, these hemodynamic changes remain incompletely understood. To address this challenge, we devised an imaging platform that enables time-resolved ultrasound microvascular imaging during the experimental induction of ischemia and reperfusion in a mouse model. The platform leverages our previous ultrasound imaging framework combined with continuous mechanical scanning, which acquires whole-brain blood-flow signals within 5 s. The experiments demonstrated that the proposed platform can visualize both local and whole-brain hemodynamic responses to the induction of ischemia and reperfusion, suggesting its potential for rapid and continuous whole-brain hemodynamic assessment in small-animal models.
Read the full abstract

All authors Yoshihisa Kaneko, Moe Kumai, Hiroyuki Igarashi, Daisuke Ando, Kuniyasu Niizuma, Hidenori Endo, Yoshifumi Saijo, Takuro Ishii

Quantitative Methods

GPCR Ligand Bioactivity Prediction with Physics-Informed Dual-State Query Learning

Predicting the bioactivity profiles of small molecules against G protein-coupled receptors (GPCRs) is a challenge in drug discovery. Although deep learning has accelerated the prediction of binding affinities, existing approaches often struggle to distinguish between functional efficacies because they neglect dynamic conformational equilibria. Furthermore, structure-based methods are frequently limited by the scarcity of high-resolution active-state crystal structures and the indistinguishability of conformational states in static representations. To bridge the gap between black-box prediction and biophysical reality, we propose Dual-State Query (DSQ), a physics-informed multimodal architecture that explicitly embeds the Monod-Wyman-Changeux (MWC) model of allostery within a neural network. Unlike conventional models that rely on explicit 3D structures, DSQ utilizes learnable orthogonal queries to extract disentangled representations of active and inactive receptor states. These latent representations are governed by a novel neural MWC gating module, which mathematically derives the probability of receptor activation from thermodynamic competition between ligand-state affinities and the receptor's intrinsic conformational energy barrier. A contrastive ranking objective is also introduced to enforce differential affinity constraints, ensuring physical consistency. Extensive experiments demonstrate that DSQ outperforms other baselines, particularly for the agonist subset. Additional homology-stratified, temperature-sensitivity, perturbation, clustering, and efficiency analyses show that DSQ provides useful thermodynamic inductive bias, while also exposing a clear limitation on low-homology receptors. The code is available at https://github.com/jiankliu/DSQ.
Read the full abstract

All authors Shuo Zhang, Huifeng Zhang, Rongqi Hong, Jian K. Liu

Monday 14 September

Quantitative Methods

Optical microelectrode arrays for differential readout of electrical and mechanical signals in cardiac cells

Simultaneous assessment of electrical excitation and mechanical contraction is essential for understanding cardiac cell function, yet these two processes are commonly measured with separate techniques or invasively. Here, changes in cellular electrical activity modulate local charge redistribution in optical microelectrodes and are converted into fluorescence signals, while cell contraction induces membrane displacement that contributes an additional mechanical component to the optical readout. By comparing recordings obtained in beating cells with those acquired after inhibition of contraction, we separate action-potential-associated electrostatic transduction from contractility-driven membrane motion. The approach offers a label-free route to support high-throughput in vitro assays for cardiotoxicity screening and electromechanical sensing. Concurrently, it unfolds the physical mechanisms governing membrane-based optical devices deployed in cardiac cell bioelectronics and mechanics.
Read the full abstract

Quantitative Methods · Cross-listed from cs.LG

Knowledge-Enriched Structured EHR Features for 30-Day Hospital Readmission Prediction on MIMIC-IV

Recent approaches to 30-day hospital readmission prediction rely on pre-trained language models applied to discharge summaries. Although these methods achieve strong performance, they depend on the availability of clinical notes, incur substantial computational costs, and yield representations that lack interpretability. We propose a knowledge-enriched feature representation that augments structured Electronic Health Record (EHR) data with four medical knowledge sources: disease ontology mapping, procedure classification, drug ingredient vocabulary, and organ system laboratory aggregation, without using clinical notes. Each feature dimension corresponds to a named clinical concept, yielding a sparse and interpretable patient representation. The approach is evaluated with six classifiers on a MIMIC-IV v2.2 cohort. Under 20-fold cross-validation, the best configuration achieves an AUROC of 0.743. This performance is comparable to that of previously reported methods on this dataset, including both those using only structured data and those incorporating clinical notes, while requiring considerably less computational cost. Interpretability analysis shows that demographics, organ system labs, drug ingredient features, and first-level ontology disease categories drive prediction, while deeper hierarchy levels contribute negligibly. These findings indicate that knowledge-enriched structured features offer a competitive and efficient alternative to embeddings from clinical notes for 30-day readmission prediction.
Read the full abstract

All authors Mohamad Najafi, Hongyun Fu, Mathias Brochhausen, Jian Wu, Yaohang Li

Quantitative Methods

Multiscale modeling of host-pathogen interactions and mucociliary clearance during non-tuberculous mycobacterial pulmonary infection

Non-tuberculous mycobacterial (NTM) infections are a clinical challenge in cystic fibrosis (CF), where impaired mucociliary clearance and altered mucus rheology promote bacterial colonization despite host immune responses. Understanding how bacterial growth, immune cell dynamics, and mucus transport regulate infection progression is difficult because these processes interact across spatial and temporal scales. We develop a computational framework bridging a mechanistic agent-based model (ABM) of NTM infection with a spatially resolved partial differential equation (PDE) model. The PDE model couples bacterial proliferation, macrophage chemotaxis, immune-mediated clearance, mucus degradation, and viscoelastic transport in a two-compartment geometry representing mucus and lung tissue. Parameters are calibrated using data from the established ABM, yielding an efficient continuum representation while preserving cellular mechanisms. The PDE model reproduces bacterial and macrophage dynamics and enables analyses of mucus-related mechanisms and therapies. Sensitivity analysis identifies mucus viscosity, bacterial diffusivity, and macrophage mobility as key regulators of bacterial persistence through mucociliary clearance and tissue colonization. Simulations reveal nonlinear effects of mucolytic therapies: enhanced clearance reduces bacterial burden in mucus, whereas excessive viscosity reduction may promote migration into lung tissue, supporting combination with antibacterial treatment. This framework provides a quantitative platform for studying pulmonary infections, evaluating therapies, and developing patient-specific digital twins.
Read the full abstract

All authors Jindong Wang, Kali Konstantinopoulos, Po-Chun Kuo, Mingchao Cai, Ning Wei, Elsje Pienaar, Wenrui Hao

Quantitative Methods

TractSpLearn: Specialized Shared-Manifold Learning for Individualized Detection of Subtle White Matter Alterations in Mild Traumatic Brain Injury

Traumatic brain injury (TBI) often leads to subtle white matter damage that remains undetected on conventional MRI. Diffusion kurtosis imaging (DKI), an extension of diffusion tensor imaging (DTI), provides complementary information on non-Gaussian water diffusion and is sensitive to complex white-matter microstructure. With the advent of ultra-high-field MRI, the spatial resolution and signal-to-noise ratios (SNR) have been significantly enhanced, enabling more precise visualization of subtle abnormalities. Building on these advances, we developed TractSpLearn, an individualized tract-based learning framework that jointly considers within-group variability and between-group differences. Unlike the original TractLearn framework, which learns a normative manifold exclusively from healthy controls, TractSpLearn incorporates both healthy controls and patients to learn a shared manifold with a healthy-anchored representation and an additional patient-related component. To assess the performance of the proposed method, we compared TractSpLearn with the original TractLearn in three cohorts: (i) healthy controls (HC), (ii) athletes with persistent post-concussive syndromes (PPCS), and (iii) athletes with repeated head injuries (RHI), with abnormalities particularly evident in axial kurtosis (AK) and mean diffusivity (MD). In RHI, TractSpLearn highlighted recurrent abnormalities across patients. In the PPCS cohort, the overall group-level differences were more modest, potentially reflecting both limited statistical power due to the small sample size and partial normalization of white-matter alterations during recovery. Still TractSpLearn identified abnormality evidence in more patients and across more affected tracts than TractLearn.
Read the full abstract

All authors Jiqing Huang, Ali Al-Husseini, Yi Chen, Anna Gard, Laurent Lamalle, Mohamed Ali Bahri, Markus Nilsson, Niklas Marklund, Christophe Phillips, Evgenios N. Kornaropoulos

Quantitative Methods

Telegraph Processes with Extrinsic Fluctuations: Burst Dynamics and an Application to Domestic Cat Activity

Animal activity often consists of intermittent bursts separated by prolonged periods of inactivity. Here, we describe this behavior using a two-state telegraph process with extrinsically fluctuating transition rates. We derived analytical expressions for the stationary behavior of the system and characterized how stationary fluctuations in the effective transition rates modified the activity and residence time statistics. In particular, we show that whereas fixed transition rates lead to exponential residence-time distributions, stationary rate fluctuations generate effective Lomax distributions with heavier tails. We also found that fluctuations can either increase or decrease the mean occupancy of the active state, and that when the fluctuations in the two transition rates are equal, their effect on the stationary occupancy becomes indistinguishable from the case without fluctuations. We assessed the parameter inference using synthetic data and applied the framework to observations of a domestic cat. For synthetic data, accounting for extrinsic fluctuations substantially improved the recovery of the kinetic parameters used to generate the data. In the experimental data, the inferred variability of the inactive-to-active transition rate was substantially larger than that of the active-to-inactive rate. These results establish a tractable framework for studying stochastic switching systems in which extrinsic variability reshapes burst and residence-time statistics.
Read the full abstract

Sunday 13 September

Quantitative Methods · Cross-listed from cs.LG

An immune world model for multiscale forecasting and therapeutic hypothesis generation

Immune therapies act across cell-intrinsic programs, tissue ecosystems, and patient-specific immune states, yet most predictors address these scales separately. We used a governed evolutionary AI Scientist to construct the Immune World Model, an action-conditioned model that learns how interventions move immune states across cellular, tissue, and individual levels. The Immune World Model–building Scientist searched candidate architectures and workflows, and the resulting world model was frozen before independent confirmation. The frozen model generalized to unseen interventions and biological contexts, recovered intervention-specific cellular programs, integrated cell and tissue information to improve ecosystem and patient-response prediction, and forecast unseen perturbation combinations. Immune World Model–guided analysis then combined measured perturbations with cross-axis inference to nominate IL-36$\gamma$ plus SIRP$\alpha$ inhibition as a complementary-axis therapeutic hypothesis, whereas a governed self-correction audit rejected every screened cytokine pair. The Immune World Model provides a framework for multiscale immune simulation that connects AI Scientist-driven model construction, intervention forecasting, and the generation of prospectively testable therapeutic hypotheses.
Read the full abstract

All authors Taoyong Cui, Xi Wang, Zonghang Li, Jinchao Ding, Lingsen You, Yuzhi Xu, Wanghan Xu, Fang Wu, Kejun Ying, Wanli Ouyang, Pheng Ann Heng, Ling Yang, Zhenfei Yin, Yingcheng Wu

Quantitative Methods

MIRAGE: Measuring Interpolation and Redundancy in Affinity GEneralization

Deep learning now underpins structure-based drug design, from complex and affinity prediction to ligand ranking and pose generation. Recent co-folding models reportedly approach free-energy-perturbation accuracy at far lower cost. Yet standard evaluation, a single held-out correlation or pooled pose-success rate, cannot separate transferable binding principles from repeated exposure to related protein families in public databases, and practical success depends on genuinely novel targets. We introduce MIRAGE (Measuring Interpolation and Redundancy in Affinity GEneralization), a plug-in benchmark treating historical public family support (through 2019) as an explicit variable, applying a family-support axis to affinity and pose prediction via matched strata, family-disjoint controls, ligand-only baselines, and temporal evaluation. Co-folder affinity accuracy rises sharply with family support, while shallow controls that cannot exploit the test family stay flat, large for co-folders and near zero for every family-disjoint or trivial control. For Nesso-1 it survives covariate, conditioning, balancing, and clustering checks; Boltz-2's endpoint is limited by coverage. It localizes to family support rather than ligand chemistry, approaching a level from family identity alone. Rankings reverse on novel families, where a family-disjoint random forest leads both co-folders, significantly vs Nesso-1. On one external low-support target, neither co-folder beats molecular weight, corroborative rather than population-level evidence. gnina shows significant support dependence in rescoring whereas smina does not; MSA-free pose engines show larger gaps than smina redocking. This redundancy-driven inflation differs from conventional leakage. We propose reporting performance across family support plus excess over a support-insensitive baseline, and release MIRAGE as an installable benchmark and dataset.
Read the full abstract

Friday 11 September

Quantitative Methods

Automated Volumetric Segmentation of Microaneurysms on OCT Using Artificial Intelligence

Purpose: To develop and validate a deep learning-based method for the automated identification and volumetric segmentation of microaneurysms (MAs) in diabetic retinopathy (DR) using OCT. Participants: A total of 125 participants were enrolled, including 20 healthy eyes, 27 with mild NPDR, 30 with moderate NPDR, 30 with severe NPDR, and 18 with PDR. Methods: We obtained multiple repeated 3x3-mm scans from each participant using a commercial 120-kHz spectral-domain OCT system (Solix; Visionix/Optovue, Inc., California, USA), which were registered and averaged to generate high-definition volumes. We developed a 3D deep learning network to segment MAs volumetrically. The input to the network consists of the original OCT volume concatenated with its reflectance-inverted counterpart. Expert graders manually delineated MAs to generate annotations. We evaluated model performance at multiple levels, including voxel-level segmentation, lesion-level detection, and eye-level diagnosis. Results: In the test dataset (20 healthy, 20 DR eyes), the algorithm demonstrated high voxel-level accuracy, with F1 scores of 79.2% on single volumes and 86.1% on averaged volumes. Lesion-level detection reached an F1 score of 96.0% (single) and 97.1% (averaged), while scan-level diagnostic accuracy for MA presence was 97.5% (single and averaged). Conclusions: A deep learning-based method can accurately identify and segment MAs volumetrically on OCT, enabling quantification and characterization, and potentially aiding in the diagnosis, monitoring, and management of DR.
Read the full abstract

All authors Min Gao, Yukun Guo, Tristan T. Hormel, Jinyi Hao, Azaz Khan, Steven T. Bailey, Thomas S. Hwang, Yali Jia

Quantitative Methods · Cross-listed from cs.CV

Same Encoder, Different Winner: A Paired-View Framework for Cell Painting Encoder Evaluation

Vision encoders for Cell Painting are typically ranked by a single evaluation, commonly replicate mean average precision (mAP). We introduce CP-BG-Bench, a paired-view evaluation framework that holds the central cell fixed across four matched views (raw crop C, segmented S, and density-augmented variants CD and SD), ablating or augmenting surrounding pixels as a controlled intervention. Instantiating the framework on three datasets (JUMP-CP, RxRx1, RxRx3-core) and three encoders (DINOv3 ViT-B/16, OpenPhenom, SubCell) under four community-standard protocols (replicate mAP, scIB batch integration, CellProfiler feature prediction, cross-batch perturbation recall), we find that the four protocols rank the same encoders systematically differently, with disagreements decomposing along three axes: cell versus background, morphology versus context, and within-study versus across-batch. The largest effect: on RxRx3-core, SubCell with segmented inputs retains 94% of crop replicate mAP but only 32% of crop R@10, so the within-study signal preserved under segmentation is largely non-transferable; density augmentation recovers 84% of the within-study C-to-S gap but only 8% of the cross-batch gap. Segmented views predict CellProfiler features as well as or better than crops on two of three datasets, inverting the replicate-mAP ranking, and the C-to-S gap varies by an order of magnitude across datasets while remaining similar across encoders, indicating that background-driven gain is set by experimental design rather than by the encoder. Single-metric ranking of Cell Painting encoders is therefore sensitive to the protocol used, and protocol disagreements are interpretable as projections onto the three axes the paired-view design exposes. We will release the paired-view datasets, reconstruction pipelines, 36 trained checkpoints, aggregated embeddings, and the full evaluation suite.
Read the full abstract

All authors Tim Treis, Nikita Moshkov, Johan Fredin Haslum, Shantanu Singh, Fabian J. Theis

Quantitative Methods

Perturbational Validity for Foundation Models of Brain Dynamics: A Controlled Proof-of-Principle Simulation

Foundation models for human brain recordings are usually evaluated by signal reconstruction, future-state prediction, and transfer to downstream tasks. However, these benchmarks do not establish whether a transferred model remains valid when the system is actively perturbed. We define perturbational validity as the preservation, after limited system-specific adaptation, of the conditional distribution of future trajectories given the current state and a controlled input, and we evaluate it at three levels: time-series accuracy, dynamical-structure similarity, and responses to perturbations excluded from calibration. We demonstrate the framework in an oracle-drift simulation of stochastic bistable systems. Two otherwise identical multilayer perceptrons were trained on drift evaluations from passive or input-driven trajectories. Next, for each held-out system the shared weights were frozen and only a three-dimensional embedding was adapted, with a correctly specified cubic model fitted from scratch as comparator. With two to five system-specific evaluations, perturbational pretraining yielded lower errors in recovering controlled flow, landscape geometry, finite-run occupancy, response distributions, and dose-transition curves. The advantage was reproduced across five independent runs, persisted under full-network adaptation of the passive model, and was attributable to input excitation rather than transition-state coverage: excitation alone lowered controlled-flow error 1.94-fold relative to coverage alone, in five of five runs. The cubic model became competitive as calibration grew, showing that the benefit is specific to few-shot transfer. This controlled demonstration does not test recovery of dynamics from noisy or partially observed brain recordings. It shows why passive prediction should be complemented by prospective evaluation under controlled inputs.
Read the full abstract

All authors José C. Garcí Alanis, Sarah Alizadeh, Marco Rothermel, Bita Shariatpanahi, Mina Kheirkhah, Stefan G. Hofmann, Tim Hahn, Hamidreza Jamalabadi

Quantitative Methods

Audiovisual diarization of overlapping click trains in sperm whale (Physeter macrocephalus) vocal sparring using a three-hydrophone array

The individual attribution of sperm whale (Physeter macrocephalus) vocalizations during surface interactions constitutes a methodological challenge due to acoustic overlaps, multipath propagation, body shadowing, and indistinguishable inter-pulse intervals among similar-sized individuals. A multimodal audio-visual workflow is presented to deinterleave and attribute uncharacterized click trains produced by immature males during “vocal sparring” using a portable three-hydrophone array coupled with synchronized video. The approach combines spectro-temporal tracking, based on inter-click interval dynamics and the Constant-Q Transform, with spatial time-difference-of-arrival modeling projected onto the image plane for optical validation. Analysis of 2,655 manually validated clicks shows that purely acoustic clustering diarization errors remain low, peaking at 25.79% only during extreme temporal superpositions. Integrating the optical modality resolves residual spatial indeterminacies; although near-planar array geometry induces vertical ambiguities, the system achieves up to 100% horizontal visual concordance for primary emitters. Crucially, this framework successfully reconstructed and assigned 11 distinct, intertwined click trains totaling 882 clicks to specific focal individuals despite near-field tactile constraints. Because ethological descriptions remain incomplete without identifying the emitter, explicitly correlating these emissions with physical kinematics provides the fine-scale resolution required to define this socio-acoustic behavior.
Read the full abstract

All authors Lara Berkenbaum, Hervé Glotin, François Sarano, Walter M X Zimmer, Véronique Sarano, Olivier Adam, Pascale Giraudet

Thursday 10 September

Quantitative Methods · Cross-listed from cs.LG

MRFlow: Mechanism-Structured Residual Flow Matching for Virtual Cell Perturbation Prediction

Predicting cellular responses to chemical and genetic perturbations is a central test for AI virtual cells, yet many single-cell models reduce each condition to an average expression shift. This target is misaligned with perturbation biology: rare states, resistant subpopulations, and context-dependent effect sizes may coexist under the same intervention. MRFlow treats virtual-cell prediction as conditional transport between empirical single-cell measures. Control and perturbed cells are embedded in a variational latent state space, greedy low-cost supervision is constructed over detached minibatch states, and a conditional flow-matching vector field transports control populations toward perturbed populations. The field decomposes into a mechanism-structured affine dictionary, which recombines shared latent response operators with a low-rank condition-specific correction, and a spectrally controlled residual field for nonlinear cell-level variation. Genetic perturbations are encoded through a regulatory-graph GCN, chemical perturbations through Morgan fingerprints, and both enter the same conditional dynamics interface. Across eight gene and drug perturbation benchmarks, MRFlow attains the highest composite score while preserving competitive raw performance on differential-expression and distributional metrics. The results support latent distribution transport as a principled and empirically effective formulation for perturbation-aware virtual cell modeling.
Read the full abstract

Quantitative Methods · Cross-listed from cs.CL

HypoKG: Evidence-Disciplined Biomedical Hypothesis Generation Beyond Endpoint Knowledge

Large language models (LLMs) can generate biomedical hypotheses, but it remains unclear whether they truly reason from scientific evidence or simply produce convincing-sounding ideas. To study this, we combine three major biological databases: the Kyoto Encyclopedia of Genes and Genomes (KEGG), Rhea, and UniProt, into a unified biochemical knowledge graph and construct a benchmark of 550 paths connecting enzyme sources to rare disease endpoints, yielding 13,200 hypotheses from six LLMs under four conditions varying the biological information each model receives: source enzyme only, full biological path, or source and disease endpoint only. Hypotheses are scored using an expert-derived five-criterion rubric on a 1-5 scale per criterion. We find that models given both the source and disease endpoint often produce the highest-scoring hypotheses, showing that LLMs can generate compelling ideas from minimal information. However, these hypotheses are less grounded in the evidence. In contrast, models given the full biological path generate hypotheses more consistent with known mechanistic relationships. We call this evidence-disciplined reasoning. To confirm this effect, we shuffled intermediate path steps while keeping endpoints fixed. Evidence grounding dropped significantly (delta = -0.793, p < 0.001), confirming models genuinely used path structure during reasoning. Our findings show that knowledge graphs support hypothesis generation in two ways: they identify biological endpoint pairs absent from the literature, and their mechanistic paths guide how LLMs reason between them.
Read the full abstract

Quantitative Methods

Robust POMDP Framework for Lung Cancer Screening Problems

Lung cancer remains a leading cause of cancer mortality because many cases are diagnosed at advanced stages. Low-dose computed tomography (LDCT) screening can reduce mortality through earlier detection. Partially observable Markov decision process (POMDP) models can personalize screening by maintaining a belief over an individual's latent cancer state. However, cancer-state transition probabilities are often generated from clinical simulations and are subject to estimation error and model misspecification. We propose a robust POMDP framework using $\ell_1$-norm ambiguity sets around the nominal transition probability. The model optimizes screening decisions against the worst-case transition probability while keeping other components fixed at nominal values. Building on the piecewise-linear and convex structure of the robust value function, we adapt point-based value iteration to compute robust screening policies. We evaluate the policies using out-of-sample simulations that perturb selected cancer-progression parameters and compare them with the nominal ENGAGE policy for representative female and male heavy-smoker cohorts at age 50. Robust POMDP policies generally outperform nominal ENGAGE in mean out-of-sample quality-adjusted life-years (QALYs), with the best performance at a moderate ambiguity radius within the tested grid. Clinical analysis shows that the robust policy reduces lung cancer deaths (LCDs) in all evaluated settings for the female cohort and in most settings for the male cohort, with additional false positives (FPs). Screening-schedule analysis shows that the robust policy recommends more LDCT screens and detects more early-stage lung cancers. These findings show that incorporating transition-model uncertainty into data-driven screening models can improve out-of-sample reliability and provide more robust decision support when clinical simulation inputs are misspecified.
Read the full abstract

Quantitative Methods

Biology-in-the-loop: Amortized Adaptive Hit Discovery in CRISPR Screens

Many biological discovery problems require experiments to be selected sequentially under constrained budgets. CRISPR screening is a prominent example, as exhaustive perturbation testing is often infeasible and candidate perturbations must instead be prioritized over multiple experimental rounds. Despite the importance of this problem, existing benchmarks for adaptive hit discovery remain limited in scale and diversity. Here, we introduce AssayBench-Loop, a large-scale benchmark for adaptive hit discovery comprising 1,389 CRISPR screens across five phenotype categories. Beyond enabling systematic evaluation, its scale makes it possible to learn acquisition strategies across historical experiments. Building on this resource, we introduce AssayLoop, a sequential experimental design framework combining AssayFormer, a transformer-based amortized acquisition policy trained across historical screens to adapt from experimental feedback, with LLM-derived biological priors through an adaptive handoff. In this view, completed experiments become training data for learning how accumulated evidence should guide what to test next, while LLMs provide prior biological knowledge to seed the search. We further introduce AssayLLM, showing that the same principle can be extended directly to an LLM through task-specific post-training. On temporally held-out screens, AssayLoop achieves a 5.67-fold enrichment over random selection and recovers 27.7% of hits after assaying approximately 5% of the candidate library, outperforming existing adaptive-design methods and standalone LLMs, and AssayFormer alone. Performance improves with increasing historical training data and transfers to phenotype categories excluded from training. These results demonstrate the value of learning acquisition policies across historical experiments and combining them with broad biological priors for efficient adaptive hit discovery.
Read the full abstract

All authors Carl Edwards, Edward De Brouwer, Xiner Li, Namkyeong Lee, Ehsan Hajiramezanali, Anne Biton, Sara Mostafavi, Gabriele Scalia

Wednesday 9 September

Quantitative Methods

PyFLI: A Python Library for Simulation, Parameter Estimation, and Benchmarking in Fluorescence Lifetime Imaging

Fluorescence lifetime imaging (FLI) measures the temporal decay of fluorescence after excitation and provides quantitative information about a fluorophore's local environment and molecular interactions. Depending on the fluorophore and experimental design, lifetime can report changes associated with pH, oxygenation, cellular metabolism, and Forster resonance energy transfer. These properties make FLI useful across microscopy, biophysics, biomedical optics, and preclinical imaging, where the same type of molecular contrast can be studied across different biological scales. FLI measurements, however, are acquired with instruments that record fluorescence in different ways. Intensified charge-coupled device (ICCD) cameras, single-photon avalanche diode (SPAD) arrays, and time-correlated single-photon counting (TCSPC) systems differ in temporal sampling, data organization, detector noise, instrument response, and file format. Lifetime-estimation methods also make different assumptions about the recorded decay, while learning-based approaches require realistic training and validation data for which the underlying parameters are known. PyFLI is an open-source framework for FLI processing and standardized data simulation. It imports measurements from acquisition systems, simulates labeled data under configurable acquisition and noise conditions, and provides complementary approaches for parameter estimation. These include nonlinear least-squares fitting (NLSF), maximum-likelihood estimation (MLE), phasor analysis, rapid lifetime determination (RLD), Laguerre-based estimation, and optional Bayesian and deep-learning inference. CPU and GPU processing support image-scale analysis, while reconstruction, visualization, statistical analysis, and cross-software comparison provide tools for evaluating results. PyFLI includes compressed-sensing reconstruction for single-pixel hyperspectral FLI.
Read the full abstract

All authors Vikas Pandey, Ismail Erbas, Margarida Barroso, Stefan Radev, Xavier Intes

Quantitative Methods · Cross-listed from eess.IV

Seamless Whole Slide Label-Free Virtual Staining

Label-free virtual staining offers a compelling, non-destructive alternative to standard histopathology; however, its clinical adoption is hindered by the computational bottlenecks inherent to processing gigapixel Whole Slide Images (WSIs). Current deep learning approaches require patch-based inference to avoid memory constraints, which disrupts global tissue continuity and introduces tiling artifacts–displaying visible seams and color shifts. To address this, we introduce the Consistency Memory Bank (COMB), a novel label-free virtual staining framework that enforces spatial and channel consistency across tiles without memory bottlenecks. COMB decouples context storage from computation, utilizing a dynamic retrieval mechanism to fetch feature representations from adjacent tiles. This enables a retrieval-based context integration strategy that adopts local padding to resolve spatial discontinuities and neighbor-aware channel attention to stabilize statistical drift. Further optimized with a sliding window schedule to ensure minimal memory overhead, our method demonstrates superior performance over state-of-the-art baselines, achieving significant improvements in both perceptual fidelity and tiling consistency, while suggesting its downstream utility in tumor segmentation. Code is available at https://github.com/dou0000/COMB.
Read the full abstract

All authors Dou Hoon Kwark, Kianoush Falahkheirkhah, Ji-hun Oh, Shirui Luo, Volodymyr Kindratenko, Rohit Bhargava