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)
A topic of intense current investigation pursues the question of how the highly crowded environment of biological cells affects the dynamic properties of passively diffusing particles. Motivated by recent experiments we report results of extensive simulations of the motion of a finite sized tracer particle in a heterogeneously crowded environment made up of quenched distributions of monodisperse crowders of varying sizes in finite circular two-dimensional domains. For given spatial distributions of monodisperse crowders we demonstrate how anomalous diffusion with strongly non-Gaussian features arises in this model system. We investigate both biologically relevant situations of particles released either at the surface of an inner domain or at the outer boundary, exhibiting distinctly different features of the observed anomalous diffusion for heterogeneous distributions of crowders. Specifically we reveal an asymmetric spreading of tracers even at moderate crowding. In addition to the mean squared displacement (MSD) and local diffusion exponent we investigate the magnitude and the amplitude scatter of the time averaged MSD of individual tracer trajectories, the non-Gaussianity parameter, and the van Hove correlation function. We also quantify how the average tracer diffusivity varies with the position in the domain with a heterogeneous radial distribution of crowders and examine the behaviour of the survival probability and the dynamics of the tracer survival probability. Inter alia, the systems we investigate are related to the passive transport of lipid molecules and proteins in two-dimensional crowded membranes or the motion in colloidal solutions or emulsions in effectively two-dimensional geometries, as well as inside supercrowded, surface adhered cells.
We develop an encounter-based approach for describing restricted diffusion with a gradient drift toward a partially reactive boundary. For this purpose, we introduce an extension of the Dirichlet-to-Neumann operator and use its eigenbasis to derive a spectral decomposition for the full propagator, i.e. the joint probability density function for the particle position and its boundary local time. This is the central quantity that determines various characteristics of diffusion-influenced reactions such as conventional propagators, survival probability, first-passage time distribution, boundary local time distribution, and reaction rate. As an illustration, we investigate the impact of a constant drift onto the boundary local time for restricted diffusion on an interval. More generally, this approach accesses how external forces may influence the statistics of encounters of a diffusing particle with the reactive boundary.
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.
We investigate a class of diffusion-controlled reactions that are initiated at the time instance when a prescribed number K among N particles independently diffusing in a solvent are simultaneously bound to a target region.
In the irreversible target-binding setting, the particles that bind to the target stay there forever, and the reaction time is the Kth fastest first-passage time to the target, whose distribution is well-known. In turn, reversible binding, which is common for most applications, renders theoretical analysis much more challenging and drastically changes the distribution of reaction times.
We develop a renewal-based approach to derive an approximate solution for the probability density of the reaction time.
This approximation turns out to be remarkably accurate for a broad range of parameters.
We also analyze the dependence of the mean reaction time or, equivalently, the inverse reaction rate, on the main parameters such as K, N, and binding/unbinding constants. Some biophysical applications and further perspectives are briefly discussed.
We study the probability density function (PDF) of the first-reaction times between a diffusive ligand and a membrane-bound, immobile imperfect target region in a restricted 'onion-shell' geometry bounded by two nested membranes of arbitrary shapes. For such a setting, encountered in diverse molecular signal transduction pathways or in the narrow escape problem with additional steric constraints, we derive an exact spectral form of the PDF, as well as present its approximate form calculated by help of the so-called self-consistent approximation. For a particular case when the nested domains are concentric spheres, we get a fully explicit form of the approximated PDF, assess the accuracy of this approximation, and discuss various facets of the obtained distributions. Our results can be straightforwardly applied to describe the PDF of the terminal reaction event in multi-stage signal transduction processes.
We study the probability density function (PDF) of the first-reaction times between a diffusive ligand and a membrane-bound, immobile imperfect target region in a restricted 'onion-shell' geometry bounded by two nested membranes of arbitrary shapes. For such a setting, encountered in diverse molecular signal transduction pathways or in the narrow escape problem with additional steric constraints, we derive an exact spectral form of the PDF, as well as present its approximate form calculated by help of the so-called self-consistent approximation. For a particular case when the nested domains are concentric spheres, we get a fully explicit form of the approximated PDF, assess the accuracy of this approximation, and discuss various facets of the obtained distributions. Our results can be straightforwardly applied to describe the PDF of the terminal reaction event in multi-stage signal transduction processes.
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.
We consider a sequential cascade of molecular first-reaction events towards a terminal reaction centre in which each reaction step is controlled by diffusive motion of the particles. The model studied here represents a typical reaction setting encountered in diverse molecular biology systems, in which, e.g. a signal transduction proceeds via a series of consecutive 'messengers': the first messenger has to find its respective immobile target site triggering a launch of the second messenger, the second messenger seeks its own target site and provokes a launch of the third messenger and so on, resembling a relay race in human competitions. For such a molecular relay race taking place in infinite one-, two- and three-dimensional systems, we find exact expressions for the probability density function of the time instant of the terminal reaction event, conditioned on preceding successful reaction events on an ordered array of target sites. The obtained expressions pertain to the most general conditions: number of intermediate stages and the corresponding diffusion coefficients, the sizes of the target sites, the distances between them, as well as their reactivities are arbitrary.
We consider the first-passage problem for N identical independent particles that are initially released uniformly in a finite domain Ω and then diffuse toward a reactive area Γ, which can be part of the outer boundary of Ω or a reaction centre in the interior of Ω. For both cases of perfect and partial reactions, we obtain the explicit formulas for the first two moments of the fastest first-passage time (fFPT), i.e., the time when the first out of the N particles reacts with Γ. Moreover, we investigate the full probability density of the fFPT. We discuss a significant role of the initial condition in the scaling of the average fFPT with the particle number N, namely, a much stronger dependence (1/N and 1/N² for partially and perfectly reactive targets, respectively), in contrast to the well known inverse-logarithmic behaviour found when all particles are released from the same fixed point. We combine analytic solutions with scaling arguments and stochastic simulations to rationalise our results, which open new perspectives for studying the relevance of multiple searchers in various situations of molecular reactions, in particular, in living cells.
Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains
(2017)
We study the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry-characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. A similar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA. We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters. We analyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations.