Structural characterization on the gel to liquid-crystal phase transition of fully hydrated DSPC and DSPE bilayers

Development of a coarse-grained lipid model, LIME 2.0, for DSPE using multistate iterative Boltzmann inversion and discontinuous molecular dynamics simulations

We suggest an improved version of the intermediate resolution implicit solvent model for lipids, LIME, that was previously developed for use with discontinuous molecular dynamics (DMD) simulations. LIME gets its geometrical and the energy parameters between bonded and nonbonded pairs of coarse-grained (CG) sites from atomistic simulations. The improved model, LIME 2.0, uses multiple square wells rather than the single square well used in original LIME to obtain intermolecular interactions that more faithfully mimic those from atomistic simulations. The multi-state iterative Boltzmann inversion (MS-IBI) scheme is used to determine the interaction parameters. This means that a single set of interaction parameters between coarse-grained sites can be used to represent the lipid bilayers at different temperatures. The physical properties of CG DSPE lipid bilayer are calculated using CG simulations and compared to atomistic simulations results to verify the improved model. The phase transition temperature of the lipid bilayer is measured accurately and the lipid translocation phenomenon, " flip-flop" is observed through CG simulation. These results suggest that CG parameterization using multiple square-well and the MS-IBI scheme is well suited to the study of lipid bilayers cross a range of temperatures with DMD simulations.

Synthetic Liposomal Mimics of Biological Viruses for the Study of Immune Responses to Infection and Vaccination

Human viruses possess very complex supramolecular structures. Both icosahedral and enveloped viruses typically display an array of viral-encoded protein antigens at varied spatial densities on the viral particle surface. The viral nucleic acid genome, on the other hand, is encapsulated inside the viral particle. Although both the surface antigen and the interior nucleic acids could independently produce immunological responses, how B cells integrate these two types of signals and respond to a typical virus particle to initiate activation is not well understood at a molecular level. The study of these fundamental biological processes would benefit from the development of viral structural mimics that are well constructed to incorporate both quantitative and qualitative viral features for presentation to B cells. These novel tools would enable researchers to systematically dissect the underlying processes. Here we report the development of such particulate antigens based on liposomes engineered to display a model protein antigen, hen egg lysozyme (HEL). We developed methods to overexpress and purify various affinity mutants of HEL from E. coli. We conjugated the purified recombinant HEL proteins onto the surface of a virion-sized liposome in an orientation-specific manner at defined spatial densities and also encapsulated nucleic acid molecules into the interior of the liposome. Both the chemical conjugation of the HEL antigen on liposome surfaces and the encapsulation of nucleic acids were stable under physiologically relevant conditions. These liposomes elicited antigen-specific B-cell responses in vitro, which validate these supramolecular structures as a novel and effective approach to mimic and systematically isolate the role of essential viral features in directing the B-cell response to particulate antigens.

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.

Impact of Additives on the Microstructural Properties of DIDMAMS Bilayers: Studied by Molecular Dynamics Simulations

To further understand the mechanism of the impact of perfume raw materials (PRMs) such as allyl heptoate (AHT) and cashmeran (CMR) on distearoyl isopropyl dimethylammonium methyl sulfate (DIDMAMS) bilayers, 90 ns molecular dynamics simulations were conducted to investigate the structure of bilayers consisting of DIDMAMS and PRMs at 350 K on the molecular scale. Structural properties such as density profiles, order parameters, radial distribution functions (RDFs), and bilayer thickness were analyzed. The bilayers appear to be the structure of the ripple phase whether PRMs are added or not. The RDF and density profiles show that CMR molecules tend to locate in the region close to head groups and AHT molecules prefer to uniformly distribute among hydrocarbon chains. The special distribution of CMR molecules results in hydrocarbon chains twining around CMR molecules. The existence of CMR molecules between bilayers and the consequent highest bilayer thickness may be the main cause of higher viscosity. We expect that this work can help to screen stable vesicular formula and understand the relationship between microstructures of the vesicles and macroscopic fluidic properties.

Lipid lateral self-diffusion drop at liquid-gel phase transition

A drop of lipid lateral self-diffusion coefficient at the liquid-gel phase transition in lipid membranes is calculated. So far this drop was missing theoretical description. Our microscopic model captures so-called subdiffusion regime, which takes place on 1 ps-100 ns timescale and reveals a jump of self-diffusion coefficient. Calculation of the jump is based on our recent study of liquid-gel phase transition. Subdiffusive regime is described within the free volume theory. Calculated values of self-diffusion coefficient are in agreement with quasielastic neutron scattering measurements. Self-diffusion coefficient is found to be composed of two factors: one is related to an area per lipid change at the phase transition, and the other one is due to an order of magnitude change in the stiffness of entropic repulsive potential.

Importance of Hydrophilic Groups on Modulating the Structural, Mechanical, and Interfacial Properties of Bilayers: A Comparative Molecular Dynamics Study of Phosphatidylcholine and Ion Pair Amphiphile Membranes

An ion pair amphiphile (IPA), a molecular complex composed of two oppositely charged amphiphiles, is a phospholipid mimic which differs from a phospholipid only in the hydrophilic compositions. Here, we utilized molecular dynamics (MD) simulations to compare the bilayer systems composed of phosphatidylcholines (PC) and alkyltrimethylammonium-alkylsulfate IPAs with various alkyl chain lengths. The membrane properties for both liquid-disordered (Ld) and gel (S) phase bilayers were examined via running simulations above and below the main transition temperatures. The electrostatic attraction between the IPA hydrophilic groups leads to a more ordered molecular packing within both S and Ld phase IPA membranes, as revealed by the molecular area, deuterium order parameter, and gauche conformation analyses. Furthermore, IPA bilayers possess a higher area compressibility modulus, molecular tilt modulus, and effective bending rigidity than PC systems. The variation of hydrophilic groups of IPA also leads to fewer hydrogen bonds on the membrane surface and smaller electrostatic potentials for the biomimetic bilayer. The non-covalently linked head groups of IPA further decouple alkyl tilting and surface water retention. The combined results reveal the importance of hydrophilic groups of amphiphiles on modulating the membrane properties, which also provides insights for designs of biomimetic membranes.

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.

Membrane phase transition during heating and cooling: molecular insight into reversible melting

With increasing temperature, lipid bilayers undergo a gel-fluid phase transition, which plays an essential role in many physiological phenomena. In the present work, this first-order phase transition was investigated for variable heating and cooling rates for a dipalmitoylphosphatidylcholine (DPPC) lipid bilayer by means of atomistic molecular dynamics simulations. Alternative methods to track the melting temperature [Formula: see text] are compared. The resulting [Formula: see text] is shown to be independent of the scan rate for small heating rates (0.05-0.3 K/ns) implying reversible melting, and increases for larger heating (0.3-4 K/ns) or cooling rates (2-0.1 K/ns). The reported dependency of the melting temperature on the heating rate is in perfect agreement with a two-state kinetic rate model as suggested previously. Expansion and shrinkage, as well as the dynamics of melting seeds is described. The simulations further exhibit a relative shift between melting seeds in opposing membrane leaflets as predicted from continuum elastic theory.

Characterizing the interactions of two lipid modifications with lipid rafts: farnesyl anchors vs. palmitoyl anchors

Farnesyl (Far) and palmitoyl (Pal) anchors play important roles in the traffic of many lipidated proteins. Herein, we show the distinctive interactions and influences of the two lipid modifications on lipid rafts (LRs) and non-raft-like membranes using molecular dynamics simulations. Palmitoyl anchors behave in a more ordered fashion, pack tighter with the lipids of LRs and diffuse at a slower rate than farnesyl anchors in LRs. When interacting with non-raft-like membranes these two types of anchors become less ordered, pack more loosely with lipids, and diffuse at a higher rate. By calculating both the number of contacts per chain and the number of contact atoms per carbon of the two anchors with the lipid components, we found that the palmitoyl chains preferred to associate with the saturated chains of lipids and cholesterol molecules in LRs, while farnesyl chains favored association with saturated chains and unsaturated chains. For non-raft-like membranes, these two lipid anchors had roughly the same preference for the three types of contact lipid chains. Additionally, palmitoyl anchors caused cholesterol to orient more perpendicular to the membrane surface, surrounding lipids to become more ordered, and lipid lateral fluidity to reduce significantly, compared to farnesyl anchors in LRs. By contrast, the POPE and DSPC became much less ordered, cholesterol became more tilted, and lipids became more fluid, when the two types lipid anchors were inserted in non-raft-like membranes. These findings are useful for understanding the traffic mechanisms of lipidated proteins with farnesyl and palmitoyl modifications in cell membranes.

Multiscale Modelling of Bionano Interface

We present a framework for coarse-grained modelling of the interface between foreign nanoparticles (NP) and biological fluids and membranes. Our model includes united-atom presentations of membrane lipids and globular proteins in implicit solvent, which are based on all-atom structures of the corresponding molecules and parameterised using experimental data or atomistic simulation results. The NPs are modelled by homogeneous spheres that interact with the beads of biomolecules via a central force that depends on the NP size. The proposed methodology is used to predict the adsorption energies for human blood plasma proteins on NPs of different sizes as well as the preferred orientation of the molecules upon adsorption. Our approach allows one to rank the proteins by their binding affinity to the NP, which can be used for predicting the composition of the NP-protein corona for the corresponding material. We also show how the model can be used for studying NP interaction with a lipid bilayer membrane and thus can provide a mechanistic insight for modelling NP toxicity.

Molecular Simulation of the DPPE Lipid Bilayer Gel Phase: Coupling between Molecular Packing Order and Tail Tilt Angle

The structural properties and thermal stability of dipalmitoylphosphatidylethanolamine (DPPE) in the ordered gel phase have been studied by molecular dynamics simulation using two force fields: the Berger united-atom model and the CHARMM C36 atomistic model. As is widely known, structural features are sensitive to the initial preparation of the gel phase structure, as some degrees of freedom are slow to equilibrate on the simulation time scale of hundreds of nanoseconds. In particular, we find that the degree of alignment of the lipids' glycerol backbones, which join the two hydrocarbon tails of each molecule, strongly affects the tilt angle of the tails in the resulting structures. Disorder in the backbone correlates with lower tilt angles: bilayer configurations initiated with aligned backbones produced tilt angles near 21° and 29° for the Berger and C36 force fields, respectively, while structures initiated with randomized backbone orientations showed average tilt angles of 7° and 18°, in closer agreement with the untilted structure observed experimentally. The transition temperature for the Berger force field gel bilayer has been determined by monitoring changes in width of gel phase stripe domains as a function of temperature and is 12 ± 5 K lower than the experimental value.

Conformation of Capping Ligands on Nanoplates: Facet-Edge-Induced Disorder and Self-Assembly-Related Ordering Revealed by Sum Frequency Generation Spectroscopy

Surface-curvature-amplified conformational disorder in alkyl capping ligands has been observed previously when the nanoparticle radii approach the ligand length. Herein, sum frequency generation studies on oleic-acid-capped nanoplates show that even on faceted surfaces with dimensions tens of times greater than the ligand length a significant proportion of gauche defects exist in the capping layer. The molecular disorder on the nanosized facets is attributed to a facet-edge effect, which is diminished when increasing the facet size or assembling the nanofacets side to side. This feature is further explored to probe the self-assembly dynamics of nanoplates.

Asymmetric lipid membranes: towards more realistic model systems

Despite the ubiquity of transbilayer asymmetry in natural cell membranes, the vast majority of existing research has utilized chemically well-defined symmetric liposomes, where the inner and outer bilayer leaflets have the same composition. Here, we review various aspects of asymmetry in nature and in model systems in anticipation for the next phase of model membrane studies.

Simulation study of the permeability of a model lipid membrane at the fluid-solid phase transition

When a range of lipid bilayers are melted to the disordered fluid phase from the (much less permeable) ordered gel phase, their permeability to a variety of permeants shows a peak at the transition temperature and drops off with increasing temperature, rather than just rising as melting proceeds. To explore this anomalous behavior, a simulated coarse-grained lipid membrane model that exhibits a phase transition upon expansion or compression was studied to determine how the permeation rate of a simple particle depends on the phase composition in the two-phase region and on particle size. The permeation rate and each phase's area fraction and area density could be directly calculated, along with the probability that the permeant would cross in either phase or in interfacial regions. For large permeants and system sizes, conditions could be found where permeability increases upon compression of the bilayer. Permeation was negligible in the gel phase and, in contrast to the predictions of the "leaky interface" hypothesis, was not enriched in interfacial regions. The anomalous effect could instead be attributed to an increase in the area per lipid of fluid-phase domains. This result motivated a model for the decrease in effective permeability barrier through fluid-phase domains arising from a decrease in the length of the gel/fluid interface at the midpoint of a permeation event.

Modeling lipid raft domains containing a mono-unsaturated phosphatidylethanolamine species

An advanced coarse-grained model for “atypical” lipid rafts was built and validated to be employed in studies of membrane-protein interactions.

Long-timescale motions in glycerol-monopalmitate lipid bilayers investigated using molecular dynamics simulation

The occurrence of long-timescale motions in glycerol-1-monopalmitate (GMP) lipid bilayers is investigated based on previously reported 600 ns molecular dynamics simulations of a 2×8×8 GMP bilayer patch in the temperature range 302-338 K, performed at three different hydration levels, or in the presence of the cosolutes methanol or trehalose at three different concentrations. The types of long-timescale motions considered are: (i) the possible phase transitions; (ii) the precession of the relative collective tilt-angle of the two leaflets in the gel phase; (iii) the trans-gauche isomerization of the dihedral angles within the lipid aliphatic tails; and (iv) the flipping of single lipids across the two leaflets. The results provide a picture of GMP bilayers involving a rich spectrum of events occurring on a wide range of timescales, from the 100-ps range isomerization of single dihedral angles, via the 100-ns range of tilt precession motions, to the multi-μs range of phase transitions and lipid-flipping events.

Physicochemical properties of liposomes as potential anticancer drugs carriers. Interaction of etoposide and cytarabine with the membrane: spectroscopic studies

The interactions between etoposide, cytarabine and 1,2-dihexadecanoyl-sn-glycerol-3-phosphocholine bilayers were studied using differential scanning calorimetry (DSC), Fourier transform infrared spectroscopy (FT-IR) and nuclear magnetic resonance (NMR). These techniques have proven to be a very powerful tool in studying the structure and dynamics of phospholipid bilayers. In particular, DSC can provide information on the phase transition temperature and cooperativity of the lipid molecules in the absence and presence of the drug. Vibrational spectroscopy is well suited to the study of drug-lipid interactions, since it allows for an investigation of the conformation of phospholipid molecules at different levels in lipid bilayers and follows structural changes that occur during the gel to liquid-crystalline phase transition. NMR supported the determination of the main phase transition temperatures (TC) of 1,2-dihexadecanoyl-sn-glycerol-3-phosphocholine (DPPC). The main phase transition temperature (TC) determined by (1)H NMR is comparable with values obtained by DSC for all studied liposomes. The location of cytarabine and etoposide in liposomes was also determined by NMR. Atomic force microscopy (AFM) images, acquired immediately after sample deposition on a mica surface, revealed the spherical shape of lipid vesicles.

Surface tension effects on the phase transition of a DPPC bilayer with and without protein: a molecular dynamics simulation

A coarse-grained molecular dynamics simulation was applied to illustrate the phase transition behavior of the pure DPPC bilayer and aquaporin-embedded DPPC bilayer under different surface tensions.

Molecular dynamics study of catanionic bilayers composed of ion pair amphiphile with double-tailed cationic surfactant

The physical stability of catanionic vesicles is important for the development of novel drug or DNA carriers. For investigating the mechanism by which catanionic vesicles are stabilized, molecular dynamics (MD) simulation is an attractive approach that provides microscopic structural information on the vesicular bilayer. In this study, MD simulation was applied to investigate the bilayer properties of catanionic vesicles composed of an ion pair amphiphile (IPA), hexadecyltrimethylammonium-dodecylsulfate (HTMA-DS), and a double-tailed cationic surfactant, ditetradecyldimethylammonium chloride (DTDAC). Structural information regarding membrane elasticity and the organization and conformation of surfactant molecules was obtained based on the resulting trajectory. Simulation results showed that a proper amount of DTDAC could be used to complement the asymmetric structure between HTMA and DS, resulting in an ordered hydrocarbon chain packing within the rigid membrane observed in the mixed HTMA-DS/DTDAC system. The coexistence of gel and fluid phases was also observed in the presence of excess DTDAC. MD simulation results agreed well with results obtained from experimental studies examining mixed HTMA-DS/DTDAB vesicles.

Using molecular dynamics to study liquid phase behavior: simulations of the ternary sodium laurate/sodium oleate/water system

The prediction of surfactant phase behavior has applications in a wide range of areas. An accurate modeling of liquid phase behavior can aid our understanding of colloidal process or be used to design phases that respond in a defined way to their environment. In this work, we use molecular dynamics to model the phase behavior of the ternary sodium laurate/sodium oleate/water system and compare the simulation results to experimental data. Simulations were performed with the GROMOS 53A6 united-atom force field and cover the entire ternary phase diagram, producing micellar, hexagonal, and lamellar phases. The aggregate simulation time for the 33 simulations performed during this study is 4.4 μs. We find that the simulations were able to model the experimentally observed liquid phase behavior accurately, showing that the carboxylate and lipid parameters of the 53A6 force field give very good quality results for the in silico prediction of liquid system phase behavior.

The role of electrostatic interactions in the membrane binding of melittin

The binding of melittin and the C-terminally truncated analogue of melittin (21Q) to a range of phospholipid bilayers was studied using surface plasmon resonance (SPR). The phospholipid model membranes included zwitterionic dimyristylphosphatidylcholine (DMPC) and dimyristylphosphatidylethanolamine (DMPE), together with mixtures DMPC/dimyristylphosphatidylglycerol (DMPG), DMPC/DMPG/cholesterol and DMPE/DMPG. Melittin bound rapidly to all membrane mixtures, whereas 21Q, which has a reduced charge, bound much more slowly on the DMPC and DMPC/DMPG mixtures reflecting the role of the initial electrostatic interaction. The loss of the cationic residues also significantly decreased the binding of 21Q with DMPC/DMPG/Cholesterol, DMPE and DMPE/DMPG. The role of electrostatics was also highlighted with NaCl in the buffer, which affected the way melittin bound to the different membranes, causing a more uniform, concentration dependant increase in response. The biosensor results were correlated with the conformation of the peptides determined by circular dichroism analysis, which indicated that high α-helicity was associated with high binding affinity. Overall, the results demonstrate that the positively charged residues at the C-terminus of melittin play an essential role in membrane binding, that modulation of peptide charge influences selectivity of binding to different phospholipids and that manipulation of the cationic regions of antimicrobial peptides can be used to modulate membrane selectivity.