Motor protein traffic regulation by supply-demand balance of resources

Availability versus carrying capacity: Phases of asymmetric exclusion processes competing for finite pools of resources

We address how the interplay between the finite availability and carrying capacity of particles at different parts of a spatially extended system can control the steady-state currents and density profiles in the one-dimensional current-carrying lanes connecting the different parts of the system. To study this, we set up a minimal model consisting of two particle reservoirs of the same finite carrying capacity connected by two equally sized antiparallel totally asymmetric simple exclusion processes (TASEPs). We focus on the steady-state currents and particle density profiles in the two TASEP lanes. The ensuing phases and the phase diagrams, which can be remarkably complex, are parametrized by the model parameters defining particle exchange between the TASEP lanes and the reservoirs and the filling fraction of the particles that determine the total resources available. These parameters may be tuned to make the densities of the two TASEP lanes globally uniform or piece-wise continuous in the form of a combination of a single localized domain wall and a spatially constant density or a pair of delocalized domain walls. Our model reveals that the two reservoirs can be preferentially populated or depopulated in the steady states.

Collective dynamics on constrained three-lane exclusion process

The motivation behind the proposed study stems from multilane traffic systems with finite availability of particles. Our investigation revolves around a totally asymmetric simple exclusion process incorporating a finite reservoir and the occurrence of lane-switching phenomena. The study delves into the system's characteristics, including phase diagrams, density profiles, phase transitions, finite-size effects, and shock positions. These analyses concern the number of particles within the system and various weak-coupling rates. The outcomes obtained from the generalized mean-field theory are cross-validated against the results derived from Monte Carlo simulations. In scrutinizing the system's dynamics, we observe several noteworthy observations. Notably, we identified critical mixed profiles featuring instances of double shocks. The system, intriguingly, demonstrates a transition known as reentrance transition. The study also reports a rare phenomenon, namely, the jumping effect within the shock profile, adding a layer of complexity to the system's behavior and proving the significance of limited resources.

Availability, storage capacity, and diffusion: Stationary states of an asymmetric exclusion process connected to two reservoirs

We explore how the interplay of finite availability, carrying capacity of particles at different parts of a spatially extended system, and particle diffusion between them control the steady-state currents and density profiles in a one-dimensional current-carrying channel connecting the different parts of the system. To study this, we construct a minimal model consisting of two particle reservoirs of finite carrying capacities connected by a totally asymmetric simple exclusion process (TASEP). In addition to particle transport via TASEP between the reservoirs, the latter can also directly exchange particles via Langmuir kinetics-like processes, modeling particle diffusion between them that can maintain a steady current in the system. We calculate the steady-state density profiles and the associated particle currents in the TASEP lane. The resulting phases and the phase diagrams are quite different from an open TASEP, and are characterized by the model parameters defining particle exchanges between the TASEP and the reservoirs, direct particle exchanges between the reservoirs, and the filling fraction of the particles that determines the total resources available. These parameters can be tuned to make the density on the TASEP lane globally uniform or piecewise continuous, and can make the two reservoirs preferentially populated or depopulated.

Competition for resources in an exclusion model with biased lane-changing mechanism

The motivation for the proposed work is drawn from the attachment-detachment observed in biological and physical transport processes that entail finite resources. We investigate the influence of limited particle availability on particle dynamics within two parallel totally asymmetric simple exclusion lanes, with one lane incorporating only particle detachment and the other considering particle attachment. We establish a theoretical framework by employing vertical mean-field theory in conjunction with singular perturbation technique. The analytical findings are supported by numerical and stochastic validation using a finite-difference scheme and the Gillespie algorithm. By utilizing these approaches, we scrutinize various stationary properties, including particle densities, phase boundaries, and particle currents for both lanes. Our analysis reveals that the complexity of the phase diagram exhibits a nonmonotonic trend in the number of stationary phases as the particle count increases. Each phase diagram is constructed with respect to the intrinsic boundary parameters, illustrating both bulk and surface transitions occurring within the lanes. The interplay between finite resources and coupling mechanisms gives rise to two phases involving upward shock in one of the lanes, while two phases exhibit synchronized downward shock in both lanes. Finally, we delve into shock dynamics to comprehend critical phase transitions occurring in the system.

Totally asymmetric simple exclusion process with local resetting in a resource-constrained environment

Inspired by the process of mRNA translation, in which the stochastic degradation of mRNA-ribosome machinery is modeled by the resetting dynamics, we study an open totally asymmetric simple exclusion process with local resetting at the entry site in a resource-constrained environment. The effect of constrained resources on the stationary properties of the system has been comprehended in the form of the filling factor. The mean-field approximations are utilized to obtain stationary state features, such as density profiles and phase diagrams. The phase diagram possesses pure phases as well as coexisting phases, including a low-density-high-density phase separation, which did not manifest under periodic boundary conditions despite the system being closed there as well. The role of the resetting rate has been investigated on the stationary properties of the system, depending on how the filling factor scales with the system size. In contrast to the resetting model for infinite resources, two distinct phase transitions are observed for the smaller values of the filling factor leading to a change in the topology of the phase diagram. The impact of the resetting rate along with the finite-size effect has also been examined on the shock dynamics. All the mean-field results are found in remarkable agreement with the Monte Carlo simulations.

Physical modeling of ribosomes along messenger RNA: Estimating kinetic parameters from ribosome profiling experiments using a ballistic model

Gene expression is the synthesis of proteins from the information encoded on DNA. One of the two main steps of gene expression is the translation of messenger RNA (mRNA) into polypeptide sequences of amino acids. Here, by taking into account mRNA degradation, we model the motion of ribosomes along mRNA with a ballistic model where particles advance along a filament without excluded volume interactions. Unidirectional models of transport have previously been used to fit the average density of ribosomes obtained by the experimental ribo-sequencing (Ribo-seq) technique in order to obtain the kinetic rates. The degradation rate is not, however, accounted for and experimental data from different experiments are needed to have enough parameters for the fit. Here, we propose an entirely novel experimental setup and theoretical framework consisting in splitting the mRNAs into categories depending on the number of ribosomes from one to four. We solve analytically the ballistic model for a fixed number of ribosomes per mRNA, study the different regimes of degradation, and propose a criterion for the quality of the inverse fit. The proposed method provides a high sensitivity to the mRNA degradation rate. The additional equations coming from using the monosome (single ribosome) and polysome (arbitrary number) ribo-seq profiles enable us to determine all the kinetic rates in terms of the experimentally accessible mRNA degradation rate.

Axonal transport during injury on a theoretical axon

Neurodevelopment, plasticity, and cognition are integral with functional directional transport in neuronal axons that occurs along a unique network of discontinuous polar microtubule (MT) bundles. Axonopathies are caused by brain trauma and genetic diseases that perturb or disrupt the axon MT infrastructure and, with it, the dynamic interplay of motor proteins and cargo essential for axonal maintenance and neuronal signaling. The inability to visualize and quantify normal and altered nanoscale spatio-temporal dynamic transport events prevents a full mechanistic understanding of injury, disease progression, and recovery. To address this gap, we generated DyNAMO, a Dynamic Nanoscale Axonal MT Organization model, which is a biologically realistic theoretical axon framework. We use DyNAMO to experimentally simulate multi-kinesin traffic response to focused or distributed tractable injury parameters, which are MT network perturbations affecting MT lengths and multi-MT staggering. We track kinesins with different motility and processivity, as well as their influx rates, in-transit dissociation and reassociation from inter-MT reservoirs, progression, and quantify and spatially represent motor output ratios. DyNAMO demonstrates, in detail, the complex interplay of mixed motor types, crowding, kinesin off/on dissociation and reassociation, and injury consequences of forced intermingling. Stalled forward progression with different injury states is seen as persistent dynamicity of kinesins transiting between MTs and inter-MT reservoirs. DyNAMO analysis provides novel insights and quantification of axonal injury scenarios, including local injury-affected ATP levels, as well as relates these to influences on signaling outputs, including patterns of gating, waves, and pattern switching. The DyNAMO model significantly expands the network of heuristic and mathematical analysis of neuronal functions relevant to axonopathies, diagnostics, and treatment strategies.

Multiple reentrance transitions in exclusion process with finite reservoir

The proposed study is motivated by the scenario of two-way vehicular traffic. We consider a totally asymmetric simple exclusion process in the presence of a finite reservoir along with the particle attachment, detachment, and lane-switching phenomena. The various system properties in terms of phase diagrams, density profiles, phase transitions, finite size effect, and shock position are analyzed, considering the available number of particles in the system and different values of coupling rate, by employing the generalized mean-field theory and the obtained results are detected to be a good match with the Monte Carlo simulation outcomes. It is discovered that the finite resources significantly affect the phase diagram for different coupling rate values, which leads to nonmonotonic changes in the number of phases in the phase plane for comparatively minor lane-changing rates and produces various exciting features. We calculate the critical value of the total number of particles in the system at which the multiple phases in the phase diagram appear or disappear. The competition between the limited particles, bidirectional motion, Langmuir kinetics, and particle lane-shifting behavior yields unanticipated and unique mixed phases, including the double shock phase, multiple reentrance and bulk-induced phase transitions, and phase segregation of the single shock phase.

Interplay of reservoirs in a bidirectional system

Motivated by the interplay of multiple species in several real world transport processes, we propose a bidirectional totally asymmetric simple exclusion process with two finite particle reservoirs regulating the inflow of oppositely directed particles corresponding to two different species. The system's stationary characteristics, such as densities, currents, etc., are investigated using a theoretical framework based on mean-field approximation and are supported by extensive Monte Carlo simulations. The impact of individual species populations, quantified by filling factor, has been comprehensively analyzed considering both equal and unequal conditions. For the equal case, the system exhibits the spontaneous symmetry-breaking phenomena and admits both symmetric as well as asymmetric phases. Moreover, the phase diagram exhibits a different asymmetric phase and displays a nonmonotonic variation in the number of phases with respect to the filling factor. For unequal filling factors, the phase schema can display at most five phases including a phase that shows maximal current for one of the species.

Exclusion process with scaled resources: Delocalized shocks and interplay of reservoirs

In this paper we study a conserved system comprised of two directed lanes having identical dynamics and two reservoirs with scaled resources that are strategically connected to the boundaries of the lanes, forming a ringlike structure. The steady-state properties of the system have been analyzed in the framework of mean-field theory. Our findings display a rich behavior, emphasizing the nontrivial effects of incorporating two reservoirs. As a consequence, two distinct phases that admit delocalized shocks emerge and occupy a significant region in the phase diagram. Moreover in one of theses phases, each lane admits a delocalized shock whose movements are perfectly synchronized. In another phase, the single shock in the system may traverse both lanes or remain restricted to a single lane, depending upon the size of the system. All the findings are validated by Monte Carlo simulations.

Asymmetric exclusion processes with fixed resources: Reservoir crowding and steady states

We study the reservoir crowding effect by considering the nonequilibrium steady states of an asymmetric exclusion process (TASEP) coupled to a reservoir with fixed available resources and dynamically coupled entry and exit rate. We elucidate how the steady states are controlled by the interplay between the coupled entry and exit rates, both being dynamically controlled by the reservoir population, and the fixed total particle number in the system. The TASEP can be in the low-density, high-density, maximal current, and shock phases. We show that such a TASEP is different from an open TASEP for all values of available resources: here the TASEP can support only localized domain walls for any (finite) amount of resources that do not tend to delocalize even for large resources, a feature attributed to the form of the dynamic coupling between the entry and exit rates. Furthermore, in the limit of infinite resources, in contrast to an open TASEP, the TASEP can be found in its high-density phase only for any finite values of the control parameters, again as a consequence of the coupling between the entry and exit rates.

Exclusion process on two intersecting lanes with constrained resources: Symmetry breaking and shock dynamics

We present a study of the exclusion process on a peculiar topology of network with two intersecting lanes, competing for the particles in a reservoir with finite capacity. To provide a theoretical ground for our findings, we exploit mean-field approximation along with domain-wall theory. The stationary properties of the system, including phase transitions, density profiles, and position of the domain wall are derived analytically. Under the similar dynamical rules, the particles of both lanes interact only at the intersected site. The symmetry of the system is maintained until the number of particles do not exceed the total number of sites. However, beyond this, the symmetry breaking phenomenon occurs, resulting in the appearance of asymmetric phases and continues to persist even for an infinite number of particles. The complexity of the phase diagram shows a nonmonotonic behavior with an increasing number of particles in the system. A bulk induced shock appears in a symmetric phase, whereas, a boundary induced shock is observed in the symmetric as well as the asymmetric phase. Monitoring the location of localized shock with increasing entry of particles, we explain the possible phase transitions. The theoretical results are supported by extensive Monte Carlo simulations and explained using simple physical arguments.

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.

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.

Active and passive transport of cargo in a corrugated channel: A lattice model study

Inside cells, cargos such as vesicles and organelles are transported by molecular motors to their correct locations via active motion on cytoskeletal tracks and passive, Brownian diffusion. During the transportation of cargos, motor-cargo complexes (MCCs) navigate the confining and crowded environment of the cytoskeletal network and other macromolecules. Motivated by this, we study a minimal two-state model of motor-driven cargo transport in confinement and predict transport properties that can be tested in experiments. We assume that the motion of the MCC is directly affected by the entropic barrier due to confinement if it is in the passive, unbound state but not in the active, bound state where it moves with a constant bound velocity. We construct a lattice model based on a Fokker Planck description of the two-state system, study it using a kinetic Monte Carlo method and compare our numerical results with analytical expressions for a mean field limit. We find that the effect of confinement strongly depends on the bound velocity and the binding kinetics of the MCC. Confinement effectively reduces the effective diffusivity and average velocity, except when it results in an enhanced average binding rate and thereby leads to a larger average velocity than when unconfined.

Limited Resources Induce Bistability in Microtubule Length Regulation

The availability of protein is an important factor for the determination of the size of the mitotic spindle. Involved in spindle-size regulation is kinesin-8, a molecular motor and microtubule (MT) depolymerase, which is known to tightly control MT length. Here, we propose and analyze a theoretical model in which kinesin-induced MT depolymerization competes with spontaneous polymerization while supplies of both tubulin and kinesin are limited. In contrast to previous studies where resources were unconstrained, we find that, for a wide range of concentrations, MT length regulation is bistable. We test our predictions by conducting in vitro experiments and find that the bistable behavior manifests in a bimodal MT length distribution.

Far-from-equilibrium bidirectional transport system with constrained entrances competing for pool of limited resources

Motivated by the wide occurrence of limited resources in many real-life systems, we investigate two-lane totally asymmetric simple exclusion process with constrained entrances under finite supply of particles. We analyze the system within the framework of mean-field theory and examine various complex phenomena, including phase separation, phase transition, and symmetry breaking. Based on the theoretical analysis, we analytically derive the phase boundaries for various symmetric as well as asymmetric phases. It has been observed that the symmetry-breaking phenomenon initiates even for very small number of particles in the system. The phases with broken symmetry originates as shock-low density phase under limited resources, which is in contrast to the scenario with infinite number of particles. As expected, the symmetry breaking continues to persist even for higher values of system particles. Seven stationary phases are observed, with three of them exhibiting symmetry-breaking phenomena. The critical values of a total number of system particles, beyond which various symmetrical and asymmetrical phases appear and disappear are identified. Theoretical outcomes are supported by extensive Monte Carlo simulations. Finally, the size-scaling effect and symmetry-breaking phenomenon on the simulation results have also been examined based on particle density histograms.

Pentagalloylglucose Blocks the Nuclear Transport and the Process of Nucleocapsid Egress to Inhibit HSV-1 Infection

Herpes simplex virus type 1 (HSV-1), a widespread virus, causes a variety of human viral diseases worldwide. The serious threat of drug-resistance highlights the extreme urgency to develop novel antiviral drugs with different mechanisms of action. Pentagalloylglucose (PGG) is a natural polyphenolic compound with significant anti-HSV activity; however, the mechanisms underlying its antiviral activity need to be defined by further studies. In this study, we found that PGG treatment delays the nuclear transport process of HSV-1 particles by inhibiting the upregulation of dynein (a cellular major motor protein) induced by HSV-1 infection. Furthermore, PGG treatment affects the nucleocapsid egress of HSV-1 by inhibiting the expression and disrupting the cellular localization of pEGFP-UL31 and pEGFP-UL34, which are indispensable for HSV-1 nucleocapsid egress from the nucleus. However, the over-expression of pEGFP-UL31 and pEGFP-UL34 could decrease the antiviral effect of PGG. In this study, for the first time, the antiviral activity of PGG against acyclovir-resistant virus was demonstrated in vitro, and the possible mechanisms of its anti-HSV activities were identified based on the inhibition of nuclear transport and nucleocapsid egress in HSV-1. It was further confirmed that PGG could be a promising candidate for HSV therapy, especially for drug-resistant strains.