Elucidating the origin of anomalous diffusion in crowded fluids

Intracellular Macromolecular Crowding within Individual Stress Fibers Analyzed by Fluorescence Correlation Spectroscopy

Introduction

The diffusion of cell components such as proteins is crucial to the function of all living cells. The abundance of macromolecules in cells is likely to cause a state of macromolecular crowding, but its effects on the extent of diffusion remain poorly understood.

Methods

Here we investigate the diffusion rate in three distinct locations in mesenchymal cell types, namely the open cytoplasm, the stress fibers in the open cytoplasm, and those below the nucleus using three kinds of biologically inert green fluorescent proteins (GFPs), namely a monomer, dimer, and trimer GFP. Fluorescence correlation spectroscopy (FCS) was used to determine the diffusion coefficients.

Results

We show that diffusion tends to be lowered on average in stress fibers and is significantly lower in those located below the nucleus. Our data suggest that the diffusive properties of GFPs, and potentially other molecules as well, are hindered by macromolecular crowding. However, although the size dependence on protein diffusion was also studied for monomer, dimer, and trimer GFPs, there was no significant difference in the diffusion rates among the GFPs of these sizes. These results could be attributed to the lack of significant change in protein size among the selected GFP multimers.

Conclusion

The data presented here would provide a basis for better understanding of the complex protein diffusion in the nonuniform cytoplasm, shedding light on cellular responses to mechanical stress, their local mechanical properties, and reduced turnover in senescent cells.

Aging and confinement in subordinated fractional Brownian motion

We study the effects of aging properties of subordinated fractional Brownian motion (FBM) with drift and in harmonic confinement, when the measurement of the stochastic process starts a time t_{a}>0 after its original initiation at t=0. Specifically, we consider the aged versions of the ensemble mean-squared displacement (MSD) and the time-averaged MSD (TAMSD), along with the aging factor. Our results are favorably compared with simulations results. The aging subordinated FBM exhibits a disparity between MSD and TAMSD and is thus weakly nonergodic, while strong aging is shown to effect a convergence of the MSD and TAMSD. The information on the aging factor with respect to the lag time exhibits an identical form to the aging behavior of subdiffusive continuous-time random walks (CTRW). The statistical properties of the MSD and TAMSD for the confined subordinated FBM are also derived. At long times, the MSD in the harmonic potential has a stationary value, that depends on the Hurst index of the parental (nonequilibrium) FBM. The TAMSD of confined subordinated FBM does not relax to a stationary value but increases sublinearly with lag time, analogously to confined CTRW. Specifically, short aging times t_{a} in confined subordinated FBM do not affect the aged MSD, while for long aging times the aged MSD has a power-law increase and is identical to the aged TAMSD.

Diffusion in Porous Rock Is Anomalous

Molecular diffusion of chemical species in subsurface environments─rock formations, soil sediments, marine, river, and lake sediments─plays a critical role in a variety of dynamic processes, many of which affect water chemistry. We investigate and demonstrate the occurrence of anomalous (non-Fickian) diffusion behavior, distinct from classically assumed Fickian diffusion. We measured molecular diffusion through a series of five chalk and dolomite rock samples over a period of about two months. We demonstrate that in all cases, diffusion behavior is significantly different than Fickian. We then analyze the results using a continuous time random walk framework that can describe anomalous diffusion in heterogeneous porous materials such as rock. This methodology shows extreme long-time tailing of tracer advance as compared to conventional Fickian diffusion processes. The finding that distinct anomalous diffusion occurs ubiquitously implies that diffusion-driven processes in subsurface zones should be analyzed using tools that account for non-Fickian diffusion.

Concentration-Dependent Diffusion of Monoclonal Antibodies: Underlying Mechanisms of Anomalous Diffusion

The development, storage, transport, and subcutaneous delivery of highly concentrated monoclonal antibody formulations pose significant challenges due to the high solution viscosity and low diffusion of the antibody molecules in crowded environments. These issues often stem from the self-associating behavior of the antibody molecules, potentially leading to aggregation. In this work, we used a dissipative particle dynamics-based coarse-grained model to investigate the diffusion behavior of IgG1 antibody molecules in aqueous solutions with 15 and 32 mM NaCl and antibody concentrations ranging from 10 to 400 mg/mL. We determined the coarse-grained interaction parameters by matching the calculated structure factor with the computational and experimental data from the literature. Our results indicate Fickian diffusion for antibody concentrations of 10 and 25 mg/mL and anomalous diffusion for concentrations exceeding 50 mg/mL. The anomalous diffusion was observed for ∼0.33 to 0.4 μs, followed by Fickian diffusion for all antibody concentrations. We observed a strong linear correlation between the diffusion behavior of the antibody molecules (diffusion coefficient D and anomalous diffusion exponent α) and the amount of aggregates present in the solution and between the amount of aggregates and the Coulomb interaction energy. The investigation of underlying mechanisms for anomalous diffusion revealed that in crowded environments at high antibody concentrations, the attractive interaction between electrostatically complementary regions of the antibody molecules could further bring the neighboring molecules closer to one another, ultimately resulting in aggregate formation. Further, the Coulomb attraction can continue to draw more molecules together, forming larger aggregates.

Ergodic criterion of a random diffusivity model

The random diffusivity, initially proposed to explain Brownian yet non-Gaussian diffusion, has garnered significant attention due to its capacity not only for elucidating the internal physical mechanism of non-Gaussian diffusion, but also for establishing an analytical framework to characterize particle motion in complex environments. In this paper, based on the correlation function C(t_{1},t_{2})=〈D(t_{1})D(t_{2})〉 of random diffusivity D(t), we quantitatively propose a general criterion of determining the ergodic property of the Langevin equation with the arbitrary random diffusivity D(t). Due to the critical role of correlation function C(t_{1},t_{2}), we derive the criterion for the two cases with stationary diffusivity or nonstationary diffusivity, respectively. By utilizing the quantitative criterion, we can directly judge the ergodic properties of the random diffusivity model based on the correlation function C(t_{1},t_{2}) of random diffusivity D(t). Several typical diffusivities, including the common square of the Brownian motion and of the (fractional) Ornstein-Uhlenbeck process, are found to contribute to different ergodic properties, which validates our proposed criterion built on the correlation function C(t_{1},t_{2}).

Molecular Crowding: The History and Development of a Scientific Paradigm

It is now generally accepted that macromolecules do not act in isolation but "live" in a crowded environment, that is, an environment populated by numerous different molecules. The field of molecular crowding has its origins in the far 80s but became accepted only by the end of the 90s. In the present issue, we discuss various aspects that are influenced by crowding and need to consider its effects. This Review is meant as an introduction to the theme and an analysis of the evolution of the crowding concept through time from colloidal and polymer physics to a more biological perspective. We introduce themes that will be more thoroughly treated in other Reviews of the present issue. In our intentions, each Review may stand by itself, but the complete collection has the aspiration to provide different but complementary perspectives to propose a more holistic view of molecular crowding.

Single-molecule imaging of stochastic interactions that drive dynein activation and cargo movement in cells

Cytoplasmic dynein 1 (dynein) is the primary minus end-directed motor protein in most eukaryotic cells. Dynein remains in an inactive conformation until the formation of a tripartite complex comprising dynein, its regulator dynactin, and a cargo adaptor. How this process of dynein activation occurs is unclear since it entails the formation of a three-protein complex inside the crowded environs of a cell. Here, we employed live-cell, single-molecule imaging to visualize and track fluorescently tagged dynein. First, we observed that only ∼30% of dynein molecules that bound to the microtubule (MT) engaged in minus end-directed movement, and that too for a short duration of ∼0.6 s. Next, using high-resolution imaging in live and fixed cells and using correlative light and electron microscopy, we discovered that dynactin and endosomal cargo remained in proximity to each other and to MTs. We then employed two-color imaging to visualize cargo movement effected by single motor binding. Finally, we performed long-term imaging to show that short movements are sufficient to drive cargo to the perinuclear region of the cell. Taken together, we discovered a search mechanism that is facilitated by dynein's frequent MT binding-unbinding kinetics: (i) in a futile event when dynein does not encounter cargo anchored in proximity to the MT, dynein dissociates and diffuses into the cytoplasm, (ii) when dynein encounters cargo and dynactin upon MT binding, it moves cargo in a short run. Several of these short runs are undertaken in succession for long-range directed movement. In conclusion, we demonstrate that dynein activation and cargo capture are coupled in a step that relies on the reduction of dimensionality to enable minus end-directed transport in cellulo and that complex cargo behavior emerges from stochastic motor-cargo interactions.

Anomalous Transport of Small Polarons Arises from Transient Lattice Relaxation or Immovable Boundaries

Elucidating transport mechanisms is crucial for advancing material design, yet state-of-the-art theory is restricted to exact simulations of small lattices with severe finite-size effects or approximate ones that assume the nature of transport. We leverage algorithmic advances to tame finite-size effects and exactly simulate small polaron formation and transport in the Holstein model. We further analyze the applicability of the ubiquitously used equilibrium-based Green-Kubo relations and nonequilibrium methods to predict charge mobility. We find that these methods can converge to different values and track this disparity to finite-size dependence and the sensitivity of Green-Kubo relations to the system's topology. Contrary to standard perturbative calculations, our results demonstrate that small polarons exhibit anomalous transport that manifests transiently due to nonequilibrium lattice relaxation or permanently as a signature of immovable boundaries. These findings can offer new interpretations of transport experiments on polymers and transition metal oxides.

Information Transmission via Molecular Communication in Astrobiological Environments

The ubiquity of information transmission via molecular communication between cells is comprehensively documented on Earth; this phenomenon might even have played a vital role in the origin(s) and early evolution of life. Motivated by these considerations, a simple model for molecular communication entailing the diffusion of signaling molecules from transmitter to receiver is elucidated. The channel capacity C (maximal rate of information transmission) and an optimistic heuristic estimate of the actual information transmission rate ℐ are derived for this communication system; the two quantities, especially the latter, are demonstrated to be broadly consistent with laboratory experiments and more sophisticated theoretical models. The channel capacity exhibits a potentially weak dependence on environmental parameters, whereas the actual information transmission rate may scale with the intercellular distance d as ℐ ∝ d-4 and could vary substantially across settings. These two variables are roughly calculated for diverse astrobiological environments, ranging from Earth's upper oceans (C ∼ 3.1 × 103 bits/s; ℐ ∼ 4.7 × 10-2 bits/s) and deep sea hydrothermal vents (C ∼ 4.2 × 103 bits/s; ℐ ∼ 1.2 × 10-1 bits/s) to the hydrocarbon lakes and seas of Titan (C ∼ 3.8 × 103 bits/s; ℐ ∼ 2.6 × 10-1 bits/s).

Multiple diffusive behaviors of the random walk in inhomogeneous environments

Anomalous diffusive behaviors are observed in highly inhomogeneous but relatively stable environments such as intracellular media and are increasingly attracting attention. In this paper we develop a coupled continuous-time random walk model in which the waiting time is power-law coupled with the local environmental diffusion coefficient. We provide two forms of the waiting time density, namely, a heavy-tailed density and an exponential density. For different waiting time densities, anomalous diffusions with the diffusion exponent between 0 and 2 and Brownian yet non-Gaussian diffusion can be realized within the present model. The diffusive behaviors are analyzed and discussed by deriving the mean-squared displacement and probability density function. In addition we derive the effective jump length density corresponding to the decoupled form to help distinguish the diffusion types. Our model unifies two kinds of anomalous diffusive behavior with different characteristics in the same inhomogeneous environment into a theoretical framework. The model interprets the random motion of particles in a complex inhomogeneous environment and reproduces the experimental results of different biological and physical systems.

Ligand-Receptor Interaction Triggers Hopping and Sliding Motions on Living Cell Membranes

Exploring the surface-capturing and releasing processes of nanocargo on the living cell membrane is critical for understanding the membrane translocation process. In this work, we achieve total internal reflection scattering (TIRS) illumination on a commercial dark-field optical microscope without the introduction of any additional optical components. By gradually reducing the diaphragm size in the excitation light path, the angle of the incident beam can be well manipulated. Under optimal conditions, the excitation light can be totally reflected at the glass/water interface, resulting in a thin layer of evanescent field for TIRS illumination. Due to the exponential decay feature of the evanescent field, the displacement of the nanocargo along the vertical direction can be directly resolved in the intensity track. With this method, we selectively monitor the dynamics of the transferrin-modified nanocargo on the living cell membrane. Transition between confined diffusion and long-range searching is involved in the binding site recognition process, which exhibits non-Gaussian and nonergodic-like behavior. More interestingly, 2D fast sliding and 3D hopping motions are also distinguished on the fluidic cell membrane, which is essentially modulated by the strength of ligand-receptor interactions, as revealed by the free-energy profiles. These heterogeneous and dynamic interactions together control the diffusion mode of the nanocargo on the lipid membrane and, thus, determine the cellular translocation efficiency.

Fractional and scaled Brownian motion on the sphere: The effects of long-time correlations on navigation strategies

We analyze fractional Brownian motion and scaled Brownian motion on the two-dimensional sphere S^{2}. We find that the intrinsic long-time correlations that characterize fractional Brownian motion collude with the specific dynamics (navigation strategies) carried out on the surface giving rise to rich transport properties. We focus our study on two classes of navigation strategies: one induced by a specific set of coordinates chosen for S^{2} (we have chosen the spherical ones in the present analysis), for which we find that contrary to what occurs in the absence of such long-time correlations, nonequilibrium stationary distributions are attained. These results resemble those reported in confined flat spaces in one and two dimensions [Guggenberger et al. New J. Phys. 21, 022002 (2019)1367-263010.1088/1367-2630/ab075f; Vojta et al. Phys. Rev. E 102, 032108 (2020)2470-004510.1103/PhysRevE.102.032108]; however, in the case analyzed here, there are no boundaries that affect the motion on the sphere. In contrast, when the navigation strategy chosen corresponds to a frame of reference moving with the particle (a Frenet-Serret reference system), then the equilibrium distribution on the sphere is recovered in the long-time limit. For both navigation strategies, the relaxation times toward the stationary distribution depend on the particular value of the Hurst parameter. We also show that on S^{2}, scaled Brownian motion, distinguished by a time-dependent diffusion coefficient with a power-scaling, is independent of the navigation strategy finding a good agreement between the analytical calculations obtained from the solution of a time-dependent diffusion equation on S^{2}, and the numerical results obtained from our numerical method to generate ensemble of trajectories.

Fluctuating landscapes and heavy tails in animal behavior

Animal behavior is shaped by a myriad of mechanisms acting on a wide range of scales. This immense variability hampers quantitative reasoning and renders the identification of universal principles elusive. Through data analysis and theory, we here show that slow non-ergodic drives generally give rise to heavy-tailed statistics in behaving animals. We leverage high-resolution recordings of C. elegans locomotion to extract a self-consistent reduced order model for an inferred reaction coordinate, bridging from sub-second chaotic dynamics to long-lived stochastic transitions among metastable states. The slow mode dynamics exhibits heavy-tailed first passage time distributions and correlation functions, and we show that such heavy tails can be explained by dynamics on a time-dependent potential landscape. Inspired by these results, we introduce a generic model in which we separate faster mixing modes that evolve on a quasi-stationary potential, from slower non-ergodic modes that drive the potential landscape, and reflect slowly varying internal states. We show that, even for simple potential landscapes, heavy tails emerge when barrier heights fluctuate slowly and strongly enough. In particular, the distribution of first passage times and the correlation function can asymptote to a power law, with related exponents that depend on the strength and nature of the fluctuations. We support our theoretical findings through direct numerical simulations.

Anomalous diffusion, non-Gaussianity, and nonergodicity for subordinated fractional Brownian motion with a drift

The stochastic motion of a particle with long-range correlated increments (the moving phase) which is intermittently interrupted by immobilizations (the trapping phase) in a disordered medium is considered in the presence of an external drift. In particular, we consider trapping events whose times follow a scale-free distribution with diverging mean trapping time. We construct this process in terms of fractional Brownian motion with constant forcing in which the trapping effect is introduced by the subordination technique, connecting "operational time" with observable "real time." We derive the statistical properties of this process such as non-Gaussianity and nonergodicity, for both ensemble and single-trajectory (time) averages. We demonstrate nice agreement with extensive simulations for the probability density function, skewness, kurtosis, as well as ensemble and time-averaged mean-squared displacements. We place a specific emphasis on the comparisons between the cases with and without drift.

Edwards-Wilkinson depinning transition in fractional Brownian motion background

There are various reports about the critical exponents associated with the depinning transition. In this study, we investigate how the disorder strength present in the support can account for this diversity. Specifically, we examine the depinning transition in the quenched Edwards-Wilkinson (QEW) model on a correlated square lattice, where the correlations are modeled using fractional Brownian motion (FBM) with a Hurst exponent of H.We identify a crossover time [Formula: see text] that separates the dynamics into two distinct regimes: for [Formula: see text], we observe the typical behavior of pinned surfaces, while for [Formula: see text], the behavior differs. We introduce a novel three-variable scaling function that governs the depinning transition for all considered H values. The associated critical exponents exhibit a continuous variation with H, displaying distinct behaviors for anti-correlated ([Formula: see text]) and correlated ([Formula: see text]) cases. The critical driving force decreases with increasing H, as the host medium becomes smoother for higher H values, facilitating fluid mobility. This fact causes the asymptotic velocity exponent [Formula: see text] to increase monotonically with H.

Normal and anomalous diffusion in a bouncing ball over an irregular surface

The problem of a bouncing ball on a nonplanar surface is investigated. We discovered that surface undulation adds a horizontal component to the impact force, which acquires a random character. Some aspects of Brownian motion are found in the horizontal distribution of the particle. On the x axis, normal and superdiffusion are observed. For the probability density's functional form, a scaling hypothesis is presented.

Extracting, quantifying, and comparing dynamical and biomechanical properties of living matter through single particle tracking

A panoply of new tools for tracking single particles and molecules has led to an explosion of experimental data, leading to novel insights into physical properties of living matter governing cellular development and function, health and disease. In this Perspective, we present tools to investigate the dynamics and mechanics of living systems from the molecular to cellular scale via single-particle techniques. In particular, we focus on methods to measure, interpret, and analyse complex data sets that are associated with forces, materials properties, transport, and emergent organisation phenomena within biological and soft-matter systems. Current approaches, challenges, and existing solutions in the associated fields are outlined in order to support the growing community of researchers at the interface of physics and the life sciences. Each section focuses not only on the general physical principles and the potential for understanding living matter, but also on details of practical data extraction and analysis, discussing limitations, interpretation, and comparison across different experimental realisations and theoretical frameworks. Particularly relevant results are introduced as examples. While this Perspective describes living matter from a physical perspective, highlighting experimental and theoretical physics techniques relevant for such systems, it is also meant to serve as a solid starting point for researchers in the life sciences interested in the implementation of biophysical methods.

Molecular Dynamics Simulations of Rhodamine B Zwitterion Diffusion in Polyelectrolyte Solutions

Polyelectrolytes continue to find wide interest and application in science and engineering, including areas such as water purification, drug delivery, and multilayer thin films. We have been interested in the dynamics of small molecules in a variety of polyelectrolyte (PE) environments; in this paper, we report simulations and analysis of the small dye molecule rhodamine B (RB) in several very simple polyelectrolyte solutions. Translational diffusion of the RB zwitterion has been measured in fully atomistic, 2 μs long molecular dynamics simulations in four different polyelectrolyte solutions. Two solutions contain the common polyanion sodium poly(styrene sulfonate) (PSS), one with a 30-mer chain and the other with 10 trimers. The other two solutions contain the common polycation poly(allyldimethylammonium) chloride (PDDA), one with two 15-mers and the other with 10 trimers. RB diffusion was also simulated in several polymer-free solutions to verify its known experimental value for the translational diffusion coefficient, D_RB, of 4.7 × 10^-6 cm²/s at 300 K. RB diffusion was slowed in all four simulated PE solutions, but to varying degrees. D_RB values of 3.07 × 10^-6 and 3.22 × 10^-6 cm²/s were found in PSS 30-mer and PSS trimer solutions, respectively, whereas PDDA 15-mer and trimer solutions yielded values of 2.19 × 10^-6 and 3.34 × 10^-6 cm²/s. Significant associations between RB and the PEs were analyzed and interpreted via a two-state diffusion model (bound and free diffusion) that describes the data well. Crowder size effects and anomalous diffusion were also analyzed. Finally, RB translation along the polyelectrolytes during association was characterized.

Occupation time statistics of the fractional Brownian motion in a finite domain

We study statistics of occupation times for a fractional Brownian motion (fBm), which is a typical model of a non-Markov process. Due to the non-Markovian nature, recurrence times to the origin depend on the history. Numerical simulations indicate that dependence on the sum of successive recurrence times becomes weak. As a result, the distribution of the occupation time in a finite domain follows the Mittag-Leffler distribution when the Hurst exponent of the fBm is close to 1/2. We show this distributional behavior of a time-averaged observable by renewal theory. This result is an extension of the distributional limit theorem known as the Darling-Kac theorem in general Markov processes to non-Markov processes.

The Fokker-Planck equation of the superstatistical fractional Brownian motion with application to passive tracers inside cytoplasm

By collecting from literature data experimental evidence of anomalous diffusion of passive tracers inside cytoplasm, and in particular of subdiffusion of mRNA molecules inside live Escherichia coli cells, we obtain the probability density function of molecules' displacement and we derive the corresponding Fokker-Planck equation. Molecules' distribution emerges to be related to the Krätzel function and its Fokker-Planck equation to be a fractional diffusion equation in the Erdélyi-Kober sense. The irreducibility of the derived Fokker-Planck equation to those of other literature models is also discussed.

Crowding and confinement act in concert to slow DNA diffusion within cell-sized droplets

Dynamics of biological macromolecules, such as DNA, in crowded and confined environments are critical to understanding cellular processes such as transcription, infection, and replication. However, the combined effects of cellular confinement and crowding on macromolecular dynamics remain poorly understood. Here, we use differential dynamic microscopy to investigate the diffusion of large DNA molecules confined in cell-sized droplets and crowded by dextran polymers. We show that confined and crowded DNA molecules exhibit universal anomalous subdiffusion with scaling that is insensitive to the degree of confinement and crowding. However, effective DNA diffusion coefficients Deff decrease up to 2 orders of magnitude as droplet size decreases—an effect that is enhanced by increased crowding. We mathematically model the coupling of crowding and confinement by combining polymer scaling theories with confinement-induced depletion effects. The generality and tunability of our system and models render them applicable to elucidating wide-ranging crowded and confined systems.

•

DNA diffusion measured in cell-sized droplets with differential dynamic microscopy

•

Combination of crowding and confinement leads to subdiffusion and slowing

•

Diffusion coefficients of DNA decrease strongly with decreasing droplet size

•

Polymer scaling theories and depletion effects predict observed dynamics

How macromolecules softness affects diffusion under crowding

Diffusion in a macromolecularly crowded environment is essential for many intracellular processes, from metabolism and catalysis to gene transcription and translation. So far, theoretical and experimental work has focused on anomalous subdiffusion, and the effects of interactions, shapes, and composition, while the compactness or softness of macromolecules has received less attention. Herein, we use Brownian dynamics simulations to study how the softness of crowders affects macromolecular diffusion. We find that in most cases, soft crowders slow down the diffusion less effectively than hard crowders like Ficoll. For instance, at a 30% occupied volume fraction, the diffusion in Ficoll70 is about 20% slower than in soft crowders of the same size. However, our simulations indicate that elongated macromolecules, such as double-stranded DNA pieces, can diffuse comparably or even faster in hard crowders. We relate these effects to the volume excluded by soft and hard crowders to different tracers. Our results show that the softness and shape of macromolecules are crucial factors determining diffusion under crowding, relevant to diverse intracellular environments.

Lévy Walk Dynamics in an External Constant Force Field in Non-Static Media

Based on the recognition of the huge change of the transport properties for diffusion particles in non-static media, we consider a Lévy walk model subjected to an external constant force in non-static media. Since the physical and comoving coordinates of non-static media are related by scale factor, we equivalently transfer the process from physical coordinate into comoving coordinate and derive the master equation governing the probability density function of the position of the particles in comoving coordinate. Utilizing the Hermite orthogonal polynomial expansions, some statistical properties are obtained, including the asymptotic behaviors of the first two moments in both coordinates and kurtosis. For some representative types of non-static media and Lévy walks, the striking and interesting phenomena originating from the interplay between non-static media, external force, and intrinsic stochastic motion are observed. The stationary distribution are also analyzed for some cases through numerical simulations.

Plasmid DNA Undergoes Two Compaction Regimes under Macromolecular Crowding

The laser light scattering experiments were performed to explore the role of dextran (size (d): 2.6, 6.9, and 17.0 nm) in compacting the plasmids (pBS: 2.9 kbps; pCMV-Tag2B: 4.3 kbps; and pET28a: 5.3 kbps) in vitro in the volume fraction (ϕ) range 0.01 to 0.15 of the macromolecular crowder. Two compaction regimes were observed in terms of the radius of gyration (R_g) for plasmid-dextran combinations, wherein the plasmid diffusivity is governed by normal diffusion and subdiffusion, respectively. Generalized scaling, R_g ∼ ϕ^-1/(1+x), where x represents the conformational geometry of plasmids, is reported. The plasmid conformation depends on the crowder's size, with larger conformational changes observed in the presence of smaller crowders. The second virial coefficient (A₂) and translational diffusion coefficient (D_t) indicate that entropically driven depletion of crowders, excluded volume, and interplasmid repulsive interactions govern plasmids' conformational changes, validated herein from the scaling of D_t with molecular weight.

Anomalous diffusion, aging, and nonergodicity of scaled Brownian motion with fractional Gaussian noise: overview of related experimental observations and models

How does a systematic time-dependence of the diffusion coefficient $D (t)$ affect the ergodic and statistical characteristics of fractional Brownian motion (FBM)? Here, we examine how the behavior of the...

Strong anomalous diffusive behaviors of the two-state random walk process

The phenomenon of the two-state process is observed in various systems and is increasingly attracting attention, such that there is a need for a theoretical model of the process. In this paper, we present a prototypal two-state random walk (TSRW) model of a renewal process alternating between the continuous-time random walk (CTRW) state and Lévy walk (LW) state. The jump length distribution of the CTRW state is assumed to be Gaussian whereas the time distributions of the two states are both considered to follow a power law. The diffusive behavior is analyzed and discussed by calculating the mean squared displacement (MSD) analytically and numerically. The results reveal that it displays strong anomalous diffusive behaviors caused by random motions of both states, i.e., two anomalous diffusion terms coexist in the expression of the MSD, and the time distribution which has the heavier tail determines their forms. Moreover, because the two diffusion terms originate from different mechanisms, we find that the diffusion can be characterized by either the term with the largest diffusion exponent or the term with the largest diffusion coefficient at long timescales, which shows very different properties from the single-state process. In addition, the two-state nature of the process of the particle moving in a velocity field makes the TSRW model applicable to describe it. Results obtained from the two-state model reveal that the diffusion can even exhibit subdiffusive behavior, which is significantly different from known results obtained using the single-state model.

Negative DNA supercoiling makes protein-mediated looping deterministic and ergodic within the bacterial doubling time

Protein-mediated DNA looping is fundamental to gene regulation and such loops occur stochastically in purified systems. Additional proteins increase the probability of looping, but these probabilities maintain a broad distribution. For example, the probability of lac repressor-mediated looping in individual molecules ranged 0–100%, and individual molecules exhibited representative behavior only in observations lasting an hour or more. Titrating with HU protein progressively compacted the DNA without narrowing the 0–100% distribution. Increased negative supercoiling produced an ensemble of molecules in which all individual molecules more closely resembled the average. Furthermore, in only 12 min of observation, well within the doubling time of the bacterium, most molecules exhibited the looping probability of the ensemble. DNA supercoiling, an inherent feature of all genomes, appears to impose time-constrained, emergent behavior on otherwise random molecular activity.

Correlated continuous-time random walk in the velocity field: the role of velocity and weak asymptotics

Within the framework of a space-time correlated continuous-time random walk model, anomalous diffusion of particles moving in the velocity field is studied in this paper. The weak asymptotic form ω(t) ∼ t^-(1+α), 1 < α < 2 for large t, is considered to be the waiting time distribution. The analytical results reveal that the diffusion in the velocity field, i.e., the mean squared displacement, can display a multi-fractional form caused by dispersive bias and space-time correlation. The numerical results indicate that the multi-fractional diffusion leads to a crossover phenomenon in-between the process at an intermediate timescale, followed by a steady state which is always determined by the largest diffusion exponent term. In addition, the role of velocity and weak asymptotics is discussed. The extremely small fluid velocity can characterize the diffusion by a diffusion coefficient instead of diffusion exponent, which is distinctly different from the former definition. In particular, for the waiting time displaying a weak asymptotic property, if the anomalous part is suppressed by the normal part, a second crossover phenomenon appears at an intermediate timescale, followed by a steady normal diffusion, which implies that the anomalies underlying the process are smoothed out at large timescales. Moreover, we discuss that the consideration of bias and correlation could help to avoid a possible not readily noticeable mistake in studying the topic concerned in this paper, which may be helpful in the relevant experimental research.

Aging power spectrum of membrane protein transport and other subordinated random walks

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.

Classification of particle trajectories in living cells: Machine learning versus statistical testing hypothesis for fractional anomalous diffusion

Single-particle tracking (SPT) has become a popular tool to study the intracellular transport of molecules in living cells. Inferring the character of their dynamics is important, because it determines the organization and functions of the cells. For this reason, one of the first steps in the analysis of SPT data is the identification of the diffusion type of the observed particles. The most popular method to identify the class of a trajectory is based on the mean-square displacement (MSD). However, due to its known limitations, several other approaches have been already proposed. With the recent advances in algorithms and the developments of modern hardware, the classification attempts rooted in machine learning (ML) are of particular interest. In this work, we adopt two ML ensemble algorithms, i.e., random forest and gradient boosting, to the problem of trajectory classification. We present a new set of features used to transform the raw trajectories data into input vectors required by the classifiers. The resulting models are then applied to real data for G protein-coupled receptors and G proteins. The classification results are compared to recent statistical methods going beyond MSD.

Single-Particle Tracking Reveals Anti-Persistent Subdiffusion in Cell Extracts

Single-particle tracking (SPT) has become a powerful tool to quantify transport phenomena in complex media with unprecedented detail. Based on the reconstruction of individual trajectories, a wealth of informative measures become available for each particle, allowing for a detailed comparison with theoretical predictions. While SPT has been used frequently to explore diffusive transport in artificial fluids and inside living cells, intermediate systems, i.e., biochemically active cell extracts, have been studied only sparsely. Extracts derived from the eggs of the clawfrog Xenopus laevis, for example, are known for their ability to support and mimic vital processes of cells, emphasizing the need to explore also the transport phenomena of nano-sized particles in such extracts. Here, we have performed extensive SPT on beads with 20 nm radius in native and chemically treated Xenopus extracts. By analyzing a variety of distinct measures, we show that these beads feature an anti-persistent subdiffusion that is consistent with fractional Brownian motion. Chemical treatments did not grossly alter this finding, suggesting that the high degree of macromolecular crowding in Xenopus extracts equips the fluid with a viscoelastic modulus, hence enforcing particles to perform random walks with a significant anti-persistent memory kernel.

Polarization-resolved single-molecule tracking reveals strange dynamics of fluorescent tracers through a deep rubbery polymer network

Tracking the movement of fluorescent single-molecule (SM) tracers has provided several new insights into the local structure and dynamics in complex environments such as soft materials and biological systems. However, SM tracking (SMT) remains unreliable at molecular length scales, as the localization error (LE) of SM trajectories (∼30-50 nm) is considerably larger than the size of molecular tracers (∼1-2 nm). Thus, instances of tracer (im)mobility in heterogeneous media, which provide indicators for underlying anomalous-transport mechanisms, remain obscured within the realms of SMT. Since the translation of passive tracers in an isotropic media is associated with fast dipolar rotation, we propose that authentic pauses within the LE can be revealed by probing the hindrance of SM reorientational dynamics. Here, we demonstrate how polarization-resolved SMT (PR-SMT) can provide emission anisotropy at each super-localized position, thereby revealing the tumbling propensity of SMs during random walks. For rhodamine 6G tracers undergoing heterogeneous transport in a hydrated polyvinylpyrrolidone (PVP) network, analysis of PR-SMT trajectories enabled us to discern instances of genuine immobility and localized motion within the LE. Our investigations on 100 SMs in (plasticized) PVP films reveal a wide distribution of dwell times and pause frequencies, demonstrating that most probes intermittently experience complete translational and rotational immobilization. This indicates that tracers serendipitously encounter compact, rigid polymer cavities during transport, implying the existence of nanoscale glass-like domains sparsely distributed in a predominantly deep-rubbery polymer network far above the glass transition.

Lévy walk dynamics in mixed potentials from the perspective of random walk theory

Lévy walk process is one of the most effective models to describe superdiffusion, which underlies some important movement patterns and has been widely observed in micro- and macrodynamics. From the perspective of random walk theory, here we investigate the dynamics of Lévy walks under the influences of the constant force field and the one combined with harmonic potential. Utilizing Hermite polynomial approximation to deal with the spatiotemporally coupled analysis challenges, some striking features are detected, including non-Gaussian stationary distribution, faster diffusion, still strongly anomalous diffusion, etc.

Single molecule characterization of anomalous transport in a thin, anisotropic film

The diffusion of small, charged molecules incorporated in an anisotropic polyelectrolyte multilayer (PEM) was tracked in three dimensions by combining single-molecule fluorescence localization (to characterize lateral diffusion) with Förster resonance energy transfer (FRET) between diffusing molecules and the supporting surface (to measure diffusion in the surface-normal direction). Analysis of the surface-normal diffusion required model-based statistical analysis to account for the inherently noisy FRET signal. Combining these distinct single-molecule methods, which are inherently sensitive to different length-scales, permitted simultaneous characterization of severely anisotropic diffusion, which was more than three orders of magnitude slower in the surface-normal direction. We hypothesize that the anomalously slow surface-normal diffusion was related to the periodic distribution of charge in the PEM, which created electrostatic barriers. The motion was strongly subdiffusive, with anomalous temporal scaling exponents in lateral and normal directions, suggesting a connection to the transient, random fractal conformation of polymer chains in the film's matrix.

Measurements and characterization of the dynamics of tracer particles in an actin network

The underlying physics governing the diffusion of a tracer particle in a viscoelastic material is a topic of some dispute. The long-term memory in the mechanical response of such materials should induce diffusive motion with a memory kernel, such as fractional Brownian motion (fBM). This is the reason that microrheology is able to provide the shear modulus of polymer networks. Surprisingly, the diffusion of a tracer particle in a network of a purified protein, actin, was found to conform to the continuous time random walk type (CTRW). We set out to resolve this discrepancy by studying the tracer particle diffusion using two different tracer particle sizes, in actin networks of different mesh sizes. We find that the ratio of tracer particle size to the characteristic length scale of a bio-polymer network plays a crucial role in determining the type of diffusion it performs. We find that the diffusion of the tracer particles has features of fBm when the particle is large compared to the mesh size, of normal diffusion when the particle is much smaller than the mesh size, and of the CTRW in between these two limits. Based on our findings, we propose and verify numerically a new model for the motion of the tracer in all regimes. Our model suggests that diffusion in actin networks consists of fBm of the tracer particle coupled with caging events with power-law distributed escape times.

Superstatistical modelling of protein diffusion dynamics in bacteria

A recent experiment (Sadoon AA, Wang Y. 2018 Phys. Rev. E98, 042411. (doi:10.1103/PhysRevE.98.042411)) has revealed that nucleoid-associated proteins (i.e. DNA-binding proteins) exhibit highly heterogeneous diffusion processes in bacteria where not only the diffusion constant but also the anomalous diffusion exponent fluctuates for the various proteins. The distribution of displacements of such proteins is observed to take a q-Gaussian form, which decays as a power law. Here, a statistical model is developed for the diffusive motion of the proteins within the bacterium, based on a superstatistics with two variables. This model hierarchically takes into account the joint fluctuations of both the anomalous diffusion exponents and the diffusion constants. A fractional Brownian motion is discussed as a possible local model. Good agreement with the experimental data is obtained.

Heterogeneous diffusion processes and nonergodicity with Gaussian colored noise in layered diffusivity landscapes

Heterogeneous diffusion processes (HDPs) with space-dependent diffusion coefficients D(x) are found in a number of real-world systems, such as for diffusion of macromolecules or submicron tracers in biological cells. Here, we examine HDPs in quenched-disorder systems with Gaussian colored noise (GCN) characterized by a diffusion coefficient with a power-law dependence on the particle position and with a spatially random scaling exponent. Typically, D(x) is considered to be centerd at the origin and the entire x axis is characterized by a single scaling exponent α. In this work we consider a spatially random scenario: in periodic intervals ("layers") in space D(x) is centerd to the midpoint of each interval. In each interval the scaling exponent α is randomly chosen from a Gaussian distribution. The effects of the variation of the scaling exponents, the periodicity of the domains ("layer thickness") of the diffusion coefficient in this stratified system, and the correlation time of the GCN are analyzed numerically in detail. We discuss the regimes of superdiffusion, subdiffusion, and normal diffusion realisable in this system. We observe and quantify the domains where nonergodic and non-Gaussian behaviors emerge in this system. Our results provide new insights into the understanding of weak ergodicity breaking for HDPs driven by colored noise, with potential applications in quenched layered systems, typical model systems for diffusion in biological cells and tissues, as well as for diffusion in geophysical systems.

Correlated continuous-time random walk in a velocity field: Anomalous bifractional crossover

The diffusion of space-time correlated continuous-time random walk moving in the velocity field, which includes the fluid flowing freely and the fluid flowing through porous media, is investigated in this paper. Results reveal that it presents anomalous diffusion merely caused by space-time correlation in the freely flowing fluid, and the bias from the velocity field only supplies a standard advection, which is verified by the corresponding generalized diffusion equation which includes a standard advection term. However, the diffusion in the fluid flowing through porous media, i.e., the mean squared displacement, can display a bifractional form of which one originates from space-time correlation and the other one originates from dispersive bias caused by sticking of the porous media. The fractional advection term emerging in the corresponding generalized diffusion equation confirms the results. Moreover, the coexistence of correlation and dispersive bias result in crossover phenomenon in-between the diffusive process at an intermediate timescale, but just as the definition of diffusion, the one owning the largest diffusion exponent always prevails at large timescales. However, since the two fractional diffusions originate from a different mechanism, even if it owns the smaller diffusion exponent, that one can dominate the diffusion if it fluctuates much stronger than the other one, which no longer obeys the previous conclusion.

Membrane Surface Modulates Slow Diffusion in Small Crowded Droplets

Membranes are ubiquitous structures in cells. The effects of membranes on various functional molecules have been reported, but their behaviors under macromolecular crowding and cell-sized confinement have not fully been understood. In this study, we model an intracellular environment by crowding micrometer-sized droplets and investigate the effects of membrane properties on molecular diffusion. The molecular diffusion inside small droplets covered with a lipid layer of phosphatidylcholine (PC) becomes slower compared with that of the corresponding bulk solutions under a crowding condition of polysaccharide dextran but not of its monomer unit, glucose. The addition of a poly(ethylene glycol) conjugated lipid (PEGylated lipid) to the PC membrane significantly alters the degree of slow diffusion observed inside small droplets of concentrated dextran. Interestingly, the change is not monotonic against dextran concentration; that is, the PEGylated membrane increases and decreases the degree of slow diffusion with increasing dextran concentration. We explain the nonmonotonic alternation from the increase in effective dextran concentration and the hindered temporal adsorption of dextran to the membrane. Because diffusion alteration by adding PEGylated lipid is observed for condensed small droplets of linear polymer PEG and hydrophilic protein bovine serum albumin, the phenomenon is general for other polymer systems as well. Furthermore, our findings may facilitate the understanding of intracellular molecular behaviors based on membrane effects as well as the development of numerous applications using polymer droplets.

Time-dependent classification of protein diffusion types: A statistical detection of mean-squared-displacement exponent transitions

In this paper, we have proposed a statistical procedure for detecting transitions of the mean-square-displacement exponent value within a single trajectory. With this procedure, we have identified three regimes of proteins dynamics on a cell membrane, namely, subdiffusion, free diffusion, and immobility. The fourth considered dynamics type, namely, superdiffusion was not detected. We show that the analyzed protein trajectories are not stationary and not ergodic. Moreover, classification of the dynamics type performed without prior detection of transitions may lead to the overestimation of the proportion of subdiffusive trajectories.

Testing of Multifractional Brownian Motion

Fractional Brownian motion (FBM) is a generalization of the classical Brownian motion. Most of its statistical properties are characterized by the self-similarity (Hurst) index 0<H<1. In nature one often observes changes in the dynamics of a system over time. For example, this is true in single-particle tracking experiments where a transient behavior is revealed. The stationarity of increments of FBM restricts substantially its applicability to model such phenomena. Several generalizations of FBM have been proposed in the literature. One of these is called multifractional Brownian motion (MFBM) where the Hurst index becomes a function of time. In this paper, we introduce a rigorous statistical test on MFBM based on its covariance function. We consider three examples of the functions of the Hurst parameter: linear, logistic, and periodic. We study the power of the test for alternatives being MFBMs with different linear, logistic, and periodic Hurst exponent functions by utilizing Monte Carlo simulations. We also analyze mean-squared displacement (MSD) for the three cases of MFBM by comparing the ensemble average MSD and ensemble average time average MSD, which is related to the notion of ergodicity breaking. We believe that the presented results will be helpful in the analysis of various anomalous diffusion phenomena.

Fractional Dynamics Identification via Intelligent Unpacking of the Sample Autocovariance Function by Neural Networks

Many single-particle tracking data related to the motion in crowded environments exhibit anomalous diffusion behavior. This phenomenon can be described by different theoretical models. In this paper, fractional Brownian motion (FBM) was examined as the exemplary Gaussian process with fractional dynamics. The autocovariance function (ACVF) is a function that determines completely the Gaussian process. In the case of experimental data with anomalous dynamics, the main problem is first to recognize the type of anomaly and then to reconstruct properly the physical rules governing such a phenomenon. The challenge is to identify the process from short trajectory inputs. Various approaches to address this problem can be found in the literature, e.g., theoretical properties of the sample ACVF for a given process. This method is effective; however, it does not utilize all of the information contained in the sample ACVF for a given trajectory, i.e., only values of statistics for selected lags are used for identification. An evolution of this approach is proposed in this paper, where the process is determined based on the knowledge extracted from the ACVF. The designed method is intuitive and it uses information directly available in a new fashion. Moreover, the knowledge retrieval from the sample ACVF vector is enhanced with a learning-based scheme operating on the most informative subset of available lags, which is proven to be an effective encoder of the properties inherited in complex data. Finally, the robustness of the proposed algorithm for FBM is demonstrated with the use of Monte Carlo simulations.

Diffusion of Anisotropic Colloids in Periodic Arrays of Obstacles

Colloidal suspensions in confined geometries exhibit rich diffusion dynamics governed by particle shapes and particle-confinement interactions. Here, we propose a colloidal system, consisting of ellipsoids in periodic array of obstacles, to investigate the confined diffusion of anisotropic colloids. From the obstacle density-dependent diffusion, we discover a decoupling of translational and rotational diffusion in which only rotational motion is localized while translational motion remains diffusive. Moreover, by evaluating the probability distributions of displacements, we found Brownian but non-Gaussian diffusion behaviors with increasing the obstacle densities, which originates from the shape anisotropy of the colloid and the multiplicity of the local configurations of the ellipsoids with respect to the obstacle. Our results suggest that the shape anisotropy and spatial confinements play a vital role in the diffusion dynamics. It is important for understanding the transportations of anisotropic objects in complex environments.

Reflected fractional Brownian motion in one and higher dimensions

Fractional Brownian motion (FBM), a non-Markovian self-similar Gaussian stochastic process with long-ranged correlations, represents a widely applied, paradigmatic mathematical model of anomalous diffusion. We report the results of large-scale computer simulations of FBM in one, two, and three dimensions in the presence of reflecting boundaries that confine the motion to finite regions in space. Generalizing earlier results for finite and semi-infinite one-dimensional intervals, we observe that the interplay between the long-time correlations of FBM and the reflecting boundaries leads to striking deviations of the stationary probability density from the uniform density found for normal diffusion. Particles accumulate at the boundaries for superdiffusive FBM while their density is depleted at the boundaries for subdiffusion. Specifically, the probability density P develops a power-law singularity, P∼r^{κ}, as a function of the distance r from the wall. We determine the exponent κ as a function of the dimensionality, the confining geometry, and the anomalous diffusion exponent α of the FBM. We also discuss implications of our results, including an application to modeling serotonergic fiber density patterns in vertebrate brains.

Macromolecular Crowding: How Shape and Interactions Affect Diffusion

A significant fraction of the cell volume is occupied by various proteins, polysaccharides, nucleic acids, etc., which considerably reduces the mobility of macromolecules. Theoretical and experimental work so far have mainly focused on the dependence of the mobility on the occupied volume, while the effect of a macromolecular shape received less attention. Herein, using fluorescence correlation spectroscopy (FCS) and Brownian dynamics (BD) simulations, we report on a dramatic slowdown of tracer diffusion by cylindrically shaped double-stranded (ds) DNAs (16 nm in length). We find, for instance, that the translational diffusion coefficient of a streptavidin tracer is reduced by about 60% for a volume fraction of dsDNA as low as just 5%. For comparison, for a spherical crowder (Ficoll70) the slowdown is only 10% at the same volume fraction and 60% reduction occurs at a volume fraction as high as 35%. BD simulations reveal that this reduction can be attributed to a larger volume excluded to a tracer by dsDNA particles, as compared with spherical Ficoll70 at the same volume fraction, and to the differences in the tracer-crowder attractive interactions. In addition, we find using BD simulations that rotational diffusion of dsDNA is less affected by the crowder shape than its translational motion. Our results show that diffusion in crowded systems is determined not merely by the occupied volume fraction, but that the shape and interactions can determine diffusion, which is relevant to the diverse intracellular environments inside living cells.

Anomalous diffusion and nonergodicity for heterogeneous diffusion processes with fractional Gaussian noise

Heterogeneous diffusion processes (HDPs) feature a space-dependent diffusivity of the form D(x)=D_{0}|x|^{α}. Such processes yield anomalous diffusion and weak ergodicity breaking, the asymptotic disparity between ensemble and time averaged observables, such as the mean-squared displacement. Fractional Brownian motion (FBM) with its long-range correlated yet Gaussian increments gives rise to anomalous and ergodic diffusion. Here, we study a combined model of HDPs and FBM to describe the particle dynamics in complex systems with position-dependent diffusivity driven by fractional Gaussian noise. This type of motion is, inter alia, relevant for tracer-particle diffusion in biological cells or heterogeneous complex fluids. We show that the long-time scaling behavior predicted theoretically and by simulations for the ensemble- and time-averaged mean-squared displacements couple the scaling exponents α of HDPs and the Hurst exponent H of FBM in a characteristic way. Our analysis of the simulated data in terms of the rescaled variable y∼|x|^{1/(2/(2-α))}/t^{H} coupling particle position x and time t yields a simple, Gaussian probability density function (PDF), P_{HDP-FBM}(y)=e^{-y^{2}}/sqrt[π]. Its universal shape agrees well with theoretical predictions for both uni- and bimodal PDF distributions.

Assessment of MRI-based anomalous diffusion changes in brain ischemic stroke with a fractional motion model

The actual diffusion process in human brain has been shown to be anomalous comparing to that predicted with traditional diffusion MRI (dMRI) theory. Recently, dMRI based on fractional motion (FM) model has demonstrated the potential to accurately describe anomalous diffusion in vivo. In this work, we explored the potential value of FM model-based dMRI in quantificational identification of ischemic stroke and compared that with the traditional apparent diffusion coefficient (ADC). We included 23 acute stroke patients, 8 of whom finished a follow-up scan, and 22 matched healthy controls. The dMRI images were acquired by using a Stejskal-Tanner single-shot spin-echo echo-planar-imaging sequence (diffusion gradients were applied in three orthogonal directions with 25 non-zero b values ranging from 248 to 4474 s/mm²) at 3.0 T MRI. We calculated the coefficient of variation (CV) for FM-related parameters in stroke lesions, and compared the mean values for FM-related parameters and ADC by using two-sample t-tests. Correlation analysis was achieved using Pearson correlation coefficient test. In acute stroke lesions, CV for FM-related parameters showed significant increase compared with normal tissues (P < 0.01), while those of ADC didn't appear statistical difference. Mean values for FM-related parameters showed significant decrease in acute lesion (P < 0.01) and their changing pattern during follow-up was positively correlated with ADC (P < 0.005). Our results initially verified the utility of the FM-model in detecting ischemic stroke compared with traditional dMRI.