Component and state separation in DMPC/DSPC lipid bilayers: a Monte Carlo simulation study

Microfluidic-assisted fabrication of phosphatidylcholine-based liposomes for controlled drug delivery of chemotherapeutics

Microfluidic enables precise control over the continuous mixing of fluid phases at the micrometre scale, aiming to optimize the processing parameters and to facilitate scale-up feasibility. The optimization of parameters to obtain monodispersed drug-loaded liposomes however is challenging. In this work, two phosphatidylcholines (PC) differing in acyl chain length were selected, and used to control the release of the chemotherapeutic agent doxorubicin hydrochloride, an effective drug used to treat breast cancer. Microfluidics was used to rapidly screen manufacturing parameters and PC formulations to obtain monodispersed unilamellar liposomal formulations with a reproducible size (i.e. < 200 nm). Cholesterol was included in all liposomal formulations; some formulations also contained DMPC(1,2-dimyristoyl-sn-glycero-3-phosphocholine) and/or DSPC(1,2-distearoyl-sn-glycero-3-phosphocholine). Systematic variations in microfluidics total flow rate (TFR) settings were performed, while keeping a constant flow rate ratio (FRR). A total of six PC-based liposomes were fabricated using the optimal manufacturing parameters (TFR 500 μL/min, FRR 0.1) for the production of reproducible, stable liposome formulations with a narrow size distribution. Liposomes actively encapsulating doxorubicin exhibited high encapsulation efficiencies (>80%) for most of the six formulations, and sustained drug release profiles in vitro over 48 h. Drug release profiles varied as a function of the DMPC/DSPC mol content in the lipid bilayer, with DMPC-based liposomes exhibiting a sustained release of doxorubicin when compared to DSPC liposomes. The PC-based liposomes, with a slower release of doxorubicin, were tested in vitro, as to investigate their cytotoxic activity against three human breast cancer cell lines: the non-metastatic ER+/PR + MCF7 cells, the triple-negative aggressive MDA-MB 231 cells, and the metastatic HER2-overexpressing/PR + BT474 cells. Similar cytotoxicity levels to that of free doxorubicin were reported for DMPC₅ and DMPC₃ binary liposomes (IC₅₀ ~ 1 μM), whereas liposomes composed of a single PC were less cytotoxic (IC₅₀ ~ 3-4 μM). These results highlight that microfluidics is suitable for the manufacture of monodispersed and size-specific PC-based liposomes in a controlled single-step; furthermore, selected PC-based liposome represent promising nanomedicines for the prolonged release of chemotherapeutics, with the aim of improving outcomes for patients.

Scaling relationships for the elastic moduli and viscosity of mixed lipid membranes

The elastic and viscous properties of biological membranes play a vital role in controlling cell functions that require local reorganization of the membrane components as well as dramatic shape changes such as endocytosis, vesicular trafficking, and cell division. These properties are widely acknowledged to depend on the unique composition of lipids within the membrane, yet the effects of lipid mixing on the membrane biophysical properties remain poorly understood. Here, we present a comprehensive characterization of the structural, elastic, and viscous properties of fluid membranes composed of binary mixtures of lipids with different tail lengths. We show that the mixed lipid membrane properties are not simply additive quantities of the single-component analogs. Instead, the mixed membranes are more dynamic than either of their constituents, quantified as a decrease in their bending modulus, area compressibility modulus, and viscosity. While the enhanced dynamics are seemingly unexpected, we show that the measured moduli and viscosity for both the mixed and single-component bilayers all scale with the area per lipid and collapse onto respective master curves. This scaling links the increase in dynamics to mixing-induced changes in the lipid packing and membrane structure. More importantly, the results show that the membrane properties can be manipulated through lipid composition the same way bimodal blends of surfactants, liquid crystals, and polymers are used to engineer the mechanical properties of soft materials, with broad implications for understanding how lipid diversity relates to biomembrane function.

Liver-Targeted siRNA Lipid Nanoparticles Treat Hepatic Cirrhosis by Dual Antifibrotic and Anti-inflammatory Activities

Previous studies on the treatment of hepatic cirrhosis have been focusing on how to inhibit liver fibrosis, while ignoring liver inflammation, a key and underlying factor that promotes cirrhosis. High mobility group box-1 (HMGB1) protein, a pro-inflammatory factor and fibroblast chemokine, can promote the proliferation of hepatic stellate cells (HSCs) and the development of hepatic inflammation and fibrosis, playing a key role in cirrhosis formation. In this study, we prepared pPB peptide (C*SRNLIDC*)-modified and HMGB1-siRNA-loaded stable nucleic acid lipid nanoparticles (HMGB1-siRNA@SNALP-pPB) to effectively treat hepatic cirrhosis by their dual antifibrotic and anti-inflammatory activities. The pPB peptide-modified and heat shock protein 47 (HSP47)-siRNA-loaded stable nucleic acid lipid nanoparticles (HSP47-siRNA@SNALP-pPB), which have only an antifibrotic effect without an anti-inflammatory effect, was used as control. The results demonstrated that HMGB1-siRNA@SNALP-pPB were actively targeted to HSCs by the mediation of pPB peptide, effectively silenced the HMGB1 gene, inhibited the activation and proliferation of HSCs, reduced the release of HMGB1 protein, inhibited collagen deposition and fibrosis formation in the liver, and significantly prolonged the survival time of cirrhotic mice models. HMGB1-siRNA@SNALP-pPB showed a stronger therapeutic effect on liver cirrhosis than HSP47-siRNA@SNALP-pPB. This study provides an actively targeted siRNA delivery system for cirrhosis treatment based on the dual antifibrotic and anti-inflammatory effects. In addition, this study clarified the role of inflammatory problems in cirrhosis treatment in addition to liver fibrosis, providing a useful idea and scientific basis for the development of cirrhosis treatment strategies in the future.

To Hop or not to Hop: Exceptions in the FCS Diffusion Law

Diffusion obstacles in membranes have not been directly visualized because of fast membrane dynamics and the occurrence of subresolution molecular complexes. To understand the obstacle characteristics, mobility-based methods are often used as an indirect way of assessing the membrane structure. Molecular movement in biological plasma membranes is often characterized by anomalous diffusion, but the exact underlying mechanisms are still elusive. Imaging total internal reflection fluorescence correlation spectroscopy (ITIR-FCS) is a well-established mobility-based method that provides spatially resolved diffusion coefficient maps and is combined with FCS diffusion law analysis to examine subresolution membrane organization. In recent years, although FCS diffusion law analysis has been instrumental in providing new insights into the membrane structure below the optical diffraction limit, there are certain exceptions and anomalies that require further clarification. To this end, we correlate the membrane structural features imaged by atomic force microscopy (AFM) with the dynamics measured using ITIR-FCS. We perform ITIR-FCS measurements on supported lipid bilayers (SLBs) of various lipid compositions to characterize the anomalous diffusion of lipid molecules in distinct obstacle configurations, along with the high-resolution imaging of the membrane structures with AFM. Furthermore, we validate our experimental results by performing simulations on image grids with experimentally determined obstacle configurations. This study demonstrates that FCS diffusion law analysis is a powerful tool to determine membrane heterogeneities implied from dynamics measurements. Our results corroborate the commonly accepted interpretations of imaging FCS diffusion law analysis, and we show that exceptions happen when domains reach the percolation threshold in a biphasic membrane and a network of domains behaves rather like a meshwork, resulting in hop diffusion.

Lipid Organization in Mixed Lipid Membranes Driven by Intrinsic Curvature Difference

Laurdan fluorescence, novel spectral fitting, and dynamic light scattering were combined to determine lateral lipid organization in mixed lipid membranes of the oxidized lipid, 1-palmitoyl-2-glutaryl-sn-glycero-3-phosphocholine (PGPC), and each of the three bilayer lipids, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and 1-palmitoyl-2-oleoylphosphatidylcholine (POPC). Second harmonic spectra were computed to determine the number of elementary emissions present. All mixtures indicated two emissions. Accordingly, spectra were fit to two log-normal distributions. Changes with PGPC mole fraction, X_PGPC, of the area of the shorter wavelength line and of dynamic light scattering-derived aggregate sizes show that: DPPC and PGPC form component-separated mixed vesicles for X_PGPC ≤ 0.2 and coexisting vesicles and micelles for X_PGPC > 0.2 in gel and liquid-ordered phases and for all X_PGPC in the liquid-disordered phase; POPC and PGPC form randomly mixed vesicles for X_PGPC ≤ 0.2 and component-separated mixed vesicles for X_PGPC > 0.2. DOPC and PGPC separate into vesicles and micelles. Component segregation is due to unstable inhomogeneous membrane curvature stemming from lipid-specific intrinsic curvature differences between mixing molecules. PGPC is inverse cone-shaped because its truncated tail with a terminal polar group points into the interface. It is similar to and mixes with POPC, also an inverse cone because of mobility of its unsaturated tail. PGPC is least similar to DOPC because mobilities of both unsaturated tails confer a cone shape to DOPC, and PGPC separates form DOPC. DPPC and PGPC do not mix in the liquid-disordered phase because mobility of both tails in this phase renders DPPC a cone. DPPC is a cylinder in the gel phase and of moderate similarity to PGPC and mixes moderately with PGPC.

Scaling of lipid membrane rigidity with domain area fraction

Biological membranes are highly heterogeneous in composition which in turn leads to local variations in the physical properties. Here we quantify how heterogeneity in stiffness determines the effective bending modulus, κ_eff, of model phase-separated membranes with coexisting soft fluid and rigid gel domains. We find that the temperature- and composition- dependent trends in membrane rigidity collapse onto a single curve, such that κ_eff directly scales with the area fraction of the rigid gel domains. Using no adjustable parameters, the measurements are found to agree with theoretical predictions for inhomogenous membranes and indicate that κ_eff is sensitive to the lateral distribution of the rigid phase within the membrane. This key finding confirms that the properties of heterogeneous membranes can be quantitatively predicted if the area fraction and properties of the individual phases are known.

Gramicidin Lateral Distribution in Phospholipid Membranes: Fluorescence Phasor Plots and Statistical Mechanical Model

When using small mole fraction increments to study gramicidins in phospholipid membranes, we found that the phasor dots of intrinsic fluorescence of gramicidin D and gramicidin A in dimyristoyl-sn-glycero-3-phosphocholine (DMPC) unilamellar and multilamellar vesicles exhibit a biphasic change with peptide content at 0.143 gramicidin mole fraction. To understand this phenomenon, we developed a statistical mechanical model of gramicidin/DMPC mixtures. Our model assumes a sludge-like mixture of fluid phase and aggregates of rigid clusters. In the fluid phase, gramicidin monomers are randomly distributed. A rigid cluster is formed by a gramicidin dimer and DMPC molecules that are condensed to the dimer, following particular stoichiometries (critical gramicidin mole fractions, Xcr including 0.143). Rigid clusters form aggregates in which gramicidin dimers are regularly distributed, in some cases, even to superlattices. At Xcr, the size of cluster aggregates and regular distributions reach a local maximum. Before a similar model was developed for cholesterol/DMPC mixtures (Sugar and Chong (2012) J. Am. Chem. Soc. 134, 1164⁻1171) and here the similarities and differences are discussed between these two models.

Neurotransmitter-Loaded Nanocapsule Triggers On-Demand Muscle Relaxation in Living Organism

This paper reports the on-demand artificial muscle relaxation using a thermosensitive liposome encapsulating γ-aminobutyric acid (GABA) inhibitory neurotransmitter. Muscle relaxation is not feasible in principle, although muscle contraction can be easily induced by electrical stimulation. Herein, thermosensitive liposomes (phase transition temperature = 40 °C) were synthesized to encapsulate GABA and were injected into a leg of a living beetle. The leg was wrapped around by a Ni-Cr wire heater integrated with a thermocouple to enable the feedback control and to manipulate the leg temperature. The injected leg was temporarily immobilized by heating it up to 45 °C. The leg did not swing even by electrically stimulating the leg muscle. Subsequently, the leg recovered to swing. The result indicates that GABA was released from liposomes and fed to the leg muscle, enabling temporal muscle relaxation.

Lipid- and Cholesterol-Mediated Time-Scale-Specific Modulation of the Outer Membrane Protein X Dynamics in Lipid Bilayers

Membrane protein function fundamentally depends on lipid-bilayer fluidity and the composition of the biological membrane. Although dynamic interdependencies between membrane proteins and the surrounding lipids are suspected, a detailed description is still missing. To uncover lipid-modulated membrane protein backbone dynamics, time-scale-specific NMR relaxation experiments with residue-resolution were recorded. The data revealed that lipid order, modified either biochemically or biophysically, changes the dynamics of the immersed membrane protein in a specific and time-scale-dependent manner. A temperature-dependent dynamics analysis furthermore suggests a direct coupling between lipid and protein dynamics in the picosecond-nanosecond, microsecond, and millisecond time scales, caused by the lipid's trans-gauche isomerization, the segmental and rotational motion of lipids, and the fluidity of the lipid phase, respectively. These observations provide evidence of a direct modulatory capability of the membrane to regulate protein function through lipid dynamics ranging from picoseconds to milliseconds.

Thermal Response Analysis of Phospholipid Bilayers Using Ellipsometric Techniques

Biomimetic planar artificial membranes have been widely studied due to their multiple applications in several research fields. Their humectation and thermal response are crucial for reaching stability; these characteristics are related to the molecular organization inside the bilayer, which is affected by the aliphatic chain length, saturations, and molecule polarity, among others. Bilayer stability becomes a fundamental factor when technological devices are developed—like biosensors—based on those systems. Thermal studies were performed for different types of phosphatidylcholine (PC) molecules: two pure PC bilayers and four binary PC mixtures. These analyses were carried out through the detection of slight changes in their optical and structural parameters via Ellipsometry and Surface Plasmon Resonance (SPR) techniques. Phospholipid bilayers were prepared by Langmuir-Blodgett technique and deposited over a hydrophilic silicon wafer. Their molecular inclination degree, mobility, and stability of the different phases were detected and analyzed through bilayer thickness changes and their optical phase-amplitude response. Results show that certain binary lipid mixtures—with differences in its aliphatic chain length—present a co-existence of two thermal responses due to non-ideal mixing.

Multiscale Modeling of Complex Formation and CD80 Depletion during Immune Synapse Development

The mechanisms that discriminate self- and foreign antigen before T cell activation are unresolved. As part of the immune system's adaptive response to specific infections or neoplasms, antigen-presenting cells (APC) and effector T cells form transcellular molecular complexes. CTLA4 expression on regulatory or effector T cells reduces T cell activation. The CTLA4 transendocytosis hypothesis proposes that CTLA4 depletes CD80 and CD86 proteins from the APC membrane, rendering the APC incapable of activating T cells. We developed a multiscale spatiotemporal model for the interaction of a T cell and APC. Formation of the immune complex between T cell and APC starts with formation of the transmembrane complexes between the major histocompatibility complex and the T cell receptor (Signal 1) and between CD80 or CD86 and CD28 (Signal 2) at the opposing membrane surfaces of the interacting cells. By 0.01 s after contact simulation, an increasing concentration gradient of the free membrane proteins develops between the opposing surfaces and spherical parts of each cell's membrane, reaching a maximum at ∼30 s. Over several hours, diffusion across the gradient equalizes the free protein concentrations. During this phase, CTLA4 surface expression and its complexation with CD80/CD86 cause internalization and degradation of CD80/CD86. The simulation results show reasonable agreement with reported experimental data and indicate that key molecular processes take place over a very broad timescale, covering five orders of magnitude. Besides the fast complexation reactions, diffusion-limited processes, especially lateral diffusion in cell membranes and geometrical constraints, considerably slow down evolution of the synapse. Our results are consistent with the CTLA4 transendocytosis hypothesis and suggest the importance of lateral diffusion of surface proteins in contributing to a gradual increase in Signal 1 and Signal 2.

The use of zeta potential as a tool to study phase transitions in binary phosphatidylcholines mixtures

Temperature dependence of the zeta potential (ZP) is proposed as a tool to analyze the thermotropic behavior of unilamellar liposomes prepared from binary mixtures of phosphatidylcholines in the absence or presence of ions in aqueous suspensions. Since the lipid phase transition influences the surface potential of the liposome reflecting a sharp change in the ZP during the transition, it is proposed as a screening method for transition temperatures in complex systems, given its high sensitivity and small amount of sample required, that is, 70% less than that required in the use of conventional calorimeters. The sensitivity is also reflected in the pre-transition detection in the presence of ions. Plots of phase boundaries for these mixed-lipid vesicles were constructed by plotting the delimiting temperatures of both main phase transition and pre-transition vs. the lipid composition of the vesicle. Differential scanning calorimetry (DSC) studies, although subject to uncertainties in interpretation due to broad bands in lipid mixtures, allowed the validation of the temperature dependence of the ZP method for determining the phase transition and pre-transition temperatures. The system chosen was dipalmitoylphosphatidylcholine/dimyristoyl phosphatidylcholine (DMPC/DPPC), the most common combination in biological membranes. This work may be considered as a starting point for further research into more complex lipid mixtures with functional biological importance.

Tuning membrane thickness fluctuations in model lipid bilayers

Membrane thickness fluctuations have been associated with a variety of critical membrane phenomena, such as cellular exchange, pore formation, and protein binding, which are intimately related to cell functionality and effective pharmaceuticals. Therefore, understanding how these fluctuations are controlled can remarkably impact medical applications involving selective macromolecule binding and efficient cellular drug intake. Interestingly, previous reports on single-component bilayers show almost identical thickness fluctuation patterns for all investigated lipid tail-lengths, with similar temperature-independent membrane thickness fluctuation amplitude in the fluid phase and a rapid suppression of fluctuations upon transition to the gel phase. Presumably, in vivo functions require a tunability of these parameters, suggesting that more complex model systems are necessary. In this study, we explore lipid tail-length mismatch as a regulator for membrane fluctuations. Unilamellar vesicles of an equimolar mixture of dimyristoylphosphatidylcholine and distearoylphosphatidylcholine molecules, with different tail-lengths and melting transition temperatures, are used as a model system for this next level of complexity. Indeed, this binary system exhibits a significant response of membrane dynamics to thermal variations. The system also suggests a decoupling of the amplitude and the relaxation time of the membrane thickness fluctuations, implying a potential for independent control of these two key parameters.

Main transition in the Pink membrane model: finite-size scaling and the influence of surface roughness

We consider the main transition in single-component membranes using computer simulations of the Pink model [D. A. Pink et al., Biochemistry 19, 349 (1980)]. We first show that the accepted parameters of the Pink model yield a main transition temperature that is systematically below experimental values. This resolves an issue that was first pointed out by Corvera and co-workers [Phys. Rev. E 47, 696 (1993)]. In order to yield the correct transition temperature, the strength of the van der Waals coupling in the Pink model must be increased; by using finite-size scaling, a set of optimal values is proposed. We also provide finite-size scaling evidence that the Pink model belongs to the universality class of the two-dimensional Ising model. This finding holds irrespective of the number of conformational states. Finally, we address the main transition in the presence of quenched disorder, which may arise in situations where the membrane is deposited on a rough support. In this case, we observe a stable multidomain structure of gel and fluid domains, and the absence of a sharp transition in the thermodynamic limit.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: theory versus experiment

Aqueous dispersions of phosphatidylglycerol (PG) lipids may present an anomalous chain-melting transition at low ionic strengths, as seen by different experimental techniques such as calorimetry or light scattering. The anomaly disappears at high ionic strengths or for longer acyl-chain lengths. In this article, we use a statistical model for the bilayer that distinguishes both lipid chain and headgroup states in order to compare model and experimental thermotropic and electrical properties. The effective van der Waals interactions among hydrophobic chains compete with the electrostatic repulsions between polar headgroups, which may be ionized (counterion dissociated) or electrically neutral (associated with counterions). Electric degrees of freedom introduce new thermotropic charge-ordered phases in which headgroup charges may be spatially ordered, depending on the electrolyte ionic strength, introducing a new rationale for experimental data on PGs. The thermal phases presented by the model for different chain lengths, at fixed ionic strength, compare well with an experimental phase diagram constructed on the basis of differential scanning calorimetry profiles. In the case of dispersions of DMPG (dimyristoyl phosphatidylglycerol) with added monovalent salt, the model properties reproduce the main features displayed by data from differential scanning calorimetry as well as the characteristic profile for the degree of ionization of the bilayer surface across the anomalous transition region, obtained from the theoretical interpretation of electrokinetic (conductivity and electrophoretic mobility) measurements.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: a statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types--each one associated, respectively, with the polar-headgroup and the acyl-chain states--which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.

Near-critical fluctuations and cytoskeleton-assisted phase separation lead to subdiffusion in cell membranes

We address the relationship between membrane microheterogeneity and anomalous subdiffusion in cell membranes by carrying out Monte Carlo simulations of two-component lipid membranes. We find that near-critical fluctuations in the membrane lead to transient subdiffusion, while membrane-cytoskeleton interaction strongly affects phase separation, enhances subdiffusion, and eventually leads to hop diffusion of lipids. Thus, we present a minimum realistic model for membrane rafts showing the features of both microscopic phase separation and subdiffusion.

Determination of local diffusion properties in heterogeneous biomaterials

The coupling between structure and diffusion properties is essential for the functionality of heterogeneous biomaterials. Structural heterogeneity is defined and its implications for time-dependent diffusion are discussed in detail. The effect of structural heterogeneity in biomaterials on diffusion and the relevance of length scales are exemplified with regard to different biomaterials such as gels, emulsions, phase separated biopolymer mixtures and chocolate. Different diffusion measurement techniques for determination of diffusion properties at different length and time scales are presented. The interplay between local and global diffusion is discussed. New measurement techniques have emerged that enable simultaneous determination of both structure and local diffusion properties. Special emphasis is given to fluorescence recovery after photobleaching (FRAP). The possibilities of FRAP at a conceptual level is presented. The method of FRAP is briefly reviewed and its use in heterogeneous biomaterials, at barriers and during dynamic changes of the structure is discussed.

Atomistic simulations of mixed-lipid bilayers in gel and fluid phases

The slow rate of diffusive mixing poses a challenge for molecular dynamics (MD) simulation studies of mixed-lipid bilayers. A mixed Monte Carlo-molecular dynamics (MC-MD) approach, which uses mutation moves to swap lipid types throughout the system within the semi-grand canonical ensemble, is here applied to a comparison of binary mixtures in the gel and liquid crystalline phases. The two lipid components modeled, distearoylphosphatidylcholine (DSPC) and dimyristoylphosphatidylcholine (DMPC), differ by four carbons in the lengths of their acyl tails and are investigated here at full hydration at a temperature (313 K) between their transition temperatures, where coexistence between a DSPC-rich gel phase and a DMPC-rich liquid crystalline phase is expected. An analysis of DSPC-DMPC mixtures in the gel phase indicates strong deviation from ideality in the thermodynamics of mixing, accompanied by a tendency of the shorter-tailed component DMPC to associate laterally and for DMPC headgroups to be displaced toward the bilayer midplane. The liquid crystal phase mixtures, in contrast, show more mild deviation from thermodynamically ideal mixing with no apparent tendency for similar lipids to cluster laterally and no difference in headgroup normal distribution profiles.

On the inner structure and topology of clusters in two-component lipid bilayers. Comparison of monomer and dimer Ising models

It has been shown on model and biological systems that membrane clusters can affect in-plane membrane reactions and can control biochemical reaction cascades. Clusters of two-component phospholipid bilayers have been simulated by two Ising-type lattice models: the monomer and the dimer model. In each model the plane of one layer of the bilayer is represented by a triangular lattice, each site of which is occupied by an acyl chain of either a component 1 or a component 2 lipid molecule. The dimer model assumes that pairs of acyl chains (lipid molecules) are permanently connected, forming dimers on the lattice, while in the case of the monomer model this covalent connection between acyl chains is ignored. Phase diagrams of two-component phospholipid bilayers were successfully calculated by both models. In this work, we use Monte Carlo techniques to calculate thermodynamic averages of global and local characteristics of the largest component 2 cluster (such as outer/inner perimeter, percolation, minimal linear size, and local density) and compare the results obtained by the two models. A new method is developed to characterize the inner structure of the clusters. Each point of a cluster is classified based on its shortest distance (or depth) from the cluster's outer perimeter. Then local cluster properties, such as density, are calculated as a function of the depth. The depth analysis reveals that toward the cluster interior the average density usually decreases in midsize clusters and remains constant in very large clusters. On the basis of the simulations the following typical cluster topologies are identified at different cluster sizes and cooperativity parameter values: (i) branch-like, (ii) circular, (iii) band-like, and (iv) planar.We did not find qualitative differences between the cluster structures in the dimer and monomer model. However, at the same cluster size and cooperativity parameter value the cluster of the dimer model is more compact. The cluster properties of the dimer model are different from that of the monomer model because of the lower mixing entropy and higher formation energy of an elementary inner island.

Phase fluctuations on the micron-submicron scale in GUVs composed of a binary lipid mixture

We used a combination of imaging and fluctuation techniques to investigate the temporal evolution of gel phase domains at the onset of phase separation, as well as the correlation between domain topology and local lipid ordering in GUVs composed of a binary mixture of DPPC/DLPC 1:1. The data acquired at temperatures immediately above the transition temperature of the two lipids suggest fluctuations in the lipid organization with a lifetime <0.1 s and a characteristic length of 1.2 microm. As the temperature is decreased below the transition temperature of one of the lipids, coupling between the two leaflets of the bilayer is observed to begin within the first five minutes after the onset of phase separation. However, domains confined to only one leaflet can be found during the first 45-50 min after the onset of phase separation. Our analysis using a two-state model (liquid and gel) indicates that for the first 45-50 min from the onset of phase separation the two lipid phases do not strongly influence the phase behavior of each other on the micron-length scale. At longer times, behavior that deviates from the two-state model is observed and appears to be correlated to domain morphology.

Lateral diffusion of molecules in two-component lipid bilayer: a Monte Carlo simulation study

Lateral diffusion of membrane components makes possible any in-plane membrane reaction and has a key role in signaling in cell membranes. In this report the equilibrium lateral diffusion of intrinsic molecules in an equimolar DMPC/DSPC mixture is simulated using a thoroughly tested two-state model of two-component phospholipid bilayers. The model has been successful in calculating the excess heat capacity function, the most frequent center-to-center distances between DSPC clusters, and the fractal dimensions of gel clusters (Sugar, I. P., Thompson, T. E., Biltonen, R. L. Biophys. J. 1999, 76, 2099-2110). In the gel/fluid mixed phase region, a diffusing intrinsic molecule may change its state from fluid to gel (or from gel to fluid) at any time. A common characterization of the diffusion of intrinsic molecules is given by the simulated average first-passage time curves. We find that these curves can be described as power functions containing two parameters, alpha and beta, except near the percolation threshold of gel/fluid or compositional clusters. We find also that the intrinsic molecules are involved in approximately normal diffusion, i.e., beta approximately 2 in the extreme gel and fluid phase regions, while in the gel/fluid and gel/gel mixed phase regions the diffusion is anomalous, i.e., beta not equal 2. In the mixed phase regions, when the initial local state of the diffusing molecule is not specified, each component is involved in sub-diffusion (beta > 2). In the gel/fluid mixed phase region molecules situated initially inside a fluid cluster are involved in sub-diffusion, but DMPC molecules situated initially inside a gel cluster are involved in super-diffusion (beta < 2). The possibility of anomalous diffusion in membranes apparently arises because the diffusing molecule visits a variety of different environments characterized by its relative proximity to various membrane components. The diffusion is actually anomalous when the components of the bilayer are nonrandomly distributed. The deviation from random distribution is strongly correlated with beta. Similar to the results of the NMR experiments, the calculated relative diffusion coefficient continuously decreases in the gel/fluid mixed phase region with decreasing temperature. In apparent contradiction, diffusion measured by fluorescence recovery after photobleaching (FRAP) demonstrates the existence of a threshold temperature, below which long-range diffusion of FRAP probe molecules is essentially blocked. This threshold temperature is highly correlated with the percolation temperature of gel clusters.

Phase behavior and dynamic heterogeneities in lipids: A coarse-grained simulation study of DPPC–DPPE mixtures

The coarse-grained Marrink-model for biomembrane simulation is used to study mixtures of dipalmitoylphosphatidylcholine (DPPC) and dipalmitoylphosphatidylethanolamine (DPPE) at various concentrations and temperatures. At high temperatures close to ideal mixing is observed. In the low temperature ordered phase dynamic heterogeneities are identified under some conditions. These are correlated with heterogeneities in the local order and define local neighborhoods.

Equilibrium distributions of dipalmitoyl phosphatidylcholine and dilauroyl phosphatidylcholine in a mixed lipid bilayer: atomistic semigrand canonical ensemble simulations

Conventional molecular dynamics (MD) simulations are seriously limited by the slow rate of diffusive mixing in their ability to predict lateral distributions of different lipid types within mixed-lipid bilayers using atomistic models. A method to overcome this limitation, using configuration-bias Monte Carlo (MC) "mutation" moves to transform lipids from one type to another in dynamic equilibrium, is demonstrated in binary fluid-phase mixtures of lipids whose tails differ in length by four carbons. The hybrid MC-MD method operates within a semigrand canonical ensemble, so that an equilibrium composition of the mixture is determined by a constant difference in chemical potential (Delta(mu)) chosen for the components. Within several nanoseconds, bilayer structures initiated as pure dipalmitoyl phosphatidylcholine (DPPC) or pure dilauroyl phosphatidylcholine (DLPC) converge to a common composition and structure in independent simulations conducted at the same Delta(mu). Trends in bilayer thickness, area per lipid, density distributions across the bilayer, and order parameters have been investigated at three mixture compositions and compared with results from the pure bilayers at 323 K. The mixtures showed a moderate increase in DPPC acyl tail sites crossing the bilayer midplane relative to pure DPPC. Correlations between lateral positions of the two lipid types within or across the bilayer were found to be weak or absent. While the lateral distribution is consistent with nearly ideal mixing, the dependence of composition on Delta(mu) indicates a positive excess free energy of mixing.

Karhunen-Loeve analysis for pattern description in phase separated lipid bilayer systems

Karhunen-Loeve analysis, a special variant of principal component analysis, is used to describe and analyze the dynamics of self-assembled pattern formation in a mixed phospholipid bilayer. The dominant modes of the evolving heterogeneities in density and dynamics are elucidated. At low temperatures the evolution of patterns can be followed by the principal modes of the systems. We find that the higher modes only evolve after the dominant modes have been established. At high temperatures no such dominant modes are found. So, a clear descriptor of an evolving self-assembled pattern can be identified and its time evolution can be monitored. This analysis suggests, additionally, a new way of determining the equilibration decision in complex systems.

Thermotropic phase behaviour of the pseudobinary mixtures of DPPC/C12E5 and DMPC/C12E5 determined by differential scanning calorimetry and ultrasonic velocimetry

The present paper reports on the phase behaviour of the pseudobinary aqueous mixtures of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC)/pentaethylene glycol monododecyl ether (C12E5) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine monohydrate (DMPC)/C12E5. Both systems exhibit a variety of mesophases, such as lamellar gel, liquid crystalline and micellar phases. The phase diagrams show peritectic and eutectic behaviours. The existence of a compound complex is established. From the phase diagrams, the temperature dependence of the solubilisation parameters is obtained. The phase diagrams, especially with respect to the solubilisation process were qualitatively explained assuming that the packing of the constituents plays a dominating role. Finally, differential scanning calorimetry and ultrasonic velocimetry are compared concerning their potentials to determine characteristics of phase transitions in pseudobinary phospholipid/surfactant mixtures.

Multi-scale modeling of phase separation in mixed lipid bilayers

An approach to bridging the phenomenological field theory description of phase separation in binary mixed lipid bilayers with coarse-grained (CG) molecular dynamics (MD) simulation is presented. CG MD simulation is carried out for a 1:1 dipalmitoylphosphatidylcholine/dipalmitoylphosphatidylethanolamine lipid mixture at the liquid-gel phase coexistence condition. The liquid-gel phase separation can be characterized by the bilayer thickness, area per lipid molecule, and orientation parameter of the lipid tails. After a local order parameter is defined using the lipid tail bond orientation parameter, the CG MD data are bridged to a mesoscopic model based on the phenomenological Landau-Ginzberg free-energy functional. All parameters in this mesoscopic model are defined from the information of the phase boundary structure and the distributions of the order parameter in the liquid and gel phases. It is found that the mesoscopic model reproduces the equilibrium properties of the system very well, including collective fluctuations in both phases, spatial correlation functions of the order parameter, and the line tension. The possibility of using a time-dependent Landau-Ginzberg model to mimic the phase-separation dynamics is also investigated, using the relaxation time constant obtained by fitting the time-dependent correlation functions of the order parameter.

Temperature and pressure dependent growth and morphology of DMPC/DSPC domains studied by Brewster angle microscopy

In this work, the temperature and pressure dependent growth of domains in DMPC/DSPC monolayers at various molar ratios was studied by Brewster angle microscopy. Upon compression, roughly discoidal domains with some branching are formed. Further compression leads to an increase in both the number and the average size of the domains, which range between ca. 5 and 20 microm. The isobaric heating of the monolayers results in a gradual decrease of the domain size until their disappearance. The size and morphology of the domains depend not only on equilibrium parameters such as temperature, pressure and composition, but appear to be also strongly dependent on non-equilibrium parameters such as the rate of perturbation. The comparison between our results and those previously published for bilayers allows us to infer that the growth behaviour in monolayers can be qualitatively but not quantitatively extrapolated to bilayers.

Simulation of domain formation in DLPC-DSPC mixed bilayers

Binary mixtures of two phosphatidylcholines of different chain lengths are simulated in the bilayer state. We find a phase transition between a liquid state and a gel state at all concentrations. This phase transition is characterized by the area per lipid headgroup, the order parameter, and a change in dynamics. At concentrations with a majority of the longer lipid, we find phase separation into a gel and a liquid state in a small temperature window. This leads to a strong dynamic heterogeneity. Experimental phase transition temperatures are reproduced semiquantitatively. We see a clear shift in the phase transition to higher temperatures with increasing concentration of the longer lipid.

Effect of acyl chain mismatch on the contact mechanics of two-component phospholipid vesicle during main phase transition

It has been recently demonstrated that acyl chain mismatch of phospholipid bilayer composed of a binary lipid mixture induces component formation on the lateral plane of the bilayer [Biophys. J. 83 (2002) 1820-1883]. In this report, the contact mechanics of unilamellar vesicles composed of binary dimyristoyl-phosphatidylcholine (DMPC)/dipalmitoyl-phosphocholine (DPPC) mixtures on fused silica and amino-modified substrates is simultaneously probed by confocal-reflectance interference contrast microscopy (C-RICM) and cross-polarized light microscopy during gel to liquid crystalline transition of the lipid bilayer. C-RICM results indicate that the average degree of vesicle deformation for DMPC-rich and DPPC-rich vesicles adhering on fused silica substrate is increased by 30% and 14%, respectively, in comparison with that in pure DMPC and DPPC vesicles. Also, lateral heterogeneity induced by acyl chain mismatch increases the average magnitude of adhesion energy in DMPC-rich and DPPC-rich vesicles of all sizes by 6.4 times and 2.3 times, respectively. Similar modulation of adhesion mechanics induced by carbon chain difference is obtained on amino-modified substrate. Most importantly, the thermotropic transition of the mixed bilayer from gel (below T(m)) to fluid phase (above T(m)) further exemplifies the effect of acyl chain mismatch on the increases of degree of vesicle deformation and adhesion energy.