Towards a full quantitative description of single-molecule reaction kinetics in biological cells

Controlling Uncertainty of Empirical First-Passage Times in the Small-Sample Regime

We derive general bounds on the probability that the empirical first-passage time τ[over ¯]_{n}≡∑_{i=1}^{n}τ_{i}/n of a reversible ergodic Markov process inferred from a sample of n independent realizations deviates from the true mean first-passage time by more than any given amount in either direction. We construct nonasymptotic confidence intervals that hold in the elusive small-sample regime and thus fill the gap between asymptotic methods and the Bayesian approach that is known to be sensitive to prior belief and tends to underestimate uncertainty in the small-sample setting. We prove sharp bounds on extreme first-passage times that control uncertainty even in cases where the mean alone does not sufficiently characterize the statistics. Our concentration-of-measure-based results allow for model-free error control and reliable error estimation in kinetic inference, and are thus important for the analysis of experimental and simulation data in the presence of limited sampling.

Diffusion-Controlled Reactions: An Overview

We review the milestones in the century-long development of the theory of diffusion-controlled reactions. Starting from the seminal work by von Smoluchowski, who recognized the importance of diffusion in chemical reactions, we discuss perfect and imperfect surface reactions, their microscopic origins, and the underlying mathematical framework. Single-molecule reaction schemes, anomalous bulk diffusions, reversible binding/unbinding kinetics, and many other extensions are presented. An alternative encounter-based approach to diffusion-controlled reactions is introduced, with emphasis on its advantages and potential applications. Some open problems and future perspectives are outlined.

Encounter-based approach to the escape problem

We revise the encounter-based approach to imperfect diffusion-controlled reactions, which employs the statistics of encounters between a diffusing particle and the reactive region to implement surface reactions. We extend this approach to deal with a more general setting, in which the reactive region is surrounded by a reflecting boundary with an escape region. We derive a spectral expansion for the full propagator and investigate the behavior and probabilistic interpretations of the associated probability flux density. In particular, we obtain the joint probability density of the escape time and the number of encounters with the reactive region before escape, and the probability density of the first-crossing time of a prescribed number of encounters. We briefly discuss generalizations of the conventional Poissonian-type surface reaction mechanism described by Robin boundary condition and potential applications of this formalism in chemistry and biophysics.

Imperfect narrow escape problem

We consider the kinetics of the imperfect narrow escape problem, i.e., the time it takes for a particle diffusing in a confined medium of generic shape to reach and to be adsorbed by a small, imperfectly reactive patch embedded in the boundary of the domain, in two or three dimensions. Imperfect reactivity is modeled by an intrinsic surface reactivity κ of the patch, giving rise to Robin boundary conditions. We present a formalism to calculate the exact asymptotics of the mean reaction time in the limit of large volume of the confining domain. We obtain exact explicit results in the two limits of large and small reactivities of the reactive patch, and a semianalytical expression in the general case. Our approach reveals an anomalous scaling of the mean reaction time as the inverse square root of the reactivity in the large-reactivity limit, valid for an initial position near the extremity of the reactive patch. We compare our exact results with those obtained within the "constant flux approximation"; we show that this approximation turns out to give exactly the next-to-leading-order term of the small-reactivity limit, and provides a good approximation of the reaction time far from the reactive patch for all reactivities, but not in the vicinity of the boundary of the reactive patch due to the above-mentioned anomalous scaling. These results thus provide a general framework to quantify the mean reaction times for the imperfect narrow escape problem.

Theoretical insights into the full description of DNA target search by subdiffusing proteins

DNA binding proteins (DBPs) diffuse in the cytoplasm to recognise and bind with their respective target sites on the DNA to initiate several biologically important processes. The first passage time distributions (FPTDs) of DBPs are useful in quantifying the timescales of the most-probable search paths in addition to the mean value of the distribution which, strikingly, are decades of order apart in time. However, extremely crowded in vivo conditions or the viscoelasticity of the cellular medium among other factors causes biomolecules to exhibit anomalous diffusion which is usually overlooked in most theoretical studies. We have obtained approximate analytical expressions of a general FPTD and the two characteristic timescales that are valid for any single subdiffusing protein searching for its target in vivo. Our results can be applied to single-particle tracking experiments of target search.

Mean first-passage time to a small absorbing target in three-dimensional elongated domains

We derive an approximate formula for the mean first-passage time (MFPT) to a small absorbing target of arbitrary shape inside an elongated domain of a slowly varying axisymmetric profile. For this purpose, the original Poisson equation in three dimensions is reduced to an effective one-dimensional problem on an interval with a semipermeable semiabsorbing membrane. The approximate formula captures correctly the dependence of the MFPT on the distance to the target, the radial profile of the domain, and the size and the shape of the target. This approximation is validated by Monte Carlo simulations.

First-passage times to anisotropic partially reactive targets

We investigate restricted diffusion in a bounded domain towards a small partially reactive target in three- and higher-dimensional spaces. We propose a simple explicit approximation for the principal eigenvalue of the Laplace operator with mixed Robin-Neumann boundary conditions. This approximation involves the harmonic capacity and the surface area of the target, the volume of the confining domain, the diffusion coefficient, and the reactivity. The accuracy of the approximation is checked by using a finite-elements method. The proposed approximation determines also the mean first-reaction time, the long-time decay of the survival probability, and the overall reaction rate on that target. We identify the relevant lengthscale of the target, which determines its trapping capacity, and we investigate its relation to the target shape. In particular, we study the effect of target anisotropy on the principal eigenvalue by computing the harmonic capacity of prolate and oblate spheroids in various space dimensions. Some implications of these results in chemical physics and biophysics are briefly discussed.

First-encounter time of two diffusing particles in two- and three-dimensional confinement

The statistics of the first-encounter time of diffusing particles changes drastically when they are placed under confinement. In the present work, we make use of Monte Carlo simulations to study the behavior of a two-particle system in two- and three-dimensional domains with reflecting boundaries. Based on the outcome of the simulations, we give a comprehensive overview of the behavior of the survival probability S(t) and the associated first-encounter time probability density H(t) over a broad time range spanning several decades. In addition, we provide numerical estimates and empirical formulas for the mean first-encounter time 〈T〉, as well as for the decay time T characterizing the monoexponential long-time decay of the survival probability. Based on the distance between the boundary and the center of mass of two particles, we obtain an empirical lower bound t_{B} for the time at which S(t) starts to significantly deviate from its counterpart for the no boundary case. Surprisingly, for small-sized particles, the dominant contribution to T depends only on the total diffusivity D=D_{1}+D_{2}, in sharp contrast to the one-dimensional case. This contribution can be related to the Wiener sausage generated by a fictitious Brownian particle with diffusivity D. In two dimensions, the first subleading contribution to T is found to depend weakly on the ratio D_{1}/D_{2}. We also investigate the slow-diffusion limit when D_{2}≪D_{1}, and we discuss the transition to the limit when one particle is a fixed target. Finally, we give some indications to anticipate when T can be expected to be a good approximation for 〈T〉.

Reversible target-binding kinetics of multiple impatient particles

Certain biochemical reactions can only be triggered after binding a sufficient number of particles to a specific target region such as an enzyme or a protein sensor. We investigate the distribution of the reaction time, i.e., the first instance when all independently diffusing particles are bound to the target. When each particle binds irreversibly, this is equivalent to the first-passage time of the slowest (last) particle. In turn, reversible binding to the target renders the problem much more challenging and drastically changes the distribution of the reaction time. We derive the exact solution of this problem and investigate the short-time and long-time asymptotic behaviors of the reaction time probability density. We also analyze how the mean reaction time depends on the unbinding rate and the number of particles. Our exact and asymptotic solutions are compared to Monte Carlo simulations.

Universal kinetics of imperfect reactions in confinement

Chemical reactions generically require that particles come into contact. In practice, reaction is often imperfect and can necessitate multiple random encounters between reactants. In confined geometries, despite notable recent advances, there is to date no general analytical treatment of such imperfect transport-limited reaction kinetics. Here, we determine the kinetics of imperfect reactions in confining domains for any diffusive or anomalously diffusive Markovian transport process, and for different models of imperfect reactivity. We show that the full distribution of reaction times is obtained in the large confining volume limit from the knowledge of the mean reaction time only, which we determine explicitly. This distribution for imperfect reactions is found to be identical to that of perfect reactions upon an appropriate rescaling of parameters, which highlights the robustness of our results. Strikingly, this holds true even in the regime of low reactivity where the mean reaction time is independent of the transport process, and can lead to large fluctuations of the reaction time - even in simple reaction schemes. We illustrate our results for normal diffusion in domains of generic shape, and for anomalous diffusion in complex environments, where our predictions are confirmed by numerical simulations.

Mesoscale Simulation of Bacterial Chromosome and Cytoplasmic Nanoparticles in Confinement

In this study we investigated, using a simple polymer model of bacterial chromosome, the subdiffusive behaviors of both cytoplasmic particles and various loci in different cell wall confinements. Non-Gaussian subdiffusion of cytoplasmic particles as well as loci were obtained in our Langevin dynamic simulations, which agrees with fluorescence microscope observations. The effects of cytoplasmic particle size, locus position, confinement geometry, and density on motions of particles and loci were examined systematically. It is demonstrated that the cytoplasmic subdiffusion can largely be attributed to the mechanical properties of bacterial chromosomes rather than the viscoelasticity of cytoplasm. Due to the randomly positioned bacterial chromosome segments, the surrounding environment for both particle and loci is heterogeneous. Therefore, the exponent characterizing the subdiffusion of cytoplasmic particle/loci as well as Laplace displacement distributions of particle/loci can be reproduced by this simple model. Nevertheless, this bacterial chromosome model cannot explain the different responses of cytoplasmic particles and loci to external compression exerted on the bacterial cell wall, which suggests that the nonequilibrium activity, e.g., metabolic reactions, play an important role in cytoplasmic subdiffusion.

Minimum time controller in a class of chemical reactors based on Lagrangian approach

The main goal of this work is the construction of a class of controller, which employs directly a Lagrangian formulation to resolve the classical brachistochrone problem, this allows to obtain an optimal controller which reaches in a minimum time the stabilization of an isothermal continuous stirred tank reactor, whose chemical kinetic model is based on the power law. The proposed methodology is compared with an input/output linearizing which achieve asymptotic and exponential closed-loop convergence, sliding-mode controller with a finite time convergence and an exact gradient optimal control to compare the time convergence performance. Numerical experiments show the satisfactory performance of the proposed controller, despite sustained disturbances in the concentration input feed.

First-encounter time of two diffusing particles in confinement

We investigate how confinement may drastically change both the probability density of the first-encounter time and the associated survival probability in the case of two diffusing particles. To obtain analytical insights into this problem, we focus on two one-dimensional settings: a half-line and an interval. We first consider the case with equal particle diffusivities, for which exact results can be obtained for the survival probability and the associated first-encounter time density valid over the full time domain. We also evaluate the moments of the first-encounter time when they exist. We then turn to the case with unequal diffusivities and focus on the long-time behavior of the survival probability. Our results highlight the great impact of boundary effects in diffusion-controlled kinetics even for simple one-dimensional settings, as well as the difficulty of obtaining analytic results as soon as the translational invariance of such systems is broken.

Statistical Tests for Force Inference in Heterogeneous Environments

We devise a method to detect and estimate forces in a heterogeneous environment based on experimentally recorded stochastic trajectories. In particular, we focus on systems modeled by the heterogeneous overdamped Langevin equation. Here, the observed drift includes a "spurious” force term when the diffusivity varies in space. We show how Bayesian inference can be leveraged to reliably infer forces by taking into account such spurious forces of unknown amplitude as well as experimental sources of error. The method is based on marginalizing the force posterior over all possible spurious force contributions. The approach is combined with a Bayes factor statistical test for the presence of forces. The performance of our method is investigated analytically, numerically and tested on experimental data sets. The main results are obtained in a closed form allowing for direct exploration of their properties and fast computation. The method is incorporated into TRamWAy, an open-source software platform for automated analysis of biomolecule trajectories.

Boundary homogenization for trapping patchy particles

Many systems in chemical and biological physics involve diffusing particles being trapped by absorbing patches on otherwise reflecting surfaces. Such systems are commonly studied by boundary homogenization, in which the heterogenous boundary condition on the patchy surface is replaced by a uniform boundary condition involving a single parameter which encapsulates the effective trapping properties of the surface. In prior works on boundary homogenization, the surface is patchy and the diffusing particles are homogeneous. In this paper, we consider the opposite scenario in which a homogenous surface traps patchy particles, which could model proteins with localized binding sites, cells with membrane receptors, or patchy colloids or nanoparticles. We derive an explicit formula for the effective trapping rate which reveals a fundamental interplay between the translational and rotational diffusivities of the patchy particle, a phenomenon not typically seen in boundary homogenization. Motivated by receptors on the cell membrane, our analysis also allows for the possibility that the patches diffuse on the surface of the particle. We formulate the system in terms of a high-dimensional, time-dependent, anisotropic diffusion equation and employ matched asymptotic analysis to derive the effective trapping rate. We confirm our results by detailed numerical simulations.

Polymer Translocation Across a Corrugated Channel: Fick⁻Jacobs Approximation Extended Beyond the Mean First-Passage Time

Polymer translocation across a corrugated channel is a paradigmatic stochastic process encountered in diverse systems. The instance of time when a polymer first arrives to some prescribed location defines an important characteristic time-scale for various phenomena, which are triggered or controlled by such an event. Here we discuss the translocation dynamics of a Gaussian polymer in a periodically-corrugated channel using an appropriately generalized Fick⁻Jacobs approach. Our main aim is to probe an effective broadness of the first-passage time distribution (FPTD), by determining the so-called coefficient of variation γ of the FPTD, defined as the ratio of the standard deviation versus the mean first-passage time (MFPT). We present a systematic analysis of γ as a function of a variety of system's parameters. We show that γ never significantly drops below 1 and, in fact, can attain very large values, implying that the MFPT alone cannot characterize the first-passage statistics of the translocation process exhaustively well.

Diffusion-limited reactions in dynamic heterogeneous media

Most biochemical reactions in living cells rely on diffusive search for target molecules or regions in a heterogeneous overcrowded cytoplasmic medium. Rapid rearrangements of the medium constantly change the effective diffusivity felt locally by a diffusing particle and thus impact the distribution of the first-passage time to a reaction event. Here, we investigate the effect of these dynamic spatiotemporal heterogeneities onto diffusion-limited reactions. We describe a general mathematical framework to translate many results for ordinary homogeneous Brownian motion to heterogeneous diffusion. In particular, we derive the probability density of the first-passage time to a reaction event and show how the dynamic disorder broadens the distribution and increases the likelihood of both short and long trajectories to reactive targets. While the disorder slows down reaction kinetics on average, its dynamic character is beneficial for a faster search and realization of an individual reaction event triggered by a single molecule.

“Diffusing diffusivity” concept has been recently put forward to account for rapid structural rearrangements in soft matter and biological systems. Here the authors propose a general mathematical framework to compute the distribution of first-passage times in a dynamically heterogeneous medium.