A data-driven alternative to the fractional Fokker–Planck equation

Abstract

Anomalous diffusion processes are ubiquitous in biology and arise in the transport of proteins, vesicles and other particles. Such anomalously diffusive behavior is attributed to a number of factors within the cell including heterogeneous environments, active transport processes and local trapping/binding. There are a number of microscopic principles?such as power law jump size and/or waiting time distributions?from which the fractional Fokker?Planck equation (FFPE) can be derived and used to provide mechanistic insight into the origins of anomalous diffusion. On the other hand, it is fair to ask if other microscopic principles could also have given rise to the evolution of an observed density profile that appears to be well fit by an FFPE. Here we discuss another possible mechanistic alternative that can give rise to densities like those generated by FFPEs. Rather than to fit a density (or concentration profile) using a solution to the spatial FFPE, we reconstruct the profile generated by an FFPE using a regular FPE with a spatial and time-dependent force. We focus on the special case of the spatial FFPE for superdiffusive processes. This special case is relevant to, for example, active transport in a biological context. We devise a prescription for extracting such forces on synthetically generated data and provide an interpretation to the forces extracted. In particular, the time-dependence of forces could tell us about ATP depletion or changes in the cell's metabolic activity. Modeling anomalous behavior with normal diffusion driven by these effective forces yields an alternative mechanistic picture that, ultimately, could help motivate future experiments.