A lattice model of ternary mixtures of lipids and cholesterol with tunable domain sizes.

Much of our understanding of the physical properties of raft domains in biological membranes, and some insight into the mechanisms underlying their formation stem from atomistic simulations of simple model systems, especially ternary mixtures consisting of saturated and unsaturated lipids, and cholesterol (Chol). To explore the properties of such systems at large spatial scales, we here present a simple ternary mixture lattice model, involving a small number of nearest neighbor interaction terms. Monte Carlo simulations of mixtures with different compositions show excellent agreement with experimental and atomistic simulation observations across multiple scales, ranging from the local distributions of lipids to the phase diagram of the system. The simplicity of the model allows us to identify the roles played by the different interactions between components, and the interplay between them. Importantly, by changing the value of one of the model parameters, we can tune the size of the liquid-ordered domains, thereby simulating both Type II mixtures exhibiting macroscopic phase separation and Type I mixtures with nanoscopic domains. The Type II mixture simulation results fit well to the experimentally determined phase diagram of mixtures containing saturated DPPC/unsaturated DOPC/Chol. When the tunable parameter is changed, we obtain the Type I version of DPPC/DOPC/Chol, i.e., a mixture not showing thermodynamic phase transitions but one that may be fitted to the same phase diagram if local measures are used to distinguish between the different states. Our model results suggest that short range packing is likely to be a key regulator of the stability and size distribution of biological rafts.

Characterizing the heterogeneity of membrane liquid-ordered domains.

We use a lattice model of a ternary mixture containing saturated and unsaturated lipids with cholesterol (Chol), to study the structural properties characterizing the coexistence between the liquid-disordered and liquid-ordered phases. Depending on the affinity of the saturated and unsaturated lipids, the system may exhibit macroscopic (thermodynamic) liquid-liquid phase separation or be divided into small-size liquid-ordered domains surrounded by a liquid-disordered matrix. In both cases, it is found that the nanoscale structure of the liquid-ordered regions is heterogeneous, and that they are partitioned into Chol-rich sub-domains and Chol-free, gel-like, nano-clusters. This emerges as a characteristic feature of the liquid-ordered state, which helps distinguishing between liquid-ordered domains in a two-phase mixture, and similar-looking domains in a one-phase mixture that are rich in saturated lipids and Chol, but are merely thermal density fluctuations. The nano-structure heterogeneity of the liquid-ordered phase can be detected by suitable experimental spectroscopic methods and is observed also in atomistic computer simulations.

Reactive Momentum Transfer Contributes to the Self-Propulsion of Janus Particles.

We report simulations of a spherical Janus particle undergoing exothermic surface reactions around one pole only. Our model excludes self-phoretic transport by design. Nevertheless, net motion occurs from direct momentum transfer between solvent and colloid, with speed scaling as the square root of the energy released during the reaction. We find that such propulsion is dominated by the system's short-time response, when neither the time dependence of the flow around the colloid nor the solvent compressibility can be ignored. Our simulations agree reasonably well with previous experiments.

A simple statistical-mechanical interpretation of Onsager reciprocal relations and Derjaguin theory of thermo-osmosis.

The application of a temperature gradient along a fluid-solid interface generates stresses in the fluid causing "thermo-osmotic" flow. Much of the understanding of this phenomenon is based on Derjaguin's work relating thermo-osmotic flows to the mechano-caloric effect, namely, the interfacial heat flow induced by a pressure gradient. This is done by using Onsager's reciprocity relationship for the equivalence of the thermo-osmotic and mechano-caloric cross-term transport coefficients. Both Derjaguin theory and Onsager framework for out-of-equilibrium systems are formulated in macroscopic thermodynamics terms and lack a clear interpretation at the molecular level. Here, we use statistical-mechanical tools to derive expressions for the transport cross-coefficients and, thereby, to directly demonstrate their equality. This is done for two basic models: i) an incopressible continuum solvent containing non-interacting solute particles, and ii) a single-component fluid without thermal expansivity. The derivation of the mechano-caloric coefficient appears to be remarkably simple, and provides a simple interpretation for the connection between interfacial heat and particle fluxes. We use this interpretation to consider yet another example, which is an electrolyte interacting with a uniformly charged surface in the strong screening (Debye-Hückel) regime.

Interplay between membrane elasticity and active cytoskeleton forces regulates the aggregation dynamics of the immunological synapse.

Adhesion between a T cell and an antigen presenting cell is achieved by TCR-pMHC and LFA1-ICAM1 protein complexes. These segregate to form a special pattern, known as the immunological synapse (IS), consisting of a central quasi-circular domain of TCR-pMHC bonds surrounded by a peripheral domain of LFA1-ICAM1 complexes. Insights gained from imaging studies had led to the conclusion that the formation of the central adhesion domain in the IS is driven by active (ATP-driven) mechanisms. Recent studies, however, suggested that passive (thermodynamic) mechanisms may also play an important role in this process. Here, we present a simple physical model, taking into account the membrane-mediated thermodynamic attraction between the TCR-pMHC bonds and the effective forces that they experience due to ATP-driven actin retrograde flow and transport by dynein motor proteins. Monte Carlo simulations of the model exhibit a good spatio-temporal agreement with the experimentally observed pattern evolution of the TCR-pMHC microclusters. The agreement is lost when one of the aggregation mechanisms is "muted", which helps to identify their respective roles in the process. We conclude that actin retrograde flow drives the centripetal motion of TCR-pMHC bonds, while the membrane-mediated interactions facilitate microcluster formation and growth. In the absence of dynein motors, the system evolves into a ring-shaped pattern, which highlights the role of dynein motors in the formation of the final concentric pattern. The interplay between the passive and active mechanisms regulates the rate of the accumulation process, which in the absence of one them proceeds either too quickly or slowly.

Formation of semi-dilute adhesion domains driven by weak elasticity-mediated interactions.

Cell-cell adhesion is established by specific binding of receptor and ligand proteins anchored in the cell membranes. The adhesion bonds attract each other and often aggregate into large clusters that are central to many biological processes. One possible origin of attractive interactions between adhesion bonds is the elastic response of the membranes to their deformation by the bonds. Here, we analyze these elasticity-mediated interactions using a novel mean-field approach. Our analysis of systems at different densities of bonds, ϕ, reveals that the phase diagram, i.e., the binodal and spinodal lines, exhibit a nearly universal behavior when the temperature T is plotted against the scaled density x = ϕξ(2), where ξ is the linear size of the membrane's region affected by the presence of a single isolated bond. The critical point (ϕc , Tc) is located at very low densities, and slightly below Tc we identify phase coexistence between two low-density phases. Dense adhesion domains are observed only when the height by which the bonds deform the membranes, h0, is much larger than their thermal roughness, Δ, which occurs at very low temperatures T≪Tc. We, thus, conclude that the elasticity-mediated interactions are weak and cannot be regarded as responsible for the formation of dense adhesion domains. The weakness of the elasticity-mediated effect and its relevance to dilute systems only can be attributed to the fact that the membrane's elastic energy saturates in the semi-dilute regime, when the typical spacing between the bonds r≳ξ, i.e., for x≲ 1. Therefore, at higher densities, only the mixing entropy of the bonds (which always favors uniform distributions) is thermodynamically relevant. We discuss the implications of our results for the question of immunological synapse formation, and demonstrate that the elasticity-mediated interactions may be involved in the aggregation of these semi-dilute membrane domains.

Correction: The thermodynamics of endosomal escape and DNA release from lipoplexes.

Correction for 'The thermodynamics of endosomal escape and DNA release from lipoplexes' by Yotam Y. Avital et al., Phys. Chem. Chem. Phys., 2016, 18, 2591-2596.

The thermodynamics of endosomal escape and DNA release from lipoplexes.

Complexes of cationic and neutral lipids and DNA (lipoplexes) are emerging as promising vectors for gene therapy applications. Their appeal stems from their non pathogenic nature and the fact that they self-assemble under conditions of thermal equilibrium. Lipoplex adhesion to the cell plasma membrane initiates a three-stage process termed transfection, consisting of (i) endocytosis, (ii) lipoplex breakdown, and (iii) DNA release followed by gene expression. As successful transfection requires lipoplex degradation, it tends to be hindered by the lipoplex thermodynamic stability; nevertheless, it is known that the transfection process may proceed spontaneously. Here, we use a simple model to study the thermodynamic driving forces governing transfection. We demonstrate that after endocytosis [stage (i)], the lipoplex becomes inherently unstable. This instability, which is triggered by interactions between the cationic lipids of the lipoplex and the anionic lipids of the enveloping plasma membrane, is entropically controlled involving both remixing of the lipids and counterions release. Our detailed calculation shows that the free energy gain during stage (ii) is approximately linear in Φ+, the mole fraction of cationic lipids in the lipoplex. This free energy gain, ΔF, reduces the barrier for fusion between the enveloping and the lipoplex bilayers, which produces a hole allowing for DNA release [stage (iii)]. The linear relationship between ΔF and the fraction of cationic lipids explains the experimentally observed exponential increase of transfection efficiency with Φ+ in lamellar lipoplexes.

On the connection between dissipative particle dynamics and the Itô-Stratonovich dilemma.

Dissipative Particle Dynamics (DPD) is a popular simulation model for investigating hydrodynamic behavior of systems with non-negligible equilibrium thermal fluctuations. DPD employs soft core repulsive interactions between the system particles, thus allowing them to overlap. This supposedly permits relatively large integration time steps, which is an important feature for simulations on large temporal scales. In practice, however, an increase in the integration time step leads to increasingly larger systematic errors in the sampling statistics. Here, we demonstrate that the prime origin of these systematic errors is the multiplicative nature of the thermal noise term in Langevin's equation, i.e., the fact that it depends on the instantaneous coordinates of the particles. This lead to an ambiguity in the interpretation of the stochastic differential Langevin equation, known as the Itô-Stratonovich dilemma. Based on insights from previous studies of the dilemma, we propose a novel algorithm for DPD simulations exhibiting almost an order of magnitude improvement in accuracy, and nearly twice the efficiency of commonly used DPD Langevin thermostats.

Formation of adhesion domains in stressed and confined membranes.

The adhesion bonds connecting a lipid bilayer to an underlying surface may undergo a condensation transition resulting from an interplay between a short range attractive potential between them, and a long range fluctuation-induced potential of mean force. Here, we use computer simulations of a coarse-grained molecular model of supported lipid bilayers to study this transition in confined membranes, and in membranes subjected to a non-vanishing surface tension. Our results show that confinement may alter significantly the condensation transition of the adhesion bonds, whereas the application of surface tension has a very minor effect on it. We also investigate domain formation in membranes under negative tension which, in free membranes, causes the enhancement of the amplitude of membrane thermal undulations. Our results indicate that in supported membranes, this effect of a negative surface tension on the fluctuation spectrum is largely eliminated by the pressure resulting from the mixing entropy of the adhesion bonds.

Small membranes under negative surface tension.

We use computer simulations and a simple free energy model to study the response of a bilayer membrane to the application of a negative (compressive) mechanical tension. Such a tension destabilizes the long wavelength undulation modes of giant vesicles, but it can be sustained when small membranes and vesicles are considered. Our negative tension simulation results reveal two regimes-(i) a weak negative tension regime characterized by stretching-dominated elasticity and (ii) a strong negative tension regime featuring bending-dominated elastic behavior. This resembles the findings of the classic Evans and Rawicz micropipette aspiration experiment in giant unilamellar vesicles (GUVs) [E. Evans and W. Rawicz, Phys, Rev. Lett. 64, 2094 (1990)]. However, in GUVs the crossover between the two elasticity regimes occurs at a small positive surface tension, while in smaller membranes it takes place at a moderate negative tension. Another interesting observation concerning the response of a small membrane to negative surface tension is related to the relationship between the mechanical and fluctuation tensions, which are equal to each other for non-negative values. When the tension decreases to negative values, the fluctuation tension γ drops somewhat faster than the mechanical tension τ in the small negative tension regime, before it saturates (and becomes larger than τ) for large negative tensions. The bending modulus exhibits an "opposite" trend. It remains almost unchanged in the stretching-dominated elastic regime, and decreases in the bending-dominated regime. Both the amplitudes of the thermal height undulations and the projected area variations diverge at the onset of mechanical instability.

Elasticity and mechanical instability of charged lipid bilayers in ionic solutions.

We use coarse-grained Monte Carlo simulations to study the elastic properties of charged membranes in solutions of monovalent and pentavalent counterions. The simulation results of the two cases reveal trends opposite to each other. The bending rigidity and projected area increase with the membrane charge density for monovalent counterions, while they decrease for the pentavalent ions. These observations can be related to the counterion screening of the lipid charges. While the monovalent counterions only weakly screen the Coulomb interactions, which implies a repulsive Coulomb system, the multivalent counterions condense on the membrane and, through spatial charge correlations, make the effective interactions due to the charged lipids attractive. The differences in the elastic properties of the charged membranes in monovalent and multivalent counterion solutions are reflected in the mechanisms leading to their mechanical instability at high charge densities. In the former case, the membranes develop pores to relieve the electrostatic tensile stresses, while in the latter case, the membrane exhibits large wavelength bending instability.

Coarse-grained molecular simulations of membrane adhesion domains.

We use a coarse-grained molecular model of supported lipid bilayers to study the formation of adhesion domains. We find that this process is a first order phase transition, triggered by a combination of pairwise short range attractive interactions between the adhesion bonds and many-body Casimir-like interactions, mediated by the membrane thermal undulations. The simulation results display an excellent agreement with the recently proposed Weil-Farago two-dimensional lattice model, in which the occupied and empty sites represent, respectively, the adhesion bonds and unbound segments of the membrane. A second phase transition, into a hexatic phase, is observed when the attraction between the adhesion bonds is further strengthened.

Constant pressure and temperature discrete-time Langevin molecular dynamics.

We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are built on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems-a one-dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three-dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb and Dünweg [J. Chem. Phys. 111, 4453 (1999)], show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation.

Confined run-and-tumble particles with non-Markovian tumbling statistics.

Confined active particles constitute simple, yet realistic, examples of systems that converge into a nonequilibrium steady state. We investigate a run-and-tumble particle in one spatial dimension, trapped by an external potential, with a given distribution g(t) of waiting times between tumbling events whose mean value is equal to τ. Unless g(t) is an exponential distribution (corresponding to a constant tumbling rate), the process is non-Markovian, which makes the analysis of the model particularly challenging. We use an analytical framework involving effective position-dependent tumbling rates to develop a numerical method that yields the full steady-state distribution (SSD) of the particle's position. The method is very efficient and requires modest computing resources, including in the large-deviation and/or small-τ regime, where the SSD can be related to the the large-deviation function, s(x), via the scaling relation P_{st}(x)∼e^{-s(x)/τ}.

300-Times-Increased Diffusive Skyrmion Dynamics and Effective Pinning Reduction by Periodic Field Excitation.

Thermally induced skyrmion dynamics, as well as skyrmion pinning effects, in thin films have attracted significant interest. While pinning poses challenges in deterministic skyrmion devices and slows down skyrmion diffusion, for applications in non-conventional computing, both pinning of an appropriate strength and skyrmion diffusion speed are key. Here, periodic field excitations are employed to realize an increase of the skyrmion diffusion by more than two orders of magnitude. Amplifying the excitation, a drastic reduction of the effective skyrmion pinning, is reported, and a transition from pinning-dominated diffusive hopping to dynamics approaching free diffusion is observed. By tailoring the field oscillation frequency and amplitude, a continuous tuning of the effective pinning and skyrmion dynamics is demonstrated, which is a key asset and enabler for non-conventional computing applications. It is found that the periodic excitations additionally allow stabilization of skyrmions at different sizes for field values that are inaccessible in static systems, opening up new approaches to ultrafast skyrmion motion by transiently exciting moving skyrmions.

Hierarchical self-organization of cytoskeletal active networks.

The structural reorganization of the actin cytoskeleton is facilitated through the action of motor proteins that crosslink the actin filaments and transport them relative to each other. Here, we present a combined experimental-computational study that probes the dynamic evolution of mixtures of actin filaments and clusters of myosin motors. While on small spatial and temporal scales the system behaves in a very noisy manner, on larger scales it evolves into several well distinct patterns such as bundles, asters and networks. These patterns are characterized by junctions with high connectivity, whose formation is possible due to the organization of the motors in 'oligoclusters' (intermediate-size aggregates). The simulations reveal that the self-organization process proceeds through a series of hierarchical steps, starting from local microscopic moves and ranging up to the macroscopic large scales where the steady-state structures are formed. Our results shed light on the mechanisms involved in processes such as cytokinesis and cellular contractility, where myosin motors organized in clusters operate cooperatively to induce the structural organization of cytoskeletal networks.

Nonequilibrium steady state for harmonically confined active particles.

We study the full nonequilibrium steady-state distribution P_{st}(X) of the position X of a damped particle confined in a harmonic trapping potential and experiencing active noise whose correlation time τ_{c} is assumed to be very short. Typical fluctuations of X are governed by a Boltzmann distribution with an effective temperature that is found by approximating the noise as white Gaussian thermal noise. However, large deviations of X are described by a non-Boltzmann steady-state distribution. We find that, in the limit τ_{c}→0, they display the scaling behavior P_{st}(X)∼e^{-s(X)/τ_{c}}, where s(X) is the large-deviation function. We obtain an expression for s(X) for a general active noise and calculate it exactly for the particular case of telegraphic (dichotomous) noise.

Thermodynamics of a Brownian particle in a nonconfining potential.

We consider the overdamped Brownian dynamics of a particle starting inside a square potential well which, upon exiting the well, experiences a flat potential where it is free to diffuse. We calculate the particle's probability distribution function (PDF) at coordinate x and time t, P(x,t), by solving the corresponding Smoluchowski equation. The solution is expressed by a multipole expansion, with each term decaying t^{1/2} faster than the previous one. At asymptotically large times, the PDF outside the well converges to the Gaussian PDF of a free Brownian particle. The average energy, which is proportional to the probability of finding the particle inside the well, diminishes as E∼1/t^{1/2}. Interestingly, we find that the free energy of the particle, F, approaches the free energy of a freely diffusing particle, F_{0}, as δF=F-F_{0}∼1/t, i.e., at a rate faster than E. We provide analytical and computational evidence that this scaling behavior of δF is a general feature of Brownian dynamics in nonconfining potential fields. Furthermore, we argue that δF represents a diminishing entropic component which is localized in the region of the potential, and which diffuses away with the spreading particle without being transferred to the heat bath.

Biased transport of elastic cytoskeletal filaments with alternating polarities by molecular motors.

We present a simple model for the bidirectional dynamics of actin bundles with alternating polarities in gliding assays with nonprocessive myosin motors. The bundle is represented as an elastic chain consisting of monomers with positive and negative polarities. The motion of the bundle is induced by the pulling forces of the underlying motors which stochastically attach to the monomers and, depending on their polarities, pull them in the right or left direction. We demonstrate that perfectly apolar chains consisting of equal numbers of monomers with positive and negative polarities may exhibit biased bidirectional motion with a nonzero drift. This effect is attributed to the elastic tension developed in the chain due to the action of the myosin motors. We also show that as a result of this tension, the attachment probability of the motors is greatly reduced and becomes strongly dependent on the length of the chain.