Subdiffusion and weak ergodicity breaking in the presence of a reactive boundary

Aging and nonergodicity beyond the Khinchin theorem

The Khinchin theorem provides the condition that a stationary process is ergodic, in terms of the behavior of the corresponding correlation function. Many physical systems are governed by nonstationary processes in which correlation functions exhibit aging. We classify the ergodic behavior of such systems and suggest a possible generalization of Khinchin's theorem. Our work also quantifies deviations from ergodicity in terms of aging correlation functions. Using the framework of the fractional Fokker-Planck equation, we obtain a simple analytical expression for the two-time correlation function of the particle displacement in a general binding potential, revealing universality in the sense that the binding potential only enters into the prefactor through the first two moments of the corresponding Boltzmann distribution. We discuss applications to experimental data from systems exhibiting anomalous dynamics.

Subdiffusive motion of a polymer composed of subdiffusive monomers

We use Brownian dynamics simulations and analytical theory to investigate the physical principles underlying subdiffusive motion of a polymer. Specifically, we examine the consequences of confinement, self-interaction, viscoelasticity, and random waiting on monomer motion, as these physical phenomena may be relevant to the behavior of biological macromolecules in vivo. We find that neither confinement nor self-interaction alter the fundamental Rouse mode relaxations of a polymer. However, viscoelasticity, modeled by fractional Langevin motion, and random waiting, modeled with a continuous time random walk, lead to significant and distinct deviations from the classic polymer-dynamics model. Our results provide diagnostic tools--the monomer mean square displacement scaling and the velocity autocorrelation function--that can be applied to experimental data to determine the underlying mechanism for subdiffusive motion of a polymer.

Subdiffusion in a bounded domain with a partially absorbing-reflecting boundary

The exit time of a subdiffusive process from a bounded domain with a partially absorbing/reflecting boundary is considered. The short-time and long-time behaviors of the exit time probability density are investigated by using a spectral decomposition on the basis of the Laplace operator eigenfunctions. Rotation-invariant domains are analyzed in depth in order to illustrate the use of theoretical formulas and to compare them to numerical simulations. The asymptotic results obtained are relevant for describing subdiffusion inside a living cell with a semipermeable membrane, in a chemical reactor filled with catalytic grains of finite reactivity, or in mineral or biological samples which are probed by nuclear magnetic resonance measurements subject to surface relaxation.

Spatio-temporal correlations can drastically change the response of a MAPK pathway

Multisite covalent modification of proteins is omnipresent in eukaryotic cells. A well-known example is the mitogen-activated protein kinase (MAPK) cascade where, in each layer of the cascade, a protein is phosphorylated at two sites. It has long been known that the response of a MAPK pathway strongly depends on whether the enzymes that modify the protein act processively or distributively. A distributive mechanism, in which the enzyme molecules have to release the substrate molecules in between the modification of the two sites, can generate an ultrasensitive response and lead to hysteresis and bistability. We study by Green's Function Reaction Dynamics (GFRD), a stochastic scheme that makes it possible to simulate biochemical networks at the particle level in time and space, a dual phosphorylation cycle in which the enzymes act according to a distributive mechanism. We find that the response of this network can differ dramatically from that predicted by a mean-field analysis based on the chemical rate equations. In particular, rapid rebindings of the enzyme molecules to the substrate molecules after modification of the first site can markedly speed up the response and lead to loss of ultrasensitivity and bistability. In essence, rapid enzyme-substrate rebindings can turn a distributive mechanism into a processive mechanism. We argue that slow ADP release by the enzymes can protect the system against these rapid rebindings, thus enabling ultrasensitivity and bistability.

Survival probability of a subdiffusive particle in a d-dimensional sea of mobile traps

We investigate the long-time behavior of the survival probability P(t) of a mobile particle in d-dimensional continuous Euclidean media doped with noninteracting mobile traps. The particle is strictly subdiffusive, implying that its mean-square displacement grows as tgamma' with 0<gamma'<1. Initially, the traps are scattered randomly and their subsequent mean-square displacement grows as tgamma with 0<gamma<or=1. Instantaneous annihilation of the particle takes place upon contact with any of the traps. The solution to this problem is obtained by deriving lower and upper asymptotic bounds of the survival probability and showing that they converge to one another for long times, thereby unambiguously determining the long-time decay of P(t). For d>or=2 we find that at late times the survival probability is that of the pure target problem (the problem where the particle remains immobile) in agreement with previous studies for the d=1 case. These decay laws remain invariant over the whole gamma range as opposed to the dynamical crossover observed for the case of a purely diffusive particle (gamma'=1) where, for gamma<2/(2+d) , the survival probability becomes that of the so-called trapping problem (the problem where the particle moves in a sea of static traps). This behavior implies that for sufficiently low values of gamma(gamma<2/(2+d)) the survival probability becomes singular in the limit gamma'-->1: trappinglike for gamma'=1 and targetlike for any gamma'<1.

How subdiffusion changes the kinetics of binding to a surface

Under molecular crowding conditions, biopolymers have been reported to subdiffuse, (r(2)(t)) approximately = t(alpha), with 0 <alpha < 1. Here we study the exchange dynamics of such a subdiffusing particle with a reactive boundary using a continuous time random walk approach. We derive the generalized boundary condition and consider the unbinding from the boundary. An ensuing weak ergodicity breaking has profound consequences for material exchange between the boundary and bulk. We discuss the effects in biological contexts such as gene regulation or membrane-bulk exchange processes. We also suggest various methods to experimentally probe the subdiffusive behavior.

Subdiffusion in time-averaged, confined random walks

Certain techniques characterizing diffusive processes, such as single-particle tracking or molecular dynamics simulation, provide time averages rather than ensemble averages. Whereas the ensemble-averaged mean-squared displacement (MSD) of an unbounded continuous time random walk (CTRW) with a broad distribution of waiting times exhibits subdiffusion, the time-averaged MSD, delta2, does not. We demonstrate that, in contrast to the unbounded CTRW, in which delta2 is linear in the lag time Delta, the time-averaged MSD of the CTRW of a walker confined to a finite volume is sublinear in Delta, i.e., for long lag times delta2 approximately Delta1-alpha. The present results permit the application of CTRW to interpret time-averaged experimental quantities.

Quantifying hopping and jumping in facilitated diffusion of DNA-binding proteins

Facilitated diffusion of DNA-binding proteins is known to speed up target site location by combining three dimensional excursions and linear diffusion along the DNA. Here we explicitly calculate the distribution of the relocation lengths of such 3D excursions, and we quantify the short-range correlated excursions, also called hops, and the long-range uncorrelated jumps. Our results substantiate recent single-molecule experiments that reported sliding and 3D excursions of the restriction enzyme EcoRV on elongated DNA molecules. We extend our analysis to the case of anomalous 3D diffusion, likely to occur in a crowded cellular medium.

Facilitated diffusion with DNA coiling

When DNA-binding proteins search for their specific binding site on a DNA molecule they alternate between linear 1-dimensional diffusion along the DNA molecule, mediated by nonspecific binding, and 3-dimensional volume excursion events between successive dissociation from and rebinding to DNA. If the DNA molecule is kept in a straight configuration, for instance, by optical tweezers, these 3-dimensional excursions may be divided into long volume excursions and short hops along the DNA. These short hops correspond to immediate rebindings after dissociation such that a rebinding event to the DNA occurs at a site that is close to the site of the preceding dissociation. When the DNA molecule is allowed to coil up, immediate rebinding may also lead to so-called intersegmental jumps, i.e., immediate rebindings to a DNA segment that is far away from the unbinding site when measured in the chemical distance along the DNA, but close by in the embedding 3-dimensional space. This effect is made possible by DNA looping. The significance of intersegmental jumps was recently demonstrated in a single DNA optical tweezers setup. Here we present a theoretical approach in which we explicitly take the effect of DNA coiling into account. By including the spatial correlations of the short hops we demonstrate how the facilitated diffusion model can be extended to account for intersegmental jumping at varying DNA densities. It is also shown that our approach provides a quantitative interpretation of the experimentally measured enhancement of the target location by DNA-binding proteins.

Random time-scale invariant diffusion and transport coefficients

Single particle tracking of mRNA molecules and lipid granules in living cells shows that the time averaged mean squared displacement delta2[over ] of individual particles remains a random variable while indicating that the particle motion is subdiffusive. We investigate this type of ergodicity breaking within the continuous time random walk model and show that delta2[over ] differs from the corresponding ensemble average. In particular we derive the distribution for the fluctuations of the random variable delta2[over ]. Similarly we quantify the response to a constant external field, revealing a generalization of the Einstein relation. Consequences for the interpretation of single molecule tracking data are discussed.

Fluorescence recovery after photobleaching: the case of anomalous diffusion

The method of FRAP (fluorescence recovery after photobleaching), which has been broadly used to measure lateral mobility of fluorescent-labeled molecules in cell membranes, is formulated here in terms of continuous time random walks (CTRWs), which offer both analytical expressions and a scheme for numerical simulations. We propose an approach based on the CTRW and the corresponding fractional diffusion equation (FDE) to analyze FRAP results in the presence of anomalous subdiffusion. The FDE generalizes the simple diffusive picture, which has been applied to FRAP when assuming regular diffusion, to account for subdiffusion. We use a subordination relationship between the solutions of the fractional and normal diffusion equations to fit FRAP recovery curves obtained from CTRW simulations, and compare the fits to the commonly used approach based on the simple diffusion equation with a time dependent diffusion coefficient (TDDC). The CTRW and TDDC describe two different dynamical schemes, and although the CTRW formalism appears to be more complicated, it provides a physical description that underlies anomalous lateral diffusion.

Probing microscopic origins of confined subdiffusion by first-passage observables

Subdiffusive motion of tracer particles in complex crowded environments, such as biological cells, has been shown to be widespread. This deviation from Brownian motion is usually characterized by a sublinear time dependence of the mean square displacement (MSD). However, subdiffusive behavior can stem from different microscopic scenarios that cannot be identified solely by the MSD data. In this article we present a theoretical framework that permits the analytical calculation of first-passage observables (mean first-passage times, splitting probabilities, and occupation times distributions) in disordered media in any dimensions. This analysis is applied to two representative microscopic models of subdiffusion: continuous-time random walks with heavy tailed waiting times and diffusion on fractals. Our results show that first-passage observables provide tools to unambiguously discriminate between the two possible microscopic scenarios of subdiffusion. Moreover, we suggest experiments based on first-passage observables that could help in determining the origin of subdiffusion in complex media, such as living cells, and discuss the implications of anomalous transport to reaction kinetics in cells.

Distribution of time-averaged observables for weak ergodicity breaking

We find a general formula for the distribution of time-averaged observables for systems modeled according to the subdiffusive continuous time random walk. For Gaussian random walks coupled to a thermal bath we recover ergodicity and Boltzmann's statistics, while for the anomalous subdiffusive case a weakly nonergodic statistical mechanical framework is constructed, which is based on Lévy's generalized central limit theorem. As an example we calculate the distribution of X, the time average of the position of the particle, for unbiased and uniformly biased particles, and show that X exhibits large fluctuations compared with the ensemble average <X>.

Subdiffusive target problem: survival probability

The asymptotic survival probability of a spherical target in the presence of a single subdiffusive trap or surrounded by a sea of subdiffusive traps in a continuous Euclidean medium is calculated. In one and two dimensions the survival probability of the target in the presence of a single trap decays to zero as a power law and as a power law with logarithmic correction, respectively. The target is thus reached with certainty, but it takes the trap an infinite time on average to do so. In dimensions higher than two a single trap may never reach the target and so the survival probability is finite and, in fact, does not depend on whether the traps move diffusively or subdiffusively. When the target is surrounded by a sea of traps, on the other hand, its survival probability decays as a stretched exponential in all dimensions (with a logarithmic correction in the exponent for d=2). A trap will therefore reach the target with certainty, and will do so in a finite time. These results may be directly related to enzyme binding kinetics on DNA in the crowded cellular environment.

Modeling diffusional transport in the interphase cell nucleus

In this paper a lattice model for the diffusional transport of particles in the interphase cell nucleus is proposed. Dense networks of chromatin fibers are created by three different methods: Randomly distributed, noninterconnected obstacles, a random walk chain model, and a self-avoiding random walk chain model with persistence length. By comparing a discrete and a continuous version of the random walk chain model, we demonstrate that lattice discretization does not alter particle diffusion. The influence of the three dimensional geometry of the fiber network on the particle diffusion is investigated in detail while varying the occupation volume, chain length, persistence length, and walker size. It is shown that adjacency of the monomers, the excluded volume effect incorporated in the self-avoiding random walk model, and, to a lesser extent, the persistence length affect particle diffusion. It is demonstrated how the introduction of the effective chain occupancy, which is a convolution of the geometric chain volume with the walker size, eliminates the conformational effects of the network on the diffusion, i.e., when plotting the diffusion coefficient as a function of the effective chain volume, the data fall onto a master curve.

Sampling the cell with anomalous diffusion - the discovery of slowness

Diffusion-mediated searching for interaction partners is an ubiquitous process in cell biology. Transcription factors, for example, search specific DNA sequences, signaling proteins aim at interacting with specific cofactors, and peripheral membrane proteins try to dock to membrane domains. Brownian motion, however, is affected by molecular crowding that induces anomalous diffusion (so-called subdiffusion) of proteins and larger structures, thereby compromising diffusive transport and the associated sampling processes. Contrary to the naive expectation that subdiffusion obstructs cellular processes, we show here by computer simulations that subdiffusion rather increases the probability of finding a nearby target. Consequently, important events like protein complex formation and signal propagation are enhanced as compared to normal diffusion. Hence, cells indeed benefit from their crowded internal state and the associated anomalous diffusion.

Fractional diffusion equation in a confined region: surface effects and exact solutions

Surface effects on a diffusion process governed by a fractional diffusion equation in a confined region with spatial and time dependent boundary conditions are investigated. First, we consider the one-dimensional case with the boundary conditions rho(0,t)=Phi0(t) and rho(a,t)=Phia(t). Subsequently, the two-dimensional case in the cylindrical symmetry with rho(a,theta,t)=Phia(theta,t) and rho(b,theta,t)=Phib(theta,t) is investigated. For these cases, we also obtain exact solutions for an arbitrary initial condition by using the Green's function approach.