Sensing membrane voltage by reorientation of dipolar transmembrane peptides.

Membrane voltage plays a vital role in the behavior and functions of the lipid bilayer membrane. For instance, it regulates the exchange of molecules across the membrane through transmembrane proteins such as ion channels. In this paper, we study the membrane voltage-sensing mechanism, which entails the reorientation of α-helices with a change in the membrane voltage. We consider a helix having a large electrical macrodipole embedded in a lipid bilayer as a model system. We performed extensive molecular dynamics simulations to study the effect of variation of membrane voltage on the tilt angle of peptides and ascertain the optimal parameters for designing such a voltage-sensing peptide. A theoretical model for the system is also developed to investigate the interplay of competing effects of hydrophobic mismatch and dipole-electric field coupling on the tilt of the peptide and further explore the parameter space. This work opens the possibility for the design and fabrication of artificial dipolar membrane voltage-sensing elements for biomedical applications.

Multilayer diffusion modeling and Coherent anti-Stokes Raman scattering microscopy for spatially resolved water diffusion measurements in human skin.

In this work using Coherent anti-Stokes Raman Scattering microscopy, it was possible to directly measure the time dependent, spatially resolved change in concentration of water (D2 O) in intact skin tissue with a spatial resolution of under 1 µm, and combined with a multilayer diffusion model, diffusion coefficients at different depths in the tissue were extracted. The results show that the diffusion varies at different layers throughout the Stratum Corneum (SC), indicating that the SC is not a homogeneous barrier but a complicated heterogeneous structure. Interestingly, averaging over the diffusion at the different depths and samples gave a relatively constant value of 0.047 ± 0.01 µm2 /second. Treating the skin with acetone or tape stripping led to an increased diffusion coefficient of 0.064 ± 0.02 µm2 /second and 0.079 ± 0.03 µm2 /second, respectively. The combined method and model presented here shows potential for wide applications for measuring spatially resolved diffusion of different substances in a variety of different samples.

Objective comparison of methods to decode anomalous diffusion.

Deviations from Brownian motion leading to anomalous diffusion are found in transport dynamics from quantum physics to life sciences. The characterization of anomalous diffusion from the measurement of an individual trajectory is a challenging task, which traditionally relies on calculating the trajectory mean squared displacement. However, this approach breaks down for cases of practical interest, e.g., short or noisy trajectories, heterogeneous behaviour, or non-ergodic processes. Recently, several new approaches have been proposed, mostly building on the ongoing machine-learning revolution. To perform an objective comparison of methods, we gathered the community and organized an open competition, the Anomalous Diffusion challenge (AnDi). Participating teams applied their algorithms to a commonly-defined dataset including diverse conditions. Although no single method performed best across all scenarios, machine-learning-based approaches achieved superior performance for all tasks. The discussion of the challenge results provides practical advice for users and a benchmark for developers.

Stochastic unfolding of nanoconfined DNA: Experiments, model and Bayesian analysis.

Nanochannels provide a means for detailed experiments on the effect of confinement on biomacromolecules, such as DNA. Here we introduce a model for the complete unfolding of DNA from the circular to linear configuration. Two main ingredients are the entropic unfolding force and the friction coefficient for the unfolding process, and we describe the associated dynamics by a non-linear Langevin equation. By analyzing experimental data where DNA molecules are photo-cut and unfolded inside a nanochannel, our model allows us to extract values for the unfolding force as well as the friction coefficient for the first time. In order to extract numerical values for these physical quantities, we employ a recently introduced Bayesian inference framework. We find that the determined unfolding force is in agreement with estimates from a simple Flory-type argument. The estimated friction coefficient is in agreement with theoretical estimates for motion of a cylinder in a channel. We further validate the estimated friction constant by extracting this parameter from DNA's center-of-mass motion before and after unfolding, yielding decent agreement. We provide publically available software for performing the required image and Bayesian analysis.

Bayesian analysis of single-particle tracking data using the nested-sampling algorithm: maximum-likelihood model selection applied to stochastic-diffusivity data.

We employ Bayesian statistics using the nested-sampling algorithm to compare and rank multiple models of ergodic diffusion (including anomalous diffusion) as well as to assess their optimal parameters for in silico-generated and real time-series. We focus on the recently-introduced model of Brownian motion with "diffusing diffusivity"-giving rise to widely-observed non-Gaussian displacement statistics-and its comparison to Brownian and fractional Brownian motion, also for the time-series with some measurement noise. We conduct this model-assessment analysis using Bayesian statistics and the nested-sampling algorithm on the level of individual particle trajectories. We evaluate relative model probabilities and compute best-parameter sets for each diffusion model, comparing the estimated parameters to the true ones. We test the performance of the nested-sampling algorithm and its predictive power both for computer-generated (idealised) trajectories as well as for real single-particle-tracking trajectories. Our approach delivers new important insight into the objective selection of the most suitable stochastic model for a given time-series. We also present first model-ranking results in application to experimental data of tracer diffusion in polymer-based hydrogels.

Fitting a function to time-dependent ensemble averaged data.

Time-dependent ensemble averages, i.e., trajectory-based averages of some observable, are of importance in many fields of science. A crucial objective when interpreting such data is to fit these averages (for instance, squared displacements) with a function and extract parameters (such as diffusion constants). A commonly overlooked challenge in such function fitting procedures is that fluctuations around mean values, by construction, exhibit temporal correlations. We show that the only available general purpose function fitting methods, correlated chi-square method and the weighted least squares method (which neglects correlation), fail at either robust parameter estimation or accurate error estimation. We remedy this by deriving a new closed-form error estimation formula for weighted least square fitting. The new formula uses the full covariance matrix, i.e., rigorously includes temporal correlations, but is free of the robustness issues, inherent to the correlated chi-square method. We demonstrate its accuracy in four examples of importance in many fields: Brownian motion, damped harmonic oscillation, fractional Brownian motion and continuous time random walks. We also successfully apply our method, weighted least squares including correlation in error estimation (WLS-ICE), to particle tracking data. The WLS-ICE method is applicable to arbitrary fit functions, and we provide a publically available WLS-ICE software.

Bayesian inference with information content model check for Langevin equations.

The Bayesian data analysis framework has been proven to be a systematic and effective method of parameter inference and model selection for stochastic processes. In this work, we introduce an information content model check that may serve as a goodness-of-fit, like the χ^{2} procedure, to complement conventional Bayesian analysis. We demonstrate this extended Bayesian framework on a system of Langevin equations, where coordinate-dependent mobilities and measurement noise hinder the normal mean-squared displacement approach.

Optimization and modeling of the remote loading of luciferin into liposomes.

We carried out a mechanistic study to characterize and optimize the remote loading of luciferin into preformed liposomes of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dipalmitoyl-sn-glycero-3-phosphoglycerol (DPPC/DPPG) 7:3 mixtures. The influence of the loading agent (acetate, propionate, butyrate), the metal counterion (Na(+), K(+), Ca(+2), Mg(+2)), and the initial extra-liposomal amount of luciferin (nL(add)) on the luciferin Loading Efficiency (LE%) and luciferin-to-lipid weight ratio, i.e., Loading Capacity (LC), in the final formulation was determined. In addition, the effect of the loading process on the colloidal stability and phase behavior of the liposomes was monitored. Based on our experimental results, a theoretical model was developed to describe the course of luciferin remote loading. It was found that the highest luciferin loading was obtained with magnesium acetate. The use of longer aliphatic carboxylates or inorganic proton donors pronouncedly reduced luciferin loading, whereas the effect of the counterion was modest. The remote-loading process barely affected the colloidal stability and drug retention of the liposomes, albeit with moderate luciferin-induced membrane perturbations. The correlation between luciferin loading, expressed as LE% and LC, and nL(add) was established, and under our conditions the maximum LC was attained using an nL(add) of around 2.6µmol. Higher amounts of luciferin tend to pronouncedly perturb the liposome stability and luciferin retention. Our theoretical model furnishes a fair quantitative description of the correlation between nL(add) and luciferin loading, and a membrane permeability coefficient for uncharged luciferin of 1×10(-8)cm/s could be determined. We believe that our study will prove very useful to optimize the remote-loading strategies of moderately polar carboxylic acid drugs in general.

Quantifying the Relationship between Curvature and Electric Potential in Lipid Bilayers.

Cellular membranes mediate vital cellular processes by being subject to curvature and transmembrane electrical potentials. Here we build upon the existing theory for flexoelectricity in liquid crystals to quantify the coupling between lipid bilayer curvature and membrane potentials. Using molecular dynamics simulations, we show that headgroup dipole moments, the lateral pressure profile across the bilayer, and spontaneous curvature all systematically change with increasing membrane potentials. In particular, there is a linear dependence between the bending moment (the product of bending rigidity and spontaneous curvature) and the applied membrane potentials. We show that biologically relevant membrane potentials can induce biologically relevant curvatures corresponding to radii of around 500 nm. The implications of flexoelectricity in lipid bilayers are thus likely to be of considerable consequence both in biology and in model lipid bilayer systems.

Imaging approaches for analysis of cholesterol distribution and dynamics in the plasma membrane.

Cholesterol is an important lipid component of the plasma membrane (PM) of mammalian cells, where it is involved in control of many physiological processes, such as endocytosis, cell migration, cell signalling and surface ruffling. In an attempt to explain these functions of cholesterol, several models have been put forward about cholesterol's lateral and transbilayer organization in the PM. In this article, we review imaging techniques developed over the last two decades for assessing the distribution and dynamics of cholesterol in the PM of mammalian cells. Particular focus is on fluorescence techniques to study the lateral and inter-leaflet distribution of suitable cholesterol analogues in the PM of living cells. We describe also several methods for determining lateral cholesterol dynamics in the PM including fluorescence recovery after photobleaching (FRAP), fluorescence correlation spectroscopy (FCS), single particle tracking (SPT) and spot variation FCS coupled to stimulated emission depletion (STED) microscopy. For proper interpretation of such measurements, we provide some background in probe photophysics and diffusion phenomena occurring in cell membranes. In particular, we show the equivalence of the reaction-diffusion approach, as used in FRAP and FCS, and continuous time random walk (CTRW) models, as often invoked in SPT studies. We also discuss mass spectrometry (MS) based imaging of cholesterol in the PM of fixed cells and compare this method with fluorescence imaging of sterols. We conclude that evidence from many experimental techniques converges towards a model of a homogeneous distribution of cholesterol with largely free and unhindered diffusion in both leaflets of the PM.

Real sequence effects on the search dynamics of transcription factors on DNA.

Recent experiments show that transcription factors (TFs) indeed use the facilitated diffusion mechanism to locate their target sequences on DNA in living bacteria cells: TFs alternate between sliding motion along DNA and relocation events through the cytoplasm. From simulations and theoretical analysis we study the TF-sliding motion for a large section of the DNA-sequence of a common E. coli strain, based on the two-state TF-model with a fast-sliding search state and a recognition state enabling target detection. For the probability to detect the target before dissociating from DNA the TF-search times self-consistently depend heavily on whether or not an auxiliary operator (an accessible sequence similar to the main operator) is present in the genome section. Importantly, within our model the extent to which the interconversion rates between search and recognition states depend on the underlying nucleotide sequence is varied. A moderate dependence maximises the capability to distinguish between the main operator and similar sequences. Moreover, these auxiliary operators serve as starting points for DNA looping with the main operator, yielding a spectrum of target detection times spanning several orders of magnitude. Auxiliary operators are shown to act as funnels facilitating target detection by TFs.

Stochastic thermodynamics of a tagged particle within a harmonic chain.

We study the stochastic thermodynamics of an overdamped harmonic chain, which can be viewed equivalently as a one-dimensional Rouse chain or as an approximate model of single file diffusion. We discuss mainly two levels of description of this system: the Markovian level for which the trajectories of all the particles of the chain are known and the non-Markovian level in which only the motion of a tagged particle is available. For each case, we analyze the energy dissipation and its dependence on initial conditions. Surprisingly, we find that the average coarse-grained entropy production rate can become transiently negative when an oscillating force is applied to the tagged particle. This occurs due to memory effects as shown in a framework based on path integrals or on a generalized Langevin equation.

Universality and nonuniversality of mobility in heterogeneous single-file systems and Rouse chains.

We study analytically the tracer particle mobility in single-file systems with distributed friction constants. Our system serves as a prototype for nonequilibrium, heterogeneous, strongly interacting Brownian systems. The long time dynamics for such a single-file setup belongs to the same universality class as the Rouse model with dissimilar beads. The friction constants are drawn from a density ϱ(ξ), and we derive an asymptotically exact solution for the mobility distribution P[µ0(s)], where µ0(s) is the Laplace-space mobility. If ϱ is light tailed (first moment exists), we find a self-averaging behavior: P[µ0(s)]=δ[µ0(s)-µ(s)], with µ(s)∝s1/2. When ϱ(ξ) is heavy tailed, ϱ(ξ)≃ξ-1-α(0<α<1) for large ξ, we obtain moments 〈[µs(0)]n〉∝sßn, where ß=1/(1+α) and there is no self-averaging. The results are corroborated by simulations.

Microscopic origin of the logarithmic time evolution of aging processes in complex systems.

There exists compelling experimental evidence in numerous systems for logarithmically slow time evolution, yet its full theoretical understanding remains elusive. We here introduce and study a generic transition process in complex systems, based on nonrenewal, aging waiting times. Each state n of the system follows a local clock initiated at t = 0. The random time τ between clock ticks follows the waiting time density ψ(τ). Transitions between states occur only at local clock ticks and are hence triggered by the local forward waiting time, rather than by ψ(τ). For power-law forms ψ(τ) ≃ τ(-1-α) (0 < α < 1) we obtain a logarithmic time evolution of the state number ⟨n(t)⟩ ≃ log(t/t(0)), while for α > 2 the process becomes normal in the sense that ⟨n(t)⟩ ≃ t. In the intermediate range 1 < α < 2 we find the power-law growth ⟨n(t)⟩ ≃ t(α-1). Our model provides a universal description for transition dynamics between aging and nonaging states.

Bulk-mediated diffusion on a planar surface: full solution.

We consider the effective surface motion of a particle that intermittently unbinds from a planar surface and performs bulk excursions. Based on a random-walk approach, we derive the diffusion equations for surface and bulk diffusion including the surface-bulk coupling. From these exact dynamic equations, we analytically obtain the propagator of the effective surface motion. This approach allows us to deduce a superdiffusive, Cauchy-type behavior on the surface, together with exact cutoffs limiting the Cauchy form. Moreover, we study the long-time dynamics for the surface motion.

Quantitative fluorescence loss in photobleaching for analysis of protein transport and aggregation.

BACKGROUND: Fluorescence loss in photobleaching (FLIP) is a widely used imaging technique, which provides information about protein dynamics in various cellular regions. In FLIP, a small cellular region is repeatedly illuminated by an intense laser pulse, while images are taken with reduced laser power with a time lag between the bleaches. Despite its popularity, tools are lacking for quantitative analysis of FLIP experiments. Typically, the user defines regions of interest (ROIs) for further analysis which is subjective and does not allow for comparing different cells and experimental settings. RESULTS: We present two complementary methods to detect and quantify protein transport and aggregation in living cells from FLIP image series. In the first approach, a stretched exponential (StrExp) function is fitted to fluorescence loss (FL) inside and outside the bleached region. We show by reaction-diffusion simulations, that the StrExp function can describe both, binding/barrier-limited and diffusion-limited FL kinetics. By pixel-wise regression of that function to FL kinetics of enhanced green fluorescent protein (eGFP), we determined in a user-unbiased manner from which cellular regions eGFP can be replenished in the bleached area. Spatial variation in the parameters calculated from the StrExp function allow for detecting diffusion barriers for eGFP in the nucleus and cytoplasm of living cells. Polyglutamine (polyQ) disease proteins like mutant huntingtin (mtHtt) can form large aggregates called inclusion bodies (IB's). The second method combines single particle tracking with multi-compartment modelling of FL kinetics in moving IB's to determine exchange rates of eGFP-tagged mtHtt protein (eGFP-mtHtt) between aggregates and the cytoplasm. This method is self-calibrating since it relates the FL inside and outside the bleached regions. It makes it therefore possible to compare release kinetics of eGFP-mtHtt between different cells and experiments. CONCLUSIONS: We present two complementary methods for quantitative analysis of FLIP experiments in living cells. They provide spatial maps of exchange dynamics and absolute binding parameters of fluorescent molecules to moving intracellular entities, respectively. Our methods should be of great value for quantitative studies of intracellular transport.

Two-photon time-lapse microscopy of BODIPY-cholesterol reveals anomalous sterol diffusion in chinese hamster ovary cells.

BACKGROUND: Cholesterol is an important membrane component, but our knowledge about its transport in cells is sparse. Previous imaging studies using dehydroergosterol (DHE), an intrinsically fluorescent sterol from yeast, have established that vesicular and non-vesicular transport modes contribute to sterol trafficking from the plasma membrane. Significant photobleaching, however, limits the possibilities for in-depth analysis of sterol dynamics using DHE. Co-trafficking studies with DHE and the recently introduced fluorescent cholesterol analog BODIPY-cholesterol (BChol) suggested that the latter probe has utility for prolonged live-cell imaging of sterol transport. RESULTS: We found that BChol is very photostable under two-photon (2P)-excitation allowing the acquisition of several hundred frames without significant photobleaching. Therefore, long-term tracking and diffusion measurements are possible. Two-photon temporal image correlation spectroscopy (2P-TICS) provided evidence for spatially heterogeneous diffusion constants of BChol varying over two orders of magnitude from the cell interior towards the plasma membrane, where D ~ 1.3 µm2/s. Number and brightness (N&B) analysis together with stochastic simulations suggest that transient partitioning of BChol into convoluted membranes slows local sterol diffusion. We observed sterol endocytosis as well as fusion and fission of sterol-containing endocytic vesicles. The mobility of endocytic vesicles, as studied by particle tracking, is well described by a model for anomalous subdiffusion on short time scales with an anomalous exponent α ~ 0.63 and an anomalous diffusion constant of Dα = 1.95 x 10-3 µm2/sα. On a longer time scale (t > ~5 s), a transition to superdiffusion consistent with slow directed transport with an average velocity of v ~ 6 x 10-3 µm/s was observed. We present an analytical model that bridges the two regimes and fit this model to vesicle trajectories from control cells and cells with disrupted microtubule or actin filaments. Both treatments reduced the anomalous diffusion constant and the velocity by ~40-50%. CONCLUSIONS: The mobility of sterol-containing vesicles on the short time scale could reflect dynamic rearrangements of the cytoskeleton, while directed transport of sterol vesicles occurs likely along both, microtubules and actin filaments. Spatially varying anomalous diffusion could contribute to fine-tuning and local regulation of intracellular sterol transport.

Effective surface motion on a reactive cylinder of particles that perform intermittent bulk diffusion.

In many biological and small scale technological applications particles may transiently bind to a cylindrical surface. In between two binding events the particles diffuse in the bulk, thus producing an effective translation on the cylindrical surface. We here derive the effective motion on the surface allowing for additional diffusion on the cylindrical surface itself. We find explicit solutions for the number of adsorbed particles at one given instant, the effective surface displacement, as well as the surface propagator. In particular sub- and superdiffusive regimes are found, as well as an effective stalling of diffusion visible as a plateau in the mean squared displacement. We also investigate the corresponding first passage problem.

Elastic moderation of intrinsically applied tension in lipid membranes.

Tension in lipid membranes is often controlled externally, by pulling on the boundary of the membrane or changing osmotic pressure across a curved membrane. But modifications of the tension can also be induced in an internal fashion, for instance as a byproduct of changing a membrane's electric potential or, as observed experimentally, by activity of membrane proteins. Here we develop a theory that demonstrates how such internal contributions to the tension are moderated through elastic stretching of the membrane when the membrane is initially in a low-tension floppy state.