Refine
Language
- English (25)
Is part of the Bibliography
- yes (25)
Keywords
- diffusion (11)
- first passage time (5)
- random diffusivity (4)
- Brownian motion (3)
- aspect ratio (3)
- cylindrical geometry (3)
- first-passage time (3)
- protein search (3)
- approximate methods (2)
- extremal values (2)
Institute
- Institut für Physik und Astronomie (23)
- Extern (1)
Stochastic models based on random diffusivities, such as the diffusing-diffusivity approach, are popular concepts for the description of non-Gaussian diffusion in heterogeneous media. Studies of these models typically focus on the moments and the displacement probability density function. Here we develop the complementary power spectral description for a broad class of random-diffusivity processes. In our approach we cater for typical single particle tracking data in which a small number of trajectories with finite duration are garnered. Apart from the diffusing-diffusivity model we study a range of previously unconsidered random-diffusivity processes, for which we obtain exact forms of the probability density function. These new processes are different versions of jump processes as well as functionals of Brownian motion. The resulting behaviour subtly depends on the specific model details. Thus, the central part of the probability density function may be Gaussian or non-Gaussian, and the tails may assume Gaussian, exponential, log-normal, or even power-law forms. For all these models we derive analytically the moment-generating function for the single-trajectory power spectral density. We establish the generic 1/f²-scaling of the power spectral density as function of frequency in all cases. Moreover, we establish the probability density for the amplitudes of the random power spectral density of individual trajectories. The latter functions reflect the very specific properties of the different random-diffusivity models considered here. Our exact results are in excellent agreement with extensive numerical simulations.
Stochastic models based on random diffusivities, such as the diffusing-diffusivity approach, are popular concepts for the description of non-Gaussian diffusion in heterogeneous media. Studies of these models typically focus on the moments and the displacement probability density function. Here we develop the complementary power spectral description for a broad class of random-diffusivity processes. In our approach we cater for typical single particle tracking data in which a small number of trajectories with finite duration are garnered. Apart from the diffusing-diffusivity model we study a range of previously unconsidered random-diffusivity processes, for which we obtain exact forms of the probability density function. These new processes are different versions of jump processes as well as functionals of Brownian motion. The resulting behaviour subtly depends on the specific model details. Thus, the central part of the probability density function may be Gaussian or non-Gaussian, and the tails may assume Gaussian, exponential, log-normal, or even power-law forms. For all these models we derive analytically the moment-generating function for the single-trajectory power spectral density. We establish the generic 1/f²-scaling of the power spectral density as function of frequency in all cases. Moreover, we establish the probability density for the amplitudes of the random power spectral density of individual trajectories. The latter functions reflect the very specific properties of the different random-diffusivity models considered here. Our exact results are in excellent agreement with extensive numerical simulations.
Towards a full quantitative description of single-molecule reaction kinetics in biological cells
(2018)
The first-passage time (FPT), i.e., the moment when a stochastic process reaches a given threshold value for the first time, is a fundamental mathematical concept with immediate applications. In particular, it quantifies the statistics of instances when biomolecules in a biological cell reach their specific binding sites and trigger cellular regulation. Typically, the first-passage properties are given in terms of mean first-passage times. However, modern experiments now monitor single-molecular binding-processes in living cells and thus provide access to the full statistics of the underlying first-passage events, in particular, inherent cell-to-cell fluctuations. We here present a robust explicit approach for obtaining the distribution of FPTs to a small partially reactive target in cylindrical-annulus domains, which represent typical bacterial and neuronal cell shapes. We investigate various asymptotic behaviours of this FPT distribution and show that it is typically very broad in many biological situations, thus, the mean FPT can differ from the most probable FPT by orders of magnitude. The most probable FPT is shown to strongly depend only on the starting position within the geometry and to be almost independent of the target size and reactivity. These findings demonstrate the dramatic relevance of knowing the full distribution of FPTs and thus open new perspectives for a more reliable description of many intracellular processes initiated by the arrival of one or few biomolecules to a small, spatially localised region inside the cell.
High diffusion-sensitizing magnetic field gradients have been more and more often applied nowadays to achieve a better characterization of the microstructure. As the resulting spin-echo signal significantly deviates from the conventional Gaussian form, various models have been employed to interpret these deviations and to relate them with the microstructural properties of a sample. In this paper, we argue that the non-Gaussian behavior of the signal is a generic universal feature of the Bloch-Torrey equation. We provide a simple yet rigorous description of the localization regime emerging at high extended gradients and identify its origin as a symmetry breaking at the reflecting boundary. We compare the consequent non-Gaussian signal decay to other diffusion NMR regimes such as slow-diffusion, motional-narrowing and diffusion-diffraction regimes. We emphasize limitations of conventional perturbative techniques and advocate for non-perturbative approaches which may pave a way to new imaging modalities in this field.
Textbook concepts of diffusion-versus kinetic-control are well-defined for reaction-kinetics involving macroscopic concentrations of diffusive reactants that are adequately described by rate-constants—the inverse of the mean-first-passage-time to the reaction-event. In contradiction, an open important question is whether the mean-first-passage-time alone is a sufficient measure for biochemical reactions that involve nanomolar reactant concentrations. Here, using a simple yet generic, exactly solvable model we study the effect of diffusion and chemical reaction-limitations on the full reaction-time distribution. We show that it has a complex structure with four distinct regimes delineated by three characteristic time scales spanning a window of several decades. Consequently, the reaction-times are defocused: no unique time-scale characterises the reaction-process, diffusion- and kinetic-control can no longer be disentangled, and it is imperative to know the full reaction-time distribution. We introduce the concepts of geometry- and reaction-control, and also quantify each regime by calculating the corresponding reaction depth.
Textbook concepts of diffusion-versus kinetic-control are well-defined for reaction-kinetics involving macroscopic concentrations of diffusive reactants that are adequately described by rate-constants—the inverse of the mean-first-passage-time to the reaction-event. In contradiction, an open important question is whether the mean-first-passage-time alone is a sufficient measure for biochemical reactions that involve nanomolar reactant concentrations. Here, using a simple yet generic, exactly solvable model we study the effect of diffusion and chemical reaction-limitations on the full reaction-time distribution. We show that it has a complex structure with four distinct regimes delineated by three characteristic time scales spanning a window of several decades. Consequently, the reaction-times are defocused: no unique time-scale characterises the reaction-process, diffusion- and kinetic-control can no longer be disentangled, and it is imperative to know the full reaction-time distribution. We introduce the concepts of geometry- and reaction-control, and also quantify each regime by calculating the corresponding reaction depth.
Diffusive search for a static target is a common problem in statistical physics with numerous applications in chemistry and biology.
We look at this problem from a different perspective and investigate the statistics of encounters between the diffusing particle and the target. While an exact solution of this problem was recently derived in the form of a spectral expansion over the eigenbasis of the Dirichlet-to-Neumann operator, the latter is generally difficult to access for an arbitrary target.
In this paper, we present three complementary approaches to approximate the probability density of the rescaled number of encounters with a small target in a bounded confining domain. In particular, we derive a simple fully explicit approximation, which depends only on a few geometric characteristics such as the surface area and the harmonic capacity of the target, and the volume of the confining domain.
We discuss the advantages and limitations of three approaches and check their accuracy.
We also deduce an explicit approximation for the distribution of the first-crossing time, at which the number of encounters exceeds a prescribed threshold. Its relations to common first-passage time problems are discussed.
The time instant-the first-passage time (FPT)-when a diffusive particle (e.g., a ligand such as oxygen or a signalling protein) for the first time reaches an immobile target located on the surface of a bounded three-dimensional domain (e.g., a hemoglobin molecule or the cellular nucleus) is a decisive characteristic time-scale in diverse biophysical and biochemical processes, as well as in intermediate stages of various inter- and intra-cellular signal transduction pathways. Adam and Delbruck put forth the reduction-of-dimensionality concept, according to which a ligand first binds non-specifically to any point of the surface on which the target is placed and then diffuses along this surface until it locates the target. In this work, we analyse the efficiency of such a scenario and confront it with the efficiency of a direct search process, in which the target is approached directly from the bulk and not aided by surface diffusion. We consider two situations: (i) a single ligand is launched from a fixed or a random position and searches for the target, and (ii) the case of 'amplified' signals when N ligands start either from the same point or from random positions, and the search terminates when the fastest of them arrives to the target. For such settings, we go beyond the conventional analyses, which compare only the mean values of the corresponding FPTs. Instead, we calculate the full probability density function of FPTs for both scenarios and study its integral characteristic-the 'survival' probability of a target up to time t. On this basis, we examine how the efficiencies of both scenarios are controlled by a variety of parameters and single out realistic conditions in which the reduction-of-dimensionality scenario outperforms the direct search.
In the scenario of the narrow escape problem (NEP) a particle diffuses in a finite container and eventually leaves it through a small 'escape window' in the otherwise impermeable boundary, once it arrives to this window and crosses an entropic barrier at the entrance to it. This generic problem is mathematically identical to that of a diffusion-mediated reaction with a partially-reactive site on the container's boundary. Considerable knowledge is available on the dependence of the mean first-reaction time (FRT) on the pertinent parameters. We here go a distinct step further and derive the full FRT distribution for the NEP. We demonstrate that typical FRTs may be orders of magnitude shorter than the mean one, thus resulting in a strong defocusing of characteristic temporal scales. We unveil the geometry-control of the typical times, emphasising the role of the initial distance to the target as a decisive parameter. A crucial finding is the further FRT defocusing due to the barrier, necessitating repeated escape or reaction attempts interspersed with bulk excursions. These results add new perspectives and offer a broad comprehension of various features of the by-now classical NEP that are relevant for numerous biological and technological systems.