Generalized Kubo relations and conditions for anomalous diffusion: Physical insights from a mathematical theorem

A microrheological examination of insulin-secreting β-cells in healthy and diabetic-like conditions

Pancreatic β-cells regulate glucose homeostasis through glucose-stimulated insulin secretion, which is hindered in type-2 diabetes. Transport of the insulin vesicles is expected to be affected by changes in the viscoelastic and transport properties of the cytoplasm. These are evaluated in situ through particle-tracking measurements using a rat insulinoma β-cell line. The use of inert probes assists in decoupling the material properties of the cytoplasm from the active transport through cellular processes. The effect of glucose-stimulated insulin secretion is examined, and the subsequent remodeling of the cytoskeleton, at constant effects of cell activity, is shown to result in reduced mobility of the tracer particles. Induction of diabetic-like conditions is identified to alter the mean-squared displacement of the passive particles in the cytoplasm and diminish its reaction to glucose stimulation.

Can the microscopic and macroscopic transport phenomena in deep eutectic solvents be reconciled?

Deep eutectic solvents (DESs) have become ubiquitous in a variety of industrial and pharmaceutical applications since their discovery. However, the fundamental understanding of their physicochemical properties and their emergence from the microscopic features is still being explored fervently. Particularly, the knowledge of transport mechanisms in DESs is essential to tune their properties, which shall aid in expanding the territory of their applications. This perspective presents the current state of understanding of the bulk/macroscopic transport properties and microscopic relaxation processes in DESs. The dependence of these properties on the components and composition of the DES is explored, highlighting the role of hydrogen bonding (H-bonding) interactions. Modulation of these interactions by water and other additives, and their subsequent effect on the transport mechanisms, is also discussed. Various models (e.g. hole theory, free volume theory, etc.) have been proposed to explain the macroscopic transport phenomena from a microscopic origin. But the formation of H-bond networks and clusters in the DES reveals the insufficiency of these models, and establishes an antecedent for dynamic heterogeneity. Even significantly above the glass transition, the microscopic relaxation processes in DESs are rife with temporal and spatial heterogeneity, which causes a substantial decoupling between the viscosity and microscopic diffusion processes. However, we propose that a thorough understanding of the structural relaxation associated to the H-bond dynamics in DESs will provide the necessary framework to interpret the emergence of bulk transport properties from their microscopic counterparts.

Dioctadecyldimethylammonium bromide, a surfactant model for the cell membrane: Importance of microscopic dynamics

Cationic lipid membranes have recently attracted huge attention both from a fundamental point of view and due to their practical applications in drug delivery and gene therapy. The dynamical behavior of the lipids in the membrane is a key parameter controlling various physiological processes and drug release kinetics. Here, we review the dynamical and thermotropic phase behavior of an archetypal cationic lipid membrane, dioctadecyldimethylammonium bromide (DODAB), as studied using neutron scattering and molecular dynamics simulation techniques. DODAB membranes exhibit interesting phase behavior, specifically showing coagel, gel, and fluid phases in addition to a large hysteresis when comparing heating and cooling cycles. The dynamics of the lipid membrane is strongly dependent on the physical state of the bilayer. Lateral diffusion of the lipids is faster, by an order of magnitude, in the fluid phase than in the ordered phase. It is not only the characteristic times but also the nature of the segmental motions that differ between the ordered and fluid phases. The effect of different membrane active molecules including drugs, stimulants, gemini surfactants, and unsaturated lipids, on the dynamical and thermotropic phase behavior of the DODAB membrane, is also discussed here. Various interesting features such as induced synchronous ordering between polar head groups and tails, sub diffusive behavior, etc., are observed. The results shed light on the interaction between these additives and the membrane, which is found to be a complex interplay between the physical state of the membrane, charge, concentration, molecular architecture of the additives, and their location within the membrane.

Generalization of the Kubo relation for confined motion and ergodicity breakdown

A time-dependent generalized Kubo relation is derived by introducing the notion of a diffusion function for a particle confined in a harmonic potential. The relation reduces to the standard Kubo relation as a special case but holds for anomalous diffusion, nonergodic processes, and bounded motion. We analyze in detail the behaviors of the diffusion and memory functions and report a generalized Stokes-Einstein relation concerning anomalous diffusion. Furthermore, we demonstrate that when a high finite-frequency cutoff is imposed on the noise spectral density, a breakdown in ergodicity accompanied by the appearance of nonstationarity in the velocity autocorrelation function occurs in forced systems. This breakdown is taken as explicit evidence for either decay-spring-memory or recovering-force effects leading to nonexponential relaxation kinematics.

Fractal analysis of lateral movement in biomembranes

Lateral movement of a molecule in a biomembrane containing small compartments (0.23-μm diameter) and large ones (0.75 μm) is analyzed using a fractal description of its walk. The early time dependence of the mean square displacement varies from linear due to the contribution of ballistic motion. In small compartments, walking molecules do not have sufficient time or space to develop an asymptotic relation and the diffusion coefficient deduced from the experimental records is lower than that measured without restrictions. The model makes it possible to deduce the molecule step parameters, namely the step length and time, from data concerning confined and unrestricted diffusion coefficients. This is also possible using experimental results for sub-diffusive transport. The transition from normal to anomalous diffusion does not affect the molecule step parameters. The experimental literature data on molecular trajectories recorded at a high time resolution appear to confirm the modeled value of the mean free path length of DOPE for Brownian and anomalous diffusion. Although the step length and time give the proper values of diffusion coefficient, the DOPE speed calculated as their quotient is several orders of magnitude lower than the thermal speed. This is interpreted as a result of intermolecular interactions, as confirmed by lateral diffusion of other molecules in different membranes. The molecule step parameters are then utilized to analyze the problem of multiple visits in small compartments. The modeling of the diffusion exponent results in a smooth transition to normal diffusion on entering a large compartment, as observed in experiments.

Non-Poisson renewal events and memory

We study two different forms of fluctuation-dissipation processes generating anomalous relaxations to equilibrium of an initial out-of-equilibrium condition, the former being based on a stationary although very slow correlation function and the latter characterized by the occurrence of crucial events, namely, non-Poisson renewal events, incompatible with the stationary condition. Both forms of regression to equilibrium have the same nonexponential Mittag-Leffler structure. We analyze the single trajectories of the two processes by recording the time distances between two consecutive origin recrossings and establishing the corresponding waiting time probability density function (PDF), ψ(t). In the former case, with no crucial events, ψ(t) is an exponential, and in the latter case, with crucial events, ψ(t) is an inverse power law PDF with a diverging first moment. We discuss the consequences that this result is expected to have for the correct interpretation of some anomalous relaxation processes.

Lévy flights versus Lévy walks in bounded domains

Lévy flights and Lévy walks serve as two paradigms of random walks resembling common features but also bearing fundamental differences. One of the main dissimilarities is the discontinuity versus continuity of their trajectories and infinite versus finite propagation velocity. As a consequence, a well-developed theory of Lévy flights is associated with their pathological physical properties, which in turn are resolved by the concept of Lévy walks. Here, we explore Lévy flight and Lévy walk models on bounded domains, examining their differences and analogies. We investigate analytically and numerically whether and under which conditions both approaches yield similar results in terms of selected statistical observables characterizing the motion: the survival probability, mean first passage time, and stationary probability density functions. It is demonstrated that the similarity of the models is affected by the type of boundary conditions and the value of the stability index defining the asymptotics of the jump length distribution.

Asymptotic neutron scattering laws for anomalously diffusing quantum particles

The paper deals with a model-free approach to the analysis of quasielastic neutron scattering intensities from anomalously diffusing quantum particles. All quantities are inferred from the asymptotic form of their time-dependent mean square displacements which grow ∝t(α), with 0 ≤ α < 2. Confined diffusion (α = 0) is here explicitly included. We discuss in particular the intermediate scattering function for long times and the Fourier spectrum of the velocity autocorrelation function for small frequencies. Quantum effects enter in both cases through the general symmetry properties of quantum time correlation functions. It is shown that the fractional diffusion constant can be expressed by a Green-Kubo type relation involving the real part of the velocity autocorrelation function. The theory is exact in the diffusive regime and at moderate momentum transfers.

Ergodicity testing for anomalous diffusion: small sample statistics

The analysis of trajectories recorded in experiments often requires calculating time averages instead of ensemble averages. According to the Boltzmann hypothesis, they are equivalent only under the assumption of ergodicity. In this paper, we implement tools that allow to study ergodic properties. This analysis is conducted in two classes of anomalous diffusion processes: fractional Brownian motion and subordinated Ornstein-Uhlenbeck process. We show that only first of them is ergodic. We demonstrate this by applying rigorous statistical methods: mean square displacement, confidence intervals, and dynamical functional test. Our methodology is universal and can be implemented for analysis of many experimental data not only if a large sample is available but also when there are only few trajectories recorded.

Non-Brownian diffusion in lipid membranes: Experiments and simulations

The dynamics of constituents and the surface response of cellular membranes-also in connection to the binding of various particles and macromolecules to the membrane-are still a matter of controversy in the membrane biophysics community, particularly with respect to crowded membranes of living biological cells. We here put into perspective recent single particle tracking experiments in the plasma membranes of living cells and supercomputing studies of lipid bilayer model membranes with and without protein crowding. Special emphasis is put on the observation of anomalous, non-Brownian diffusion of both lipid molecules and proteins embedded in the lipid bilayer. While single component, pure lipid bilayers in simulations exhibit only transient anomalous diffusion of lipid molecules on nanosecond time scales, the persistence of anomalous diffusion becomes significantly longer ranged on the addition of disorder-through the addition of cholesterol or proteins-and on passing of the membrane lipids to the gel phase. Concurrently, experiments demonstrate the anomalous diffusion of membrane embedded proteins up to macroscopic time scales in the minute time range. Particular emphasis will be put on the physical character of the anomalous diffusion, in particular, the occurrence of ageing observed in the experiments-the effective diffusivity of the measured particles is a decreasing function of time. Moreover, we present results for the time dependent local scaling exponent of the mean squared displacement of the monitored particles. Recent results finding deviations from the commonly assumed Gaussian diffusion patterns in protein crowded membranes are reported. The properties of the displacement autocorrelation function of the lipid molecules are discussed in the light of their appropriate physical anomalous diffusion models, both for non-crowded and crowded membranes. In the last part of this review we address the upcoming field of membrane distortion by elongated membrane-binding particles. We discuss how membrane compartmentalisation and the particle-membrane binding energy may impact the dynamics and response of lipid membranes. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Communication: A scaling approach to anomalous diffusion

The paper presents a rigorous derivation of the velocity autocorrelation function for an anomalously diffusing slow solute particle in a bath of fast solvent molecules. The result is obtained within the framework of the generalized Langevin equation and uses only scaling arguments and identities which are based on asymptotic analysis. It agrees with the velocity autocorrelation function of an anomalously diffusing Rayleigh particle whose dynamics is described by a fractional Ornstein-Uhlenbeck process in velocity space. A simple semi-analytical example illustrates under which conditions the latter model is appropriate.

Fractal model of anomalous diffusion

An equation of motion is derived from fractal analysis of the Brownian particle trajectory in which the asymptotic fractal dimension of the trajectory has a required value. The formula makes it possible to calculate the time dependence of the mean square displacement for both short and long periods when the molecule diffuses anomalously. The anomalous diffusion which occurs after long periods is characterized by two variables, the transport coefficient and the anomalous diffusion exponent. An explicit formula is derived for the transport coefficient, which is related to the diffusion constant, as dependent on the Brownian step time, and the anomalous diffusion exponent. The model makes it possible to deduce anomalous diffusion properties from experimental data obtained even for short time periods and to estimate the transport coefficient in systems for which the diffusion behavior has been investigated. The results were confirmed for both sub and super-diffusion.

Temperature Dependence of Logarithmic-like Relaxational Dynamics of Hydrated tRNA

The dynamics of RNA within the β-relaxation region of 10 ps to 1 ns is crucial to its biological function. Because of its simpler chemical building blocks and the lack of the side methyl groups, faster relaxational dynamics of RNA compared to proteins can be expected. However, the situation is actually opposite. In this work, the relaxational dynamics of tRNA is measured by quasielastic neutron scattering and analyzed using the mode coupling theory, originally developed for glass-forming liquids. Our results reveal that the dynamics of tRNA follows a log-decay within the β-relaxation region, which is an important trait demonstrated by the dynamics of proteins. The dynamics of hydrated tRNA and lysozyme compared in the time domain further demonstrate that the slower dynamics of tRNA relative to proteins originates from the difference in the folded states of tRNA and proteins, as well as the influence of their hydration water.

Kinetic lattice Monte Carlo simulation of viscoelastic subdiffusion

We propose a kinetic Monte Carlo method for the simulation of subdiffusive random walks on a cartesian lattice. The random walkers are subject to viscoelastic forces which we compute from their individual trajectories via the fractional Langevin equation. At every step the walkers move by one lattice unit, which makes them differ essentially from continuous time random walks, where the subdiffusive behavior is induced by random waiting. To enable computationally inexpensive simulations with n-step memories, we use an approximation of the memory and the memory kernel functions with a complexity O(log n). Eventual discretization and approximation artifacts are compensated with numerical adjustments of the memory kernel functions. We verify with a number of analyses that this new method provides binary fractional random walks that are fully consistent with the theory of fractional brownian motion.

Individual degrees of freedom and the solvation properties of water

Using molecular dynamics simulations in conjunction with home-developed Split Integration Symplectic Method we effectively decouple individual degrees of freedom of water molecules and connect them to corresponding thermostats. In this way, we facilitate elucidation of structural, dynamical, spectral, and hydration properties of bulk water at any given combination of rotational, translational, and vibrational temperatures. Elevated rotational temperature of the water medium is found to severely hinder hydration of polar molecules, to affect hydration of ionic species in a nonmonotonous way and to somewhat improve hydration of nonpolar species. As proteins consist of charged, polar, and nonpolar amino-acid residues, the developed methodology is also applied to critically evaluate the hypothesis that the overall decrease in protein hydration and the change in the subtle balance between hydration of various types of amino-acid residues provide a plausible physical mechanism through which microwaves enhance aberrant protein folding and aggregation.