Anomalous diffusion in lipid bilayers is usually attributed to viscoelastic behavior. We compute the scaling exponent of relative fluctuations of the time-averaged mean square displacement in a lipid bilayer, by using a molecular dynamics simulation. According to the continuous time random walk theory, this exponent indicates non-Gaussian behavior caused by a power-law trapping time. Our results provide the first evidence that a lipid bilayer has not only viscoelastic properties but also trapping times distributed according to a power law.

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.

SummaryA molecular simulation pipeline for determining the mode of interaction of pleckstrin homology (PH) domains with phosphatidylinositol phosphate (PIP)-containing lipid bilayers is presented. We evaluate our methodology for the GRP1 PH domain via comparison with structural and biophysical data. Coarse-grained simulations yield a 2D density landscape for PH/membrane interactions alongside residue contact profiles. Predictions of the membrane localization and interactions of 13 PH domains reveal canonical, non-canonical, and dual PIP-binding sites on the proteins. Thus, the PH domains associate with the PIP molecules in the membrane via a highly positively charged loop. Some PH domains exhibit modes of interaction with PIP-containing membranes additional to this canonical binding mode. All 13 PH domains cause a degree of local clustering of PIP molecules upon binding to the membrane. This provides a global picture of PH domain interactions with membranes. The high-throughput approach could be extended to other families of peripheral membrane proteins.

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.

No abstract