Distinct speed and direction memories of migrating dendritic cells diversify their search strategies

Migrating cells exhibit various motility patterns, resulting from different migration mechanisms, cell properties, or cell-environment interactions. The complexity of cell dynamics is reflected, e.g., in the diversity of the observed forms of velocity autocorrelation function-which has been widely served as a measure of diffusivity and spreading. By analyzing the dynamics of migrating dendritic cells in vitro, we disentangle the contributions of direction θ and speed v to the velocity autocorrelation. We find that the ability of cells to maintain their speed or direction of motion is unequal, reflected in different temporal decays of speed and direction autocorrelation functions, ACv(t)∼t-1.2 and ACθ(t)∼t-0.5, respectively. The larger power-law exponent of ACv(t) indicates that the cells lose their speed memory considerably faster than the direction memory. Using numerical simulations, we investigate the influence of ACθ and ACv as well as the direction-speed cross correlation Cθ-v on the search time of a persistent random walker to find a randomly located target in confinement. Although ACθ and Cθ-v play the major roles, we find that the speed autocorrelation ACv can be also tuned to minimize the search time. Adopting an optimal ACv can reduce the search time even up to 10% compared with uncorrelated spontaneous speeds. Our results suggest that migrating cells can improve their search efficiency, especially in crowded environments, through the directional or speed persistence or the speed-direction correlation.

Effects of vimentin on the migration, search efficiency, and mechanical resilience of dendritic cells

Dendritic cells use amoeboid migration to pass through narrow passages in the extracellular matrix and confined tissue in search for pathogens and to reach the lymph nodes and alert the immune system. Amoeboid migration is a migration mode that, instead of relying on cell adhesion, is based on mechanical resilience and friction. To better understand the role of intermediate filaments in ameboid migration, we studied the effects of vimentin on the migration of dendritic cells. We show that the lymph node homing of vimentin-deficient cells is reduced in our in vivo experiments in mice. Lack of vimentin also reduces the cell stiffness, the number of migrating cells, and the migration speed in vitro in both 1D and 2D confined environments. Moreover, we find that lack of vimentin weakens the correlation between directional persistence and migration speed. Thus, vimentin-expressing dendritic cells move faster in straighter lines. Our numerical simulations of persistent random search in confined geometries verify that the reduced migration speed and the weaker correlation between the speed and direction of motion result in longer search times to find regularly located targets. Together, these observations show that vimentin enhances the ameboid migration of dendritic cells, which is relevant for the efficiency of their random search for pathogens.

Kinematics of persistent random walkers with two distinct modes of motion

We study the stochastic motion of active particles that undergo spontaneous transitions between two distinct modes of motion. Each mode is characterized by a velocity distribution and an arbitrary (anti)persistence. We present an analytical formalism to provide a quantitative link between these two microscopic statistical properties of the trajectory and macroscopically observable transport quantities of interest. For exponentially distributed residence times in each state, we derive analytical expressions for the initial anomalous exponent, the characteristic crossover time to the asymptotic diffusive dynamics, and the long-term diffusion constant. We also obtain an exact expression for the time evolution of the mean square displacement over all timescales and provide a recipe to obtain higher displacement moments. Our approach enables us to disentangle the combined effects of velocity, persistence, and switching probabilities between the two states on the kinematics of particles in a wide range of stochastic active or passive processes and to optimize the transport quantities of interest with respect to any of the particle dynamics properties.

Gravity Governs Shear Localization in Confined Dense Granular Flows

The prediction of flow profiles of slowly sheared granular materials is a major geophysical and industrial challenge. Understanding the role of gravity is particularly important for future planetary exploration in varying gravitational environments. Using the principle of minimization of energy dissipation, and combining experiments and variational analysis, we disentangle the contributions of the gravitational acceleration, confining pressure, and layer thickness on shear strain localization induced by moving fault boundaries at the bottom of a granular layer. The flow profile is independent of the gravity for geometries with a free top surface. However, under a confining pressure or if the sheared layer withstands the weight of the upper layers, increasing gravity promotes the transition from closed shear zones buried in the bulk to open ones that intersect the top surface. We show that the center position and width of the shear zone and the axial angular velocity at the top surface follow universal scaling laws when properly scaled by the gravity, applied pressure, and layer thickness. Our finding that the flow profiles lie on a universal master curve opens the possibility to predict the quasistatic shear flow of granular materials in varying gravitational environments.

Orientational memory of active particles in multistate non-Markovian processes

The orientational memory of particles can serve as an effective measure of diffusivity, spreading, and search efficiency in complex stochastic processes. We develop a theoretical framework to describe the decay of directional correlations in a generic class of stochastic active processes consisting of distinct states of motion characterized by their persistence and switching probabilities between the states. For exponentially distributed sojourn times, the orientation autocorrelation is analytically derived and the characteristic times of its crossovers are obtained in terms of the persistence of each state and the switching probabilities. We show how nonexponential sojourn-time distributions of interest, such as Gaussian and power-law distributions, can result from history-dependent transitions between the states. The relaxation behavior of the correlation function in such non-Markovian processes is governed by the history dependence of the switching probabilities and cannot be solely determined by the mean sojourn times of the states.

Evolution of shear zones in granular packings under pressure

Stress transmission in realistic granular media often occurs under external load and in the presence of boundary slip. We investigate shear localization in a split-bottom Couette cell with smooth walls subject to a confining pressure experimentally and by means of numerical simulations. We demonstrate how the characteristics of the shear zone, such as its center position and width, evolve as the confining pressure and wall slip modify the local effective friction coefficient of the material. For increasing applied pressure, the shear zone evolves toward the center of the cylinder and grows wider and the angular velocity reduces compared to the driving rate of the bottom disk. Moreover, the presence of slip promotes the transition from open shear zones at the top surface to closed shear zones inside the bulk. We also systematically vary the ratio of the effective friction near the bottom plate and in the bulk in simulations and observe the resulting impact on the surface flow profile. Besides the boundary conditions and external load, material properties such as grain size are also known to influence the effective friction coefficient. However, our numerical results reveal that the center position and width of the shear zone are insignificantly affected by the choice of the grain size as far as it remains small compared to the radius of the rotating bottom disk.

Shearing of granular materials in a confined split-bottom Couette cell

Formation of shear bands is one of the most remarkable phenomena in the dynamics of granular matter. Several parameters have been so far identified to influence the behavior of the shear bands. We carried out experiments to investigate the evolution of the shear bands in the split-bottom Couette cell in the presence of confining pressure. We employed the Particle Image Velocimetry (PIV) to characterize the shear band both in the absence and presence of external pressure. Our results show that the location and width of the shear band are affected by both the confining pressure and the filling height. The shear zone evolves towards the middle of the cylinder and expands to a broader region with increasing applied pressure or filling height; also the angular velocity decreases relative to the rotation rate of the bottom disk. Our findings are consistent with prior empirical observations on the formation of wide shear bands at free surfaces.

Persistence-Speed Coupling Enhances the Search Efficiency of Migrating Immune Cells

Migration of immune cells within the human body allows them to fulfill their main function of detecting pathogens. We present experimental evidence showing the optimality of the search strategy of these cells, which is of crucial importance to achieve an efficient immune response. We find that the speed and directional persistence of migrating dendritic cells in our in vitro experiments are highly correlated, which enables them to reduce their search time. We introduce theoretically a new class of random search optimization problems by minimizing the mean first-passage time (MFPT) with respect to the strength of the coupling between influential parameters. We derive an analytical expression for the MFPT in a confined geometry and verify that the correlated motion enhances the search efficiency if the mean persistence length is sufficiently shorter than the confinement size. Our correlated search optimization approach provides an efficient searching recipe and predictive power in a broad range of correlated stochastic processes.

Trapping in and Escape from Branched Structures of Neuronal Dendrites

We present a coarse-grained model for stochastic transport of noninteracting chemical signals inside neuronal dendrites and show how first-passage properties depend on the key structural factors affected by neurodegenerative disorders or aging: the extent of the tree, the topological bias induced by segmental decrease of dendrite diameter, and the trapping probabilities in biochemical cages and growth cones. We derive an exact expression for the distribution of first-passage times, which follows a universal exponential decay in the long-time limit. The asymptotic mean first-passage time exhibits a crossover from power-law to exponential scaling upon reducing the topological bias. We calibrate the coarse-grained model parameters and obtain the variation range of the mean first-passage time when the geometrical characteristics of the dendritic structure evolve during the course of aging or neurodegenerative disease progression (a few disorders for which clear trends for the pathological changes of dendritic structure have been reported in the literature are chosen and studied). We prove the validity of our analytical approach under realistic fluctuations of structural parameters by comparison to the results of Monte Carlo simulations. Moreover, by constructing local structural irregularities, we analyze the resulting influence on transport of chemical signals and formation of heterogeneous density patterns. Because neural functions rely on chemical signal transmission to a large extent, our results open the possibility of establishing a direct link between the disease progression and neural functions.

Run-and-pause dynamics of cytoskeletal motor proteins

Cytoskeletal motor proteins are involved in major intracellular transport processes which are vital for maintaining appropriate cellular function. When attached to cytoskeletal filaments, the motor exhibits distinct states of motility: active motion along the filaments, and pause phase in which it remains stationary for a finite time interval. The transition probabilities between motion and pause phases are asymmetric in general, and considerably affected by changes in environmental conditions which influences the efficiency of cargo delivery to specific targets. By considering the motion of individual non-interacting molecular motors on a single filament as well as a dynamic filamentous network, we present an analytical model for the dynamics of self-propelled particles which undergo frequent pause phases. The interplay between motor processivity, structural properties of filamentous network, and transition probabilities between the two states of motility drastically changes the dynamics: multiple transitions between different types of anomalous diffusive dynamics occur and the crossover time to the asymptotic diffusive or ballistic motion varies by several orders of magnitude. We map out the phase diagrams in the space of transition probabilities, and address the role of initial conditions of motion on the resulting dynamics.

Tracking of plus-ends reveals microtubule functional diversity in different cell types

Many cellular processes are tightly connected to the dynamics of microtubules (MTs). While in neuronal axons MTs mainly regulate intracellular trafficking, they participate in cytoskeleton reorganization in many other eukaryotic cells, enabling the cell to efficiently adapt to changes in the environment. We show that the functional differences of MTs in different cell types and regions is reflected in the dynamic properties of MT tips. Using plus-end tracking proteins EB1 to monitor growing MT plus-ends, we show that MT dynamics and life cycle in axons of human neurons significantly differ from that of fibroblast cells. The density of plus-ends, as well as the rescue and catastrophe frequencies increase while the growth rate decreases toward the fibroblast cell margin. This results in a rather stable filamentous network structure and maintains the connection between nucleus and membrane. In contrast, plus-ends are uniformly distributed along the axons and exhibit diverse polymerization run times and spatially homogeneous rescue and catastrophe frequencies, leading to MT segments of various lengths. The probability distributions of the excursion length of polymerization and the MT length both follow nearly exponential tails, in agreement with the analytical predictions of a two-state model of MT dynamics.

Cargo binding promotes KDEL receptor clustering at the mammalian cell surface

Transmembrane receptor clustering is a ubiquitous phenomenon in pro- and eukaryotic cells to physically sense receptor/ligand interactions and subsequently translate an exogenous signal into a cellular response. Despite that receptor cluster formation has been described for a wide variety of receptors, ranging from chemotactic receptors in bacteria to growth factor and neurotransmitter receptors in mammalian cells, a mechanistic understanding of the underlying molecular processes is still puzzling. In an attempt to fill this gap we followed a combined experimental and theoretical approach by dissecting and modulating cargo binding, internalization and cellular response mediated by KDEL receptors (KDELRs) at the mammalian cell surface after interaction with a model cargo/ligand. Using a fluorescent variant of ricin toxin A chain as KDELR-ligand (eGFP-RTA^H/KDEL), we demonstrate that cargo binding induces dose-dependent receptor cluster formation at and subsequent internalization from the membrane which is associated and counteracted by anterograde and microtubule-assisted receptor transport to preferred docking sites at the plasma membrane. By means of analytical arguments and extensive numerical simulations we show that cargo-synchronized receptor transport from and to the membrane is causative for KDELR/cargo cluster formation at the mammalian cell surface.

Enhanced diffusion and anomalous transport of magnetic colloids driven above a two-state flashing potential

We combine experiments and theory to investigate the diffusive and the subdiffusive dynamics of paramagnetic colloids driven above a two-state flashing potential. The magnetic potential was realized by periodically modulating the stray field of a magnetic bubble lattice in a uniaxial ferrite garnet film. At large amplitudes H0 of the driving field, the dynamics of the particle resemble an ordinary random walk with a frequency-dependent diffusion coefficient. However, subdiffusive and oscillatory dynamics at short time scales are observed when decreasing H0. We present a persistent random walk model to elucidate the underlying mechanism of motion, and perform numerical simulations to demonstrate that the anomalous motion originates from the dynamic disorder in the structure of the magnetic lattice, induced by the slightly irregular shape of bubbles.

Persistent-random-walk approach to anomalous transport of self-propelled particles

The motion of self-propelled particles is modeled as a persistent random walk. An analytical framework is developed that allows the derivation of exact expressions for the time evolution of arbitrary moments of the persistent walk's displacement. It is shown that the interplay of step length and turning angle distributions and self-propulsion produces various signs of anomalous diffusion at short time scales and asymptotically a normal diffusion behavior with a broad range of diffusion coefficients. The crossover from the anomalous short-time behavior to the asymptotic diffusion regime is studied and the parameter dependencies of the crossover time are discussed. Higher moments of the displacement distribution are calculated and analytical expressions for the time evolution of the skewness and the kurtosis of the distribution are presented.

Anomalous diffusion of self-propelled particles in directed random environments

We theoretically study the transport properties of self-propelled particles on complex structures, such as motor proteins on filament networks. A general master equation formalism is developed to investigate the persistent motion of individual random walkers, which enables us to identify the contributions of key parameters: the motor processivity, and the anisotropy and heterogeneity of the underlying network. We prove the existence of different dynamical regimes of anomalous motion, and that the crossover times between these regimes as well as the asymptotic diffusion coefficient can be increased by several orders of magnitude within biologically relevant control parameter ranges. In terms of motion in continuous space, the interplay between stepping strategy and persistency of the walker is established as a source of anomalous diffusion at short and intermediate time scales.

Evolution of the force distributions in jammed packings of soft particles

The evolution of the force distributions during the isotropic compression of two-dimensional packings of soft frictional particles is investigated numerically. Regardless of the applied deformation, the normal contact force distribution P(f(n)) can be fitted by the product of a power law, and a stretched exponential, while the tangential force distribution P(f(t)) is fitted well by a Gaussian. With increasing strain, the asymptotic behavior at large forces does not change, but both P(f(n)) and P(f(t)) exhibit a broadening, even though, when scaled with the average forces, their widths decrease. Furthermore, the distribution of friction mobilization P(η) is a decreasing function of η=|f(t)|/(μf(n)), except for an increased probability of fully mobilized contacts (η=1). The excess coordination number of the packings increases with the applied strain, indicating that the more a packing is compressed the more stable it becomes.

Coexistence and transition between shear zones in slow granular flows

We report experiments on slow granular flows in a split-bottom Couette cell that show novel strain localization features. Nontrivial flow profiles have been observed which are shown to be the consequence of simultaneous formation of shear zones in the bulk and at the boundaries. The fluctuating band model based on a minimization principle can be fitted to the experiments over a large variation of morphology and filling height with one single fit parameter, the relative friction coefficient μ(rel) between wall and bulk. The possibility of multiple shear zone formation is controlled by μ(rel). Moreover, we observe that the symmetry of an initial state, with coexisting shear zones at both side walls, breaks spontaneously below a threshold value of the shear velocity. A dynamical transition between two asymmetric flow states happens over a characteristic time scale which depends on the shear strength.

Long-range interactions in randomly driven granular fluids

We study the long-range spatial correlations in the nonequilibrium steady state of a randomly driven granular fluid with the emphasis on obtaining the explicit form of the static structure factors. The presence of immobile particles immersed in such a fluidized bed of fine particles leads to the confinement of the fluctuation spectrum of the hydrodynamic fields, which results in effective long-range interactions between the intruders. The analytical predictions are in agreement with the results of discrete element method simulations. By changing the shape and orientation of the intruders, we address how the effective force is affected by small changes in the boundary conditions.

Nonadditivity of fluctuation-induced forces in fluidized granular media

We investigate the effective long-range interactions between intruder particles immersed in a randomly driven granular fluid. The effective Casimir-like force between two intruders, induced by the fluctuations of the hydrodynamic fields, can change its sign when varying the control parameters: the volume fraction, the distance between the intruders, and the restitution coefficient. More interestingly, by inserting more intruders, we verify that the fluctuation-induced interaction is not pairwise additive. The simulation results are qualitatively consistent with the theoretical predictions based on mode coupling calculations. These results shed new light on the underlying mechanisms of collective behaviors in fluidized granular media.

Influence of polydispersity on micromechanics of granular materials

We study the effect of polydispersity on the macroscopic physical properties of granular packings in two and three dimensions. A mean-field approach is developed to approximate the macroscale quantities as functions of the microscopic ones. We show that the trace of the fabric and stress tensors are proportional to the mean packing properties (e.g., packing fraction, average coordination number, and average normal force) and dimensionless correction factors, which depend only on the moments of the particle-size distribution. Similar results are obtained for the elements of the stiffness tensor of isotropic packings in the linear affine response regime. Our theoretical predictions are in good agreement with the simulation results.

Extent of force indeterminacy in packings of frictional rigid disks

Static packings of frictional rigid particles are investigated by means of discrete element simulations. We explore the ensemble of allowed force realizations in the space of contact forces for a given packing structure. We estimate the extent of force indeterminacy with different methods. The indeterminacy exhibits a nonmonotonic dependence on the interparticle friction coefficient. We verify directly that larger force indeterminacy is accompanied by a more robust behavior against local perturbations. We also investigate the local indeterminacy of individual contact forces. The probability distribution of local indeterminacy changes its shape depending on friction. We find that local indeterminacy tends to be larger on force chains for intermediate friction. This correlation disappears in the large friction limit.

Diffusive transport of light in a two-dimensional disordered packing of disks: analytical approach to transport mean free path

We study photon diffusion in a two-dimensional random packing of monodisperse disks as a simple model of granular media and wet foams. We assume that the intensity reflectance of disks is a constant r . We present an analytic expression for the transport mean free path l;{*} in terms of the velocity of light in the disks and host medium, radius R and packing fraction of the disks, and the intensity reflectance. For glass beads immersed in air or water, we estimate transport mean free paths about half the experimental ones. For air bubbles immersed in water, l;{*}R is a linear function of 1epsilon , where epsilon is the liquid volume fraction of the model wet foam. This throws light on the empirical law of Vera [Appl. Opt. 40, 4210 (2001)] and promotes more realistic models.

Unjamming due to local perturbations in granular packings with and without gravity

We investigate the unjamming response of disordered packings of frictional hard disks with the help of computer simulations. First, we generate jammed configurations of the disks and then force them to move again by local perturbations. We study the spatial distribution of the stress and displacement response and find long range effects of the perturbation in both cases. We record the penetration depth of the displacements and the critical force that is needed to make the system yield. These quantities are tested in two types of systems: in ideal homogeneous packings in zero gravity and in packings settled under gravity. The penetration depth and the critical force are sensitive to the interparticle friction coefficient. Qualitatively, the same nonmonotonic friction dependence is found both with and without gravity, however the location of the extrema are at different friction values. We discuss the role of the connectivity of the contact network and of the pressure gradient in the unjamming response.

Unjamming of granular packings due to local perturbations: stability and decay of displacements

We study the mechanical response generated by local deformations in jammed packings of rigid disks. Based on discrete element simulations we determine the critical force of the local perturbation that is needed to break the mechanical equilibrium and examine the generated displacement field. Displacements decay as a power law of the distance from the perturbation point. The decay exponent and the critical force exhibit a nontrivial dependence on the friction: Both quantities are nonmonotonic and have a sharp maximum at the friction coefficient 0.1. We find that the mechanical response properties are closely related to the problem of force indeterminate where similar nonmonotonic behavior was observed previously. We establish a direct connection between the critical force and the ensemble of static force networks.

Characteristics of vehicular traffic flow at a roundabout

We construct a stochastic cellular automata model for the description of vehicular traffic at a roundabout designed at the intersection of two perpendicular streets. The vehicular traffic is controlled by a self-organized scheme in which traffic lights are absent. This controlling method incorporates a yield-at-entry strategy for the approaching vehicles to the circulating traffic flow in the roundabout. Vehicular dynamics is simulated and the delay experienced by the traffic at each individual street is evaluated. We discuss the impact of the geometrical properties of the roundabout on the total delay. We compare our results with traffic-light signalization schemes, and obtain the critical traffic volume over which the intersection is optimally controlled through traffic-light signalization schemes.