Phase behavior of model lipid bilayers

Protonation State of a Bioactive Compound Regulates Its Release from Lamellar Gel-Phase Bilayers

Lamellar gel networks (LGNs) in personal care or pharmaceutical lotions and creams provide an opaque cream appearance and a creamy texture to these products. Within the LGNs, the lamellar gel (Lβ) phase composed of regularly spaced bilayers of surfactants and long-chain fatty alcohols is predominately responsible for the unique rheological properties of the LGNs. To extend the shelf life of LGN-containing products, bioactive compounds with antimicrobial properties are often incorporated into the formulation. However, how the protonation state of the bioactive compounds regulates their release from the Lβ-phase bilayers is currently unknown. Using molecular dynamics simulations, we found that the protonated (neutral) form of cinnamic acid, a common antimicrobial food additive, has a retention ratio higher than that of its deprotonated (charged) counterpart in the Lβ-phase bilayer. From free energy calculations, we determined that not only is the protonated molecule more stable in the hydrophobic interior of the bilayer but also the formation of hydrogen-bonded dimers significantly enhances its stability within the bilayer. Thus, the protonation state has a profound impact on bioavailability of the compounds. Our results also highlight the importance of considering possible oligomeric states of molecules when performing calculations to estimate the permeability of molecules within various bilayers.

Thermal control of photothermal implants inspired by polar bear skin for the treatment of infected bone defects

Photothermal therapy (PTT) encounters challenges in addressing deep tissue infections, characterized by limited penetration or potential hyperthermal damage to surrounding tissues, initiating undesirable inflammatory cascades. Inspired by polar bear thermal regulation, we present a "bio-based endogenic thermal-adaptive booster" implant coating. This coating integrates a photothermal poly(tannic acid) (pTA) layer, mimicking the "polar bear dark skin", securely linked with anti-inflammatory dexamethasone (Dex), resembling the "secretion", and a red blood cell membrane (RBCM) layer, forming the insulating "transparent fur". The RBCM "fur" demonstrates unexpectedly superior local heat storage, amplifying the photothermal effect of the pTA "skin" by 1.30 times and boosting localized photothermal antibacterial efficiency by 1.30-fold (approximately 99%) compared to those without RBCM. Furthermore, RBCM sustains Dex release and offers additional protection against thermal inflammation, releasing Dex 1.90 times more under NIR irradiation than under non-photothermal conditions. In a rat infectious bone model, the photothermal-boosting implant coating provides a favorable biological interface and achieves a 99.97% photothermal antibacterial ratio, enhancing osseointegration without evident tissue harm, evidenced by a 2.47-fold increase in bone volume fraction and a 2.24-fold reduction in pro-inflammatory cytokines compared to those lacking a RBCM. Insights derived from cell membrane-based thermal-adaptive coatings herald a paradigm shift in efficient and safe PTT.

Temperature response of sucrose palmitate solutions: Role of ratio between monoesters and diesters

HYPOTHESIS

Aqueous solutions of long-chain water-soluble sucrose ester surfactants exhibit non-trivial response to temperature variations, revealing a peak in viscosity around 40-50 °C. While previous investigations have explored the structures within sucrose stearate systems at various constant temperatures, a comprehensive understanding of the entire temperature dependence and the underlying molecular factors, contributing to this phenomenon is currently missing.

EXPERIMENTS

Temperature dependent properties and supramolecular structures formed in aqueous solutions of commercial sucrose palmitate were examined using SAXS/WAXS, DSC, optical microscopy, rheological measurements, NMR, and cryo-TEM.

FINDINGS

The underlying mechanism governing this unusual behavior is revealed and is shown to relate to the mono- to di-esters ratio in the solutions. Solutions primarily containing sucrose monoesters (monoesters molecules ≳ 98% of all surfactant molecules) exhibit behavior typical of nonionic surfactants, with minimal changes with temperature. In contrast, the coexistence of mono- and di-esters results in the formation of discrete monodisperse diester particles and a network of partially fused diester particles at low temperature. As the temperature approaches the diesters' melting point, wormlike mixed micelles form, causing a viscosity peak. The height of this peak increases significantly with the diester concentration. Further temperature increase leads to fluidization of surfactant tails and formation of branched micelles, while excess diester molecules phase separate into distinct droplets.

Design Principles for Smart Linear Polymer Ligand Carriers with Efficient Transcellular Transport Capabilities

The surface functionalization of polymer-mediated drug/gene delivery holds immense potential for disease therapy. However, the design principles underlying the surface functionalization of polymers remain elusive. In this study, we employed computer simulations to demonstrate how the stiffness, length, density, and distribution of polymer ligands influence their penetration ability across the cell membrane. Our simulations revealed that the stiffness of polymer ligands affects their ability to transport cargo across the membrane. Increasing the stiffness of polymer ligands can promote their delivery across the membrane, particularly for larger cargoes. Furthermore, appropriately increasing the length of polymer ligands can be more conducive to assisting cargo to enter the lower layer of the membrane. Additionally, the distribution of polymer ligands on the surface of the cargo also plays a crucial role in its transport. Specifically, the one-fourth mode and stripy mode distributions of polymer ligands exhibited higher penetration ability, assisting cargoes in penetrating the membrane. These findings provide biomimetic inspiration for designing high-efficiency functionalization polymer ligands for drug/gene delivery.

Effects of shear-induced crystallization on the complex viscosity of lamellar-structured concentrated surfactant solutions

Material relationships at low temperatures were determined for concentrated surfactant solutions using a combination of rheological experiments, cross-polarized microscopy, calorimetry, and small angle X-ray scattering. A lamellar structured 70 wt% solution of sodium laureth sulfate in water was used as a model system. At cold temperatures (5 °C and 10 °C), the formation of surfactant crystals resulted in extremely high viscosity. The bulk flow behavior of multi-lamellar vesicles (20 °C) and focal conic defects (90 °C) in the lamellar phase was similar. Shear-induced crystallization at temperatures higher than the equilibrium crystallization temperature range resulted in an unusual complex viscosity peak. The effects of processing-relevant parameters including temperature, cooling time, and applied shear were investigated. Knowledge of key low-temperature structure-property-processing relationships for concentrated feedstocks is essential for the sustainable design and manufacturing of surfactant-based consumer products for applications such as cold-water laundry.

Exploring the Physical Properties of Lipid Membranes with Polyhydroxy Oxanorbornane Head Group Using NBD-Conjugated and DPH Fluorescent Probes

The present study focuses on exploring the physical properties of lipid membranes based on the polyhydroxy oxanorbornane (PH-ONB) headgroup, designed as synthetic analogues of naturally occurring archaeal lipid membranes. Specifically, we study two variants of PH-ONB headgroup-based lipids differing in the number of hydroxy groups present in the headgroup, with one having two hydroxy groups (ONB-2OH) and the other having three (ONB-3OH). These lipids form stable bilayer membranes. The study begins with a comprehensive analysis of the fluorescence characteristics of nitrobenzoxadiazole (NBD)-tagged ONB-based lipids in different solvent environments and within a model lipid membrane 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Subsequently, the physical properties of the ONB-based membranes were examined by using an NBD-tagged ONB-based probe and a commonly used extrinsic 1,6-diphenyl-1,3,5-hexatriene (DPH) fluorescent probe. The steady-state and time-resolved fluorescence properties of the NBD-tagged ONB-based probe and DPH were used to compare the physical properties of the ONB-based membranes, including polarity, fluidity, phase transition, order, hydration, location, heterogeneity, and rotational diffusion. The solid gel to liquid crystalline phase transition temperatures of ONB-2OH and ONB-3OH lipid membranes are found to be (68 ± 1) °C and (74 ± 1) °C, respectively. The variation in organization (size), fluidity, and phase transition temperature of ONB-based lipid membranes is explained by the extent of hydrogen bonding interactions between lipid head groups. ONB-based membranes exhibit characteristics similar to those of phospholipid membranes and possess a notably high phase transition temperature. These properties make them a promising and cost-effective synthetic alternative to archaeal lipid membranes with a wide range of potential applications.

Glycan-Driven Formation of Raft-Like Domains with Hierarchical Periodic Nanoarrays on Dendrimersome Synthetic Cells

The accurate spatial segregation into distinct phases within cell membranes coordinates vital biochemical processes and functionalities in living organisms. One of nature's strategies to localize reactivity is the formation of dynamic raft domains. Most raft models rely on liquid-ordered L0 phases in a liquid-disordered Ld phase lacking correlation and remaining static, often necessitating external agents for phase separation. Here, we introduce a synthetic system of bicomponent glycodendrimersomes coassembled from Janus dendrimers and Janus glycodendrimers (JGDs), where lactose-lactose interactions exclusively drive lateral organization. This mechanism results in modulated phases across two length scales, yielding raft-like microdomains featuring nanoarrays at the nanoscale. By varying the density of lactose and molecular architecture of JGDs, the nanoarray type and size, shape, and spacing of the domains were controlled. Our findings offer insight into the potential primordial origins of rudimentary raft domains and highlight the crucial role of glycans within the glycocalyx.

Brownian motion of soft particles near a fluctuating lipid bilayer

The dynamics of a soft particle suspended in a viscous fluid can be changed by the presence of an elastic boundary. Understanding the mechanisms and dynamics of soft-soft surface interactions can provide valuable insights into many important research fields, including biomedical engineering, soft robotics development, and materials science. This work investigates the anomalous transport properties of a soft nanoparticle near a visco-elastic interface, where the particle consists of a polymer assembly in the form of a micelle and the interface is represented by a lipid bilayer membrane. Mesoscopic simulations using a dissipative particle dynamics model are performed to examine the impact of micelle's proximity to the membrane on its Brownian motion. Two different sizes are considered, which correspond to ≈10-20nm in physical units. The wavelengths typically seen by the largest micelle fall within the range of wavenumbers where the Helfrich model captures fairly well the bilayer mechanical properties. Several independent simulations allowed us to compute the micelle trajectories during an observation time smaller than the diffusive time scale (whose order of magnitude is similar to the membrane relaxation time of the largest wavelengths), this time scale being hardly accessible by experiments. From the probability density function of the micelle normal position with respect to the membrane, it is observed that the position remains close to the starting position during ≈0.05τd (where τd corresponds to the diffusion time), which allowed us to compare the negative excess of mean-square displacement (MSD) to existing theories. In that time range, the MSD exhibits different behaviors along parallel and perpendicular directions. When the micelle is sufficiently close to the bilayer (its initial distance from the bilayer equals approximately twice its gyration radius), the micelle motion becomes quickly subdiffusive in the normal direction. Moreover, the temporal evolution of the micelle MSD excess in the perpendicular direction follows that of a nanoparticle near an elastic membrane. However, in the parallel direction, the MSD excess is rather similar to that of a nanoparticle near a liquid interface.

Anomalous lateral diffusion of lipids during the fluid/gel phase transition of a lipid membrane

A lipid membrane undergoes a phase transition from fluid to gel phase upon changing external thermodynamic conditions, such as decreasing temperature and increasing pressure. Extremophilic organisms face the challenge of preventing this deleterious phase transition. The main focus of their adaptive strategy is to facilitate effective temperature sensing through sensor proteins, relying on the drastic changes in packing density and membrane fluidity during the phase transition. Although the changes in packing density parameters due to the fluid/gel phase transition are studied in detail, the impact on membrane fluidity is less explored in the literature. Understanding the lateral diffusive dynamics of lipids in response to temperature, particularly during the fluid/gel phase transition, is albeit crucial. Here we have simulated the phase transition of a single component lipid membrane composed of dipalmitoylphosphatidylcholine (DPPC) lipids using a coarse-grained (CG) model and studied the changes of the structural and dynamical properties. It is observed that near the phase transition point, both fluid and gel phase domains coexist together. The dynamics remains highly non-Gaussian for a long time even when the mean square displacement reaches the Fickian regime at a much earlier time. This Fickian yet non-Gaussian diffusion (FnGD) is a characteristic of a highly heterogeneous system, previously observed for the lateral diffusion of lipids in raft mimetic membranes having liquid-ordered and liquid-disordered phases co-existing together. We have analyzed the molecular trajectories and calculated the jump-diffusion of the lipids, stemming from sudden jump translations, using a translational jump-diffusion (TJD) approach. An overwhelming contribution of the jump-diffusion of the lipids is observed suggesting anomalous diffusion of lipids during fluid/gel phase transition of the membrane. These results are important in unravelling the intricate nature of lipid diffusion during the phase transition of the membrane and open up a new possibility of investigating the most significant change of membrane properties during phase transition, which can be effectively sensed by proteins.

Understanding the phase behavior of a protobiomembrane

The rich thermotropic behavior of lipid bilayers is addressed using phenomenological theory informed by many experiments. The most recent experiment not yet addressed by theory has shown that the tilt modulus in DMPC lipid bilayers decreases dramatically as the temperature is lowered toward the main transition temperature T_{M}. It is shown that this behavior can be understood by introducing a simple free energy functional for tilt that couples to the area per molecule. This is combined with a chain melting free energy functional in which the area is the primary order parameter that is the driver of the main transition. Satisfactory agreement with experiment is achieved with values of the model parameters determined by experiments, but the transition is directly into the gel phase. The theory is then extended to include the enigmatic ripple phase by making contact with the most recent experimentally determined ripple structure.

Dynamic Processes and Mechanical Properties of Lipid-Nanoparticle Mixtures

In this study, we investigate the dynamic processes and mechanical properties of lipid nanoparticle mixtures in a melt via dissipation particle dynamic simulation. By investigating the distribution of nanoparticles in lamellar and hexagonal lipid matrices in equilibrium state and dynamic processes, we observe that the morphology of such composites depends not only on the geometric features of the lipid matrix but also on the concentration of nanoparticles. The dynamic processes are also demonstrated by calculating the average radius of gyration, which indicates the isotropic conformation of lipid molecules in the x-y plane and that the lipid chains are stretched in the z direction with the addition of nanoparticles. Meanwhile, we predict the mechanical properties of lipid-nanoparticle mixtures in lamellar structures by analyzing the interfacial tensions. Results show that the interfacial tension decreased with the increase in nanoparticle concentration. These results provide molecular-level information for the rational and a priori design of new lipid nanocomposites with ad hoc tailored properties.

The role of hydrophobic patches of de novo designed MSI-78 and VG16KRKP antimicrobial peptides on fragmenting model bilayer membranes

Antimicrobial peptides (AMPs) with cell membrane lysing capability are considered potential candidates for the development of the next generation of antibiotics. Designing novel AMPs requires an in-depth understanding of the mechanism of action of the peptides. In this work, we used various biophysical techniques including ³¹P solid-state NMR to examine the interaction of model membranes with amphipathic de novo-designed peptides. Two such peptides, MSI-78 and VG16KRKP, were designed with different hydrophobicity and positive charges. The model lipid membranes were constituted by mixing lipids of varying degrees of 'area per lipid' (APL), which directly affected the packing properties of the membrane. The observed emergence of the isotropic peak in ³¹P NMR spectra as a function of time is a consequence of the fragmentation of the membrane mediated by the peptide interaction. The factors such as the charges, overall hydrophilicity of the AMPs, as well as lipid membrane packing, contributed to the kinetics of membrane fragmentation. Furthermore, we anticipate the designed AMPs follow the carpet and toroidal pore mechanisms when lysing the cell membrane. This study highlights the significance of the effect of the overall charges and the hydrophobicity of the novel AMPs designed for antimicrobial activity.

Effects of Shape on Interaction Dynamics of Tetrahedral Nanoplastics and the Cell Membrane

Cellular uptake of nanoplastics is instrumental in their environmental accumulation and transfer to humans through the food chain. Despite extensive studies using spherical plastic nanoparticles, the influence of the morphological characteristics of environmentally released nanoplastics is understudied. Using dissipative particle dynamics simulations, we modeled the interactions between a cell membrane and hydrophobic nanotetrahedra, which feature high shape anisotropy and large surface curvature seen for environmental nanoplastics. We observe robust uptake of nanotetrahedra with sharp vertices and edges by the lipid membrane. Two local energy minimum configurations of nanotetrahedra embedded in the membrane bilayer were identified for particles of large sizes. Further analysis of particle dynamics within the membrane shows that the two interaction states exhibit distinct translational and rotational dynamics in the directions normal and parallel to the plane of the membrane. The membrane confinement significantly arrests the out-of-plane motion, resulting in caged translation and subdiffusive rotation. While the in-plane diffusion remains Brownian, we find that the translational and rotational modes decouple from each other as the particle size increases. The rotational diffusion decreases by a greater extent compared to the translational diffusion, deviating from the continuum theory predictions. These results provide fundamental insights into the shape effect on the nanoparticle dynamics in crowded lipid membranes.

The potential of antifungal peptide Sesquin as natural food preservative

Sesquin is a wide spectrum antimicrobial peptide displaying a remarkable activity on fungi. Contrarily to most antimicrobial peptides, it presents an overall negative charge. In the present study, we elucidate the molecular basis of its mode of action towards biomimetic membranes by NMR and MD experiments. While a specific recognition of phosphatidylethanolamine (PE) might explain its activity in a variety of different organisms (including bacteria), a further interaction with ergosterol accounts for its strong antifungal activity. NMR data reveal a charge gradient along its amide protons allowing the peptide to reach the membrane phosphate groups despite its negative charge. Subsequently, the peptide gets structured inside the bilayer, reducing its order. MD simulations predict that its activity is retained in conditions commonly used for food preservation: low temperatures, high pressure, or the presence of electric field pulses, making Sesquin a good candidate as food preservative.

Effect of polyols on membrane structures of liposomes: A study using small-angle X-ray scattering data and generalized indirect Fourier transformation

This study aimed to evaluate the membrane structure of distearoylphosphatidylcholine (DSPC) liposomes dispersed in water containing various types of polyols with low molecular weight such as glycerin (Gly), 1,3-butandiol (BG), and propylene glycol (PG). To clarify the detailed membrane structure, generalized indirect Fourier transformation (GIFT) analysis, which provides information about the bilayer spacing, bilayer thickness, number of lamellar layers, and membrane flexibility, was applied to small-angle X-ray scattering (SAXS) data of the present system. The GIFT results showed that the bilayer thickness of the DSPC liposomes followed the order Gly>>BG>PG. In addition, the membrane flexibility estimated by the Caille parameter was in the order Gly>>BG>PG; this result was supported by the gel-liquid crystal phase transition temperature (T_c) obtained by differential scanning calorimetry (DSC). These results, together with the Raman spectra, suggest that BG and PG incorporated into the bilayers of DSPC liposomes result in the formation of an interdigitated lamellar structure.

Biomembrane lipids: When physics and chemistry join to shape biological activity

Biomembranes constitute the first lines of defense of cells. While small molecules can often permeate cell walls in bacteria and plants, they are generally unable to penetrate the barrier constituted by the double layer of phospholipids, unless specific receptors or channels are present. Antimicrobial or cell-penetrating peptides are in fact highly specialized molecules able to bypass this barrier and even discriminate among different cell types. This capacity is made possible by the intrinsic properties of its phospholipids, their distribution between the internal and external leaflet, and their ability to mutually interact, modulating the membrane fluidity and the exposition of key headgroups. Although common phospholipids can be found in the membranes of most organisms, some are characteristic of specific cell types. Here, we review the properties of the most common lipids and describe how they interact with each other in biomembrane. We then discuss how their assembly in bilayers determines some key physical-chemical properties such as permeability, potential and phase status. Finally, we describe how the exposition of specific phospholipids determines the recognition of cell types by membrane-targeting molecules.

Selective membrane wrapping on differently sized nanoparticles regulated by clathrin assembly: A computational model

Nanoparticles (NPs) enter cells via multiple pathways, all of which are NP size dependent. Previous studies indicated that the clathrin-mediated endocytosis has different selectivity for the NP size, but the regulatory mechanism remains unclear and difficult to study at the molecular scale in vivo. By means of computer simulation, here we design membrane wrapping on differently sized NPs with mimic clathrin assembly at the opposite membrane side. With relatively large NPs readily wrapped by a pure membrane as manifested, clathrin modulates the process and tunes the size selectivity. Although finite curvature can be generated by cage-like clathrin assembly to facilitate membrane wrapping on relatively small NPs, the clathrin assemblage has a certain range of size, which is mismatched with larger NPs. Besides, the local membrane patch is rigidified by clathrin to increase the barrier of membrane wrapping on larger NPs. Competition of these items determines whether membrane wrapping on NPs is promoted or suppressed, and can be tuned by the NP-membrane adhesion strength, clathrin concentration, and inter-NP distance. Our results highlight the significance of complex environment in altering the nature of NP interaction with cell membranes, and are expected to help design NPs for biomedical applications requiring precise control of NP uptake or cell membrane attachment.

A generic force field for simulating native protein structures using dissipative particle dynamics

A coarse-grained force field for molecular dynamics simulations of the native structures of proteins in a dissipative particle dynamics (DPD) framework is developed. The parameters for bonded interactions are derived by mapping the bonds and angles for 20 amino acids onto target distributions obtained from fully atomistic simulations in explicit solvent. A dual-basin potential is introduced for stabilizing backbone angles, to cover a wide spectrum of protein secondary structures. The backbone dihedral potential enables folding of the protein from an unfolded initial state to the folded native structure. The proposed force field is validated by evaluating the structural properties of several model peptides and proteins including the SARS-CoV-2 fusion peptide, consisting of α-helices, β-sheets, loops and turns. Detailed comparisons with fully atomistic simulations are carried out to assess the ability of the proposed force field to stabilize the different secondary structures present in proteins. The compact conformations of the native states were evident from the radius of gyration and the high intensity peaks of the root mean square deviation histograms, which were found to be within 0.4 nm. The Ramachandran-like energy landscape on the phase space of backbone angles (θ) and dihedrals (ϕ) effectively captured the conformational phase space of α-helices at ∼(ϕ = 50°,θ = 90°) and β-strands at ∼(ϕ = ±180°,θ = 90-120°). Furthermore, the residue-residue native contacts were also well reproduced by the proposed DPD model. The applicability of the model to multidomain complexes was assessed using lysozyme and a large α-helical bacterial pore-forming toxin, cytolysin A. Our study illustrates that the proposed force field is generic, and can potentially be extended for efficient in silico investigations of membrane bound polypeptides and proteins using DPD simulations.

Nanocomplex made up of antimicrobial metallo-supramolecules and model biomembranes - characterization and enhanced fluorescence

Antimicrobial pentatopic 2,2':6',2''-terpyridines that form 3-D supramolecular hexagonal prisms with Cd2+ through coordination driven self-assembly can be entrapped by lipid discoidal bicelles, composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine, 1,2-dihexanoyl-sn-glycero-3-phosphocholine and 1,2-dipalmitoyl-sn-glycero-3-phospho-(1'-rac-glycerol) lipid, forming a well-defined nanocomplex. Structural characterization performed by very small angle neutron scattering, small angle X-ray scattering and transmission electron microscopy suggests that the hexagonal prisms are preferably located at the rim of bicellar discs with the hexagonal face in parallel with the bilayers, instead of face-to-face stacking. Such a configuration reduces the π-π interaction and consequently enhances the fluorescence emission. Since novel supramolecules were reported to have antibiotic functions, this study provides insight into the interactions of antimicrobial supermolecules with lipid membranes, leading to potential theranostic applications.

Double-Chain Cationic Surfactants: Swelling, Structure, Phase Transitions and Additive Effects

Double-chain amphiphilic compounds, including surfactants and lipids, have broad significance in applications like personal care and biology. A study on the phase structures and their transitions focusing on dioctadecyldimethylammonium chloride (DODAC), used inter alia in hair conditioners, is presented. The phase behaviour is dominated by two bilayer lamellar phases, Lβ and Lα, with "solid" and "melted" alkyl chains, respectively. In particular, the study is focused on the effect of additives of different polarity on the phase transitions and structures. The main techniques used for investigation were differential scanning calorimetry (DSC) and small- and wide-angle X-ray scattering (SAXS and WAXS). From the WAXS reflections, the distance between the alkyl chains in the bilayers was obtained, and from SAXS, the thicknesses of the surfactant and water layers. The Lα phase was found to have a bilayer structure, generally found for most surfactants; a Lβ phase made up of bilayers with considerable chain tilting and interdigitation was also identified. Depending mainly on the polarity of the additives, their effects on the phase stabilities and structure vary. Compounds like urea have no significant effect, while fatty acids and fatty alcohols have significant effects, but which are quite different depending on the nonpolar part. In most cases, Lβ and Lα phases exist over wide composition ranges; certain additives induce transitions to other phases, which include cubic, reversed hexagonal liquid crystals and bicontinuous liquid phases. For a system containing additives, which induce a significant lowering of the Lβ-Lα transition, we identified the possibility of a triggered phase transition via dilution with water.

Evaluating Coarse-Grained MARTINI Force-Fields for Capturing the Ripple Phase of Lipid Membranes

Phospholipids, which are an integral component of cell membranes, exhibit a rich variety of lamellar phases modulated by temperature and composition. Molecular dynamics (MD) simulations have greatly enhanced our understanding of phospholipid membranes by capturing experimentally observed phases and phase transitions at molecular resolution. However, the ripple (P_β') membrane phase, observed as an intermediate phase below the main gel-to-liquid crystalline transition with some lipids, has been challenging to capture with MD simulations, both at all-atom and coarse-grained (CG) resolutions. Here, with an aggregate ∼2.5 μs all-atom and ∼122 μs CGMD simulations, we systematically assess the ability of six CG MARTINI 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid and water force-field (FF) variants, parametrized to capture the DPPC gel and fluid phases, for their ability to capture the P_β' phase, and compared observations with those from an all-atom FF. Upon cooling from the fluid phase to below the phase transition temperature with smaller (380-lipid) and larger (>2200-lipid) MARTINI and all-atom (CHARMM36 FF) DPPC lipid bilayers, we observed that smaller bilayers with both all-atom and MARTINI FFs sampled interdigitated P_β' and ripple-like states, respectively. However, while all-atom simulations of the larger DPPC membranes exhibited the formation of the P_β' phase, MARTINI membranes did not sample interdigitated ripple-like states at larger system sizes. We then demonstrated that the ripple-like states in smaller MARTINI membranes were kinetically trapped structures caused by finite size effects rather than being representative of true P_β' phases. We showed that a MARTINI FF variant that could capture the tilted L_β' gel phase, a prerequisite for stabilizing the P_β' phase, was unable to capture the rippled phase upon cooling. Our study reveals that the current MARTINI FFs (including MARTINI3) may require specific reparametrization of the interaction potentials to stabilize lipid interdigitation, a characteristic of the ripple phase.

Effects of Trehalose on Lipid Membranes under Rapid Cooling using All-Atom and Coarse-Grained Molecular Simulations

We investigate the effect of the cryopreservative α-α-trehalose on a model 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid membrane undergoing cooling from 350 to 250 K using all-atom (AA) and coarse-grained (CG) molecular dynamics simulation. In the AA simulations, we find that the addition of trehalose alters the L_α (liquid crystalline) to P_β (ripple) phase transition, suppressing the major domain of the P_β phase and increasing the degree of leaflet interdigitation (the minor domain) which yields a thinner membrane with a higher area per lipid. Calculation of dihedral angle distributions for the lipid tails shows a greater fraction of gauche angles in the P_β phase as trehalose concentration is increased, indicating that trehalose increases lipid disorder in the membrane. In contrast, the CG simulations transition directly from the L_α to the L_β (gel) phase upon cooling without exhibiting the P_β phase (likely due to increased lipid mobility in the CG system). Even so, the CG simulations show that the addition of trehalose clearly suppresses the L_α to L_β phase transition, demonstrating that trehalose increases lipid disorder at low temperatures for the CG system, similar to the AA. Analysis using a two-state binding model provides net affinity coefficients between trehalose and the membrane as well as trehalose partition coefficients between the membrane interface and the bulk solution for both the AA and CG systems.

Membrane Permeability and Its Activation Energies in Dependence on Analyte, Lipid, and Phase Type Obtained by the Fluorescent Artificial Receptor Membrane Assay

Time-resolved monitoring of the permeability of analytes is of utmost importance in membrane research. Existing methods are restricted to single-point determinations or flat synthetic membranes, limiting access to biologically relevant kinetic parameters (permeation rate constant, permeation coefficients). We now use the recently introduced fluorescent artificial receptor membrane assay (FARMA) as a method to monitor, in real time, the permeation of indole derivatives through liposomal membranes of different lipid compositions. This method is based on the liposomal encapsulation of a chemosensing ensemble or "fluorescent artificial receptor", consisting of 2,7-dimethyldiazapyrenium as a fluorescent dye and cucurbit[8]uril as the macrocyclic receptor, that responds to the complexation of a permeating aromatic analyte by fluorescence quenching. FARMA does not require a fluorescent labeling of the analytes and allows access to permeability coefficients in the range from 10^-8 to 10^-4 cm s^-1. The effect of temperature on the permeation rate of a series of indole derivatives across the phospholipid membranes was studied. The activation energies for permeation through POPC/POPS phospholipid membranes were in the range of 28-96 kJ mol^-1. To study the effect of different lipid phases on the membrane permeability, we performed experiments with DPPC/DOPS vesicles, which showed a phase transition from a gel phase to a liquid-crystalline phase, where the activation energies for the permeation process were expected to show a dramatic change. Accordingly, for the permeation of the indole derivatives into the DPPC/DOPS liposomes, discontinuities were observed in the Arrhenius plots, from which the permeation activation energies for the distinct phases could be determined, for example, for tryptamine 245 kJ mol^-1 in the gel phase and 47 kJ mol^-1 in the liquid-crystalline phase.

Coarse-grained implicit solvent lipid force field with a compatible resolution to the Cα protein representation

Biological membranes have been prominent targets for coarse-grained (CG) molecular dynamics simulations. While minimal CG lipid models with three beads per lipid and quantitative CG lipid models with >10 beads per lipid have been well studied, in between them, CG lipid models with a compatible resolution to residue-level CG protein models are much less developed. Here, we extended a previously developed three-bead lipid model into a five-bead model and parameterized it for two phospholipids, POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine). The developed model, iSoLF, reproduced the area per lipid, hydrophobic thickness, and phase behaviors of the target phospholipid bilayer membranes at the physiological temperature. The model POPC and DPPC membranes were in liquid and gel phases, respectively, in accordance with experiments. We further examined the spontaneous formation of a membrane bilayer, the temperature dependence of physical properties, the vesicle dynamics, and the POPC/DPPC two-component membrane dynamics of the CG lipid model, showing some promise. Once combined with standard Cα protein models, the iSoLF model will be a powerful tool to simulate large biological membrane systems made of lipids and proteins.

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.

The ENUF method-Ewald summation based on nonuniform fast Fourier transform: Implementation, parallelization, and application

Computer simulations of model systems are widely used to explore striking phenomena in promising applications spanning from physics, chemistry, biology, to materials science and engineering. The long range electrostatic interactions between charged particles constitute a prominent factor in determining structures and states of model systems. How to efficiently calculate electrostatic interactions in simulation systems subjected to partial or full periodic boundary conditions has been a grand challenging task. In the past decades, a large variety of computational schemes has been proposed, among which the Ewald summation method is the most reliable route to accurately deal with electrostatic interactions between charged particles in simulation systems. In addition, extensive efforts have been done to improve computational efficiencies of the Ewald summation based methods. Representative examples are approaches based on cutoffs, reaction fields, multi-poles, multi-grids, and particle-mesh schemes. We sketched an ENUF method, an abbreviation for the Ewald summation method based on the nonuniform fast Fourier transform technique, and have implemented this method in particle-based simulation packages to calculate electrostatic energies and forces at micro- and mesoscopic levels. Extensive computational studies of conformational properties of polyelectrolytes, dendrimer-membrane complexes, and ionic fluids demonstrated that the ENUF method and its derivatives conserve both energy and momentum to floating point accuracy, and exhibit a computational complexity of O N log N with optimal physical parameters. These ENUF based methods are attractive alternatives in molecular simulations where high accuracy and efficiency of simulation methods are needed to accelerate calculations of electrostatic interactions at extended spatiotemporal scales.

Role of the lipid bilayer in outer membrane protein folding in Gram-negative bacteria

β-Barrel outer membrane proteins (OMPs) represent the major proteinaceous component of the outer membrane (OM) of Gram-negative bacteria. These proteins perform key roles in cell structure and morphology, nutrient acquisition, colonization and invasion, and protection against external toxic threats such as antibiotics. To become functional, OMPs must fold and insert into a crowded and asymmetric OM that lacks much freely accessible lipid. This feat is accomplished in the absence of an external energy source and is thought to be driven by the high thermodynamic stability of folded OMPs in the OM. With such a stable fold, the challenge that bacteria face in assembling OMPs into the OM is how to overcome the initial energy barrier of membrane insertion. In this review, we highlight the roles of the lipid environment and the OM in modulating the OMP-folding landscape and discuss the factors that guide folding in vitro and in vivo We particularly focus on the composition, architecture, and physical properties of the OM and how an understanding of the folding properties of OMPs in vitro can help explain the challenges they encounter during folding in vivo Current models of OMP biogenesis in the cellular environment are still in flux, but the stakes for improving the accuracy of these models are high. OMP folding is an essential process in all Gram-negative bacteria, and considering the looming crisis of widespread microbial drug resistance it is an attractive target. To bring down this vital OMP-supported barrier to antibiotics, we must first understand how bacterial cells build it.

Hybrid-Supported Bilayers Formed with Mixed-Charge Surfactants on C18-Functionalized Silica Surfaces

Mixtures of cationic-anionic surfactants have been shown to spontaneously form ordered monolayers at hydrophobic-hydrophilic boundaries, including air-water and oil-water interfaces. In this work, confocal Raman microscopy is used to investigate the structure of hybrid-supported surfactant bilayers (HSSBs) formed by deposition of a distal leaflet of mixed cationic-anionic surfactants onto a proximal leaflet of n-alkane (C₁₈) chains on the interior surfaces of chromatographic silica particles. The surface coverage of the two surfactants in a hybrid bilayer was determined from carbon analysis and the relative Raman scattering of their respective head-groups. Within the measurement uncertainty, the stoichiometric ratio of the two surfactants is one-to-one, equivalent to mixed-charge-surfactant monolayers at air-water and oil-water interfaces and consistent with the role of the head-group electrostatic interactions in their formation. When self-assembled on the hydrophobic surface, pairs of oppositely charged n-alkyl chain surfactants resemble a phospholipid (phosphatidylcholine) molecule, with its zwitterionic head-group and two hydrophobic acyl chain tails. Indeed, the structure of these hybrid-supported surfactant bilayers on C₁₈-modified silica surfaces is similar to that of hybrid-supported lipid bilayers (HSLBs) on the same supports, but with denser and more-ordered n-alkyl chains. Hybrid-supported surfactant bilayers exhibit a melting phase transition (gel to liquid-crystalline phase) with structural and energetic characteristics similar to those of hybrid-supported bilayers prepared from a zwitterionic phospholipid of the same alkyl chain length. These mixed-charge surfactants on n-alkane-modified silica are stable in water over time (months), results that suggest the potential use of these hybrid bilayers for generating supported lipid-bilayer-like surfaces or for separation applications.

Ultrafast Dynamics at the Lipid-Water Interface: DMSO Modulates H-Bond Lifetimes

Dimethyl sulfoxide (DMSO) is a common cosolvent and cryopreservation agent used to freeze cells and tissues. DMSO alters the H-bond structure of water, but its interactions with biomolecules and, specifically, with biological interfaces remain poorly understood. Here we investigate the effects of DMSO on the H-bond dynamics at the lipid-water interface using a combination of ultrafast two-dimensional infrared (2D IR) spectroscopy and molecular dynamics simulations. Ester carbonyl absorption spectra show that DMSO dehydrates the interface, and simulations show that the area per lipid is decreased. Ultrafast 2D IR spectra measure the time scales of frequency fluctuations at the ester carbonyl positions located precisely between the hydrophobic and hydrophilic regions of the membrane. 2D IR measurements show that low DMSO concentrations (<10 mol %) induce ∼40% faster H-bond dynamics compared with pure water, whereas increased concentrations (>10-20 mol %) once again slow down the dynamics. This slow-fast-slow trend is described in terms of two different solvation regimes. Below 10 mol %, DMSO weakens the interfacial H bond, leading to faster "bulk-like" dynamics, whereas above 10 mol %, water molecules become "relatively immobilized" as the H-bond networks becoming disrupted by the H-bond donor/acceptor imbalance at the interface. These studies are an important step toward characterizing the environments around lipid membranes, which are essential to numerous biological processes.

Physicochemical properties of L-alpha dipalmitoyl phosphatidylcholine large unilamellar vesicles: Effect of hydrophobic block (PLA/PCL) of amphipathic diblock copolymers

In the present work, we show how amphipathic diblock copolymers affect the physicochemical properties of the lipid bilayer of DPPC liposome. Diblock copolymers proposed for this study are focused in the difference between PLA and PCL hydrophobic block, because PLA and PCL differ in their glass transition temperature, where a higher ratio of PLA, lowers the flexibility of the diblock copolymer. On the contrary, a greater proportion of PCL makes the diblock copolymer more flexible. This flexibility difference between hydrophobic block would affect the physicochemical properties of lipid bilayer of DPPC. The difference of rigidity or flexibility of hydrophobic block and their interaction with DPPC large unilamellar vesicles (LUVs) was evaluated at low and high copolymers concentration. The copolymer concentrations used were chosen based on their respective cmc. We measure (a) Thermotropic behavior from GP of Laurdan and fluorescence anisotropy of DPH; (b) Relation between wavelength excitation and generalized polarization of Laurdan; (c) Time-resolved fluorescence anisotropy of DPH; (d) Water outflow through the lipid bilayer and (e) calcein release from DPPC LUVs. Furthermore, large unilamellar vesicles in the absence and in the presence of different copolymers were characterized by size and zeta-potential. The results show that the diblock copolymer at high PLA/PCL ratio, that is, greater rigidity of hydrophobic block produces an increase of the phase transition temperature (Tm). For DPPC LUVs, Tm increase 3.5 °C at low and about 4.5 °C at high copolymers concentration, sensed by Laurdan and DPH fluorescent probes, although the DPPC/copolymers molar ratio for Cop4 is higher than Cop3, Cop2 and Cop1. In addition, we observed a decrease in the polarity of microenvironments in the bilayer and an increase in the order of the acyl chains in the bilayer to a high proportion of PLA. Furthermore, the presence of diblock copolymer with high proportion of PLA, decreases water outflow from DPPC liposome and water efflux is slower; leading to a decrease in calcein release from DPPC liposomes. Our results clearly show that the greater the stiffness of the hydrophobic block, greater degree of packaging of the lipid bilayer, greater the order of the acyl chains, and greater retention of water and calcein inside the liposome. Therefore, the presence of AB-type diblock copolymers with a more rigid hydrophobic block, stabilizes the lipid bilayer and would allow a more controlled release of water, and encapsulated molecules inside of the DPPC liposome.

Establishing an Electrostatics Paradigm for Membrane Electroporation in the Framework of Dissipative Particle Dynamics

With an exclusive aim to looking into a mechanism of membrane electroporation on mesoscopic length and time scales, we report the dissipative particle dynamics (DPD) simulation results for systems with and without electrolytes. A polarizable DPD model of water is employed for accurate modeling of long-range electrostatics near the water-lipid interfaces. A great deal of discussion on field induced change in dipole moments of water and lipids together with the special variation of electric field is made in order to understand the dielectrophoretic movement of water, initiating a pore formation via an intrusion through the bilayer core. The presence of salt alters the dipolar arrangement of lipids and water, and thereby it reduces the external field required to create a pore in the membrane. The species fluxes through the pore, distributions for bead density, electrostatic potential, stresses across the membrane, etc. are used to answer some of the key questions pertaining to mechanism of electroporation. The findings are compared with the molecular dynamics simulation results found in the literature, and the comparison successfully establishes an electrostatics paradigm for biomembrane studies using DPD simulations.

Amphiphilic copolymers change the nature of the ordered-to-disordered phase transition of lipid membranes from discontinuous to continuous

The phase behaviour as a function of temperature is explored for pure phospholipid (DPPC) and hybrid lipid-polymer (DPPC/Pluronic L64) bilayers with the aid of atomistic MD simulations. The range of the fixed-temperature simulations includes temperatures below and above the known melting temperature (Tm) for DPPC membranes. For the pure lipid bilayer, the main phase transition is discontinuous, as verified by the abrupt changes observed in the membrane structure, elasticity and the lipid diffusivity near the critical temperature Tm, which lies in the region 298.15-303.15 K. A pre-transition step is detected at 298.15 K which has been identified as the ripple phase (Pβ'), where ordered and disordered lipids coexist, causing thickness fluctuations. In the ordered gel phase, the positional ordering as assessed by the lipid radial distribution functions is long-range and some degree of hexagonal packing is measured. The hybrid bilayers on the other hand, transform from a more ordered to a disordered phase in a continuous manner, without finite jumps in their properties. No signs of the ripple phase are identified and the ordered phase exhibits very limited hexagonal packing and some positional ordering that decays fast. The effect of the inserted polymers in the two phases is reversed; at low temperatures, they render the membrane thinner, less cohesive and less ordered compared to the pure one, with the lipids assuming faster diffusion rates, whereas at high temperatures, the polymer interaction with the lipids acts reducing their diffusivity, but also increasing the lipid tail ordering and the membrane stiffness. The ability of the amphiphilic L64 copolymers to change the nature of the main phase transition of lipid membranes and their properties both in the ordered and the disordered phase is of vital importance for the prediction of both the efficacy of hybrid lipid/polymer nanoparticles as drug delivery vehicles as well as their potential adverse implications during interactions with healthy cell membranes.

Effect of Imidazolium-Based Ionic Liquids on the Structure and Phase Behavior of Palmitoyl-oleoyl-phosphatidylethanolamine

Among various applications, ionic liquids (ILs) have been used as antimicrobial agents in laboratories, possibly through induction of the leakage of bacteria. A molecular-level understanding of the mechanism that describes how ILs enhance the permeation of membranes is still lacking. In this study, the effects of imidazolium-based ILs with different alky chain lengths on the structure and phase behavior of 1-palmitoyl-2-oleoyl-phosphatidylethanolamine (POPE), which is a representative bacteria-membrane-rich lipid, have been investigated. By employing differential scanning calorimetry and synchrotron small- and wide-angle X-ray scattering techniques, we found that ILs with longer alkyl chains influenced the phase behavior more effectively, and lower IL concentrations are needed to induce phase separation for both lamellar liquid crystalline phase and nonlamellar inverted hexagonal phase of POPE. Interestingly, the IL with an alkyl chain of 12 carbon atoms ([C₁₂mim]Cl) shows a difference. It exhibits a stronger disturbing effect on the POPE bilayer structure than [C₁₆mim]Cl, indicating that the ability of ILs to influence the membrane structures is dependent not only on the alkyl chain length of ILs, but also on the degree of matching of the alkyl chain lengths of ILs and lipids. The new lamellar and nonlamellar structures induced by ILs both have smaller repeat distances than that of pure POPE, implying thinner membrane structures. Data of the fluorescence-based vesicle dye leakage assay are consistent with these results, particularly the defects caused by IL-induced phase separation can enhance the membrane permeability markedly. The present work may shed light on our understanding of the antimicrobial mechanism of ILs.

Coupling of Membrane Nanodomain Formation and Enhanced Electroporation near Phase Transition

Biological cells are enveloped by a heterogeneous lipid bilayer that prevents the uncontrolled exchange of substances between the cell interior and its environment. In particular, membranes act as a continuous barrier for salt and macromolecules to ensure proper physiological functions within the cell. However, it has been shown that membrane permeability strongly depends on temperature and, for phospholipid bilayers, displays a maximum at the transition between the gel and fluid phase. Here, extensive molecular dynamics simulations of dipalmitoylphosphatidylcholine bilayers were employed to characterize the membrane structure and dynamics close to phase transition, as well as its stability with respect to an external electric field. Atomistic simulations revealed the dynamic appearance and disappearance of spatially related nanometer-sized thick ordered and thin interdigitating domains in a fluid-like bilayer close to the phase transition temperature (Tm). These structures likely represent metastable precursors of the ripple phase that vanished at increased temperatures. Similarly, a two-phase bilayer with coexisting gel and fluid domains featured a thickness minimum at the interface because of splaying and interdigitating lipids. For all systems, application of an external electric field revealed a reduced bilayer stability with respect to pore formation for temperatures close to Tm. Pore formation occurred exclusively in thin interdigitating membrane nanodomains. These findings provide a link between the increased membrane permeability and the structural heterogeneity close to phase transition.

Extracting lipid vesicles from plasma membranes via self-assembly of clathrin-inspired scaffolding nanoparticles

Single-cell analysis is a new and rapidly expanding field, the goal of which is obtaining fresh information from individual cells to understand the regulatory mechanisms of cell development and diseases. Conventional approaches generally rely on the cell lysis which, however, is destructive to cells and against multiple sampling from the living cell. Here, we propose and design a scaffolding nanoparticle (NP) system that enables us to sample cytoplasmic contents without rupturing the cellular membrane, by mimicking the unusual features of clathrin. Our simulation results reveal the design principles, following which scaffolding NPs can extract lipid vesicles from plasma membranes, with both the pathway and the mechanism resembling the clathrin-mediated endocytosis, i.e. multiple NPs deposit at the membrane, assembling into cage-like structures to deform the membrane into a vesicle shape. As important design parameters, the interaction between different NPs should be properly stronger than that between each NP and the membrane to ensure the cage formation, and optimal NP concentration and the membrane surface tension are also requisite for extracting lipid vesicles. Our results provide useful guidelines for design of bio-inspired scaffolding NPs as an intelligent machine for practical use in but not limited to the single-cell analysis.

Effects of phosphatidylcholine membrane fluidity on the conformation and aggregation of N-terminally acetylated α-synuclein

Membrane association of α-synuclein (α-syn), a neuronal protein associated with Parkinson's disease (PD), is involved in α-syn function and pathology. Most previous studies on α-syn-membrane interactions have not used the physiologically relevant N-terminally acetylated (N-acetyl) α-syn form nor the most naturally abundant cellular lipid, i.e. phosphatidylcholine (PC). Here, we report on how PC membrane fluidity affects the conformation and aggregation propensity of N-acetyl α-syn. It is well established that upon membrane binding, α-syn adopts an α-helical structure. Using CD spectroscopy, we show that N-acetyl α-syn transitions from α-helical to disordered at the lipid melting temperature (T_m ). We found that this fluidity sensing is a robust characteristic, unaffected by acyl chain length (T_m = 34-55 °C) and preserved in its homologs β- and γ-syn. Interestingly, both N-acetyl α-syn membrane binding and amyloid formation trended with lipid order (1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) > 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)/sphingomyelin/cholesterol (2:2:1) ≥ DOPC), with gel-phase vesicles shortening aggregation kinetics and promoting fibril formation compared to fluid membranes. Furthermore, we found that acetylation enhances binding to PC micelles and small unilamellar vesicles with high curvature (r ∼16-20 nm) and that DPPC binding is reduced in the presence of cholesterol. These results confirmed that the exposure of hydrocarbon chains (i.e. packing defects) is essential for binding to zwitterionic gel membranes. Collectively, our in vitro results suggest that N-acetyl α-syn localizes to highly curved, ordered membranes inside a cell. We propose that age-related changes in membrane fluidity can promote the formation of amyloid fibrils, insoluble materials associated with PD.

Intramolecular structural parameters are key modulators of the gel-liquid transition in coarse grained simulations of DPPC and DOPC lipid bilayers

The capability of coarse-grained models based on the MARTINI mapping to reproduce the gel-liquid phase transition in saturated and unsaturated model lipids was investigated. We found that the model is able to reproduce a lower critical temperature for 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) with respect to 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). Nonetheless, the appearance of a gel phase for DOPC is strictly dependent on the intramolecular parameters chosen to model its molecular structure. In particular, we show that the bending angle at the coarse-grained bead corresponding to the unsaturated carbon-carbon bond acts as an order parameter determining the temperature of the phase transition. Structural analysis of the molecular dynamics simulations runs evidences that in the gel phase, the packing of the lipophilic tails of DOPC assume a different conformation than in the liquid phase. In the latter phase, the DOPC geometry resembles that of the relaxed free molecule. DPPC:DOPC mixtures show a single phase transition temperature, indicating that the observation of a phase separation between the two lipids requires the simulation of systems with sizes much larger than the ones used here.