Optimizing intermittent reaction paths

Numerical simulations of confined Brownian-yet-non-Gaussian motion

Brownian motion is a central scientific paradigm. Recently, due to increasing efforts and interests towards miniaturization and small-scale physics or biology, the effects of confinement on such a motion have become a key topic of investigation. Essentially, when confined near a wall, a particle moves much slower than in the bulk due to friction at the boundaries. The mobility is therefore locally hindered and space-dependent, which in turn leads to the apparition of so-called multiplicative noises, and associated non-Gaussianities which remain difficult to resolve at all times. Here, we exploit simple, optimized and efficient numerical simulations to address Brownian motion in confinement in a broadrange and quantitative way. To do so, we integrate the overdamped Langevin equation governing the thermal dynamics of a negatively-buoyant single spherical colloid within a viscous fluid confined by two rigid walls, including surface charges. From the produced large set of long random trajectories, we perform a complete statistical analysis and extract all the key quantities, such as the probability distributions in displacements and their main moments. In particular, we propose a novel method to compute high-order cumulants by reducing convergence problems, and employ it to efficiently characterize the inherent non-Gaussianity of the confined process.

Distribution of the span of one-dimensional confined random processes before hitting a target

We derive the distribution of the number of distinct sites visited by a random walker before hitting a target site of a finite one-dimensional (1D) domain. Our approach holds for the general class of Markovian processes with connected span-i.e., whose trajectories have no "holes." We show that the distribution can be simply expressed in terms of splitting probabilities only. We provide explicit results for classical examples of random processes with relevance to target search problems, such as simple symmetric random walks, biased random walks, persistent random walks, and resetting random walks. As a by-product, explicit expressions for the splitting probabilities of all these processes are given. Extensions to reflecting boundary conditions, continuous processes, and an example of a random process with a nonconnected span are discussed.

Three-Dimensional Tracking of Interfacial Hopping Diffusion

Theoretical predictions have suggested that molecular motion at interfaces-which influences processes including heterogeneous catalysis, (bio)chemical sensing, lubrication and adhesion, and nanomaterial self-assembly-may be dominated by hypothetical "hops" through the adjacent liquid phase, where a diffusing molecule readsorbs after a given hop according to a probabilistic "sticking coefficient." Here, we use three-dimensional (3D) single-molecule tracking to explicitly visualize this process for human serum albumin at solid-liquid interfaces that exert varying electrostatic interactions on the biomacromolecule. Following desorption from the interface, a molecule experiences multiple unproductive surface encounters before readsorption. An average of approximately seven surface collisions is required for the repulsive surfaces, decreasing to approximately two and a half for surfaces that are more attractive. The hops themselves are also influenced by long-range interactions, with increased electrostatic repulsion causing hops of longer duration and distance. These findings explicitly demonstrate that interfacial diffusion is dominated by biased 3D Brownian motion involving bulk-surface coupling and that it can be controlled by influencing short- and long-range adsorbate-surface interactions.

Model of chromosomal loci dynamics in bacteria as fractional diffusion with intermittent transport

The short-time dynamics of bacterial chromosomal loci is a mixture of subdiffusive and active motion, in the form of rapid relocations with near-ballistic dynamics. While previous work has shown that such rapid motions are ubiquitous, we still have little grasp on their physical nature, and no positive model is available that describes them. Here, we propose a minimal theoretical model for loci movements as a fractional Brownian motion subject to a constant but intermittent driving force, and compare simulations and analytical calculations to data from high-resolution dynamic tracking in E. coli. This analysis yields the characteristic time scales for intermittency. Finally, we discuss the possible shortcomings of this model, and show that an increase in the effective local noise felt by the chromosome associates to the active relocations.

Interfacial Molecular Searching Using Forager Dynamics

Many biological and technological systems employ efficient non-Brownian intermittent search strategies where localized searches alternate with long flights. Coincidentally, molecular species exhibit intermittent behavior at the solid-liquid interface, where periods of slow motion are punctuated by fast flights through the liquid phase. Single-molecule tracking was used here to observe the interfacial search process of DNA for complementary DNA. Measured search times were qualitatively consistent with an intermittent-flight model, and ∼10 times faster than equivalent Brownian searches, suggesting that molecular searches for reactive sites benefit from similar efficiencies as biological organisms.

Correlated Electrochemical and Optical Detection Reveals the Chemical Reactivity of Individual Silver Nanoparticles

Electrochemical (EC) impacts of single nanoparticles (NPs) on an ultramicroelectrode are coupled with optics to identify chemical processes at the level of individual NPs. While the EC signals characterize the charge transfer process, the optical monitoring gives a complementary picture of the transport and chemical transformation of the NPs. This is illustrated in the case of electrodissolution of Ag NPs. In the simplest case, the optically monitored dissolution of individual NPs is synchronized with individual EC spikes. Optics then validates in situ the concept of EC nanoimpacts for sizing and counting of NPs. Chemical complexity is introduced by using a precipitating agent, SCN(-), which tunes the overall electrodissolution kinetics. Particularly, the charge transfer and dissolution steps occur sequentially as the synchronicity between the EC and optical signals is lost. This demonstrates the level of complexity that can be revealed from such electrochemistry/optics coupling.

Mixed random walks with a trap in scale-free networks including nearest-neighbor and next-nearest-neighbor jumps

Random walks including non-nearest-neighbor jumps appear in many real situations such as the diffusion of adatoms and have found numerous applications including PageRank search algorithm; however, related theoretical results are much less for this dynamical process. In this paper, we present a study of mixed random walks in a family of fractal scale-free networks, where both nearest-neighbor and next-nearest-neighbor jumps are included. We focus on trapping problem in the network family, which is a particular case of random walks with a perfect trap fixed at the central high-degree node. We derive analytical expressions for the average trapping time (ATT), a quantitative indicator measuring the efficiency of the trapping process, by using two different methods, the results of which are consistent with each other. Furthermore, we analytically determine all the eigenvalues and their multiplicities for the fundamental matrix characterizing the dynamical process. Our results show that although next-nearest-neighbor jumps have no effect on the leading scaling of the trapping efficiency, they can strongly affect the prefactor of ATT, providing insight into better understanding of random-walk process in complex systems.

Cell migration and antigen capture are antagonistic processes coupled by myosin II in dendritic cells

The immune response relies on the migration of leukocytes and on their ability to stop in precise anatomical locations to fulfil their task. How leukocyte migration and function are coordinated is unknown. Here we show that in immature dendritic cells, which patrol their environment by engulfing extracellular material, cell migration and antigen capture are antagonistic. This antagonism results from transient enrichment of myosin IIA at the cell front, which disrupts the back-to-front gradient of the motor protein, slowing down locomotion but promoting antigen capture. We further highlight that myosin IIA enrichment at the cell front requires the MHC class II-associated invariant chain (Ii). Thus, by controlling myosin IIA localization, Ii imposes on dendritic cells an intermittent antigen capture behaviour that might facilitate environment patrolling. We propose that the requirement for myosin II in both cell migration and specific cell functions may provide a general mechanism for their coordination in time and space.

Tuning the Flight Length of Molecules Diffusing on a Hydrophobic Surface

Transport at solid-liquid interfaces is critical to self-assembly, biosensing, and heterogeneous catalysis, but surface diffusion remains difficult to characterize and rationally manipulate, due to the inherent heterogeneity of adsorption on solid surfaces. Using single-molecule tracking, we characterized the diffusion of a fluorescent long-chain surfactant on a hydrophobic surface, which involved periods of confinement alternating with bulk-mediated "flights". The concentration of methanol in solution was varied to tune the strength of the hydrophobic surface-molecule interaction. The frequency of confinement had a nonmonotonic dependence on methanol concentration that reflected the relative influence of anomalously strong adsorption sites. By carefully accounting for the effect of this surface heterogeneity, we demonstrated that flight lengths increased monotonically as the hydrophobic attraction decreased, in agreement with theoretical predictions for bulk-mediated surface diffusion. The theory provided an accurate description of surface diffusion, despite the system being heterogeneous, and can be leveraged to optimize molecular search and assembly processes.

Invariance property of wave scattering through disordered media

A fundamental insight in the theory of diffusive random walks is that the mean length of trajectories traversing a finite open system is independent of the details of the diffusion process. Instead, the mean trajectory length depends only on the system's boundary geometry and is thus unaffected by the value of the mean free path. Here we show that this result is rooted on a much deeper level than that of a random walk, which allows us to extend the reach of this universal invariance property beyond the diffusion approximation. Specifically, we demonstrate that an equivalent invariance relation also holds for the scattering of waves in resonant structures as well as in ballistic, chaotic or in Anderson localized systems. Our work unifies a number of specific observations made in quite diverse fields of science ranging from the movement of ants to nuclear scattering theory. Potential experimental realizations using light fields in disordered media are discussed.

Accelerating search kinetics by following boundaries

We derive exact expressions of the mean first-passage time to a bulk target for a random searcher that performs boundary-mediated diffusion in a circular domain. Although nonintuitive for bulk targets, it is found that boundary excursions, if fast enough, can minimize the search time. A scaling analysis generalizes these findings to domains of arbitrary shapes and underlines their robustness. Overall, these results provide a generic mechanism of optimization of search kinetics in interfacial systems, which could have important implications in chemical physics. In the context of animal behavior sciences, it shows that following the boundaries of a domain can accelerate a search process, and therefore suggests that thigmotactism could be a kinetically efficient behavior.

Mean cover time of one-dimensional persistent random walks

The cover time is defined as the time needed for a random walker to visit every site of a confined domain. Here, we focus on persistent random walks, which provide a minimal model of random walks with short-range memory. We derive the exact expression of the mean cover time of a one-dimensional lattice by such a persistent random walk, both for periodic and reflecting boundary conditions.

Mechanisms of surface-mediated DNA hybridization

Single-molecule total internal reflection fluorescence microscopy was employed in conjunction with resonance energy transfer (RET) to observe the dynamic behavior of donor-labeled ssDNA at the interface between aqueous solution and a solid surface decorated with complementary acceptor-labeled ssDNA. At least 100,000 molecular trajectories were determined for both complementary strands and negative control ssDNA. RET was used to identify trajectory segments corresponding to the hybridized state. The vast majority of molecules from solution adsorbed nonspecifically to the surface, where a brief two-dimensional search was performed with a 7% chance of hybridization. Successful hybridization events occurred with a characteristic search time of ∼0.1 s, and unsuccessful searches resulted in desorption from the surface, ultimately repeating the adsorption and search process. Hybridization was reversible, and two distinct modes of melting (i.e., dehybridization) were observed, corresponding to long-lived (∼15 s) and short-lived (∼1.4 s) hybridized time intervals. A strand that melted back onto the surface could rehybridize after a brief search or desorb from the interface. These mechanistic observations provide guidance for technologies that involve DNA interactions in the near-surface region, suggesting a need to design surfaces that both enhance the complex multidimensional search process and stabilize the hybridized state.

Single-molecule diffusion in a periodic potential at a solid-liquid interface

We used single-molecule tracking experiments to observe the motion of small hydrophobic fluorescent molecules at the interface between water and a solid surface that exhibited periodic chemical patterns. The dynamics were characterized by non-ergodic, continuous time random walk statistics. The step-size distributions displayed enhanced probability of steps to periodic distances, consistent with theoretical predictions for diffusion in an atomic/molecular scale periodic potential. Surprisingly, this general behavior was observed here for surfaces exhibiting characteristic length scales three orders of magnitude larger than atomic/molecular dimensions, and may provide a new way to understand and control solid-liquid interfacial diffusion for molecular targeting applications.

Physical constraints determine the logic of bacterial promoter architectures

Site-specific transcription factors (TFs) bind to their target sites on the DNA, where they regulate the rate at which genes are transcribed. Bacterial TFs undergo facilitated diffusion (a combination of 3D diffusion around and 1D random walk on the DNA) when searching for their target sites. Using computer simulations of this search process, we show that the organization of the binding sites, in conjunction with TF copy number and binding site affinity, plays an important role in determining not only the steady state of promoter occupancy, but also the order at which TFs bind. These effects can be captured by facilitated diffusion-based models, but not by standard thermodynamics. We show that the spacing of binding sites encodes complex logic, which can be derived from combinations of three basic building blocks: switches, barriers and clusters, whose response alone and in higher orders of organization we characterize in detail. Effective promoter organizations are commonly found in the E. coli genome and are highly conserved between strains. This will allow studies of gene regulation at a previously unprecedented level of detail, where our framework can create testable hypothesis of promoter logic.

Splitting probabilities and interfacial territory covered by two-dimensional and three-dimensional surface-mediated diffusion

We consider the mean territory covered by a particle that performs surface-mediated diffusion inside a spherical confining domain (in two and three dimensions) before exit through an opening on the surface. This quantity can be expressed in terms of the splitting probability between two targets on the surface. We derive a general formula that relates this splitting probability to the mean first passage time to a single target that has been recently calculated for such a surface-mediated diffusion process. This formula is exact for pointlike targets and is shown to be accurate for extended targets. The mean covered territory is then found and analyzed for an arbitrary extension of the exit region in both two- and three-dimensional spherical domains.

Simulation of mRNA diffusion in the nuclear environment

A mathematical model is devised to study the diffusion of mRNA in the nucleus from the site of synthesis to a nuclear pore where it is exported to the cytoplasm. This study examines the role that nuclear structure can play in determining the kinetics of export by considering models in which elements of the nuclear skeleton and confinement by chromatin direct the mRNA movement. As a rule, a dense chromatin layer favours rapid export by reducing the effective volume for diffusion. However, it may also result in a heavy tail in the export time distribution because of the low mobility of molecules that accidentally find their way deep into the dense layer. An anisotropic solid-state transport system can also assist export. There exist both an optimal ratio of the anisotropy and an optimal depth of the solid-state transport layer that favour rapid export.

Computational models for large-scale simulations of facilitated diffusion

The binding of site-specific transcription factors to their genomic target sites is a key step in gene regulation. While the genome is huge, transcription factors belong to the least abundant protein classes in the cell. It is therefore fascinating how short the time frame is that they require to home in on their target sites. The underlying search mechanism is called facilitated diffusion and assumes a combination of three-dimensional diffusion in the space around the DNA combined with one-dimensional random walk on it. In this review, we present the current understanding of the facilitated diffusion mechanism and identify questions that lack a clear or detailed answer. One way to investigate these questions is through stochastic simulation and, in this manuscript, we support the idea that such simulations are able to address them. Finally, we review which biological parameters need to be included in such computational models in order to obtain a detailed representation of the actual process.

Geometry-induced bursting dynamics in gene expression

In prokaryotes and eukaryotes, genes are transcribed stochastically according to various temporal patterns that range from simple first-order kinetics to marked bursts, resulting in temporal and cell-to-cell variations of mRNA and protein levels. Here, we consider the effect of the transport of regulatory molecules on the noise in gene expression by taking into account explicitly the dynamics of a finite number of transcription factors confined in the cell. We calculate analytically time-dependent correlation functions of mRNA levels for a wide range of transport mechanisms and find that in the limit of small-transcription-factor copy number, the results differ significantly from standard approaches, which ignore confinement. It is shown how such dynamical quantities, which can now be obtained experimentally, can be used to identify the underlying mechanisms of transcription. Of particular importance, it is demonstrated that the geometry of transcription-factor trajectories in the cellular environment plays a key role in transcription kinetics, and can intrinsically generate the observed various transcription patterns ranging from simple first-order kinetics to bursts.

Interfacial territory covered by surface-mediated diffusion

We consider a minimal model of heterogeneous catalysis in which a molecule performs surface-mediated diffusion inside a confining domain whose boundary contains catalytic sites. We explicitly take into account the combination of surface and bulk diffusion, and we obtain exact results for the mean and variance of the territory covered on the boundary by the particle before its exit in the case of a two-dimensional spherical domain. Depending on the relative positions of the entrance and exit points, very different behaviors with respect to the mean adsorption time of the molecule on the surface are found. We also determine both exact lower and upper bounds and an approximate expression of the probability of reacting with catalytic sites before exiting the domain. These results provide a quantitative measure of the efficiency of an idealized catalyst.

Optimizing persistent random searches

We consider a minimal model of persistent random searcher with a short range memory. We calculate exactly for such a searcher the mean first-passage time to a target in a bounded domain and find that it admits a nontrivial minimum as function of the persistence length. This reveals an optimal search strategy which differs markedly from the simple ballistic motion obtained in the case of Poisson distributed targets. Our results show that the distribution of targets plays a crucial role in the random search problem. In particular, in the biologically relevant cases of either a single target or regular patterns of targets, we find that, in strong contrast to repeated statements in the literature, persistent random walks with exponential distribution of excursion lengths can minimize the search time, and in that sense perform better than any Levy walk.

Exact calculations of first-passage quantities on recursive networks

We present general methods to exactly calculate mean first-passage quantities on self-similar networks defined recursively. In particular, we calculate the mean first-passage time and the splitting probabilities associated to a source and one or several targets; averaged quantities over a given set of sources (e.g., same-connectivity nodes) are also derived. The exact estimate of such quantities highlights the dependency of first-passage processes with respect to the source-target distance, which has recently revealed to be a key parameter in characterizing transport in complex media. We explicitly perform calculations for different classes of recursive networks [finitely ramified fractals, scale-free (trans)fractals, nonfractals, mixtures between fractals and nonfractals, nondecimable hierarchical graphs] of arbitrary size. Our approach unifies and significantly extends the available results in the field.

Classes of fast and specific search mechanisms for proteins on DNA

Problems of search and recognition appear over different scales in biological systems. In this review we focus on the challenges posed by interactions between proteins, in particular transcription factors, and DNA and possible mechanisms which allow for fast and selective target location. Initially we argue that DNA-binding proteins can be classified, broadly, into three distinct classes which we illustrate using experimental data. Each class calls for a different search process and we discuss the possible application of different search mechanisms proposed over the years to each class. The main thrust of this review is a new mechanism which is based on barrier discrimination. We introduce the model and analyze in detail its consequences. It is shown that this mechanism applies to all classes of transcription factors and can lead to a fast and specific search. Moreover, it is shown that the mechanism has interesting transient features which allow for stability at the target despite rapid binding and unbinding of the transcription factor from the target.

Universality classes of first-passage-time distribution in confined media

We study the first-passage time (FPT) distribution to a target site for a random walker evolving in a bounded domain. We show that in the limit of large volume of the confining domain, this distribution falls into universality classes indexed by the walk dimension d(w) and the fractal dimension d(f) of the medium, which have been recently identified previously [Bénichou et al., Nat. Chem. 2, 472 (2010)]. We present in this paper a complete derivation of these universal distributions, discuss extensively the range of applicability of the results, and extend the method to continuous-time random walks. This analysis puts forward the importance of the geometry, and in particular the position of the starting point, in first-passage statistics. Analytical results are validated by numerical simulations, applied to various models of transport in disordered media, which illustrate the universality classes of the FPT distribution.

Facilitated diffusion of proteins on chromatin

We present a theoretical model of facilitated diffusion of proteins in the cell nucleus. This model, which takes into account the successive binding and unbinding events of proteins to DNA, relies on a fractal description of the chromatin which has been recently evidenced experimentally. Facilitated diffusion is shown quantitatively to be favorable for a fast localization of a target locus by a transcription factor and even to enable the minimization of the search time by tuning the affinity of the transcription factor with DNA. This study shows the robustness of the facilitated diffusion mechanism, invoked so far only for linear conformations of DNA.

Response to targeted perturbations for random walks on networks

We introduce a general framework, applicable to a broad class of random walks on networks, that quantifies the response of the mean first-passage time to a target node to a local perturbation of the network, both in the context of attacks (damaged link) or strategies of transport enhancement (added link). This approach enables to determine explicitly the dependence of this response on geometric parameters (such as the network size and the localization of the perturbation) and on the intensity of the perturbation. In particular, it is showed that the relative variation of the mean first-passage time is independent of the network size, and remains significant in the large size limit. Furthermore, in the case of noncompact exploration of the network, it is found that a targeted perturbation keeps a substantial impact on transport properties for any localization of the damaged link.

Optimal reaction time for surface-mediated diffusion

We present an exact calculation of the mean first-passage time to a small target on the surface of a 2D or 3D spherical domain, for a molecule performing surface-mediated diffusion. This minimal model of interfacial reactions, which explicitly takes into account the combination of surface and bulk diffusion, shows the importance of correlations induced by the coupling of the switching dynamics to the geometry of the confinement, ignored so far. Our results show that, in the context of interfacial systems in confinement, the reaction time can be minimized as a function of the desorption rate from the surface, which puts forward a general mechanism of enhancement and regulation of chemical and biological reactivity.

Geometry-controlled kinetics

It has long been appreciated that the transport properties of molecules can control reaction kinetics. This effect can be characterized by the time it takes a diffusing molecule to reach a target-the first-passage time (FPT). Determining the FPT distribution in realistic confined geometries has until now, however, seemed intractable. Here, we calculate this FPT distribution analytically and show that transport processes as varied as regular diffusion, anomalous diffusion, and diffusion in disordered media and fractals, fall into the same universality classes. Beyond the theoretical aspect, this result changes our views on standard reaction kinetics and we introduce the concept of 'geometry-controlled kinetics'. More precisely, we argue that geometry-and in particular the initial distance between reactants in 'compact' systems-can become a key parameter. These findings could help explain the crucial role that the spatial organization of genes has in transcription kinetics, and more generally the impact of geometry on diffusion-limited reactions.

Optimization of the residence time of a Brownian particle in a spherical subdomain

In this communication, we show that the residence time of a Brownian particle, defined as the cumulative time spent in a given region of space, can be optimized as a function of the diffusion coefficient. We discuss the relevance of this effect to several schematic experimental situations classified in nature--random or deterministic--both of the observation time and of the starting position of the Brownian particle.

The first-passage problem for diffusion through a cylindrical pore with sticky walls

We calculate the first-passage time distribution for diffusion through a cylindrical pore with sticky walls. A particle diffusively explores the interior of the pore through a series of binding and unbinding events with the cylinder wall. Through a diagrammatic expansion we obtain first-passage time statistics for the particle's exit from the pore. Connections between the model and nucleocytoplasmic transport in cells are discussed.

Searching fast for a target on DNA without falling to traps

Genomic expression depends critically on both the ability of regulatory proteins to locate specific target sites on DNA within seconds and on the formation of long-lived (many minutes) complexes between these proteins and the DNA. Equilibrium experiments show that indeed regulatory proteins bind tightly to their target site. However, they also find strong binding to other nonspecific sites which act as traps that can dramatically increase the time needed to locate the target. This gives rise to a conflict between the speed and stability requirements. Here we suggest a simple mechanism which can resolve this long-standing paradox.

Robustness of optimal intermittent search strategies in one, two, and three dimensions

Search problems at various scales involve a searcher, be it a molecule before reaction or a foraging animal, which performs an intermittent motion. Here we analyze a generic model based on such type of intermittent motion, in which the searcher alternates phases of slow motion allowing detection and phases of fast motion without detection. We present full and systematic results for different modeling hypotheses of the detection mechanism in space in one, two, and three dimensions. Our study completes and extends the results of our recent letter [Loverdo, Nat. Phys. 4, 134 (2008)] and gives the necessary calculation details. In addition, another modeling of the detection case is presented. We show that the mean target detection time can be minimized as a function of the mean duration of each phase in one, two, and three dimensions. Importantly, this optimal strategy does not depend on the details of the modeling of the slow detection phase, which shows the robustness of our results. We believe that this systematic analysis can be used as a basis to study quantitatively various real search problems involving intermittent behaviors.

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.