FreeBridge: Variational Schrödinger Bridges for Cellular Transition Dynamics

Modeling single-cell perturbation response as endpoint-only Schrödinger Bridge transport, constrained to the manifold of real cell morphologies.

MICCAI 2026
Xurui Wang1,2 ★    Qin Ren1 ★    Jun Ma3    Haibin Ling1    Chenyu You1 ✉
1 Stony Brook University       2 University of Toronto       3 University Health Network
Equal Contribution       Corresponding Author
A Schrödinger bridge transports control cells (left) to perturbed cells (right), constrained to the observed-morphology manifold.
15.35
FID ↓
overall · BBBC021
0.11
Rviol
vs 0.31 best prior
72.1%
MoA ↑
mechanism retention

News

Overview

Modeling how cells respond to perturbations is central to phenotypic drug discovery — yet high-content imaging fixes cells at acquisition, so a single cell's trajectory is never observed, only control and perturbed populations as separate marginals.

Matching endpoints alone leaves the interior undetermined: many stochastic processes share the same control and perturbed marginals while passing through morphologies that never occur in real cells.

FreeBridge is a Schrödinger Bridge for single-cell transition modeling under endpoint-only supervision — stochastic transport regularized to stay on the manifold of real single-cell morphologies.

Keywords: High-Content Imaging · Cellular Morphology Modeling · Schrödinger Bridges · Perturbation Dynamics

Method

FreeBridge casts endpoint-only supervision as a variational Schrödinger Bridge: among all stochastic processes matching the observed control and perturbed marginals, it seeks the one of minimal control energy whose intermediate states stay on the empirical morphology manifold, trading transport energy against a latent-support penalty $V(z_t)$:

$$\min_{f_{\theta}}\;\mathbb{E}\left[\int_{0}^{1}\left(\frac{1}{2\sigma^{2}}\,\|f_{\theta}(z_t, t, c)\|^{2}+\lambda\,V(z_{t})\right)dt\right]$$

The Cell Engine

The model is realized as a Cell Engine that separates state specification from stochastic transport. Microscopy fields are instance-segmented into single-cell crops, and a frozen encoder embeds each crop as $z=\phi(x)$; the collection of embeddings forms an empirical latent support set $\mathcal{B}$ that pins down the manifold of real morphologies. Transport is a time-conditioned drift driving the SDE $dz_t = f_{\theta}(z_t,t,c)\,dt + \sigma\,dW_t$ from control cells $z_0\sim\nu_{z,c}$ at $t=0$ to perturbed cells $z_1\sim\mu_z$ at $t=1$. The drift is trained variationally to minimize control energy while the support cost $V(z)=\min_{b\in\mathcal{B}}\|z-b\|^{2}$ penalizes excursions off the manifold; the weight $\lambda_{\text{bank}}$ trades transport energy against support feasibility. Intra-batch pairing within each plate supplies the coupling, and instance segmentation ensures the bridge operates on single cells rather than whole fields.

FreeBridge framework overview
Figure 1. Overview of FreeBridge. (a) Microscopy fields are segmented into single-cell states and embedded as $z=\phi(x)$; a drift $dz_t=f_\theta(z_t,t,c)dt+\sigma dW_t$ transports Source ($z_0\sim\nu_{z,c}$, control) to Target ($z_1\sim\mu_z$, perturbed) cells under support cost $V(z)=\min_{b\in\mathcal{B}}\|z-b\|^2$. (b) On BBBC021, FreeBridge reaches lower overall and condition-level FID and KID than PhenDiff, IMPA, and CellFlux.

Why a Schrödinger Bridge?

GANs map noise directly to data and train unstably; diffusion models discretize generation into many denoising steps; both match the perturbed endpoint but say nothing about the path taken to reach it. A Schrödinger Bridge instead selects the minimal-energy stochastic process between two given marginals — exactly the control-to-perturbed setting. Adding the empirical support cost ties that process to observed morphology, so intermediate states remain biologically plausible rather than drifting through regions of latent space that no real cell occupies.

Datasets

FreeBridge is evaluated on three high-content imaging benchmarks under a single unified single-cell protocol.

BenchmarkSourceAssayPerturbations
BBBC021Broad BBBCMCF-7 breast-cancer cells, 3-channel imaging38 compounds, 12 mechanism-of-action classes
RxRx1RecursionFluorescent cells across batches / cell types1,108 siRNA genetic knockdowns
JUMPJUMP-CP ConsortiumCell Painting, multi-channel morphologyLarge-scale chemical & genetic perturbations

Results

We evaluate FreeBridge on BBBC021, RxRx1, and JUMP under a unified protocol.

Quantitative comparison

FreeBridge improves overall and condition-level FID and KID over PhenDiff, IMPA, and CellFlux across all three benchmarks, and cuts the support-violation rate on BBBC021 by roughly 3×.

BBBC021RxRx1JUMP
Method FIDoFIDcKIDoKIDcRviol FIDoFIDcKIDoKIDc FIDoFIDcKIDoKIDc
PhenDiff48.25107.503.083.260.3867.12175.625.285.3448.57126.525.135.23
IMPA33.9175.762.742.790.3439.94163.912.892.8714.81101.031.101.00
CellFlux18.7256.801.621.590.3133.00162.852.402.419.0083.210.650.67
FreeBridge (Ours)15.3547.281.101.020.1128.94147.331.982.038.5179.050.590.53

5k generated samples per condition, averaged over three seeds. Rviol (Support Violation Rate): fraction of intermediate states outside the empirical morphology support — a support-feasibility measure, not verified temporal ordering.

Does FreeBridge generate realistic perturbed cells?

We compare endpoint morphology on BBBC021 against IMPA, PhenDiff, and CellFlux for two perturbations with distinct phenotypes.

Qualitative comparison of generated single cells
Figure 2. Endpoint morphology on BBBC021 (Demecolcine; Mevinolin / Lovastatin): real target vs. baselines and FreeBridge.
1.Compound-specific structure. FreeBridge preserves nuclear compaction and cytoplasmic shrinkage characteristic of each perturbation.
2.Sharper than diffusion baselines. PhenDiff and IMPA smooth texture or under-perturb, blurring phenotype-defining detail.
3.Best endpoint fidelity. Lowest FIDo/FIDc across all three datasets in the table above.
Do intermediate states stay on the cell manifold?

Endpoint alignment leaves the interior underdetermined. We visualize trajectories across $t = 0 \to 1$ and sweep the support weight $\lambda_{\text{bank}}$.

Intermediate trajectories under endpoint-constrained transport
Figure 3. Trajectories across $t=0.00\to1.00$. Top (unconstrained): endpoints align but intermediate states cross low-support regions. Bottom (FreeBridge): support regularization keeps samples on the manifold.

1. Endpoints align regardless. Both constrained and unconstrained drifts reproduce the control and perturbed marginals.

2. Unconstrained drifts off-manifold. Without the support cost, intermediate samples traverse low-support regions and produce artifacts.

3. Support cost pulls trajectories back. Raising $\lambda_{\text{bank}}$ lowers feasibility energy and violation rate together.

VariantEfeasRviol
$\lambda_{\text{bank}}=0$ (≡ CellFlux)2.480.31
$\lambda_{\text{bank}}=0.3$1.850.19
$\lambda_{\text{bank}}=0.5$1.320.11

Setting $\lambda_{\text{bank}}=0$ removes the support cost and recovers the CellFlux baseline — matching its Rviol of 0.31 in the main table.

What does each component contribute?

We ablate the three ingredients of the Cell Engine, reporting endpoint fidelity and mechanism-of-action retention.

VariantFIDoMoA%↑KIDo
w/o intra-batch pairing19.9463.271.81
w/o support cost $V(\cdot)$17.1865.081.55
w/o instance segmentation28.5451.342.45
Full model15.3772.081.13

1. Instance segmentation matters most. Removing it collapses MoA to 51% and nearly doubles FID — single-cell geometry is essential.

2. Support cost lifts mechanism retention. Dropping $V(\cdot)$ costs ~7 MoA points and raises feasibility violations.

3. Intra-batch pairing. Removing the coupling degrades FID from 15.37 to 19.94.

Full-model numbers are from an independent ablation run and differ marginally from the main table.

How sensitive is FreeBridge to its hyperparameters?

We sweep support weight $\lambda_{\text{bank}}$, base noise $\sigma$, and inference steps (NFE) on BBBC021.

Hyperparameter sensitivity
Figure 4. Sensitivity to $\lambda_{\text{bank}}$, $\sigma$, and NFE (20, 50, 100) across $\text{FID}_o$, $\text{FID}_c$, $\text{MoA}\%$, and $\text{KID}_o$.
1.Endpoint fidelity is stable. FID and KID vary little across the full range of support weight, noise, and step count.
2.$\lambda_{\text{bank}}$ trades feasibility for energy. Mechanism retention and support feasibility peak around the default 0.5.
3.Few inference steps suffice. NFE from 20 to 100 yields near-identical FID — cheap sampling at test time.

Quick Start

The single-cell pipeline runs control/perturbed endpoints through a frozen latent encoder, then trains the bridge with Hydra experiment profiles.

Environment
# PyTorch-Lightning 1.x required (pinned to 1.8.5.post0) — 2.x is not supported
conda env create -f environment.yml
conda activate freebridge
pip install -e .
Data preparation · BBBC021
# crop single cells → train VAE → export latents → build endpoints
python data_prep/crop_bbbc021_singlecell.py --mask_root <masks>
python data_prep/train_vae_full.py
python data_prep/export_latents_from_vae.py
python data_prep/build_bbbc021_endpoints.py \
    --control_regex DMSO --target_regex taxol \
    --out_dir data/bbbc021
Training
# default FreeBridge, support weight lambda_bank = 0.5
python train.py experiment=bbbc021

# support-weight sweep
python train.py experiment=bbbc021 state_cost.support_weight=0.0
python train.py experiment=bbbc021 state_cost.support_weight=1.0

Epoch 0 uses a straight-line warm start; manifold-constrained paths activate from epoch 1, so set max_epochs ≥ 2. Full setup and configs are in the GitHub repository.

Citation

@article{wang2026freebridge,
    title={FreeBridge: Variational Schr{\"o}dinger Bridges for Cellular Transition Dynamics},
    author={Wang, Xurui and Ren, Qin and Ma, Jun and Ling, Haibin and You, Chenyu},
    journal={arXiv preprint arXiv:2606.11286},
    year={2026},
}

Acknowledgements

Built on GSBM (Liu et al., ICLR 2024); the single-cell task and unified evaluation protocol follow CellFlux (Zhang et al., 2025). We thank the authors for releasing their code.