Non-Gaussian fluctuations resulting from power-law trapping in a lipid bilayer

Multiple diffusive behaviors of the random walk in inhomogeneous environments

Anomalous diffusive behaviors are observed in highly inhomogeneous but relatively stable environments such as intracellular media and are increasingly attracting attention. In this paper we develop a coupled continuous-time random walk model in which the waiting time is power-law coupled with the local environmental diffusion coefficient. We provide two forms of the waiting time density, namely, a heavy-tailed density and an exponential density. For different waiting time densities, anomalous diffusions with the diffusion exponent between 0 and 2 and Brownian yet non-Gaussian diffusion can be realized within the present model. The diffusive behaviors are analyzed and discussed by deriving the mean-squared displacement and probability density function. In addition we derive the effective jump length density corresponding to the decoupled form to help distinguish the diffusion types. Our model unifies two kinds of anomalous diffusive behavior with different characteristics in the same inhomogeneous environment into a theoretical framework. The model interprets the random motion of particles in a complex inhomogeneous environment and reproduces the experimental results of different biological and physical systems.

Statistics of the number of renewals, occupation times, and correlation in ordinary, equilibrium, and aging alternating renewal processes

The renewal process is a point process where an interevent time between successive renewals is an independent and identically distributed random variable. Alternating renewal process is a dichotomous process and a slight generalization of the renewal process, where the interevent time distribution alternates between two distributions. We investigate statistical properties of the number of renewals and occupation times for one of the two states in alternating renewal processes. When both means of the interevent times are finite, the alternating renewal process can reach an equilibrium. On the other hand, an alternating renewal process shows aging when one of the means diverges. We provide analytical calculations for the moments of the number of renewals, occupation time statistics, and the correlation function for several case studies in the interevent-time distributions. We show anomalous fluctuations for the number of renewals and occupation times when the second moment of interevent time diverges. When the mean interevent time diverges, distributional limit theorems for the number of events and occupation times are shown analytically. These are known as the Mittag-Leffler distribution and the generalized arcsine law in probability theory.

Heterogeneous biological membranes regulate protein partitioning via fluctuating diffusivity

Cell membranes phase separate into ordered L o and disordered L d domains depending on their compositions. This membrane compartmentalization is heterogeneous and regulates the localization of specific proteins related to cell signaling and trafficking. However, it is unclear how the heterogeneity of the membranes affects the diffusion and localization of proteins in L o and L d domains. Here, using Langevin dynamics simulations coupled with the phase-field (LDPF) method, we investigate several tens of milliseconds-scale diffusion and localization of proteins in heterogeneous biological membrane models showing phase separation into L o and L d domains. The diffusivity of proteins exhibits temporal fluctuations depending on the field composition. Increases in molecular concentrations and domain preference of the molecule induce subdiffusive behavior due to molecular collisions by crowding and confinement effects, respectively. Moreover, we quantitatively demonstrate that the protein partitioning into the L o domain is determined by the difference in molecular diffusivity between domains, molecular preference of domain, and molecular concentration. These results pave the way for understanding how biological reactions caused by molecular partitioning may be controlled in heterogeneous media. Moreover, the methodology proposed here is applicable not only to biological membrane systems but also to the study of diffusion and localization phenomena of molecules in various heterogeneous systems.

Occupation time statistics of the fractional Brownian motion in a finite domain

We study statistics of occupation times for a fractional Brownian motion (fBm), which is a typical model of a non-Markov process. Due to the non-Markovian nature, recurrence times to the origin depend on the history. Numerical simulations indicate that dependence on the sum of successive recurrence times becomes weak. As a result, the distribution of the occupation time in a finite domain follows the Mittag-Leffler distribution when the Hurst exponent of the fBm is close to 1/2. We show this distributional behavior of a time-averaged observable by renewal theory. This result is an extension of the distributional limit theorem known as the Darling-Kac theorem in general Markov processes to non-Markov processes.

Role of ergodicity, aging, and Gaussianity in resolving the origins of biomolecule subdiffusion

The internal motions of biomolecules are essential to their function. Although biological macromolecules conventionally show subdiffusive dynamics, only recently has subdiffusion been associated with non-ergodicity. These findings have stimulated new questions in biophysics and statistical mechanics. Is non-ergodic subdiffusion a general strategy shared by biomolecules? What underlying mechanisms are responsible for it? Here, we performed extensive molecular dynamics (MD) simulations to characterize the internal dynamics of six different biomolecules, ranging from single or double-stranded DNA, a single domain protein (KRAS), two globular proteins (PGK and SHP2), to an intrinsically disordered protein (SNAP-25). We found that the subdiffusive behavior of these biomolecules falls into two classes. The internal motion of the first three cases is ergodic subdiffusion and can be interpreted by fractional Brownian motion (FBM), while the latter three cases involve non-ergodic subdiffusion and can be modeled by mixed origins of continuous-time random walk (CTRW) and FBM.

Rheology of sliding leaflets in coarse-grained DSPC lipid bilayers

Amphiphilic lipid bilayers modify the friction properties of the surfaces on top of which they are deposited. In particular, the measured sliding friction coefficient can be significantly reduced compared with the native surface. We investigate in this work the friction properties of a numerical coarse-grained model of DSPC (1,2-distearoyl-sn-glycero-3-phosphocholine) lipid bilayer subject to longitudinal shear. The interleaflet friction coefficient is obtained from out-of-equilibrium pulling or from relaxation simulations. In particular, we gain access to the transient viscoelastic response of a sheared bilayer. The bilayer mechanical response is found to depend significantly on the membrane physical state, with evidence in favor of a linear response regime in the fluid but not in the gel region. The linear response validity domain is established, and the timescales appearing in the membrane response discussed.

Interdomain exchange and the flip-flop of cholesterol in ternary component lipid membranes and their effects on heterogeneous cholesterol diffusion

Cell membranes are heterogeneous with a variety of lipids, cholesterol, and proteins and are composed of domains of different compositions. Such heterogeneous environments make the transport of cholesterol complicated: cholesterol not only diffuses within a particular domain but also travels between domains. Cholesterol also flip-flops between upper and lower leaflets such that cholesterol may reside both within leaflets and in the central region between two leaflets. How the presence of multiple domains and the interdomain exchange of cholesterol would affect the cholesterol transport, however, remains elusive. In this study, therefore, we perform molecular dynamics simulations up to 100μs for ternary component lipid membranes, which consist of saturated lipids (dipalmitoylphosphatidylcholine, DPPC), unsaturated lipids (dilinoleylphosphatidylcholine, DIPC), and cholesterol. The ternary component membranes in our simulations form two domains readily: DPPC and DIPC domains. We find that the diffusion of cholesterol molecules is much more heterogeneous and non-Gaussian than expected for binary component lipid membranes of lipids and cholesterol. The non-Gaussian parameter of the cholesterol molecules is about four times larger in the ternary component lipid membranes than in the binary component lipid membranes. Such non-Gaussian and heterogeneous transport of cholesterol arises from the interplay among the interdomain kinetics, the different diffusivity of cholesterol in different domains, and the flip-flop of cholesterol. This suggests that in cell membranes that consist of various domains and proteins, the cholesterol transport can be very heterogeneous. We also find that the mechanism of the interdomain exchange differs for different domains: cholesterol tends to exit the DIPC domain along the central region of the membrane for the DIPC-to-DPPC transition, while the cholesterol is likely to exit the DPPC domain within the membrane leaflet for the DPPC-to-DIPC transition. Also, the interdomain exchange kinetics of cholesterol for the DPPC-to-DIPC transition is up to 7.9 times slower than the DIPC-to-DPPC transition.

Subdiffusive-Brownian crossover in membrane proteins: a generalized Langevin equation-based approach

In this work, we propose a generalized Langevin equation-based model to describe the lateral diffusion of a protein in a lipid bilayer. The memory kernel is represented in terms of a viscous (instantaneous) and an elastic (noninstantaneous) component modeled through a Dirac δ function and a three-parameter Mittag-Leffler type function, respectively. By imposing a specific relationship between the parameters of the three-parameter Mittag-Leffler function, the different dynamical regimes-namely ballistic, subdiffusive, and Brownian, as well as the crossover from one regime to another-are retrieved. Within this approach, the transition time from the ballistic to the subdiffusive regime and the spectrum of relaxation times underlying the transition from the subdiffusive to the Brownian regime are given. The reliability of the model is tested by comparing the mean-square displacement derived in the framework of this model and the mean-square displacement of a protein diffusing in a membrane calculated through molecular dynamics simulations.

Tracking of Nanoparticle Diffusion at a Liquid-Liquid Interface Adsorbed by Nonionic Surfactants

Emulsions stabilized by both nanoparticles and surfactants often display longer shelf life than those stabilized by nanoparticles or surfactants alone. Although numerous works have been conducted to understand the effect of nanoparticles and surfactants on the variation of interfacial tension, little is known about interfacial diffusion when both nanoparticles and surfactants are present at interfaces. In this work, we used single-particle fluorescence tracking to study the lateral diffusion of individual hydrophobic nanoparticles at hexane-glycerol interfaces adsorbed by different amounts of nonionic surfactants. When the surfactant concentration is over a threshold, we found that the nanoparticle diffusion exhibits a two-regime behavior involving short-time Brownian and the emergence of subdiffusive, non-Gaussian, and dynamically anticorrelated diffusion in the long lag time regime. A stepwise analysis rationalized diffusion in different lag time regimes, leading to a mechanistic interpretation regarding the two-regime behavior. These results could provide insight into the understanding of the synergistic effect for the surfactant-assistant Pickering emulsion.

The effects of the lipid type on the spatial arrangement and dynamics of cholesterol in binary component lipid membranes

Intermolecular interactions between cholesterol and lipids in cell membranes, which play critical roles in cellular processes such as the formation of nano-domains, depend on the molecular structure of the lipids. The diffusion and the spatial arrangement of cholesterol within the lipid membranes also change with the type of lipids. For example, the flip-flop, an important transport mechanism for cholesterol in the membranes, can be facilitated significantly by the presence of unsaturated lipids. However, how the structure of lipids affects the spatial arrangement and the dynamics of cholesterol remains elusive at a molecular level. In this study, we investigate the effects of lipid-cholesterol interactions on the spatial arrangement and the dynamics of cholesterol. We perform molecular dynamics simulations for the binary component membranes of lipids and cholesterol. We employ seven different kinds of lipids by changing either the degree of a saturation level or the length of lipid tails. We find from our simulations that the rate of cholesterol flip-flop is enhanced as the lipids are either less saturated or shorter, which is consistent with previous studies. Interestingly, when the lipid tails are fully saturated and sufficiently long, the center in between two leaflets becomes metastable for cholesterol to stay at. Because the cholesterol at the membrane center diffuses faster than that within leaflets, regardless of the lipid type, such an emergence of the metastable state (in terms of the cholesterol position) complicates the cholesterol diffusion significantly.

Anomalous nanoparticle surface diffusion in LCTEM is revealed by deep learning-assisted analysis

The motion of nanoparticles near surfaces is of fundamental importance in physics, biology, and chemistry. Liquid cell transmission electron microscopy (LCTEM) is a promising technique for studying motion of nanoparticles with high spatial resolution. Yet, the lack of understanding of how the electron beam of the microscope affects the particle motion has held back advancement in using LCTEM for in situ single nanoparticle and macromolecule tracking at interfaces. Here, we experimentally studied the motion of a model system of gold nanoparticles dispersed in water and moving adjacent to the silicon nitride membrane of a commercial LC in a broad range of electron beam dose rates. We find that the nanoparticles exhibit anomalous diffusive behavior modulated by the electron beam dose rate. We characterized the anomalous diffusion of nanoparticles in LCTEM using a convolutional deep neural-network model and canonical statistical tests. The results demonstrate that the nanoparticle motion is governed by fractional Brownian motion at low dose rates, resembling diffusion in a viscoelastic medium, and continuous-time random walk at high dose rates, resembling diffusion on an energy landscape with pinning sites. Both behaviors can be explained by the presence of silanol molecular species on the surface of the silicon nitride membrane and the ionic species in solution formed by radiolysis of water in presence of the electron beam.

Reaction-diffusion and reaction-subdiffusion equations on arbitrarily evolving domains

Reaction-diffusion equations are widely used as the governing evolution equations for modeling many physical, chemical, and biological processes. Here we derive reaction-diffusion equations to model transport with reactions on a one-dimensional domain that is evolving. The model equations, which have been derived from generalized continuous time random walks, can incorporate complexities such as subdiffusive transport and inhomogeneous domain stretching and shrinking. Inhomogeneously growing domains are frequently encountered in biological phenomena involving stochastic transport, such as tumor growth and morphogen gradient formation. A method for constructing analytic expressions for short-time moments of the position of the particles is developed and moments calculated from this approach are shown to compare favorably with results from random walk simulations and numerical integration of the reaction transport equation. The results show the important role played by the initial condition. In particular, it strongly affects the time dependence of the moments in the short-time regime by introducing additional drift and diffusion terms. We also discuss how our reaction transport equation could be applied to study the spreading of a population on an evolving interface. From a more general perspective, our findings help to mitigate the scarcity of analytic results for reaction-diffusion problems in geometries displaying nonuniform growth. They are also expected to pave the way for further results, including the treatment of first-passage problems associated with encounter-controlled reactions in such domains.

Influence of External NaCl Salt on Membrane Rigidity of Neutral DOPC Vesicles

Sodium chloride (NaCl) is a very common molecule in biotic and abiotic aqueous environments. In both cases, variation of ionic strength is inevitable. In addition to the osmotic variation posed by such perturbations, the question of whether the interactions of monovalent ions Na⁺ and Cl^-, especially with the neutral head groups of phospholipid membranes are impactful enough to change the membrane rigidity, is still not entirely understood. We investigated the dynamics of 1,2-di-(octadecenoyl)-sn-glycero-3-phosphocholine (DOPC) vesicles with zwitterionic neutral head groups in the fluid phase with increasing external salt concentration. At higher salt concentrations, we observe an increase in bending rigidity from neutron spin echo (NSE) spectroscopy and an increase in bilayer thickness from small-angle X-ray scattering (SAXS). We compared different models to distinguish membrane undulations, lipid tail motions, and the translational diffusion of the vesicles. All of the models indicate an increase in bending rigidity by a factor of 1.3-3.6. We demonstrate that even down to t > 10 ns and for Q > 0.07 Å^-1, the observed NSE relaxation spectra are influenced by translational diffusion of the vesicles. For t < 5 ns, the lipid tail motion dominates the intermediate dynamic structure factor. As the salt concentration increases, this contribution diminishes. We introduced a time-dependent analysis for the bending rigidity that highlights only a limited Zilman-Granek time window in which the rigidity is physically meaningful.

Diffusive transport of nanoscale objects through cell membranes: a computational perspective

Diffusion is an essential and fundamental means of transport of substances on cell membranes, and the dynamics of biomembranes plays a crucial role in the regulation of numerous cellular processes. The understanding of the complex mechanisms and the nature of particle diffusion have a bearing on establishing guidelines for the design of efficient transport materials and unique therapeutic approaches. Herein, this review article highlights the most recent advances in investigating diffusion dynamics of nanoscale objects on biological membranes, focusing on the approaches of tailored computer simulations and theoretical analysis. Due to the presence of the complicated and heterogeneous environment on native cell membranes, the diffusive transport behaviors of nanoparticles exhibit unique and variable characteristics. The general aspects and basic theories of normal diffusion and anomalous diffusion have been introduced. In addition, the influence of a series of external and internal factors on the diffusion behaviors is discussed, including particle size, membrane curvature, particle-membrane interactions or particle-inclusion, and the crowding degree of membranes. Finally, we seek to identify open problems in the existing experimental, simulation, and theoretical research studies, and to propose challenges for future development.

Modeling the dynamics of phospholipids in the fluid phase of liposomes

We present the derivation of a new model to describe neutron spin echo spectroscopy and quasi-elastic neutron scattering data on liposomes. We compare the new model with existing approaches and benchmark it with experimental data. The analysis indicates the importance of including all major contributions in the modeling of the intermediate scattering function. Simultaneous analysis of the experimental data on lipids with full contrast and tail contrast matched samples reveals highly confined lipid tail motion. A comparison of their dynamics demonstrates the statistical independence of tail-motion and height-height correlation of the membrane. A more detailed analysis indicates that the lipid tail relaxation is confined to a potential with cylindrical symmetry, in addition to the undulation and diffusive motion of the liposome. Despite substantial differences in the chemistry of the fatty acid tails, the observation indicates a universal behavior. The analysis of partially deuterated systems confirms the strong contribution of the lipid tail to the intermediate scattering function. Within the time range from 5 to 100 ns, the intermediate scattering function can be described by the height-height correlation function. The existence of the fast-localized tail motion and the contribution of slow translational diffusion of liposomes determine the intermediate scattering function for t < 5 ns and t > 100 ns, respectively. Taking into account the limited time window lowers the bending moduli by a factor of 1.3 (DOPC) to 2 (DMPC) compared to the full range.

Confinement and activity regulate bacterial motion in porous media

Understanding how bacteria move in porous media is critical to applications in healthcare, agriculture, environmental remediation, and chemical sensing. Recent work has demonstrated that E. coli, which moves by run-and-tumble dynamics in a homogeneous medium, exhibits a new form of motility when confined in a disordered porous medium: hopping-and-trapping motility, in which cells perform rapid, directed hops punctuated by intervals of slow, undirected trapping. Here, we use direct visualization to shed light on how these processes depend on pore-scale confinement and cellular activity. We find that hopping is determined by pore-scale confinement, and is independent of cellular activity; by contrast, trapping is determined by the competition between pore-scale confinement and cellular activity, as predicted by an entropic trapping model. These results thus help to elucidate the factors that regulate bacterial motion in porous media, and could help aid the development of new models of motility in heterogeneous environments.

Bacterial hopping and trapping in porous media

Diverse processes-e.g. bioremediation, biofertilization, and microbial drug delivery-rely on bacterial migration in disordered, three-dimensional (3D) porous media. However, how pore-scale confinement alters bacterial motility is unknown due to the opacity of typical 3D media. As a result, models of migration are limited and often employ ad hoc assumptions. Here we reveal that the paradigm of run-and-tumble motility is dramatically altered in a porous medium. By directly visualizing individual Escherichia coli, we find that the cells are intermittently and transiently trapped as they navigate the pore space, exhibiting diffusive behavior at long time scales. The trapping durations and the lengths of "hops" between traps are broadly distributed, reminiscent of transport in diverse other disordered systems; nevertheless, we show that these quantities can together predict the long-time bacterial translational diffusivity. Our work thus provides a revised picture of bacterial motility in complex media and yields principles for predicting cellular migration.

Facilitated and Non-Gaussian Diffusion of Cholesterol in Liquid Ordered Phase Bilayers Depends on the Flip-Flop and Spatial Arrangement of Cholesterol

The diffusion of cholesterol in biological membranes is critical to cellular processes such as the formation of cholesterol-enriched domains. The cholesterol diffusion may be complicated especially when cholesterol flip-flops and/or stays at the membrane center. Understanding the diffusion mechanism of cholesterol at a molecular level should be, therefore, a topic of interest. We perform molecular dynamics simulations up to 100 μs for lipid bilayers with various concentrations of cholesterol. We find that cholesterol diffusion in the liquid ordered phase depends on whether it is within leaflets or at the bilayer center, is non-Gaussian for several microseconds, and is enhanced significantly compared to that of lipids. Cholesterol at the bilayer center diffuses fast, while cholesterol in the hydrocarbon region with upright orientation diffuses relatively slowly. Such position-dependent dynamics of cholesterol leads to facilitated and non-Gaussian diffusion.

Biased continuous-time random walks for ordinary and equilibrium cases: facilitation of diffusion, ergodicity breaking and ageing

We examine renewal processes with power-law waiting time distributions (WTDs) and non-zero drift via computing analytically and by computer simulations their ensemble and time averaged spreading characteristics. All possible values of the scaling exponent α are considered for the WTD ψ(t) ∼ 1/t1+α. We treat continuous-time random walks (CTRWs) with 0 < α < 1 for which the mean waiting time diverges, and investigate the behaviour of the process for both ordinary and equilibrium CTRWs for 1 < α < 2 and α > 2. We demonstrate that in the presence of a drift CTRWs with α < 1 are ageing and non-ergodic in the sense of the non-equivalence of their ensemble and time averaged displacement characteristics in the limit of lag times much shorter than the trajectory length. In the sense of the equivalence of ensemble and time averages, CTRW processes with 1 < α < 2 are ergodic for the equilibrium and non-ergodic for the ordinary situation. Lastly, CTRW renewal processes with α > 2-both for the equilibrium and ordinary situation-are always ergodic. For the situations 1 < α < 2 and α > 2 the variance of the diffusion process, however, depends on the initial ensemble. For biased CTRWs with α > 1 we also investigate the behaviour of the ergodicity breaking parameter. In addition, we demonstrate that for biased CTRWs the Einstein relation is valid on the level of the ensemble and time averaged displacements, in the entire range of the WTD exponent α.

Ergodicity, rejuvenation, enhancement, and slow relaxation of diffusion in biased continuous-time random walks

Bias plays an important role in the enhancement of diffusion in periodic potentials. Using the continuous-time random walk in the presence of a bias, we report on an interesting phenomenon for the enhancement of diffusion by the start of the measurement in a random energy landscape. When the variance of the waiting time diverges, in contrast to the bias-free case, the dynamics with bias becomes superdiffusive. In the superdiffusive regime, we find a distinct initial ensemble dependence of the diffusivity. Moreover, the diffusivity can be increased by the aging time when the initial ensemble is not in equilibrium. We show that the time-averaged variance converges to the corresponding ensemble-averaged variance; i.e., ergodicity is preserved. However, trajectory-to-trajectory fluctuations of the time-averaged variance decay unexpectedly slowly. Our findings provide a rejuvenation phenomenon in the superdiffusive regime, that is, the diffusivity for a nonequilibrium initial ensemble gradually increases to that for an equilibrium ensemble when the start of the measurement is delayed.

Impact of the green tea ingredient epigallocatechin gallate and a short pentapeptide (Ile-Ile-Ala-Glu-Lys) on the structural organization of mixed micelles and the related uptake of cholesterol

BACKGROUND

High levels of blood cholesterol are conventionally linked to an increased risk of developing cardiovascular disease (Grundy, 1986). Here we examine the molecular mode of action of natural products with known cholesterol-lowering activity, such as for example the green tea ingredient epigallocatechin gallate and a short pentapeptide, Ile-Ile-Ala-Glu-Lys.

METHODS

Molecular Dynamics simulations are used to gain insight into the formation process of mixed micelles and, correspondingly, how active agents epigallocatechin gallate and Ile-Ile-Ala-Glu-Lys could possibly interfere with it.

RESULTS

Self-assembly of physiological micelles occurs on the order of 35-50 ns; most of the structural properties of mixed micelles are unaffected by epigallocatechin gallate or Ile-Ile-Ala-Glu-Lys which integrate into the micellar surface; the diffusive motion of constituting lipids palmitoyl-oleoyl-phosphatidylcholine and cholesterol is significantly down-regulated by both epigallocatechin gallate and Ile-Ile-Ala-Glu-Lys; CONCLUSIONS: The molecular mode of action of natural compounds epigallocatechin gallate and Ile-Ile-Ala-Glu-Lys is a significant down-regulation of the diffusive motion of micellar lipids.

GENERAL SIGNIFICANCE

Natural compounds like the green tea ingredient epigallocatechin gallate and a short pentapeptide, Ile-Ile-Ala-Glu-Lys, lead to a significant down-regulation of the diffusive motion of micellar lipids thereby modulating cholesterol absorption into physiological micelles.

Dynamics of Phospholipid Membranes beyond Thermal Undulations

We investigated the molecular dynamics of unilamellar liposomes by neutron spin echo spectroscopy. We report the first experimental evidence of a short-range motion at the length scale of the size of the headgroup of a lipid. The associated mean squared displacement shows a t^0.26 dependence in the pico- to nanosecond region that indicates another process beyond the predictions of the Zilman-Granek (ZG) model ( t^0.66) and translational diffusion ( t¹). A comparison with theory shows that the observed low exponent is associated with a non-Gaussian transient trapping of lipid molecules in a local area and supports the continuous time random walk model. The analysis of the mean squared displacement leads to the important conclusion that the friction at the interface between water and liposomes plays a minor role. Center of mass diffusion of liposomes and transient trapping of lipids define the range in which the ZG model can be applied to analyze membrane fluctuations.

Phase behaviors of deeply supercooled bilayer water unseen in bulk water

Akin to bulk water, water confined to an isolated nanoslit can show a wealth of new 2D phases of ice and amorphous ice, as well as unusual phase behavior. Indeed, 2D water phases, such as bilayer hexagonal ice and monolayer square ice, have been detected in the laboratory, confirming earlier computational predictions. Herein, we report theoretical evidence of a hitherto unreported state, namely, bilayer very low density amorphous ice (BL-VLDA), as well as evidence of a strong first-order transition between BL-VLDA and the BL amorphous ice (BL-A), and a weak first-order transition between BL-VLDA and the BL very low density liquid (BL-VLDL) water. The diffusivity of BL-VLDA is typically in the range of 10-9 cm2/s to 10-10 cm2/s. Similar to bulk (3D) water, 2D water can exhibit two forms of liquid in the deeply supercooled state. However, unlike supercooled bulk water, for which the two forms of liquid can coexist and merge into one at a critical point, the 2D BL-VLDL and BL high-density liquid (BL-HDL) phases are separated by the highly stable solid phase of BL-A whose melting line exhibits the isochore end point (IEP) near 220 K in the temperature-pressure diagram. Above the IEP temperature, BL-VLDL and BL-HDL are indistinguishable. At negative pressures, the metastable BL-VLDL exhibits a spatially and temporally heterogeneous structure induced by dynamic changes in the nanodomains, a feature much less pronounced in the BL-HDL.

Non-self-averaging behaviors and ergodicity in quenched trap models with finite system sizes

Tracking tracer particles in heterogeneous environments plays an important role in unraveling material properties. These heterogeneous structures are often static and depend on the sample realizations. Sample-to-sample fluctuations of such disorder realizations sometimes become considerably large. When we investigate the sample-to-sample fluctuations, fundamental averaging procedures are a thermal average for a single disorder realization and the disorder average for different disorder realizations. Here we report on non-self-averaging phenomena in quenched trap models with finite system sizes, where we consider the periodic and the reflecting boundary conditions. Sample-to-sample fluctuations of diffusivity greatly exceed trajectory-to-trajectory fluctuations of diffusivity in the corresponding annealed model. For a single disorder realization, the time-averaged mean square displacement and position-dependent observables converge to constants because of the existence of the equilibrium distribution. This is a manifestation of ergodicity. As a result, the time-averaged quantities depend neither on the initial condition nor on the thermal histories but depend crucially on the disorder realization.

Non-Gaussianity, population heterogeneity, and transient superdiffusion in the spreading dynamics of amoeboid cells

What is the underlying diffusion process governing the spreading dynamics and search strategies employed by amoeboid cells?

Detection of transition times from single-particle-tracking trajectories

In heterogeneous environments, the diffusivity is not constant but changes with time. It is important to detect changes in the diffusivity from single-particle-tracking trajectories in experiments. Here, we devise a novel method for detecting the transition times of the diffusivity from trajectory data. A key idea of this method is the introduction of a characteristic time scale of the diffusive states, which is obtained by a fluctuation analysis of the time-averaged mean square displacements. We test our method in silico by using the Langevin equation with a fluctuating diffusivity. We show that our method can successfully detect the transition times of diffusive states and obtain the diffusion coefficient as a function of time. This method will provide a quantitative description of the fluctuating diffusivity in heterogeneous environments and can be applied to time series with transitions of states.

Computational and theoretical approaches for studies of a lipid recognition protein on biological membranes

Many cellular functions, including cell signaling and related events, are regulated by the association of peripheral membrane proteins (PMPs) with biological membranes containing anionic lipids, e.g., phosphatidylinositol phosphate (PIP). This association is often mediated by lipid recognition modules present in many PMPs. Here, I summarize computational and theoretical approaches to investigate the molecular details of the interactions and dynamics of a lipid recognition module, the pleckstrin homology (PH) domain, on biological membranes. Multiscale molecular dynamics simulations using combinations of atomistic and coarse-grained models yielded results comparable to those of actual experiments and could be used to elucidate the molecular mechanisms of the formation of protein/lipid complexes on membrane surfaces, which are often difficult to obtain using experimental techniques. Simulations revealed some modes of membrane localization and interactions of PH domains with membranes in addition to the canonical binding mode. In the last part of this review, I address the dynamics of PH domains on the membrane surface. Local PIP clusters formed around the proteins exhibit anomalous fluctuations. This dynamic change in protein-lipid interactions cause temporally fluctuating diffusivity of proteins, i.e., the short-term diffusivity of the bound protein changes substantially with time, and may in turn contribute to the formation/dissolution of protein complexes in membranes.

Non-Markovian dynamics of reaction coordinate in polymer folding

We develop a theoretical description of the critical zipping dynamics of a self-folding polymer. We use tension propagation theory and the formalism of the generalized Langevin equation applied to a polymer that contains two complementary parts which can bind to each other. At the critical temperature, the (un)zipping is unbiased and the two strands open and close as a zipper. The number of broken base pairs n(t) displays a subdiffusive motion characterized by a variance growing as 〈Δn²(t)〉 ∼ t^α with α < 1 at long times. Our theory provides an estimate of both the asymptotic anomalous exponent α and of the subleading correction term, which are both in excellent agreement with numerical simulations. The results indicate that the tension propagation theory captures the relevant features of the dynamics and shed some new insights on related polymer problems characterized by anomalous dynamical behavior.

Complementary mode analyses between sub- and superdiffusion

Several subdiffusive stochastic processes in nature, e.g., the motion of a tagged monomer in polymers, the height fluctuation of interfaces, particle dynamics in single-file diffusion, etc., can be described rigorously or approximately by the superposition of various modes whose relaxation times are broadly distributed. In this paper, we propose a mode analysis generating superdiffusion, which is paired with or complementary to subdiffusion. The key point in our discussion lies in the identification of a pair of conjugated variables, which undergo sub- and superdiffusion, respectively. We provide a simple interpretation for the sub- and superdiffusion duality for these variables using the language of polymer physics. The analysis also suggests the usefulness of looking at the force fluctuation in experiments, where a polymer is driven by a constant velocity.

Manipulation and Motion of Organelles and Single Molecules in Living Cells

The biomolecule is among the most important building blocks of biological systems, and a full understanding of its function forms the scaffold for describing the mechanisms of higher order structures as organelles and cells. Force is a fundamental regulatory mechanism of biomolecular interactions driving many cellular processes. The forces on a molecular scale are exactly in the range that can be manipulated and probed with single molecule force spectroscopy. The natural environment of a biomolecule is inside a living cell, hence, this is the most relevant environment for probing their function. In vivo studies are, however, challenged by the complexity of the cell. In this review, we start with presenting relevant theoretical tools for analyzing single molecule data obtained in intracellular environments followed by a description of state-of-the art visualization techniques. The most commonly used force spectroscopy techniques, namely optical tweezers, magnetic tweezers, and atomic force microscopy, are described in detail, and their strength and limitations related to in vivo experiments are discussed. Finally, recent exciting discoveries within the field of in vivo manipulation and dynamics of single molecule and organelles are reviewed.

Dynamic interactions between a membrane binding protein and lipids induce fluctuating diffusivity

Pleckstrin homology (PH) domains are membrane-binding lipid recognition proteins that interact with phosphatidylinositol phosphate (PIP) molecules in eukaryotic cell membranes. Diffusion of PH domains plays a critical role in biological reactions on membrane surfaces. Although diffusivity can be estimated by long-time measurements, it lacks information on the short-time diffusive nature. We reveal two diffusive properties of a PH domain bound to the surface of a PIP-containing membrane using molecular dynamics simulations. One is fractional Brownian motion, attributed to the motion of the lipids with which the PH domain interacts. The other is temporally fluctuating diffusivity; that is, the short-time diffusivity of the bound protein changes substantially with time. Moreover, the diffusivity for short-time measurements is intrinsically different from that for long-time measurements. This fluctuating diffusivity results from dynamic changes in interactions between the PH domain and PIP molecules. Our results provide evidence that the complexity of protein-lipid interactions plays a crucial role in the diffusion of proteins on biological membrane surfaces. Changes in the diffusivity of PH domains and related membrane-bound proteins may in turn contribute to the formation/dissolution of protein complexes in membranes.

Langevin equation with fluctuating diffusivity: A two-state model

Recently, anomalous subdiffusion, aging, and scatter of the diffusion coefficient have been reported in many single-particle-tracking experiments, though the origins of these behaviors are still elusive. Here, as a model to describe such phenomena, we investigate a Langevin equation with diffusivity fluctuating between a fast and a slow state. Namely, the diffusivity follows a dichotomous stochastic process. We assume that the sojourn time distributions of these two states are given by power laws. It is shown that, for a nonequilibrium ensemble, the ensemble-averaged mean-square displacement (MSD) shows transient subdiffusion. In contrast, the time-averaged MSD shows normal diffusion, but an effective diffusion coefficient transiently shows aging behavior. The propagator is non-Gaussian for short time and converges to a Gaussian distribution in a long-time limit; this convergence to Gaussian is extremely slow for some parameter values. For equilibrium ensembles, both ensemble-averaged and time-averaged MSDs show only normal diffusion and thus we cannot detect any traces of the fluctuating diffusivity with these MSDs. Therefore, as an alternative approach to characterizing the fluctuating diffusivity, the relative standard deviation (RSD) of the time-averaged MSD is utilized and it is shown that the RSD exhibits slow relaxation as a signature of the long-time correlation in the fluctuating diffusivity. Furthermore, it is shown that the RSD is related to a non-Gaussian parameter of the propagator. To obtain these theoretical results, we develop a two-state renewal theory as an analytical tool.

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.

Anomalous diffusion in time-fluctuating non-stationary diffusivity landscapes

We investigate the diffusive and ergodic properties of massive and confined particles in a model disordered medium, in which the local diffusivity fluctuates in time while its mean has a power law dependence on the diffusion time.

Anomalous Dynamics of a Lipid Recognition Protein on a Membrane Surface

Pleckstrin homology (PH) domains are lipid-binding modules present in peripheral membrane proteins which interact with phosphatidyl-inositol phosphates (PIPs) in cell membranes. We use multiscale molecular dynamics simulations to characterize the localization and anomalous dynamics of the DAPP1 PH domain on the surface of a PIP-containing lipid bilayer. Both translational and rotational diffusion of the PH domain on the lipid membrane surface exhibit transient subdiffusion, with an exponent α ≈ 0.5 for times of less than 10 ns. In addition to a PIP₃ molecule at the canonical binding site of the PH domain, we observe additional PIP molecules in contact with the protein. Fluctuations in the number of PIPs associated with the PH domain exhibit 1/f noise. We suggest that the anomalous diffusion and long-term correlated interaction of the PH domain with the membrane may contribute to an enhanced probability of encounter with target complexes on cell membrane surfaces.

Geometry controlled anomalous diffusion in random fractal geometries: looking beyond the infinite cluster

We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law ∼T(-h) with h < 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided.

Fluctuation analysis of time-averaged mean-square displacement for the Langevin equation with time-dependent and fluctuating diffusivity

The mean-square displacement (MSD) is widely utilized to study the dynamical properties of stochastic processes. The time-averaged MSD (TAMSD) provides some information on the dynamics which cannot be extracted from the ensemble-averaged MSD. In particular, the relative standard deviation (RSD) of the TAMSD can be utilized to study the long-time relaxation behavior. In this work, we consider a class of Langevin equations which are multiplicatively coupled to time-dependent and fluctuating diffusivities. Various interesting dynamics models such as entangled polymers and supercooled liquids can be interpreted as the Langevin equations with time-dependent and fluctuating diffusivities. We derive a general formula for the RSD of the TAMSD for the Langevin equation with the time-dependent and fluctuating diffusivity. We show that the RSD can be expressed in terms of the correlation function of the diffusivity. The RSD exhibits the crossover at the long time region. The crossover time is related to a weighted average relaxation time for the diffusivity. Thus the crossover time gives some information on the relaxation time of fluctuating diffusivity which cannot be extracted from the ensemble-averaged MSD. We discuss the universality and possible applications of the formula via some simple examples.

Transition from distributional to ergodic behavior in an inhomogeneous diffusion process: Method revealing an unknown surface diffusivity

Diffusion of molecules in cells plays an important role in providing a biological reaction on the surface by finding a target on the membrane surface. The water retardation (slow diffusion) near the target assists the searching molecules to recognize the target. Here, we consider effects of the surface diffusivity on the effective diffusivity, where diffusion on the surface is slower than that in bulk. We show that the ensemble-averaged mean-square displacements increase linearly with time when the desorption rate from the surface is finite, which is valid even when the diffusion on the surface is anomalous (subdiffusion). Moreover, this slow diffusion on the surface affects the fluctuations of the time-averaged mean-square displacements (TAMSDs). We find that fluctuations of the TAMSDs remain large when the measurement time is smaller than a characteristic time, and decays according to an increase of the measurement time for a relatively large measurement time. Therefore, we find a transition from nonergodic (distributional) to ergodic diffusivity in a target search process. Moreover, this fluctuation analysis provides a method to estimate an unknown surface diffusivity.

Driven anomalous diffusion: An example from polymer stretching

The way tension propagates along a chain is a key to govern many anomalous dynamics in macromolecular systems. After introducing the weak and the strong force regimes of the tension propagation, we focus on the latter, in which the dynamical fluctuations of a segment in a long polymer during its stretching process is investigated. We show that the response, i.e., average drift, is anomalous, which is characterized by the nonlinear memory kernel, and its relation to the fluctuation is nontrivial. These features are discussed on the basis of the generalized Langevin equation, in which the role of the temporal change in spring constant due to the stress hardening is pinpointed. We carried out the molecular dynamics simulation, which supports our theory.

Origin of 1/f noise in hydration dynamics on lipid membrane surfaces

Water molecules on lipid membrane surfaces are known to contribute to membrane stability by connecting lipid molecules and acting as a water bridge. Although water structures and diffusivities near the membrane surfaces have been extensively studied, hydration dynamics on the surfaces has remained an open question. Here we investigate residence time statistics of water molecules on the surface of lipid membranes using all-atom molecular dynamics simulations. We show that hydration dynamics on the lipid membranes exhibits 1/f noise. Constructing a dichotomous process for the hydration dynamics, we find that residence times in each state follow a power-law with exponential cutoff and that the process can be regarded as a correlated renewal process where interoccurrence times are correlated. The results imply that the origin of the 1/f noise in hydration dynamics on the membrane surfaces is a combination of a power-law distribution with cutoff of interoccurrence times of switching events and a long-term correlation between the interoccurrence times. These results suggest that the 1/f noise attributed to the correlated renewal process may contribute to the stability of the hydration layers and lipid membranes.

Anomalous diffusion in a quenched-trap model on fractal lattices

Models with mixed origins of anomalous subdiffusion have been considered important for understanding transport in biological systems. Here one such mixed model, the quenched-trap model (QTM) on fractal lattices, is investigated. It is shown that both ensemble- and time-averaged mean-square displacements (MSDs) show subdiffusion with different scaling exponents, i.e., this system shows weak ergodicity breaking. Moreover, time-averaged MSD exhibits aging and converges to a random variable following the modified Mittag-Leffler distribution. It is also shown that the QTM on a fractal lattice cannot be reduced to the continuous-time random walks if the spectral dimension of the fractal lattice is less than 2.

Weak ergodicity breaking and ageing in anomalous diffusion

Modern single particle tracking techniques and many large scale simulations produce time series r(t) of the position of a tracer particle. Standardly these are evaluated in terms of the time averaged mean squared displacement. For ergodic processes such as Brownian motion, one can interpret the results of such an analysis in terms of the known theories for the corresponding ensemble averaged mean squared displacement, if only the measurement time is sufficiently long. In anomalous diffusion processes, that are widely observed over many orders of magnitude, the equivalence between (long) time and ensemble averages may be broken (weak ergodicity breaking). In such cases the time averages may no longer be interpreted in terms of ensemble theories. Here we collect some recent results on weakly non-ergodic systems with respect to the time averaged mean squared displacement and the inherent irreproducibility of individual measurements. We also address the phenomenon of ageing, the dependence of physical observables on the time span between initial preparation of the system and the start of the measurement.

Poisson property of the occurrence of flip-flops in a model membrane

How do lipid molecules in membranes perform a flip-flop? The flip-flops of lipid molecules play a crucial role in the formation and flexibility of membranes. However, little has been determined about the behavior of flip-flops, either experimentally, or in molecular dynamics simulations. Here, we provide numerical results of the flip-flops of model lipid molecules in a model membrane and investigate the statistical properties, using millisecond-order coarse-grained molecular simulations (dissipative particle dynamics). We find that there are three different ways of flip-flops, which can be clearly characterized by their paths on the free energy surface. Furthermore, we found that the probability of the number of the flip-flops is well fitted by the Poisson distribution, and the probability density function for the inter-occurrence times of flip-flops coincides with that of the forward recurrence times. These results indicate that the occurrence of flip-flops is a Poisson process, which will play an important role in the flexibilities of membranes.

Origin of subdiffusion of water molecules on cell membrane surfaces

Water molecules play an important role in providing unique environments for biological reactions on cell membranes. It is widely believed that water molecules form bridges that connect lipid molecules and stabilize cell membranes. Using all-atom molecular dynamics simulations, we show that translational and rotational diffusion of water molecules on lipid membrane surfaces exhibit subdiffusion and aging. Moreover, we provide evidence that both divergent mean trapping time (continuous-time random walk) and long-correlated noise (fractional Brownian motion) contribute to this subdiffusion. These results suggest that subdiffusion on cell membranes causes the water retardation, an enhancement of cell membrane stability, and a higher reaction efficiency.

Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity

We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells.