Queueing induced by bidirectional motor motion near the end of a microtubule

Modelling the motion of organelles in an elongated cell via the coordination of heterogeneous drift-diffusion and long-range transport

 Cellular distribution of organelles in living cells is achieved via a variety of transport mechanisms, including directed motion, mediated by molecular motors along microtubules (MTs), and diffusion which is predominantly heterogeneous in space. In this paper, we introduce a model for particle transport in elongated cells that couples poleward drift, long-range bidirectional transport and diffusion with spatial heterogeneity in a three-dimensional space. Using stochastic simulations and analysis of a related population model, we find parameter regions where the three-dimensional model can be reduced to a coupled one-dimensional model or even a one-dimensional scalar model. We explore the efficiency with which individual model components can overcome drift towards one of the cell poles to reach an approximately even distribution. In particular, we find that if lateral movement is well mixed, then increasing the binding ability of particles to MTs is an efficient way to overcome a poleward drift, whereas if lateral motion is not well mixed, then increasing the axial diffusivity away from MTs becomes an efficient way to overcome the poleward drift. Our three-dimensional model provides a new tool that will help to understand the mechanisms by which eukaryotic cells organize their organelles in an elongated cell, and in particular when the one-dimensional models are applicable.

Getting around the cell: physical transport in the intracellular world

Eukaryotic cells face the challenging task of transporting a variety of particles through the complex intracellular milieu in order to deliver, distribute, and mix the many components that support cell function. In this review, we explore the biological objectives and physical mechanisms of intracellular transport. Our focus is on cytoplasmic and intra-organelle transport at the whole-cell scale. We outline several key biological functions that depend on physically transporting components across the cell, including the delivery of secreted proteins, support of cell growth and repair, propagation of intracellular signals, establishment of organelle contacts, and spatial organization of metabolic gradients. We then review the three primary physical modes of transport in eukaryotic cells: diffusive motion, motor-driven transport, and advection by cytoplasmic flow. For each mechanism, we identify the main factors that determine speed and directionality. We also highlight the efficiency of each transport mode in fulfilling various key objectives of transport, such as particle mixing, directed delivery, and rapid target search. Taken together, the interplay of diffusion, molecular motors, and flows supports the intracellular transport needs that underlie a broad variety of biological phenomena.

A mathematical understanding of how cytoplasmic dynein walks on microtubules

Cytoplasmic dynein 1 (hereafter referred to simply as dynein) is a dimeric motor protein that walks and transports intracellular cargos towards the minus end of microtubules. In this article, we formulate, based on physical principles, a mechanical model to describe the stepping behaviour of cytoplasmic dynein walking on microtubules from the cell membrane towards the nucleus. Unlike previous studies on physical models of this nature, we base our formulation on the whole structure of dynein to include the temporal dynamics of the individual subunits such as the cargo (for example, an endosome, vesicle or bead), two rings of six ATPase domains associated with diverse cellular activities (AAA+ rings) and the microtubule-binding domains which allow dynein to bind to microtubules. This mathematical framework allows us to examine experimental observations on dynein across a wide range of different species, as well as being able to make predictions on the temporal behaviour of the individual components of dynein not currently experimentally measured. Furthermore, we extend the model framework to include backward stepping, variable step size and dwelling. The power of our model is in its predictive nature; first it reflects recent experimental observations that dynein walks on microtubules using a weakly coordinated stepping pattern with predominantly not passing steps. Second, the model predicts that interhead coordination in the ATP cycle of cytoplasmic dynein is important in order to obtain the alternating stepping patterns and long run lengths seen in experiments.

Motor Protein Accumulation on Antiparallel Microtubule Overlaps

Biopolymers serve as one-dimensional tracks on which motor proteins move to perform their biological roles. Motor protein phenomena have inspired theoretical models of one-dimensional transport, crowding, and jamming. Experiments studying the motion of Xklp1 motors on reconstituted antiparallel microtubule overlaps demonstrated that motors recruited to the overlap walk toward the plus end of individual microtubules and frequently switch between filaments. We study a model of this system that couples the totally asymmetric simple exclusion process for motor motion with switches between antiparallel filaments and binding kinetics. We determine steady-state motor density profiles for fixed-length overlaps using exact and approximate solutions of the continuum differential equations and compare to kinetic Monte Carlo simulations. Overlap motor density profiles and motor trajectories resemble experimental measurements. The phase diagram of the model is similar to the single-filament case for low switching rate, while for high switching rate we find a new (to our knowledge) low density-high density-low density-high density phase. The overlap center region, far from the overlap ends, has a constant motor density as one would naïvely expect. However, rather than following a simple binding equilibrium, the center motor density depends on total overlap length, motor speed, and motor switching rate. The size of the crowded boundary layer near the overlap ends is also dependent on the overlap length and switching rate in addition to the motor speed and bulk concentration. The antiparallel microtubule overlap geometry may offer a previously unrecognized mechanism for biological regulation of protein concentration and consequent activity.

Spatial organization of organelles in fungi: Insights from mathematical modelling

•

Modelling of dynein motility reveals a stochastic role in dynein comet formation.

•

Modelling helps to elucidate mechanisms in spatial organization of early endosomes.

•

A combination of diffusion and directed motion distributes ribosomes and peroxisomes.

Mathematical modelling in cellular systems aims to describe biological processes in a quantitative manner. Most accurate modelling is based on robust experimental data. Here we review recent progress in the theoretical description of motor behaviour, early endosome motility, ribosome distribution and peroxisome transport in the fungal model system Ustilago maydis and illustrate the power of modelling in our quest to understand molecular details and cellular roles of membrane trafficking in filamentous fungi.

Phase-plane analysis of the totally asymmetric simple exclusion process with binding kinetics and switching between antiparallel lanes

Motor protein motion on biopolymers can be described by models related to the totally asymmetric simple exclusion process (TASEP). Inspired by experiments on the motion of kinesin-4 motors on antiparallel microtubule overlaps, we analyze a model incorporating the TASEP on two antiparallel lanes with binding kinetics and lane switching. We determine the steady-state motor density profiles using phase-plane analysis of the steady-state mean field equations and kinetic Monte Carlo simulations. We focus on the density-density phase plane, where we find an analytic solution to the mean field model. By studying the phase-space flows, we determine the model's fixed points and their changes with parameters. Phases previously identified for the single-lane model occur for low switching rate between lanes. We predict a multiple coexistence phase due to additional fixed points that appear as the switching rate increases: switching moves motors from the higher-density to the lower-density lane, causing local jamming and creating multiple domain walls. We determine the phase diagram of the model for both symmetric and general boundary conditions.

Application of quasi-steady state methods to molecular motor transport on microtubules in fungal hyphae

We consider bidirectional transport of cargo by molecular motors dynein and kinesin that walk along microtubules, and/or diffuse in the cell. The motors compete to transport cargo in opposite directions with respect to microtubule polarity (towards the plus or minus end of the microtubule). In recent work, Gou et al. (2014) used a hierarchical set of models, each consisting of continuum transport equations to track the evolution of motors and their cargo (early endosomes) in the specific case of the fungus Ustilago maydis. We complement their work using a framework of quasi-steady state analysis developed by Newby and Bressloff (2010) and Bressloff and Newby (2013) to reduce the models to an approximating steady state Fokker-Plank equation. This analysis allows us to find analytic approximations to the steady state solutions in many cases where the full models are not easily solved. Consequently, we can make predictions about parameter dependence of the resulting spatial distributions. We also characterize the overall rates of bulk transport and diffusion, and how these are related to state transition parameters, motor speeds, microtubule polarity distribution, and specific assumptions made.

Mathematical model with spatially uniform regulation explains long-range bidirectional transport of early endosomes in fungal hyphae

Cellular cargo transported bidirectionally along microtubules by dynein and kinesin can be organized by spatially nonuniform upstream regulation or can self-organize. A mathematical model of early endosome transport in fungal hyphae demonstrates that spatiotemporally uniform regulation results in cargo dynamics consistent with experiment.

In many cellular contexts, cargo is transported bidirectionally along microtubule bundles by dynein and kinesin-family motors. Upstream factors influence how individual cargoes are locally regulated, as well as how long-range transport is regulated at the whole-cell scale. Although the details of local, single-cargo bidirectional switching have been extensively studied, it remains to be elucidated how this results in cell-scale spatial organization. Here we develop a mathematical model of early endosome transport in Ustilago maydis. We demonstrate that spatiotemporally uniform regulation, with constant transition rates, results in cargo dynamics that is consistent with experimental data, including data from motor mutants. We find that microtubule arrays can be symmetric in plus-end distribution but asymmetric in binding-site distribution in a manner that affects cargo dynamics and that cargo can travel past microtubule ends in microtubule bundles. Our model makes several testable predictions, including secondary features of dynein and cargo distributions.

Occupation probabilities and fluctuations in the asymmetric simple inclusion process

The asymmetric simple inclusion process (ASIP), a lattice-gas model of unidirectional transport and aggregation, was recently proposed as an "inclusion" counterpart of the asymmetric simple exclusion process. In this paper we present an exact closed-form expression for the probability that a given number of particles occupies a given set of consecutive lattice sites. Our results are expressed in terms of the entries of Catalan's trapezoids-number arrays which generalize Catalan's numbers and Catalan's triangle. We further prove that the ASIP is asymptotically governed by the following: (i) an inverse square-root law of occupation, (ii) a square-root law of fluctuation, and (iii) a Rayleigh law for the distribution of interexit times. The universality of these results is discussed.

Motor-mediated bidirectional transport along an antipolar microtubule bundle: a mathematical model

Long-distance bidirectional transport of organelles depends on the coordinated motion of various motor proteins on the cytoskeleton. Recent quantitative live cell imaging in the elongated hyphal cells of Ustilago maydis has demonstrated that long-range motility of motors and their endosomal cargo occurs on unipolar microtubules (MTs) near the extremities of the cell. These MTs are bundled into antipolar bundles within the central part of the cell. Dynein and kinesin-3 motors coordinate their activity to move early endosomes (EEs) in a bidirectional fashion where dynein drives motility towards MT minus ends and kinesin towards MT plus ends. Although this means that one can easily assign the drivers of bidirectional motion in the unipolar section, the bipolar orientations in the bundle mean that it is possible for either motor to drive motion in either direction. In this paper we use a multilane asymmetric simple exclusion process modeling approach to simulate and investigate phases of bidirectional motility in a minimal model of an antipolar MT bundle. In our model, EE cargos (particles) change direction on each MT with a turning rate Ω and there is switching between MTs in the bundle at the minus ends. At these ends, particles can hop between MTs with rate q(1) on passing from a unipolar to a bipolar section (the obstacle-induced switching rate) or q(2) on passing in the other direction (the end-induced switching rate). By a combination of numerical simulations and mean-field approximations, we investigate the distribution of particles along the MTs for different values of these parameters and of Θ, the overall density of particles within this closed system. We find that even if Θ is low, the system can exhibit a variety of phases with shocks in the density profiles near plus and minus ends caused by queuing of particles. We discuss how the parameters influence the type of particle that dominates active transport in the bundle.

Limit laws for the asymmetric inclusion process

The Asymmetric Inclusion Process (ASIP) is a unidirectional lattice-gas flow model which was recently introduced as an exactly solvable 'Bosonic' counterpart of the 'Fermionic' asymmetric exclusion process. An iterative algorithm that allows the computation of the probability generating function (PGF) of the ASIP's steady state exists but practical considerations limit its applicability to small ASIP lattices. Large lattices, on the other hand, have been studied primarily via Monte Carlo simulations and were shown to display a wide spectrum of intriguing statistical phenomena. In this paper we bypass the need for direct computation of the PGF and explore the ASIP's asymptotic statistical behavior. We consider three different limiting regimes: heavy-traffic regime, large-system regime, and balanced-system regime. In each of these regimes we obtain-analytically and in closed form-stochastic limit laws for five key ASIP observables: traversal time, overall load, busy period, first occupied site, and draining time. The results obtained yield a detailed limit-laws perspective of the ASIP, numerical simulations demonstrate the applicability of these laws as useful approximations.

From the cell membrane to the nucleus: unearthing transport mechanisms for dynein

Mutations in the motor protein cytoplasmic dynein have been found to cause Charcot-Marie-Tooth disease, spinal muscular atrophy, and severe intellectual disabilities in humans. In mouse models, neurodegeneration is observed. We sought to develop a novel model which could incorporate the effects of mutations on distance travelled and velocity. A mechanical model for the dynein mediated transport of endosomes is derived from first principles and solved numerically. The effects of variations in model parameter values are analysed to find those that have a significant impact on velocity and distance travelled. The model successfully describes the processivity of dynein and matches qualitatively the velocity profiles observed in experiments.

Bidirectional transport in a multispecies totally asymmetric exclusion-process model

We study a minimal lattice model which describes bidirectional transport of "particles" driven along a one-dimensional track, as is observed in microtubule based, motor protein driven bidirectional transport of cargo vesicles, lipid bodies, and organelles such as mitochondria. This minimal model, a multispecies totally asymmetric exclusion process (TASEP) with directional switching, can provide a framework for understanding the interplay between the switching dynamics of individual particles and the collective movement of particles in one dimension. When switching is much faster than translocation, the steady-state density and current profiles of the particles are homogeneous in the bulk and are well described by mean-field (MF) theory, as determined by comparison to a Monte Carlo simulation. In this limit, we can map this model to the exactly solvable partially asymmetric exclusion-process (PASEP) model. Away from this fast switching regime the MF theory fails, although the average bulk density profile still remains homogeneous. We study the steady-state behavior as a function of the ratio of the translocation and net switching rates Q and find a unique first-order phase transition at a finite Q associated with a discontinuous change of the bulk density. When the switching rate is decreased further (keeping translocation rate fixed), the system approaches a jammed phase with a net current that tends to zero as J~1/Q. We numerically construct the phase diagram for finite Q.

Fungal development of the plant pathogen Ustilago maydis

The maize pathogen Ustilago maydis has to undergo various morphological transitions for the completion of its sexual life cycle. For example, haploid cells respond to pheromone by forming conjugation tubes that fuse at their tips. The resulting dikaryon grows filamentously, expanding rapidly at the apex and inserting retraction septa at the basal pole. In this review, we present progress on the underlying mechanisms regulating such defined developmental programmes. The key findings of the postgenomic era are as follows: (1) endosomes function not only during receptor recycling, but also as multifunctional transport platforms; (2) a new transcriptional master regulator for pathogenicity is part of an intricate transcriptional network; (3) determinants for uniparental mitochondrial inheritance are encoded at the a2 mating-type locus; (4) microtubule-dependent mRNA transport is important in determining the axis of polarity; and (5) a battery of fungal effectors encoded in gene clusters is crucial for plant infection. Importantly, most processes are tightly controlled at the transcriptional, post-transcriptional and post-translational levels, resulting in a complex regulatory network. This intricate system is crucial for the timing of the correct order of developmental phases. Thus, new insights from all layers of regulation have substantially advanced our understanding of fungal development.

Controlled and stochastic retention concentrates dynein at microtubule ends to keep endosomes on track

Bidirectional transport of early endosomes (EEs) involves microtubules (MTs) and associated motors. In fungi, the dynein/dynactin motor complex concentrates in a comet-like accumulation at MT plus-ends to receive kinesin-3-delivered EEs for retrograde transport. Here, we analyse the loading of endosomes onto dynein by combining live imaging of photoactivated endosomes and fluorescent dynein with mathematical modelling. Using nuclear pores as an internal calibration standard, we show that the dynein comet consists of ∼55 dynein motors. About half of the motors are slowly turned over (T(1/2): ∼98 s) and they are kept at the plus-ends by an active retention mechanism involving an interaction between dynactin and EB1. The other half is more dynamic (T(1/2): ∼10 s) and mathematical modelling suggests that they concentrate at MT ends because of stochastic motor behaviour. When the active retention is impaired by inhibitory peptides, dynein numbers in the comet are reduced to half and ∼10% of the EEs fall off the MT plus-ends. Thus, a combination of stochastic accumulation and active retention forms the dynein comet to ensure capturing of arriving organelles by retrograde motors.