Refine
Has Fulltext
- no (16)
Document Type
- Article (16) (remove)
Language
- English (16)
Is part of the Bibliography
- yes (16)
Keywords
- critical avalanche dynamics (2)
- memory and delay (2)
- neuronal networks (2)
- stochastic models (2)
- 1/f noise (1)
- Anomalous transport (1)
- Brownian motors (1)
- Flashing ratchets (1)
- Generalized Langevin equation (1)
- Memory effects (1)
Institute
Viscoelastic subdiffusion governed by a fractional Langevin equation is studied numerically in a random Gaussian environment modeled by stationary Gaussian potentials with decaying spatial correlations. This anomalous diffusion is archetypal for living cells, where cytoplasm is known to be viscoelastic and a spatial disorder also naturally emerges. We obtain some first important insights into it within a model one-dimensional study. Two basic types of potential correlations are studied: short-range exponentially decaying and algebraically slow decaying with an infinite correlation length, both for a moderate (several kBT, in the units of thermal energy), and strong (5–10kBT) disorder. For a moderate disorder, it is shown that on the ensemble level viscoelastic subdiffusion can easily overcome the medium's disorder. Asymptotically, it is not distinguishable from the disorder-free subdiffusion. However, a strong scatter in single-trajectory averages is nevertheless seen even for a moderate disorder. It features a weak ergodicity breaking, which occurs on a very long yet transient time scale. Furthermore, for a strong disorder, a very long transient regime of logarithmic, Sinai-type diffusion emerges. It can last longer and be faster in the absolute terms for weakly decaying correlations as compared with the short-range correlations. Residence time distributions in a finite spatial domain are of a generalized log-normal type and are reminiscent also of a stretched exponential distribution. They can be easily confused for power-law distributions in view of the observed weak ergodicity breaking. This suggests a revision of some experimental data and their interpretation.
We study subdiffusive overdamped Brownian ratchets periodically rocked by an external zero-mean force in viscoelastic media within the framework of a non-Markovian generalized Langevin equation approach and associated multidimensional Markovian embedding dynamics. Viscoelastic deformations of the medium caused by the transport particle are modeled by a set of auxiliary Brownian quasiparticles elastically coupled to the transport particle and characterized by a hierarchy of relaxation times which obey a fractal scaling. The most slowly relaxing deformations which cannot immediately follow to the moving particle imprint long-range memory about its previous positions and cause subdiffusion and anomalous transport on a sufficiently long time scale. This anomalous behavior is combined with normal diffusion and transport on an initial time scale of overdamped motion. Anomalously slow directed transport in a periodic ratchet potential with broken space inversion symmetry emerges due to a violation of the thermal detailed balance by a zero-mean periodic driving and is optimized with frequency of driving, its amplitude, and temperature. Such optimized anomalous transport can be low dispersive and characterized by a large generalized Peclet number. Moreover, we show that overdamped subdiffusive ratchets can sustain a substantial load and do useful work. The corresponding thermodynamic efficiency decays algebraically in time since the useful work done against a load scales sublinearly with time following to the transport particle position, but the energy pumped by an external force scales with time linearly. Nevertheless, it can be transiently appreciably high and compare well with the thermodynamical efficiency of the normal diffusion overdamped ratchets on sufficiently long temporal and spatial scales.
We consider a simple Markovian class of the stochastic Wilson-Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around -1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.
Stochastic Wilson
(2015)
We consider a simple Markovian class of the stochastic Wilson–Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around −1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.
Magnetic nanoparticles are met across many biological species ranging from magnetosensitive bacteria, fishes, bees, bats, rats, birds, to humans. They can be both of biogenetic origin and due to environmental contamination, being either in paramagnetic or ferromagnetic state. The energy of such naturally occurring single-domain magnetic nanoparticles can reach up to 10-20 room k(B)T in the magnetic field of the Earth, which naturally led to supposition that they can serve as sensory elements in various animals. This work explores within a stochastic modeling framework a fascinating hypothesis of magnetosensitive ion channels with magnetic nanoparticles serving as sensory elements, especially, how realistic it is given a highly dissipative viscoelastic interior of living cells and typical sizes of nanoparticles possibly involved.
Can the statistical properties of single-electron transfer events be correctly predicted within a common equilibrium ensemble description? This fundamental in nanoworld question of ergodic behavior is scrutinized within a very basic semi-classical curve-crossing problem. It is shown that in the limit of non-adiabatic electron transfer (weak tunneling) well-described by the Marcus-Levich-Dogonadze (MLD) rate the answer is yes. However, in the limit of the so-called solvent-controlled adiabatic electron transfer, a profound breaking of ergodicity occurs. Namely, a common description based on the ensemble reduced density matrix with an initial equilibrium distribution of the reaction coordinate is not able to reproduce the statistics of single-trajectory events in this seemingly classical regime. For sufficiently large activation barriers, the ensemble survival probability in a state remains nearly exponential with the inverse rate given by the sum of the adiabatic curve crossing (Kramers) time and the inverse MLD rate. In contrast, near to the adiabatic regime, the single-electron survival probability is clearly non-exponential, even though it possesses an exponential tail which agrees well with the ensemble description. Initially, it is well described by a Mittag-Leffler distribution with a fractional rate. Paradoxically, the mean transfer time in this classical on the ensemble level regime is well described by the inverse of the nonadiabatic quantum tunneling rate on a single particle level. An analytical theory is developed which perfectly agrees with stochastic simulations and explains our findings.
Perfect anomalous transport of subdiffusive cargos by molecular motors in viscoelastic cytosol
(2019)
Multiple experiments show that various submicron particles such as magnetosomes, RNA messengers, viruses, and even much smaller nanoparticles such as globular proteins diffuse anomalously slow in viscoelastic cytosol of living cells. Hence, their sufficiently fast directional transport by molecular motors such as kinesins is crucial for the cell operation. It has been shown recently that the traditional flashing Brownian ratchet models of molecular motors are capable to describe both normal and anomalous transport of such subdiffusing cargos by molecular motors with a very high efficiency. This work elucidates further an important role of mechanochemical coupling in such an anomalous transport. It shows a natural emergence of a perfect subdiffusive ratchet regime due to allosteric effects, where the random rotations of a "catalytic wheel" at the heart of the motor operation become perfectly synchronized with the random stepping of a heavily loaded motor, so that only one ATP molecule is consumed on average at each motor step along microtubule. However, the number of rotations made by the catalytic engine and the traveling distance both scale sublinearly in time. Nevertheless, this anomalous transport can be very fast in absolute terms.
Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion
(2014)
The discovery of anomalous diffusion of larger biopolymers and submicron tracers such as endogenous granules, organelles, or virus capsids in living cells, attributed to the viscoelastic nature of the cytoplasm, provokes the question whether this complex environment equally impacts the active intracellular transport of submicron cargos by molecular motors such as kinesins: does the passive anomalous diffusion of free cargo always imply its anomalously slow active transport by motors, the mean transport distance along microtubule growing sublinearly rather than linearly in time? Here we analyze this question within the widely used two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubules driven by a flashing binding potential, where the cargo particle is elastically attached to the motor. Depending on the cargo size, the loading force, the amplitude of the binding potential, the turnover frequency of the molecular motor enzyme, and the linker stiffness we demonstrate that the motor transport may turn out either normal or anomalous, as indeed measured experimentally. We show how a highly efficient normal active transport mediated by motors may emerge despite the passive anomalous diffusion of the cargo, and study the intricate effects of the elastic linker. Under different, well specified conditions the microtubule-based motor transport becomes anomalously slow and thus significantly less efficient.
We propose and study a model of hypothetical magnetosensitive ionic channels which are long thought to be a possible candidate to explain the influence of weak magnetic fields on living organisms ranging from magnetotactic bacteria to fishes, birds, rats, bats, and other mammals including humans. The core of the model is provided by a short chain of magnetosomes serving as a sensor, which is coupled by elastic linkers to the gating elements of ion channels forming a small cluster in the cell membrane. The magnetic sensor is fixed by one end on cytoskeleton elements attached to the membrane and is exposed to viscoelastic cytosol. Its free end can reorient stochastically and subdiffusively in viscoelastic cytosol responding to external magnetic field changes and can open the gates of coupled ion channels. The sensor dynamics is generally bistable due to bistability of the gates which can be in two states with probabilities which depend on the sensor orientation. For realistic parameters, it is shown that this model channel can operate in the magnetic field of Earth for a small number (five to seven) of single-domain magnetosomes constituting the sensor rod, each of which has a typical size found in magnetotactic bacteria and other organisms or even just one sufficiently large nanoparticle of a characteristic size also found in nature. It is shown that, due to the viscoelasticity of the medium, the bistable gating dynamics generally exhibits power law and stretched exponential distributions of the residence times of the channels in their open and closed states. This provides a generic physical mechanism for the explanation of the origin of such anomalous kinetics for other ionic channels whose sensors move in a viscoelastic environment provided by either cytosol or biological membrane, in a quite general context, beyond the fascinating hypothesis of magnetosensitive ionic channels we explore.
Life and death of stationary linear response in anomalous continuous time random walk dynamics
(2014)
Linear theory of stationary response in systems at thermal equilibrium requires to find equilibrium correlation function of unperturbed responding system. Studies of the response of the systems exhibiting anomalously slow dynamics are often based on the continuous time random walk description (CTRW) with divergent mean waiting times. The bulk of the literature on anomalous response contains linear response functions like one by Cole-Cole calculated from such a CTRW theory and applied to systems at thermal equilibrium. Here we show within a fairly simple and general model that for the systems with divergent mean waiting times the stationary response at thermal equilibrium is absent, in accordance with some recent studies. The absence of such stationary response (or dying to zero non-stationary response in aging experiments) would confirm CTRW with divergent mean waiting times as underlying physical relaxation mechanism, but reject it otherwise. We show that the absence of stationary response is closely related to the breaking of ergodicity of the corresponding dynamical variable. As an important new result, we derive a generalized Cole-Cole response within ergodic CTRW dynamics with finite waiting time. Moreover, we provide a physically reasonable explanation of the origin and wide presence of 1/f noise in condensed matter for ergodic dynamics close to normal, rather than strongly deviating.