Microtubules interacting with a boundary: mean length and mean first-passage times

Close encounters of the sticky kind: Brownian motion at absorbing boundaries

Encounter-based models of diffusion provide a probabilistic framework for analyzing the effects of a partially absorbing reactive surface, in which the probability of absorption depends upon the amount of surface-particle contact time. In this paper we develop a class of encounter-based models that deal with absorption at sticky boundaries. Sticky boundaries occur in a diverse range of applications, including cell biology, colloidal physics, finance, and human crowd dynamics. They also naturally arise in active matter, where confined active particles tend to spontaneously accumulate at boundaries even in the absence of any particle-particle interactions. We begin by constructing a one-dimensional encounter-based model of sticky Brownian motion (BM), which is based on the zero-range limit of nonsticky BM with a short-range attractive potential well near the origin. In this limit, the boundary-contact time is given by the amount of time (occupation time) that the particle spends at the origin. We calculate the joint probability density or propagator for the particle position and the occupation time, and then identify an absorption event as the first time that the occupation time crosses a randomly generated threshold. We illustrate the theory by considering diffusion in a finite interval with a partially absorbing sticky boundary at one end. We show how various quantities, such as the mean first passage time (MFPT) for single-particle absorption and the relaxation to steady state at the multiparticle level, depend on moments of the random threshold distribution. Finally, we determine how sticky BM can be obtained by taking a particular diffusion limit of a sticky run-and-tumble particle (RTP).

Critical threshold for microtubule amplification through templated severing

The cortical microtubule array of dark-grown hypocotyl cells of plant seedlings undergoes a striking, and developmentally significant, reorientation on exposure to light. This process is driven by the exponential amplification of a population of longitudinal microtubules, created by severing events localized at crossovers with the microtubules of the pre-existing transverse array. We present a dynamic one-dimensional model for microtubule amplification through this type of templated severing. We focus on the role of the probability of immediate stabilization-after-severing of the newly created lagging microtubule, observed to be a characteristic feature of the reorientation process. Employing stochastic simulations, we show that in the dynamic regime of unbounded microtubule growth, a finite value of this probability is not required for amplification to occur but does strongly influence the degree of amplification and hence the speed of the reorientation process. In contrast, in the regime of bounded microtubule growth, we show that amplification only occurs above a critical threshold. We construct an approximate analytical theory, based on a priori limiting the number of crossover events considered, which allows us to predict the observed critical value of the stabilization-after-severing probability with reasonable accuracy.

Impulsive signaling model of cytoneme-based morphogen gradient formation

Morphogen protein gradients play a vital role in regulating spatial pattern formation during development. The most commonly accepted mechanism of protein gradient formation involves the diffusion and degradation of morphogens from a localized source. However, there is growing experimental evidence for a direct cell-to-cell signaling mechanism via thin actin-rich cellular extensions known as cytonemes. Recent modeling studies of cytoneme-based morphogenesis in invertebrates ignore the discrete nature of vesicular transport along cytonemes, focusing on deterministic continuum models. In this paper, we develop an impulsive signaling model of morphogen gradient formation in invertebrates, which takes into account the discrete and stochastic nature of vesicular transport along cytonemes. We begin by solving a first passage time problem with sticky boundaries to determine the expected time to deliver a vesicle to a target cell, assuming that there is a 'nucleation' time for injecting the vesicle into the cytoneme. We then use queuing theory to analyze the impulsive model of morphogen gradient formation in the case of multiple cytonemes and multiple targets. In particular, we determine the steady-state mean and variance of the morphogen distribution across a one-dimensional array of target cells. The mean distribution recovers the spatially decaying morphogen gradient of previous deterministic models. However, the burst-like nature of morphogen transport can lead to Fano factors greater than unity across the array of cells, resulting in significant fluctuations at more distant target sites.

Search-and-capture model of cytoneme-mediated morphogen gradient formation

Morphogen protein gradients play an essential role in the spatial regulation of patterning during embryonic development. The most commonly accepted mechanism of protein gradient formation involves the diffusion and degradation of morphogens from a localized source. Recently, an alternative mechanism has been proposed, which is based on cell-to-cell transport via thin actin-rich cellular extensions known as cytonemes. Very little is currently known about the precise nature of the contacts between cytonemes and their target cells. Important unresolved issues include how cytoneme tips find their targets, how they are stabilized at their contact sites, and how vesicles are transferred to a receiving cell and subsequently internalized. It has been hypothesized that cytonemes find their targets via a random search process based on alternating periods of retraction and growth, perhaps mediated by some chemoattractant. This is an actin-based analog of the search-and-capture model of microtubules of the mitotic spindle searching for cytochrome binding sites (kinetochores) prior to separation of cytochrome pairs. In this paper we develop a search-and-capture model of cytoneme-based morphogenesis, in which nucleating cytonemes from a source cell dynamically grow and shrink along the surface of a one-dimensional array of target cells until making contact with one of the target cells. We analyze the first-passage-time problem for making contact and then use this to explore the formation of morphogen gradients under the mechanism proposed for Wnt in vertebrates. That is, we assume that morphogen is localized at the tip of a growing cytoneme, which is delivered as a "morphogen burst" to a target cell when the cytoneme makes temporary contact with a target cell before subsequently retracting. We show how multiple rounds of search-and-capture, morphogen delivery, cytoneme retraction, and nucleation events lead to the formation of a morphogen gradient. We proceed by formulating the morphogen bursting model as a queuing process, analogous to the study of translational bursting in gene networks. In order to analyze the expected times for cytoneme contact, we introduce an efficient method for solving first-passage-time problems in the presence of sticky boundaries, which exploits some classical concepts from probability theory, namely, stopping times and the strong Markov property. We end the paper by demonstrating how this method simplifies previous analyses of a well-studied problem in cell biology, namely, the search-and-capture model of microtubule-kinetochore attachment. Although the latter is completely unrelated to cytoneme-based morphogenesis from a biological perspective, it shares many of the same mathematical elements.

Search and Capture Efficiency of Dynamic Microtubules for Centrosome Relocation during IS Formation

Upon contact with antigen-presenting cells, cytotoxic T lymphocytes (T cells) establish a highly organized contact zone denoted as the immunological synapse (IS). The formation of the IS implies relocation of the microtubule organizing center (MTOC) toward the contact zone, which necessitates a proper connection between the MTOC and the IS via dynamic microtubules (MTs). The efficiency of the MTs finding the IS within the relevant timescale is, however, still illusive. We investigate how MTs search the three-dimensional constrained cellular volume for the IS and bind upon encounter to dynein anchored at the IS cortex. The search efficiency is estimated by calculating the time required for the MTs to reach the dynein-enriched region of the IS. In this study, we develop simple mathematical and numerical models incorporating relevant components of a cell and propose an optimal search strategy. Using the mathematical model, we have quantified the average search time for a wide range of model parameters and proposed an optimized set of values leading to the minimal capture time. Our results show that search times are minimal when the IS formed at the nearest or at the farthest sites on the cell surface with respect to the perinuclear MTOC. The search time increases monotonically away from these two specific sites and is maximal at an intermediate position near the equator of the cell. We observed that search time strongly depends on the number of searching MTs and distance of the MTOC from the nuclear surface.

Effects of aging in catastrophe on the steady state and dynamics of a microtubule population

Several independent observations have suggested that the catastrophe transition in microtubules is not a first-order process, as is usually assumed. Recent in vitro observations by Gardner et al. [M. K. Gardner et al., Cell 147, 1092 (2011)] showed that microtubule catastrophe takes place via multiple steps and the frequency increases with the age of the filament. Here we investigate, via numerical simulations and mathematical calculations, some of the consequences of the age dependence of catastrophe on the dynamics of microtubules as a function of the aging rate, for two different models of aging: exponential growth, but saturating asymptotically, and purely linear growth. The boundary demarcating the steady-state and non-steady-state regimes in the dynamics is derived analytically in both cases. Numerical simulations, supported by analytical calculations in the linear model, show that aging leads to nonexponential length distributions in steady state. More importantly, oscillations ensue in microtubule length and velocity. The regularity of oscillations, as characterized by the negative dip in the autocorrelation function, is reduced by increasing the frequency of rescue events. Our study shows that the age dependence of catastrophe could function as an intrinsic mechanism to generate oscillatory dynamics in a microtubule population, distinct from hitherto identified ones.

Dynamics and length distribution of microtubules under force and confinement

We investigate the microtubule polymerization dynamics with catastrophe and rescue events for three different confinement scenarios, which mimic typical cellular environments: (i) The microtubule is confined by rigid and fixed walls, (ii) it grows under constant force, and (iii) it grows against an elastic obstacle with a linearly increasing force. We use realistic catastrophe models and analyze the microtubule dynamics, the resulting microtubule length distributions, and force generation by stochastic and mean field calculations; in addition, we perform stochastic simulations. Freely growing microtubules exhibit a phase of bounded growth with finite microtubule length and a phase of unbounded growth. The main results for the three confinement scenarios are as follows: (i) In confinement by fixed rigid walls, we find exponentially decreasing or increasing stationary microtubule length distributions instead of bounded or unbounded phases, respectively. We introduce a realistic model for wall-induced catastrophes and investigate the behavior of the average length as a function of microtubule growth parameters. (ii) Under a constant force, the boundary between bounded and unbounded growth is shifted to higher tubulin concentrations and rescue rates. The critical force f(c) for the transition from unbounded to bounded growth increases logarithmically with tubulin concentration and the rescue rate, and it is smaller than the stall force. (iii) For microtubule growth against an elastic obstacle, the microtubule length and polymerization force can be regulated by microtubule growth parameters. For zero rescue rate, we find that the average polymerization force depends logarithmically on the tubulin concentration and is always smaller than the stall force in the absence of catastrophes and rescues. For a nonzero rescue rate, we find a sharply peaked steady-state length distribution, which is tightly controlled by microtubule growth parameters. The corresponding average microtubule length self-organizes such that the average polymerization force equals the critical force f(c) for the transition from unbounded to bounded growth. We also investigate the force dynamics if growth parameters are perturbed in dilution experiments. Finally, we show the robustness of our results against changes of catastrophe models and load distribution factors.