Active-Matter models commonly consider particles with overdamped dynamics subject to a force (speed) with constant modulus and random direction. Some models also include random noise in particle displacement (a Wiener process), resulting in diffusive motion at short time scales. On the other hand, Ornstein–Uhlenbeck processes apply Langevin dynamics to the particles’ velocity and predict motion that is not diffusive at short time scales. Experiments show that migrating cells have gradually varying speeds at intermediate and long time scales, with short-time diffusive behavior. While Ornstein–Uhlenbeck processes can describe the moderate-and long-time speed variation, Active-Matter models for over-damped particles can explain the short-time diffusive behavior. Isotropic models cannot explain both regimes, because short-time diffusion renders instantaneous velocity ill-defined, and prevents the use of dynamical equations that require velocity time-derivatives. On the other hand, both models correctly describe some of the different temporal regimes seen in migrating biological cells and must, in the appropriate limit, yield the same observable predictions. Here we propose and solve analytically an Anisotropic Ornstein–Uhlenbeck process for polarized particles, with Langevin dynamics governing the particle’s movement in the polarization direction and a Wiener process governing displacement in the orthogonal direction. Our characterization provides a theoretically robust way to compare movement in dimensionless simulations to movement in experiments in which measurements have meaningful space and time units. We also propose an approach to deal with inevitable finite-precision effects in experiments and simulations.
more »
« less
Aging power spectrum of membrane protein transport and other subordinated random walks
Abstract Single-particle tracking offers detailed information about the motion of molecules in complex environments such as those encountered in live cells, but the interpretation of experimental data is challenging. One of the most powerful tools in the characterization of random processes is the power spectral density. However, because anomalous diffusion processes in complex systems are usually not stationary, the traditional Wiener-Khinchin theorem for the analysis of power spectral densities is invalid. Here, we employ a recently developed tool named aging Wiener-Khinchin theorem to derive the power spectral density of fractional Brownian motion coexisting with a scale-free continuous time random walk, the two most typical anomalous diffusion processes. Using this analysis, we characterize the motion of voltage-gated sodium channels on the surface of hippocampal neurons. Our results show aging where the power spectral density can either increase or decrease with observation time depending on the specific parameters of both underlying processes.
more »
« less
- Award ID(s):
- 2102832
- PAR ID:
- 10307941
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The stochastic trajectories of molecules in living cells, as well as the dynamics in many other complex systems, often exhibit memory in their path over long periods of time. In addition, these systems can show dynamic heterogeneities due to which the motion changes along the trajectories. Such effects manifest themselves as spatiotemporal correlations. Despite the broad occurrence of heterogeneous complex systems in nature, their analysis is still quite poorly understood and tools to model them are largely missing. We contribute to tackling this problem by employing an integral representation of Mandelbrot’s fractional Brownian motion that is compliant with varying motion parameters while maintaining long memory. Two types of switching fractional Brownian motion are analysed, with transitions arising from a Markovian stochastic process and scale-free intermittent processes. We obtain simple formulas for classical statistics of the processes, namely the mean squared displacement and the power spectral density. Further, a method to identify switching fractional Brownian motion based on the distribution of displacements is described. A validation of the model is given for experimental measurements of the motion of quantum dots in the cytoplasm of live mammalian cells that were obtained by single-particle tracking.more » « less
-
{"Abstract":["Datasets generated in the report "Aging power spectrum of membrane protein transport and other subordinated random walks". Included data are:<\/p>\n\nNumerical simulations <\/strong>\nRWdata1.mat: 10,000 realizations, subordinated random walk with Hurst exponent, H<\/em>=0.3 and \\(\\alpha\\)=0.4.\nRWdata3.mat: 10,000 realizations, subordinated random walk with Hurst exponent, H<\/em>=0.7 and \\(\\alpha\\)=0.4.\nRWdata8.mat: 5,000 realizations, subordinated random walk with Hurst exponent, H<\/em>=0.75 and \\(\\alpha\\)=0.8.\nRWdataCTRW.mat: 10,000 realizations, continuous time random walk (CTRW), \\(\\alpha\\)=0.7.<\/p>\n\nSpectra of simulations<\/strong>\nPSDdata1.mat: Power spectral density (PSD) of a subordinated random walk with Hurst exponent, H<\/em>=0.3 and \\(\\alpha\\)=0.4. Five different realization times are used to compute the PDS: 2^8, 2^10, 2^12, 2^14, and 2^16.\nPSDdata3.mat: PSD of a subordinated random walk with Hurst exponent, H<\/em>=0.7 and \\(\\alpha\\)=0.4. Five different realization times are used to compute the PDS: 2^8, 2^10, 2^12, 2^14, and 2^16.\nPSDdata8.mat: PSD of a subordinated random walk with Hurst exponent, H<\/em>=0.75 and \\(\\alpha\\)=0.8. Four different realization times are used to compute the PDS: 2^15, 2^16, 2^17, and 2^18.\nPSDs_CTRW.mat: PSD of a continuous-time random walk (CTRW), \\(\\alpha\\)=0.7. Five different realization times are used to compute the PDS: 2^8, 2^10, 2^12, 2^14, and 2^16.<\/p>\n\nExperimental data of Nav1.6 channels in the soma of hippocampal neurons<\/strong>\nNavMSDtimes.csv: ensemble-averaged (EA) MSD and time-averaged (TA) MSD. The TA-MSD is measured for three observation times, 64, 128, and 256 frames (3.2, 6.4, and 12.8 s).\nNavPSD.csv: Power spectral density (PSD) measured for three observation times, 64, 128, and 256 frames.<\/p>"],"Other":["We acknowledge the support of the National Science Foundation grant 2102832 (to DK) and Israel Science Foundation grant 1898/17 (to EB).","{"references": ["Fox, Z.R., Barkai, E. & Krapf, D. Aging power spectrum of membrane protein transport and other subordinated random walks. Nat Commun 12, 6162 (2021). https://doi.org/10.1038/s41467-021-26465-8"]}"]}more » « less
-
The motion of nanoparticles near surfaces is of fundamental importance in physics, biology, and chemistry. Liquid cell transmission electron microscopy (LCTEM) is a promising technique for studying motion of nanoparticles with high spatial resolution. Yet, the lack of understanding of how the electron beam of the microscope affects the particle motion has held back advancement in using LCTEM for in situ single nanoparticle and macromolecule tracking at interfaces. Here, we experimentally studied the motion of a model system of gold nanoparticles dispersed in water and moving adjacent to the silicon nitride membrane of a commercial LC in a broad range of electron beam dose rates. We find that the nanoparticles exhibit anomalous diffusive behavior modulated by the electron beam dose rate. We characterized the anomalous diffusion of nanoparticles in LCTEM using a convolutional deep neural-network model and canonical statistical tests. The results demonstrate that the nanoparticle motion is governed by fractional Brownian motion at low dose rates, resembling diffusion in a viscoelastic medium, and continuous-time random walk at high dose rates, resembling diffusion on an energy landscape with pinning sites. Both behaviors can be explained by the presence of silanol molecular species on the surface of the silicon nitride membrane and the ionic species in solution formed by radiolysis of water in presence of the electron beam.more » « less
-
Mucus is a complex fluid that coats multiple organs in animals. Various physicochemical properties can alter the diffusion of microscopic particles in mucus, impacting drug delivery, virus infection, and disease development. The simultaneous effect of these physicochemical properties in particle diffusion, however, remains elusive. Here, we analyzed 106 published experiments to identify the most dominant factors controlling particle diffusion in mucus. The effective diffusion—defined using a one-second sampling time window across experiments—spanned seven orders of magnitude, from 10 –5 to 10 2 μm 2 /s. Univariate and multivariate statistical analyses identified the anomalous exponent (the logarithmic slope of the mean-squared displacement) as the strongest predictor of effective diffusion, revealing an exponential relationship that explained 89 % of the variance. A theoretical scaling analysis revealed that a stronger correlation of the anomalous exponent over the generalized diffusion constant occurs for sampling times two orders of magnitude larger than the characteristic molecular (or local) displacement time. This result predicts that at these timescales, the molecular properties controlling the anomalous exponent, like particle–mucus unbinding times or the particle to mesh size ratio, would be the most relevant physicochemical factors involved in passive microrheology of particles in mucus. Our findings contrast with the fact that only one-third of the studies measured the anomalous exponent, and most experiments did not report the associated molecular properties predicted to dominate the motion of particles in mucus. The theoretical foundation of our work can be extrapolated to other systems, providing a guide to identify dominant molecular mechanisms regulating the mobility of particles in mucus and other polymeric fluids.more » « less
An official website of the United States government
