Weakly coupled, antiparallel, totally asymmetric simple exclusion processes

Exclusion in junction geometries

We investigate the dynamics of the asymmetric exclusion process at a junction. When two input roads are initially fully occupied and a single output road is initially empty, the ensuing rarefaction wave has a rich spatial structure. The density profile also changes dramatically as the initial densities are varied. Related phenomenology arises when one road feeds into two. Finally, we determine the phase diagram of the open system, where particles are fed into two roads at rate α for each road, the two roads merge into one, and particles are extracted from the single output road at rate β.

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.

Driven transport on open filaments with interfilament switching processes

We study a two-filament driven lattice gas model with oppositely directed species of particles moving on two parallel filaments with filament-switching processes and particle inflow and outflow at filament ends. The filament-switching process is correlated with the occupation number of the adjacent site such that particles switch filaments with finite probability only when oppositely directed particles meet on the same filament. This model mimics some of the coarse-grained features observed in context of microtubule-(MT) based intracellular transport, wherein cellular cargo loaded and off-loaded at filament ends are transported on multiple parallel MT filaments and can switch between the parallel microtubule filaments. We focus on a regime where the filaments are weakly coupled, such that filament-switching rate of particles scale inversely as the length of the filament. We find that the interplay of (off-) loading processes at the boundaries and the filament-switching process of particles leads to some distinctive features of the system. These features includes occurrence of a variety of phases in the system with inhomogeneous density profiles including localized density shocks, density difference across the filaments, and bidirectional current flows in the system. We analyze the system by developing a mean field (MF) theory and comparing the results obtained from the MF theory with the Monte Carlo (MC) simulations of the dynamics of the system. We find that the steady-state density and current profiles of particles and the phase diagram obtained within the MF picture matches quite well with MC simulation results. These findings maybe useful for studying multifilament intracellular transport.

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.

Theoretical Analysis of Dynamic Processes for Interacting Molecular Motors

Biological transport is supported by collective dynamics of enzymatic molecules that are called motor proteins or molecular motors. Experiments suggest that motor proteins interact locally via short-range potentials. We investigate the fundamental role of these interactions by analyzing a new class of totally asymmetric exclusion processes where interactions are accounted for in a thermodynamically consistent fashion. It allows us to connect explicitly microscopic features of motor proteins with their collective dynamic properties. Theoretical analysis that combines various mean-field calculations and computer simulations suggests that dynamic properties of molecular motors strongly depend on interactions, and correlations are stronger for interacting motor proteins. Surprisingly, it is found that there is an optimal strength of interactions (weak repulsion) that leads to a maximal particle flux. It is also argued that molecular motors transport is more sensitive to attractive interactions. Applications of these results for kinesin motor proteins are discussed.

Bulk induced phase transition in driven diffusive systems

This paper studies a weakly and asymmetrically coupled three-lane driven diffusive system. A non-monotonically changing density profile in the middle lane has been observed. When the extreme value of the density profile reaches ρ = 0.5, a bulk induced phase transition occurs which exhibits a shock and a continuously and smoothly decreasing density profile which crosses ρ = 0.5 upstream or downstream of the shock. The existence of double shocks has also been observed. A mean-field approach has been used to interpret the numerical results obtained by Monte Carlo simulations. The current minimization principle has excluded the occurrence of two or more bulk induced shocks in the general case of nonzero lane changing rates.

Phase diagrams of the Katz-Lebowitz-Spohn process on lattices with a junction

This paper studies the Katz-Lebowitz-Spohn (KLS) process on lattices with a junction, where particles move on parallel lattice branches that combine into a single lattice at the junction. It is shown that 11 kinds of phase diagrams could be observed, depending on the two parameters ε and δ in the KLS process. We have investigated the phase diagrams as well as bulk density analytically based on flow rate conservation and the extremal current principle. Extensive Monte Carlo computer simulations are performed, and it is found that they are in excellent agreement with theoretical prediction.

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.

Multispecies model with interconversion, chipping, and injection

Motivated by the phenomenology of transport through the Golgi apparatus of cells, we study a multispecies model with boundary injection of one species of particle, interconversion between the different species of particle, and driven diffusive movement of particles through the system by chipping of a single particle from a site. The model is analyzed in one dimension using equations for particle currents. It is found that, depending on the rates of various processes and the asymmetry in the hopping, the system may exist either in a steady phase, in which the average mass at each site attains a time-independent value, or in a "growing" phase, in which the total mass grows indefinitely in time, even in a finite system. The growing phases have interesting spatial structure. In particular, we find phases in which some spatial regions of the system have a constant average mass, while other regions show unbounded growth.

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.

Driven transport on parallel lanes with particle exclusion and obstruction

We investigate a driven two-channel system where particles on different lanes mutually obstruct each other's motion, extending an earlier model by Popkov and Peschel [Phys. Rev. E 64, 026126 (2001)]. This obstruction may occur in biological contexts due to steric hinderance where motor proteins carry cargos by "walking" on microtubules. Similarly, the model serves as a description for classical spin transport where charged particles with internal states move unidirectionally on a lattice. Three regimes of qualitatively different behavior are identified, depending on the strength of coupling between the lanes. For small and large coupling strengths the model can be mapped to a one-channel problem, whereas a rich phase behavior emerges for intermediate ones. We derive an approximate but quantitatively accurate theoretical description in terms of a one-site cluster approximation, and obtain insight into the phase behavior through the current-density relations combined with an extremal-current principle. Our results are confirmed by stochastic simulations.

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

Recent live observations of motors in long-range microtubule (MT) dependent transport in the fungus Ustilago maydis have reported bidirectional motion of dynein and an accumulation of the motors at the polymerization-active (the plus-end) of the microtubule. Quantitative data derived from in vivo observation of dynein has enabled us to develop an accurate, quantitatively-valid asymmetric simple exclusion process (ASEP) model that describes the coordinated motion of anterograde and retrograde motors sharing a single oriented microtubule. We give approximate expressions for the size and distribution of the accumulation, and discuss queueing properties for motors entering this accumulation. We show for this ASEP model, that the mean accumulation can be modeled as an M/M/∞ queue that is Poisson distributed with mean F(arr)/p(d), where F(arr) is the flux of motors that arrives at the tip and p(d) is the rate at which individual motors change direction from anterograde to retrograde motion. Deviations from this can in principle be used to gain information about other processes at work in the accumulation. Furthermore, our work is a significant step toward a mathematical description of the complex interactions of motors in cellular long-range transport of organelles.

Domain wall delocalization, dynamics and fluctuations in an exclusion process with two internal states

We investigate the delocalization transition appearing in an exclusion process with two internal states, respectively on two parallel lanes. At the transition, delocalized domain walls form in the density profiles of both internal states, in agreement with a mean-field approach. Remarkably, the topology of the system's phase diagram allows for the delocalization of a (localized) domain wall when approaching the transition. We quantify the domain wall's delocalization close to the transition by analytic results obtained within the framework of the domain wall picture. Power law dependences of the domain wall width on the distance to the delocalization transition as well as on the system size are uncovered, they agree with numerical results.

Weak and strong coupling in a two-lane asymmetric exclusion process

This paper studies a two-lane totally asymmetric simple exclusion process, in which particles could jump between the two lanes with asymmetric rates. In the weak coupling situation, the rates are inversely proportional to system size L . The appearance of localized shock in one lane and the discontinuous phase transition as revealed by Juhász [Phys. Rev. E 76, 021117 (2007)] are also reproduced. The density profiles and phase diagrams are constructed in the hydrodynamic limit, by numerically solving the steady state equations. The phase diagram in our model exhibits asymmetry, which is different from the symmetric one in Juhász's model. We have studied the phase boundary and discontinuous line analytically. The analytical results are in good agreement with that obtained from numerical integration. We also study the strong coupling situation, in which the lane changing rates are independent of L . Results completely different from that arising from weak coupling are presented. Furthermore, features different from that of the model presented by Pronina and Kolomeisky [Physica A 372, 12 (2006)] are revealed.