BLINK: Behavioral Latent Modeling of NK Cell Cytotoxicity
Authors: Iman Nematollahi, Jose Francisco Villena-Ossa, Alina Moter, Kiana Farhadyar, Gabriel Kalweit, Abhinav Valada, Toni Cathomen, Evelyn Ullrich, Maria Kalweit
Organizations: Department of Computer Science, University of Freiburg, Freiburg, Germany · Institute for Transfusion Medicine and Gene Therapy, University Medical Center Freiburg, Freiburg, Germany · Goethe University, Department of Pediatrics, Experimental Immunology and Cell Therapy, Frankfurt am Main, Germany · IMBIT//BrainLinks-BrainTools, University of Freiburg, Freiburg, Germany · Collaborative Research Institute Intelligent Oncology (CRIION), Freiburg, Germany
Machine learning models of cellular interaction dynamics hold promise for understanding cell behavior. Natural killer (NK) cell cytotoxicity is a prominent example of such interaction dynamics and is commonly studied using time-resolved multi-channel fluorescence microscopy. Although tumor cell death events can be annotated at single frames, NK cytotoxic outcome emerges over time from cellular interactions and cannot be reliably inferred from frame-wise classification alone. We introduce BLINK, a trajectory-based recurrent state-space model that serves as a cell world model for NK-tumor interactions. BLINK learns latent interaction dynamics from partially observed NK-tumor interaction sequences and predicts apoptosis increments that accumulate into cytotoxic outcomes. Experiments on long-term time-lapse NK-tumor recordings show improved cytotoxic outcome detection and enable forecasting of future outcomes, together with an interpretable latent representation that organizes NK trajectories into coherent behavioral modes and temporally structured interaction phases. BLINK provides a unified framework for quantitative evaluation and structured modeling of NK cytotoxic behavior at the single-cell level.
Figures & tables
Figure 1: Overview of BLINK : (a) Multi-channel fluorescence microscopy captures NK cells, tumor cells, apoptosis, and morphology. (b) Segmentation and tracking yield NK-centered interaction trajectories. (c) BLINK encodes these sequences into a recurrent latent state-space model that captures interaction dynamics under partial observability, supports latent rollouts for future cytotoxic outcomes, and predicts NK-induced apoptosis increments that accumulate monotonically. The learned latent space organizes NK behaviors into coherent modes.
Model
MAE ↓
RMSE ↓
Corr ↑
Within ±1 (%) ↑
F-MAE 30 ↓
Zero
1.28 ± 0.16
1.74 ± 0.15
0 ± 0.0
54.3% ± 7.0%
0.12 ± 0.06
Mean
1.04 ± 0.07
1.18 ± 0.08
0 ± 0.0
49.6% ± 6.3%
0.24 ± 0.05
FrameAE
0.95 ± 0.11
1.14 ± 0.13
0.32 ± 0.07
64.9% ± 6.8%
X
GRU-regress
1.25 ± 0.14
1.72 ± 0.14
0 ± 0.0
55.7% ± 6.9%
0.12 ± 0.06
GRU-monotone
0.74 ± 0.09
1.04 ± 0.11
0.57 ± 0.04
71.9% ± 3.3%
0.22 ± 0.04
BLINK-no-action
0.80 ± 0.06
1.14 ± 0.09
0.61 ± 0.04
69.4% ± 7.3%
0.09 ± 0.01
Table 1: Track-level cumulative cytotoxic outcome prediction on the held-out test set, showing improvements of BLINK across error and forecasting metrics.
Figure 2: Real interaction trajectory (top) and world model-decoded latent trajectory (bottom). Predicted and ground truth apoptosis increments align.
Figure 3: Latent behavioral structure of NK trajectories. (a) UMAP of training window embeddings clustered into four modes. (b) Test tracks projected into the embedding. (c) State transition matrix showing temporal mode progression.
Understanding non-genetic determinants of cell fate is critical for developing and improving cancer therapies, as genetically identical cells can exhibit divergent outcomes under the same treatment conditions. In this work, we present a deep learning approach for cell fate prediction from raw long-term live-cell recordings of cancer cell populations under chemotherapeutic treatment. Our Transformer model is trained to predict cell fate directly from raw image sequences, without relying on predefined morphological or molecular features. Beyond classification, we introduce a comprehensive explainability framework for interpreting the temporal and morphological cues guiding the model's predictions. We demonstrate that prediction of cell outcomes is possible based on the video only, our model achieves balanced accuracy of 0.94 and an F1-score of 0.93. Attention and masking experiments further indicate that the signal predictive of the cell fate is not uniquely located in the final frames of a cell trajectory, as reliable predictions are possible up to 10 h before the event. Our analysis reveals distinct temporal distribution of predictive information in the mitotic and apoptotic sequences, as well as the role of cell morphology and p53 signaling in determining cell outcomes. Together, these findings demonstrate that attention-based temporal models enable accurate cell fate prediction while providing biologically interpretable insights into non-genetic determinants of cellular decision-making. The code is available at https://github.com/bozeklab/Cell-Fate-Prediction.
Florian Bürger, Martim Dias Gomes, Adrián E. Granada +2
Institute for Biomedical Informatics, Faculty of Medicine and University Hospital Cologne, University of Cologne, Cologne, Germany · Faculty of Mathematics and Natural Sciences, University of Cologne, Cologne, Germany · Cologne Excellence Cluster on Cellular Stress Responses in Aging-Associated Diseases (CECAD), University of Cologne, Cologne, Germany +2
Dead cell detection in fluorescence microscopy is essential for phenotypic screening and drug efficacy assessment. However, automated identification remains challenging due to signal variability, optical artifacts, and colocalization ambiguities across multiple channels. Markers that define nuclear boundaries do not neatly coincide spatially or temporally with apoptotic markers. Furthermore, marker accumulation kinetics vary across individual cells and apoptosis stages, making universal thresholds unreliable. We present a systematic comparison of several computational approaches for detecting apoptotic cells in nucleus-stained specimens: (1) heuristic methods using intensity thresholds, distributions, correlations, and domain guidance; (2) classification using numerical features (Logistic Regression, Random Forest, XGBoost and LightGBM) and image features (ResNet-18 and EfficientNet-B0); and (3) anomaly detection, also using numerical features (Isolation Forest, one-class SVM and PCA reconstruction) and image features (convolutional autoencoders and Deep SVDD). Using a dataset of 258 images containing 65,527 cells/nuclei across multiple scenes with 243 manually annotated ground truth dead cells, we evaluate each method on metrics such as PR-AUC, Precision, Recall, F1-score, and macro-F1. Our results indicate that classification using image features (deep learning) and numerical features (traditional) achieved the best performance, with ResNet-18 and XGBoost significantly outperforming manually designed heuristics and anomaly detection methods. Our findings demonstrate that such biomedical problems remain highly challenging to solve even with the state-of-the-art computational and learning techniques. We aim to provide an effective decision-support tool to researchers while tackling such problems either fully-automatically or in a human-in-the-loop setting as an aide to domain scientists.
Prateek Verma, Benjamin Butter, Jian Zhang +1
Department of Electrical Engineering and Computer Science University of Arkansas, Fayetteville, Arkansas, USA · Department of Biomedical Engineering University of Arkansas, Fayetteville, Arkansas, USA
Background: Molecular testing for NPM1 and FLT3-ITD mutations guides critical early treatment decisions in acute myeloid leukemia (AML), but results can take weeks, long after these decisions must be made. Flow cytometry, already performed within hours of admission as part of routine care, may carry enough signal to predict these mutations directly, without added cost or delay. Methods: We developed an interpretable multi-instance learning classifier based on a decision tree, in which each patient sample is modeled as a collection of individual cells and mutation status is inferred from cell-level predictions. The model was benchmarked against a random forest trained on clinical variables and a deep convolutional neural network adapted for multitube flow cytometry data. Performance was assessed by cross-validation on a discovery cohort of 197 patients and tested on an independent cohort of 161 patients, using the area under the receiver operating characteristic curve (AUROC) and positive predictive value. Results: In cross-validation on the discovery cohort, the MIL model achieved mean AUROCs of 0.96 (SD=0.05) for NPM1 and 0.86 (SD=0.10) for FLT3-ITD, outperforming the clinical baseline and matching deep learning approaches. The model then successfully generalized to the independent test cohort of 161 patients, reaching AUROCs of 0.90 (NPM1) and 0.82 (FLT3-ITD), with positive predictive values of 0.87 and 0.68, respectively. Cell-level interpretation recovered established immunophenotypic signatures (CD33+ /CD34___ for NPM1-mutated cases, CD33+ /low side-scatter for FLT3-ITD), directly linking model predictions to known biology. Conclusions: These results show that an interpretable model applied to data already collected in routine care can predict AML molecular status within hours, offering a practical route to earlier, biology-informed treatment decisions.
Jonathan Legrand, Aguirre Mimoun, Baudouin Denis de Senneville +3