Search papers, labs, and topics across Lattice.
KoopCell addresses the failure of flow matching and neural differential equations to extrapolate population trajectories and capture cell-fate branching from temporally sparse, unpaired distribution snapshots. Leveraging Koopman-Mori-Zwanzig operator theory, the framework derives a closed-form least-squares estimator with convergence guarantees for the Koopman generator via the weak continuity equation, while incorporating non-Markovian memory through Markovian embeddings to explicitly capture branching. Across synthetic benchmarks and three scRNA-seq datasets, the approach accurately reconstructs latent Koopman spectra and outperforms continuous generative baselines on high-dimensional distribution forecasting.
Complex developmental branching and population extrapolation can be resolved with closed-form linear latent operators rather than brittle neural ODEs by framing cell memory through Koopman-Mori-Zwanzig theory.
Learning population dynamics from temporally sparse, unpaired distribution snapshots is a fundamental challenge in developmental biology. Recent approaches based on neural differential equations and flow matching can interpolate between observed population snapshots, but may struggle to extrapolate beyond the training horizon and often lack an explicit mechanism for modeling developmental branching. We propose KoopCell, a unified generative framework based on Koopman-Mori-Zwanzig theory that jointly learns representations and predictive linear latent dynamics. Theoretically, using the weak continuity equation, we derive a closed-form least-squares estimator for the Koopman generator from distribution snapshots and establish convergence guarantees under suitable assumptions. To model branching dynamics, we further develop KoopCell-M, which incorporates non-Markovian memory into the latent Koopman dynamics through a Markovian embedding. Experiments on synthetic systems and three scRNA-seq datasets demonstrate the ability of our framework to recover Koopman spectra, model branching through memory, and scale to predicting high-dimensional gene expression distributions, achieving state-of-the-art performance among the evaluated methods.