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

Interdependence of cholesterol distribution and conformational order in lipid bilayers

We show, via molecular simulations, that not only does cholesterol induce a lipid order, but the lipid order also enhances cholesterol localization within the lipid leaflets. Therefore, there is a strong interdependence between these two phenomena. In the ordered phase, cholesterol molecules are predominantly present in the bilayer leaflets and orient themselves parallel to the bilayer normal. In the disordered phase, cholesterol molecules are mainly present near the center of the bilayer at the midplane region and are oriented orthogonal to the bilayer normal. At the melting temperature of the lipid bilayers, cholesterol concentration in the leaflets and the bilayer midplane is equal. This result suggests that the localization of cholesterol in the lipid bilayers is mainly dictated by the degree of ordering of the lipid bilayer. We validate our findings on 18 different lipid bilayer systems, obtained from three different phospholipid bilayers with varying concentrations of cholesterol. To cover a large temperature range in simulations, we employ the Dry Martini force field. We demonstrate that the Dry and the Wet Martini (with polarizable water) force fields produce comparable results.

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.

LiPyphilic: A Python Toolkit for the Analysis of Lipid Membrane Simulations

Molecular dynamics simulations are now widely used to study emergent phenomena in lipid membranes with complex compositions. Here, we present LiPyphilic-a fast, fully tested, and easy-to-install Python package for analyzing such simulations. Analysis tools in LiPyphilic include the identification of cholesterol flip-flop events, the classification of local lipid environments, and the degree of interleaflet registration. LiPyphilic is both force field- and resolution-agnostic, and by using the powerful atom selection language of MDAnalysis, it can handle membranes with highly complex compositions. LiPyphilic also offers two on-the-fly trajectory transformations to (i) fix membranes split across periodic boundaries and (ii) perform nojump coordinate unwrapping. Our implementation of nojump unwrapping accounts for fluctuations in the box volume under the NPT ensemble-an issue that most current implementations have overlooked. The full documentation of LiPyphilic, including installation instructions and links to interactive online tutorials, is available at https://lipyphilic.readthedocs.io/en/latest.

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.

Quantitative relationship between cholesterol distribution and ordering of lipids in asymmetric lipid bilayers

The plasma membrane of eukaryotic cells is known to be compositionally asymmetric. Certain phospholipids, such as sphingomyelin and phosphatidylcholine species, are predominantly localized in the outer leaflet, while phosphatidylethanolamine and phosphatidylserine species primarily reside in the inner leaflet. While phospholipid asymmetry between the membrane leaflets is well established, there is no consensus about cholesterol distribution between the two leaflets. We have performed a systematic study, via molecular simulations, of how the spatial distribution of cholesterol molecules in different "asymmetric" lipid bilayers are affected by the lipids' backbone, head-type, unsaturation, and chain-length by considering an asymmetric bilayer mimicking the plasma membrane lipids of red blood cells, as well as seventeen other asymmetric bilayers comprising of different lipid types. Our results reveal that the distribution of cholesterol in the leaflets is solely a function of the extent of ordering of the lipids within the leaflets. The ratio of the amount of cholesterol matches the ratio of lipid order in the two leaflets, thus providing a quantitative relationship between the two. These results are understood by the observation that asymmetric bilayers with equimolar amount of lipids in the two leaflets develop tensile and compressive stresses due to differences in the extent of lipid order. These stresses are alleviated by the transfer of cholesterol from the leaflet in compressive stress to the one in tensile stress. These findings are important in understanding the biology of the cell membrane, especially with regard to the composition of the membrane leaflets.

Interleaflet Decoupling in a Lipid Bilayer at Excess Cholesterol Probed by Spectroscopic Ellipsometry and Simulations

Artificial lipid membranes are often investigated as a replica of the cell membrane in the form of supported lipid bilayers (SLBs). In SLBs, the phase state of a lipid bilayer strongly depends on the presence of molecules such as cholesterol, ceramide, and physical parameters such as temperature. Cholesterol is a key molecule of biological membranes and it exerts condensing effect on lipid bilayers. In this paper, we demonstrate the influence of excess cholesterol content on a supported lipid bilayer of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) (fluid-phase) using spectroscopic ellipsometry (SE) and coarse-grained (CG) molecular dynamics (MD) simulations. The results show the condensation effect due to cholesterol addition up to 30% and interleaflet decoupling at excess cholesterol beyond 30%. SE results show the separation of individual leaflets of the bilayer and influence of cholesterol on the biophysical properties such as thickness and optical index. CG simulations were performed at different ratios of DOPC:cholesterol mixtures to explore cholesterol-driven bilayer properties and stability. The simulations displayed the accumulation of cholesterol molecules at the interface of the lower and upper leaflets of the bilayer, thus leading to undulations in the bilayer. This work reports the successful application of SE technique to study lipid-cholesterol interactions for the first time.

Nanoscale Probing of Cholesterol-Rich Domains in Single Bilayer Dimyristoyl-Phosphocholine Membranes Using Near-Field Spectroscopic Imaging

Cholesterol is believed to induce the formation of membrane domains, "rafts", which are implicated in a range of natural and pathologic membrane processes. Therefore, it is important to understand the role that cholesterol plays in the formation of these structures. Here, we use label-free spectroscopic imaging to investigate cholesterol fractioning in supported bilayer membranes at nanoscale. Scattering-type scanning near-field optical microscopy (s-SNOM) was used to visualize the formation of cholesterol-induced domains in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes. Our results revealed the coexistence of phase separated domains in DMPC lipids with 10 mol % cholesterol content, whereas a mostly homogeneous bilayer was found at low (5 mol %) and high (15 mol %) cholesterol content. Near-field nano-FTIR spectroscopy was used to identify the cholesterol-rich domains based on their qualitative chemical compositions. It was determined that cholesterol binds to phosphodiester and alkyl glycerol ester moieties, likely via hydrogen bonding of the alcohol to either of the ester oxygens. The results also confirm the existence of an ideal cholesterol-lipid mixture ratio (∼15:85) with a geometrically defined packing. At lower cholesterol content there is phase separation between liquid ordered and almost neat DMPC domains. Thus, the liquid ordered phase exists at an energy minimum at a given lipid-cholesterol ratio.

Spatial Dependence of Non-Gaussian Diffusion of Nanoparticles in Free-Standing Thin Polymer Films

The addition of nanoparticles (NPs) to a free-standing polymer film affects the properties of the film such as viscosity and glass transition temperature. Recent experiments, for example, showed that the glass transition temperature of thin polymer films was dependent on how NPs were distributed within the polymer films. However, the spatial arrangement of NPs in free-standing polymer films and its effect on the diffusion of NPs and polymers remain elusive at a molecular level. In this study, we employ generic coarse-grained models for polymers and NPs and perform extensive molecular dynamics simulations to investigate the diffusion of polymers and NPs in free-standing thin polymer films. We find that small NPs are likely to stay at the interfacial region of the polymer film, while large NPs tend to stay at the center of the film. On the other hand, as the interaction between a NP and a monomer becomes more attractive, the NP is more likely to be placed at the film center. The diffusion of monomers slows down slightly as more NPs are added to the film. Interestingly, the NP diffusion is dependent strongly on the spatial arrangement of the NPs: NPs at the interfacial region diffuse faster and undergo more non-Gaussian diffusion than NPs at the film center, which implies that the interfacial region would be more mobile and dynamically heterogeneous than the film center. We also find that the mechanism for non-Gaussian diffusion of NPs at the film center differs from that at the interfacial region and that the NP diffusion would reflect the local viscosity of the polymer films.

Anomalous Dynamics of in Vivo Cargo Delivery by Motor Protein Multiplexes

Vesicle transport conducted by motor protein multiplexes (MPMs), which is ubiquitous among eukaryotes, shows anomalous and stochastic dynamics qualitatively different from the dynamics of thermal motion and artificial active matter; the relationship between in vivo vesicle-delivery dynamics and the underlying physicochemical processes is not yet quantitatively understood. Addressing this issue, we perform accurate tracking of individual vesicles, containing upconverting nanoparticles, transported by kinesin-dynein-multiplexes along axonal microtubules. The mean-square-displacement of vesicles along the microtubule exhibits unusual dynamic phase transitions that are seemingly inconsistent with the scaling behavior of the mean-first-passage time over the travel length. These paradoxical results and the vesicle displacement distribution are quantitatively explained and predicted by a multimode MPM model, developed in the current work, where ATP-hydrolysis-coupled motion of MPM has both unidirectional and bidirectional modes.

Rapid diffusion of cholesterol along polyunsaturated membranes via deep dives

Cholesterol regulates the function of membrane proteins either via direct or membrane-mediated effects. Therefore, its ready availability is crucial for many protein-governed cellular processes. Recent studies suggest that cholesterol can partition to the core of polyunsaturated membranes, where cholesterol binding sites of many membrane proteins are also located. This core region is characterized by a lower viscosity. Therefore, we hypothesized that cholesterol partitioning into the membrane interior increases the rate of its diffusion in polyunsaturated membrane environments. We studied the behavior of cholesterol in membranes with increasing level of lipid chain unsaturation using a combination of atomistic and coarse-grained molecular dynamics simulations. Our simulations suggest a strong correlation between entropy-driven enhanced cholesterol partitioning to the membrane core and its faster lateral diffusion, which indicates that the less viscous membrane core indeed provides an efficient means for cholesterol movement in polyunsaturated membrane environments.