Abstract:Numerous models have been developed to account for the complex properties of the random walks of biomolecules. However, when analysing experimental data, conditions are rarely met to ensure model identification. The dynamics may simultaneously be influenced by spatial and temporal heterogeneities of the environment, out-of-equilibrium fluxes and conformal changes of the tracked molecules. Recorded trajectories are often too short to reliably discern such multi-scale dynamics, which precludes unambiguous assessment of the type of random walk and its parameters. Furthermore, the motion of biomolecules may not be well described by a single, canonical random walk model. Here, we develop a two-step statistical testing scheme for comparing biomolecule dynamics observed in different experimental conditions without having to identify or make strong prior assumptions about the model generating the recorded random walks. We first train a graph neural network to perform simulation-based inference and thus learn a rich summary statistics vector describing individual trajectories. We then compare trajectories obtained in different biological conditions using a non-parametric maximum mean discrepancy (MMD) statistical test on their so-obtained summary statistics. This procedure allows us to characterise sets of random walks regardless of their generating models, without resorting to model-specific physical quantities or estimators. We first validate the relevance of our approach on numerically simulated trajectories. This demonstrates both the statistical power of the MMD test and the descriptive power of the learnt summary statistics compared to estimates of physical quantities. We then illustrate the ability of our framework to detect changes in α -synuclein dynamics at synapses in cultured cortical neurons, in response to membrane depolarisation, and show that detected differences are largely driven by increased protein mobility in the depolarised state, in agreement with previous findings. The method provides a means of interpreting the differences it detects in terms of single trajectory characteristics. Finally, we emphasise the interest of performing various comparisons to probe the heterogeneity of experimentally acquired datasets at different levels of granularity (e.g., biological replicates, fields of view, and organelles). The continuous improvement of methods for single molecule tracking in live cells are driving our understanding of how biomolecules move inside cells. Analysing trajectories of single molecules is complicated by their highly erratic and noisy nature and thus requires the use of statistical models of their motion. However, it is often not possible to unambiguously determine a model from a set of short and noisy trajectories. Furthermore, the heterogeneous nature of the cellular environment means that the molecules' motion is often not properly described by a single model. In this paper we develop a new statistical testing scheme to detect differences in biomolecule dynamics within organelles without needing to identify a model of their motion. We train a graph neural network on large-scale simulations of random walks to learn a latent representation that captures relevant physical properties of a trajectory. We use a kernel-based statistical test within that latent space to compare the properties of two sets of trajectories recorded under different biological conditions. We apply our approach to detect differences in the dynamics of α -synuclein, a presynaptic protein, in axons and boutons during synaptic stimulation. This represents an important step towards automated single-molecule-based read-out of pharmacological action.

Simulation-based inference of single-molecule experiments

Simulation-based inference of single-molecule force spectroscopy

A Comprehensive Guide to Simulation-based Inference in Computational Biology

Amortized template-matching of molecular conformations from cryo-electron microscopy images using simulation-based inference

Application of Molecular Simulation in Biosensor Development

Simulation-based inference for efficient identification of generative models in computational connectomics

InferenceMAP: Mapping of Single-Molecule Dynamics with Bayesian Inference

Fast, accurate and lightweight sequential simulation-based inference using Gaussian locally linear mappings

sbi reloaded: a toolkit for simulation-based inference workflows

Simulation-based inference for non-parametric statistical comparison of biomolecule dynamics

Compositional simulation-based inference for time series

Methods and considerations for estimating parameters in biophysically detailed neural models with simulation based inference

Simultaneous identification of models and parameters of scientific simulators

Cost-aware Simulation-based Inference

Generalized Bayesian Inference for Scientific Simulators via Amortized Cost Estimation

Comparing AI versus Optimization Workflows for Simulation-Based Inference of Spatial-Stochastic Systems

Simulation-based Inference for Cardiovascular Models

Biomolecular Simulation: A Computational Microscope for Molecular Biology

Accelerating All-Atom Simulations and Gaining Mechanistic Understanding of Biophysical Systems through State Predictive Information Bottleneck

The confluence of machine learning and multiscale simulations

Optimal Design of Experiments for Simulation-Based Inference of Mechanistic Acyclic Biological Networks