Molecular motor traffic in a half-open tube

The interplay of active and passive mechanisms in slow axonal transport

A combination of intermittent active movement of transient aggregates and a paused state that intervenes between periods of active transport has been proposed to underlie the slow, directed transport of soluble proteins in axons. A component of passive diffusion in the axoplasm may also contribute to slow axonal transport, although quantitative estimates of the relative contributions of diffusive and active movement in the slow transport of a soluble protein, and in particular how they might vary across developmental stages, are lacking. Here, we propose and study a model for slow axonal transport, addressing data from bleach recovery measurements on a small, soluble, protein, choline acetyltransferase, in thin axons of the lateral chordotonal (lch5) sensory neurons of Drosophila. Choline acetyltransferase is mainly present in soluble form in the axon and catalyzes the acetylation of choline at the synapse. It does not form particulate structures in axons and moves at rates characteristic of slow component b (≈ 1-10 mm/day or 0.01-0.1 μm/s). Using our model, which incorporates active transport with paused and/or diffusive states, we predict bleach recovery, transport rates, and cargo trajectories obtained through kymographs, comparing these with experimental observations at different developmental stages. We show that changes in the diffusive fraction of cargo during these developmental stages dominate bleach recovery and that a combination of active motion with a paused state alone cannot reproduce the data. We compared predictions of the model with results from photoactivation experiments. The importance of the diffusive state in reproducing the bleach recovery signal in the slow axonal transport of small soluble proteins is our central result.

Modelling the effect of ribosome mobility on the rate of protein synthesis

Translation is one of the main steps in the synthesis of proteins. It consists of ribosomes that translate sequences of nucleotides encoded on mRNA into polypeptide sequences of amino acids. Ribosomes bound to mRNA move unidirectionally, while unbound ribosomes diffuse in the cytoplasm. It has been hypothesized that finite diffusion of ribosomes plays an important role in ribosome recycling and that mRNA circularization enhances the efficiency of translation, see e.g. Lodish et al. (Molecular cell biology, 8th edn, W.H. Freeman and Company, San Francisco, 2016). In order to estimate the effect of cytoplasmic diffusion on the rate of translation, we consider a totally asymmetric simple exclusion process coupled to a finite diffusive reservoir, which we call the ribosome transport model with diffusion. In this model, we derive an analytical expression for the rate of protein synthesis as a function of the diffusion constant of ribosomes, which is corroborated with results from continuous-time Monte Carlo simulations. Using a wide range of biological relevant parameters, we conclude that diffusion is not a rate limiting factor in translation initiation because diffusion is fast enough in biological cells.

Motor Protein Transport Along Inhomogeneous Microtubules

Many cellular processes rely on the cell's ability to transport material to and from the nucleus. Networks consisting of many microtubules and actin filaments are key to this transport. Recently, the inhibition of intracellular transport has been implicated in neurodegenerative diseases such as Alzheimer's disease and Amyotrophic Lateral Sclerosis. Furthermore, microtubules may contain so-called defective regions where motor protein velocity is reduced due to accumulation of other motors and microtubule-associated proteins. In this work, we propose a new mathematical model describing the motion of motor proteins on microtubules which incorporate a defective region. We take a mean-field approach derived from a first principle lattice model to study motor protein dynamics and density profiles. In particular, given a set of model parameters we obtain a closed-form expression for the equilibrium density profile along a given microtubule. We then verify the analytic results using mathematical analysis on the discrete model and Monte Carlo simulations. This work will contribute to the fundamental understanding of inhomogeneous microtubules providing insight into microscopic interactions that may result in the onset of neurodegenerative diseases. Our results for inhomogeneous microtubules are consistent with prior work studying the homogeneous case.

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.

Active cargo positioning in antiparallel transport networks

Cytoskeletal filaments assemble into dense parallel, antiparallel, or disordered networks, providing a complex environment for active cargo transport and positioning by molecular motors. The interplay between the network architecture and intrinsic motor properties clearly affects transport properties but remains poorly understood. Here, by using surface micropatterns of actin polymerization, we investigate stochastic transport properties of colloidal beads in antiparallel networks of overlapping actin filaments. We found that 200-nm beads coated with myosin Va motors displayed directed movements toward positions where the net polarity of the actin network vanished, accumulating there. The bead distribution was dictated by the spatial profiles of local bead velocity and diffusion coefficient, indicating that a diffusion-drift process was at work. Remarkably, beads coated with heavy-mero-myosin II motors showed a similar behavior. However, although velocity gradients were steeper with myosin II, the much larger bead diffusion observed with this motor resulted in less precise positioning. Our observations are well described by a 3-state model, in which active beads locally sense the net polarity of the network by frequently detaching from and reattaching to the filaments. A stochastic sequence of processive runs and diffusive searches results in a biased random walk. The precision of bead positioning is set by the gradient of net actin polarity in the network and by the run length of the cargo in an attached state. Our results unveiled physical rules for cargo transport and positioning in networks of mixed polarity.

A stochastic model that explains axonal organelle pileups induced by a reduction of molecular motors

Nerve cells are critically dependent on the transport of intracellular cargoes, which are moved by motor proteins along microtubule tracks. Impairments in this movement are thought to explain the focal accumulations of axonal cargoes and axonal swellings observed in many neurodegenerative diseases. In some cases, these diseases are caused by mutations that impair motor protein function, and genetic depletion of functional molecular motors has been shown to lead to cargo accumulations in axons. The evolution of these accumulations has been compared to the formation of traffic jams on a highway, but this idea remains largely untested. In this paper, we investigated the underlying mechanism of local axonal cargo accumulation induced by a global reduction of functional molecular motors in axons. We hypothesized that (i) a reduction in motor number leads to a reduction in the number of active motors on each cargo which in turn leads to less persistent movement, more frequent stops and thus shorter runs; (ii) as cargoes stop more frequently, they impede the passage of other cargoes, leading to local 'traffic jams'; and (iii) collisions between moving and stopping cargoes can push stopping cargoes further away from their microtubule tracks, preventing them from reattaching and leading to the evolution of local cargo accumulations. We used a lattice-based stochastic model to test whether this mechanism can lead to the cargo accumulation patterns observed in experiments. Simulation results of the model support the hypothesis and identify key questions that must be tested experimentally.

Self-organized system-size oscillation of a stochastic lattice-gas model

The totally asymmetric simple exclusion process (TASEP) is a paradigmatic stochastic model for nonequilibrium physics, and has been successfully applied to describe active transport of molecular motors along cytoskeletal filaments. Building on this simple model, we consider a two-lane lattice-gas model that couples directed transport (TASEP) to diffusive motion in a semiclosed geometry, and simultaneously accounts for spontaneous growth and particle-induced shrinkage of the system's size. This particular extension of the TASEP is motivated by the question of how active transport and diffusion might influence length regulation in confined systems. Surprisingly, we find that the size of our intrinsically stochastic system exhibits robust temporal patterns over a broad range of growth rates. More specifically, when particle diffusion is slow relative to the shrinkage dynamics, we observe quasiperiodic changes in length. We provide an intuitive explanation for the occurrence of these self-organized temporal patterns, which is based on the imbalance between the diffusion and shrinkage speed in the confined geometry. Finally, we formulate an effective theory for the oscillatory regime, which explains the origin of the oscillations and correctly predicts the dependence of key quantities, such as the oscillation frequency, on the growth rate.

Generic Transport Mechanisms for Molecular Traffic in Cellular Protrusions

Transport of molecular motors along protein filaments in a half-closed geometry is a common feature of biologically relevant processes in cellular protrusions. Using a lattice-gas model we study how the interplay between active and diffusive transport and mass conservation leads to localized domain walls and tip localization of the motors. We identify a mechanism for task sharing between the active motors (maintaining a gradient) and the diffusive motion (transport to the tip), which ensures that energy consumption is low and motor exchange mostly happens at the tip. These features are attributed to strong nearest-neighbor correlations that lead to a strong reduction of active currents, which we calculate analytically using an exact moment identity, and might prove useful for the understanding of correlations and active transport also in more elaborate systems.

Exclusion and Hierarchy of Time Scales Lead to Spatial Segregation of Molecular Motors in Cellular Protrusions

Molecular motors that carry cargo along biopolymer filaments within cells play a crucial role in the functioning of the cell. In particular, these motors are essential for the formation and maintenance of the cellular protrusions that play key roles in motility and specific functionalities, such as the stereocilia in hair cells. Typically, there are several species of motors, carrying different cargos, that share the same track. Furthermore, it was observed that in the mature stereocilia, the different motors occupy well-segregated bands as a function of distance from the tip. We use a totally asymmetric exclusion process model with two- and three-motor species, to study the conditions that give rise to such spatial patterns. We find that the well-segregated bands appear for motors with a strong hierarchy of attachment or detachment rates. This is a striking example of pattern formation in nonequilibrium, low-dimensional systems.

Motor protein traffic regulation by supply-demand balance of resources

In cells and in in vitro assays the number of motor proteins involved in biological transport processes is far from being unlimited. The cytoskeletal binding sites are in contact with the same finite reservoir of motors (either the cytosol or the flow chamber) and hence compete for recruiting the available motors, potentially depleting the reservoir and affecting cytoskeletal transport. In this work we provide a theoretical framework in which to study, analytically and numerically, how motor density profiles and crowding along cytoskeletal filaments depend on the competition of motors for their binding sites. We propose two models in which finite processive motor proteins actively advance along cytoskeletal filaments and are continuously exchanged with the motor pool. We first look at homogeneous reservoirs and then examine the effects of free motor diffusion in the surrounding medium. We consider as a reference situation recent in vitro experimental setups of kinesin-8 motors binding and moving along microtubule filaments in a flow chamber. We investigate how the crowding of linear motor proteins moving on a filament can be regulated by the balance between supply (concentration of motor proteins in the flow chamber) and demand (total number of polymerized tubulin heterodimers). We present analytical results for the density profiles of bound motors and the reservoir depletion, and propose novel phase diagrams that present the formation of jams of motor proteins on the filament as a function of two tuneable experimental parameters: the motor protein concentration and the concentration of tubulins polymerized into cytoskeletal filaments. Extensive numerical simulations corroborate the analytical results for parameters in the experimental range and also address the effects of diffusion of motor proteins in the reservoir. We then propose experiments for validating our models and discuss how the 'supply-demand' effects can regulate motor traffic also in in vivo conditions.

Traffic jams and shocks of molecular motors inside cellular protrusions

Molecular motors are involved in key transport processes inside actin-based cellular protrusions. The motors carry cargo proteins to the protrusion tip which participate in regulating the actin polymerization and play a key role in facilitating the growth and formation of such protrusions. It is observed that the motors accumulate at the tips of cellular protrusions and form aggregates that are found to drift towards the protrusion base at the rate of actin treadmilling. We present a one-dimensional driven lattice model, where motors become inactive after delivering their cargo at the tip, or by loosing their cargo to a cargoless neighbor. The results suggest that the experimental observations may be explained by the formation of traffic jams that form at the tip. The model is solved using a novel application of mean-field and shock analysis. We find a new class of shocks that undergo intermittent collapses. Extensions with attachment and detachment events and relevance to experiments are briefly described.

Theory of active transport in filopodia and stereocilia

The biological processes in elongated organelles of living cells are often regulated by molecular motor transport. We determined spatial distributions of motors in such organelles, corresponding to a basic scenario when motors only walk along the substrate, bind, unbind, and diffuse. We developed a mean-field model, which quantitatively reproduces elaborate stochastic simulation results as well as provides a physical interpretation of experimentally observed distributions of Myosin IIIa in stereocilia and filopodia. The mean-field model showed that the jamming of the walking motors is conspicuous, and therefore damps the active motor flux. However, when the motor distributions are coupled to the delivery of actin monomers toward the tip, even the concentration bump of G actin that they create before they jam is enough to speed up the diffusion to allow for severalfold longer filopodia. We found that the concentration profile of G actin along the filopodium is rather nontrivial, containing a narrow minimum near the base followed by a broad maximum. For efficient enough actin transport, this nonmonotonous shape is expected to occur under a broad set of conditions. We also find that the stationary motor distribution is universal for the given set of model parameters regardless of the organelle length, which follows from the form of the kinetic equations and the boundary conditions.

ACTIVE BIO-SYSTEMS: FROM SINGLE MOTOR MOLECULES TO COOPERATIVE CARGO TRANSPORT

Living cells contain a large number of molecular motors that convert the chemical energy released from nucleotide hydrolysis into mechanical work. This review focusses on stepping motors that move along cytoskeletal filaments. The behavior of these motors involves three distinct nonequilibrium processes that cover a wide range of length and time scales: (i) Directed stepping of single motors bound to a filament; (ii) Composite motor walks of single motors consisting of directed stepping interrupted by diffusive motion; and (iii) Cooperative transport by teams of several motors. On the molecular scale, the energy conversion of these motors leads to single steps along the filaments with a step size of about 10 nm. The corresponding chemomechanical coupling is governed by several distinct motor cycles, which represent the dominant pathways for different values of nucleotide concentrations and load force. For the kinesin motor, the competition of two such cycles determines the stall force, at which the motor velocity vanishes and the motor reverses the direction of its motion. Because of thermal noise, the stepping motors unbind from the filaments after a certain run time and run length. For kinesin, the run time is about 1 s and the run length is about 1 μm for high ATP concentration and low load force. On length scales that are large compared to the run length, a single motor undergoes composite walks consisting of directed stepping interrupted by diffusive motion. The relative importance of bound and unbound motor states depends on the binding and unbinding rates of the motors. The effective transport velocity and diffusion coefficient of the motors are determined by the geometry of the compartments, in which the motors move. The effective diffusion coefficient can be enhanced by several orders of magnitude if the motors undergo active diffusion by interacting with certain filament patterns. In vivo, stepping motors are responsible for the transport of vesicles and other types of intracellular cargo particles that shuttle between the different cell compartments. This cargo transport is usually performed by teams of motors. If all motors belong to the same molecular species, the cooperative action of the motors leads to uni-directional transport with a strongly increased run length and to a characteristic force dependence of the velocity distributions. If two antagonistic species of motors pull on the cargo, they perform a stochastic tug-of-war, which is characterized by a subtle force balance between the two motor teams and leads to seven distinct patterns of uni- and bi-directional transport. So far, all experimental observations on bi-directional transport are consistent with such a tug-of-war. Finally, the traffic of interacting motors is also briefly discussed. Depending on their mutual interactions and the compartment geometry, the motors form various spatio-temporal patterns such as traffic jams, and undergo nonequilibrium phase transitions between such transport patterns.

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.

Traffic of cytoskeletal motors with disordered attachment rates

Motivated by experimental results on the interplay between molecular motors and tau proteins, we extend lattice-based models of intracellular transport to include a second species of particle which locally influences the motor-filament attachment rate. We consider various exactly solvable limits of a stochastic multiparticle model before focusing on the low-motor-density regime. Here, an approximate treatment based on the random-walk behavior of single motors gives good quantitative agreement with simulation results for the tau dependence of the motor current. Finally, we discuss the possible physiological implications of our results.

Dynamic boundaries in asymmetric exclusion processes

We investigate the dynamics of a one-dimensional asymmetric exclusion process with Langmuir kinetics and a fluctuating wall. At the left-hand boundary, particles are injected onto the lattice; from there, the particles hop to the right. Along the lattice, particles can adsorb or desorb, and the right-hand boundary is defined by a wall particle. The confining wall particle has intrinsic forward and backward hopping, a net leftward drift, and cannot desorb. Performing Monte Carlo simulations and using a moving-frame finite segment approach coupled to mean field theory, we find the parameter regimes in which the wall acquires a steady-state position. In other regimes, the wall will either drift to the left and fall off the lattice at the injection site, or drift indefinitely to the right. Our results are discussed in the context of nonequilibrium phases of the system, fluctuating boundary layers, and particle densities in the laboratory frame versus the frame of the fluctuating wall.

Weakly coupled, antiparallel, totally asymmetric simple exclusion processes

We study a system composed of two parallel totally asymmetric simple exclusion processes with open boundaries, where the particles move in the two lanes in opposite directions and are allowed to jump to the other lane with rates inversely proportional to the length of the system. Stationary density profiles are determined and the phase diagram of the model is constructed in the hydrodynamic limit, by solving the differential equations describing the steady state of the system, analytically for vanishing total current and numerically for nonzero total current. The system possesses phases with a localized shock in the density profile in one of the lanes, similarly to exclusion processes endowed with nonconserving kinetics in the bulk. Besides, the system undergoes a discontinuous phase transition, where coherently moving delocalized shocks emerge in both lanes and the fluctuation of the global density is described by an unbiased random walk. This phenomenon is analogous to the phase coexistence observed at the coexistence line of the totally asymmetric simple exclusion process, however, as a consequence of the interaction between lanes, the density profiles are deformed and in the case of asymmetric lane change, the motion of the shocks is confined to a limited domain.

Directing Brownian motion by oscillating barriers

We consider the Brownian motion of a colloidal particle in a symmetric, periodic potential, whose potential barriers are subjected to temporal oscillations. Experimentally, the potential is generated by two arrays of trapped, negatively charged particles whose positions are periodically modulated with light forces. This results in a structured channel geometry of locally variable width. If all potential barriers are oscillating in synchrony, a resonance-like peak of the effective diffusion coefficient upon variation of the oscillation period is observed. For asynchronously oscillating barriers, the particle can be steered with great reliability into one or the other direction by properly choosing the oscillation periods of the different barriers along the channel.