Mixed Bilayers Research Articles

Mixed lipid bilayers with locally varying spontaneous curvature and bending.

A model of lipid bilayers made of a mixture of two lipids with different average compositions on both leaflets, is developed. A Landau Hamiltonian describing the lipid-lipid interactions on each leaflet, with two lipidic fields ψ 1 and ψ 2, is coupled to a Helfrich one, accounting for the membrane elasticity, via both a local spontaneous curvature, which varies as C 0 + C 1(ψ 1 - ψ 2/2), and a bending modulus equal to κ 0 + κ 1(ψ 1 + ψ 2)/2. This model allows us to define curved patches as membrane domains where the asymmetry in composition, ψ 1 - ψ 2, is large, and thick and stiff patches where ψ 1 + ψ 2 is large. These thick patches are good candidates for being lipidic rafts, as observed in cell membranes, which are composed primarily of saturated lipids forming a liquid-ordered domain and are known to be thick and flat nano-domains. The lipid-lipid structure factors and correlation functions are computed for globally spherical membranes and planar ones and for a whole set of parameters including the surface tension and the coupling in the two leaflet compositions. Phase diagrams are established, within a Gaussian approximation, showing the occurrence of two types of Structure Disordered phases, with correlations between either curved or thick patches, and an Ordered phase, corresponding to the divergence of the structure factor at a finite wave vector. The varying bending modulus plays a central role for curved membranes, where the driving force κ 1 C 0 (2) is balanced by the line tension, to form raft domains of size ranging from 10 to 100 nm. For planar membranes, raft domains emerge via the cross-correlation with curved domains. A global picture emerges from curvature-induced mechanisms, described in the literature for planar membranes, to coupled curvature- and bending-induced mechanisms in curved membranes forming a closed vesicle.

Molecular dynamics study of changes in physico-chemical properties of DMPC lipid bilayers by addition of nonionic surfactant C12E10

Changes in physico-chemical properties of dimyristoyl phosphatidylcholine (DMPC) lipid bilayers caused by the addition of 9.4 mol% nonionic surfactant decaoxyethylene monododecyl ethers (C12E10) have been investigated by molecular dynamics calculations. In spite of addition of single chain C12E10, the lipid bilayers showed an increase of membrane area. Isothermal area compressibility, which is a measure of membrane softness in lateral direction, also increased by 50% for DMPC/C12E10 mixed bilayers. Furthermore, the order parameter of C–H vector for DMPC acyl tails decreased. We found that these changes are caused by the hydrophilic head groups of C12E10 which are located near the glycerol backbone of the DMPC molecules and have bulky random coil conformation without any preferential ordered structures.

Vancomycin-eluting niosomes: a new approach to the inhibition of staphylococcal biofilm on abiotic surfaces.

A new vancomycin (VCM)-eluting mixed bilayer niosome formulation was evaluated for the control of staphylococcal colonization and biofilm formation on abiotic surfaces, a niosome application not explored to date. Cosurfactant niosomes were prepared using a Span 60/Tween 40/cholesterol blend (1: 1: 2). Tween 40, a polyethoxylated amphiphile, was included to enhance VCM entrapment and confer niosomal surface properties precluding bacterial adhesion. VCM-eluting niosomes showed good quality attributes including relatively high entrapment efficiency (∼50%), association of Tween 40 with vesicles in a constant proportion (∼87%), biphasic release profile suitable for inhibiting early bacterial colonization, and long-term stability at 4°C for a 12-month study period. Niosomes significantly enhanced VCM activity against planktonic bacteria of nine staphylococcal strains. Using microtiter plates as abiotic surface, VCM-eluting niosomes proved superior to VCM in inhibiting biofilm formation, eradicating surface-borne biofilms, inhibiting biofilm growth, and interfering with biofilm induction by VCM subminimal inhibitory concentrations. Data suggest dual functionality of cosurfactant VCM-eluting niosomes as passive colonization inhibiting barrier and active antimicrobial-controlled delivery system, two functions recognized in infection control of abiotic surfaces and medical devices.

Self-assembly of phosphate amphiphiles in mixtures of prebiotically plausible surfactants.

The spontaneous formation of closed bilayer structures from prebiotically plausible amphiphiles is an essential requirement for the emergence of early cells on prebiotic Earth. The sources of amphiphiles could have been both endo- and exogenous (accretion of meteorite carbonaceous material or interstellar dust particles). Among all prebiotic possible amphiphile candidates, those containing phosphate are the least investigated species because their self-assembly occurs in a seemingly too narrow range of conditions. The self-assembly of simple phosphate amphiphiles should, however, be of great interest, as contemporary membranes predominantly contain phospholipids. In contrast to common expectations, we show that these amphiphiles can be easily synthesized under prebiotically plausible environmental conditions and can efficiently form bilayer structures in the presence of various co-surfactants across a large range of pH values. Vesiculation was even observed in crude reaction mixtures that contained 1-decanol as the amphiphile precursor. The two best co-surfactants promoted vesicle formation over the entire pH range in aqueous solutions. Expanding the pH range where bilayer membranes self-assemble and remain intact is a prerequisite for the emergence of early cell-like compartments and their preservation under fluctuating environmental conditions. These mixed bilayers also retained small charged solutes, such as dyes. These results demonstrate that alkyl phosphate amphiphiles might have played a significant role as early compartment building blocks.

Preferred Orientations of Phosphoinositides in Bilayers and Their Implications in Protein Recognition Mechanisms

Phosphoinositides (PIPs), phosphorylated derivatives of phosphatidylinositol (PI), are essential regulatory lipids involved in various cellular processes, including signal transduction, membrane trafficking, and cytoskeletal remodeling. To gain insight into the protein-PIPs recognition process, it is necessary to study the inositol ring orientation (with respect to the membrane) of PIPs with different phosphorylation states. In this study, 8 PIPs (3 PIP, 2 PIP2, and 3 PIP3) with different phosphorylation and protonation sites have been separately simulated in two mixed bilayers (one with 20% phosphatidylserine (PS) lipids and another with PS lipids switched to phosphatidylcholine (PC) lipids), which roughly correspond to yeast membranes. Uniformity of the bilayer properties including hydrophobic thickness, acyl chain order parameters, and heavy atom density profiles is observed in both PS-contained and PC-enriched membranes due to the same hydrophobic core composition. The relationship between the inositol ring orientation (tilt and rotation angles) and its solvent-accessible surface area indicates that the orientation is mainly determined by its solvation energy. Different PIPs exhibit a clear preference in the inositol ring rotation angle. Surprisingly, a larger proportion of PIPs inositol rings stay closer to the surface of PS-contained membranes compared to PC-enriched ones. Such a difference is rationalized with the formation of more hydrogen bonds between the PS/PI headgroups and the PIPs inositol rings in PS-contained membranes. This hydrogen bond network could be functionally important; thus, the present results can potentially add important and detailed features into the existing protein-PIPs recognition mechanism.

Correlation of band gap position with composition in high-k films

Two series of ultrathin high-k samples (mixed layers HfxAlyOz and bilayers HfO2 on Al2O3) prepared by atomic layer deposition were investigated using spectroscopic ellipsometry in the energy range of 0.7–10.0 eV. The (effective) optical gap of both mixed layer and bilayer structures can be tuned by the film composition. The optical gap of mixed layers is linearly dependent on the Hf fraction from 5.77 (±0.02) eV for pure HfO2 to 6.71 (±0.02) eV for pure Al2O3. The effective absorption gap of bilayers measured in reflection geometry is lower than that of mixed layers with a comparable Hf fraction due to the dominant effect of the top layer. An increase in film thickness as well as a decrease in refractive index and gap energy was observed after sample storage for two months in atmosphere. The aging effect is likely due to further oxidation of the oxygen deficient high-k films caused by the oxygen diffusion from air into the films.

Synthetic Multivalent Antifungal Peptides Effective against Fungi

Taking advantage of the cluster effect observed in multivalent peptides, this work describes antifungal activity and possible mechanism of action of tetravalent peptide (B4010) which carries 4 copies of the sequence RGRKVVRR through a branched lysine core. B4010 displayed better antifungal properties than natamycin and amphotericin B. The peptide retained significant activity in the presence of monovalent/divalent cations, trypsin and serum and tear fluid. Moreover, B4010 is non-haemolytic and non-toxic to mice by intraperitoneal (200 mg/kg) or intravenous (100 mg/kg) routes. S. cerevisiae mutant strains with altered membrane sterol structures and composition showed hyper senstivity to B4010. The peptide had no affinity for cell wall polysaccharides and caused rapid dissipation of membrane potential and release of vital ions and ATP when treated with C. albicans. We demonstrate that additives which alter the membrane potential or membrane rigidity protect C. albicans from B4010-induced lethality. Calcein release assay and molecular dynamics simulations showed that the peptide preferentially binds to mixed bilayer containing ergosterol over phophotidylcholine-cholesterol bilayers. The studies further suggested that the first arginine is important for mediating peptide-bilayer interactions. Replacing the first arginine led to a 2–4 fold decrease in antifungal activities and reduced membrane disruption properties. The combined in silico and in vitro approach should facilitate rational design of new tetravalent antifungal peptides.

Cellular Scale Biophysical Models of Membrane Sculpting by the Proteins During Endocytosis and Exocytosis

We employ coarse-grained molecular dynamics (CGMD) simulations and continuum mechanics methods to quantify the membrane-bending action of exo70. CGMD simulations show that exo70 produces the strongest curvature when bound to a key lipid PIP2, and when dimerized. To relate these simulations to experiments, we first quantify the deformation field generated by exo70 using CGMD simulation, and then investigate the resulting cellular morphologies using a continuum mechanics model for membrane bending. In the first step, we characterize the deformation field by directly analyzing curved membrane conformations generated by CGMD simulations of mutant and wild-type exo70 on mixed lipid bilayers. Monte Carlo simulations yield membrane morphologies which we relate to cryo-transmission electron microscopy images of exo70 remodeling a GUV. This model agrees with experiments which show that PIP2 binding and dimerization are necessary for generating membrane tubulues and protrusions. We show that this multi-scale model allows us to relate the structural as well as thermodynamic details of exo70 membrane-binding and self-association to the induced deformation field and resulting cellular morphologies. We also extend our studies to ENTH, BAR, and ESCRT family proteins. We show that thermal undulations in the membrane and cooperativity in the curvature fields, collectively drive the membrane into different morphological states (buds, tubules, etc.) that resemble those in cellular experiments in vivo and vesicle experiments in vitro. We determine the relative stability of the above mentioned shapes based on the free energy of these membrane configurations, determined using the methods of Thermodynamic Integration (TI) and Bennett Acceptance (BA), and Widom Insertion techniques. Results are shown for the case of Exo70 protein and ENTH, and N-BAR domains and also compared against measurements determined from experiments.

Surface-Tension Replica-Exchange Molecular Dynamics Method for Efficient Conformational Sampling of Biological Membrane Systems

Conformational sampling is fundamentally important for simulating complex bio-molecular systems. Generalized-ensemble algorithm, especially the temperature replica-exchange molecular dynamics method (REMD), is one of the widely used methods to explore structures of bio-molecules. Most temperature REMD simulations have focused on soluble proteins rather than membrane proteins or lipid bilayers, because explicit membranes do not keep their structural integrity at high temperature. Here, we propose a new generalized-ensemble algorithm for membrane systems, which we call the surface-tension REMD method. Each replica is simulated in the NPγT ensemble, where surface tensions in a pair of replicas are exchanged at certain intervals to enhance conformational sampling of the target membrane system. We tested our method on a fully hydrated DPPC lipid bilayer. During the simulation, a random walk in surface tension space is realized. Large-scale lateral deformation of DPPC membranes takes place in all of the replicas without collapse of the lipid bilayer structure. Our method could be applicable to a wide variety of biological membrane systems including mixed lipid bilayers and membrane-protein systems.

Computational Design of Multi-component Bio-Inspired Bilayer Membranes

Our investigation is motivated by the need to design bilayer membranes with tunable interfacial and mechanical properties for use in a range of applications, such as targeted drug delivery, sensing and imaging. We draw inspiration from biological cell membranes and focus on their principal constituents. In this paper, we present our results on the role of molecular architecture on the interfacial, structural and dynamical properties of bio-inspired membranes. We focus on four lipid architectures with variations in the head group shape and the hydrocarbon tail length. Each lipid species is composed of a hydrophilic head group and two hydrophobic tails. In addition, we study a model of the Cholesterol molecule to understand the interfacial properties of a bilayer membrane composed of rigid, single-tail molecular species. We demonstrate the properties of the bilayer membranes to be determined by the molecular architecture and rigidity of the constituent species. Finally, we demonstrate the formation of a stable mixed bilayer membrane composed of Cholesterol and one of the phospholipid species. Our approach can be adopted to design multi-component bilayer membranes with tunable interfacial and mechanical properties. We use a Molecular Dynamics-based mesoscopic simulation technique called Dissipative Particle Dynamics that resolves the molecular details of the components through soft-sphere coarse-grained models and reproduces the hydrodynamic behavior of the system over extended time scales.

Computational Design of Multi-component Bio-Inspired Bilayer Membranes

Our investigation is motivated by the need to design bilayer membranes with tunable interfacial and mechanical properties for use in a range of applications, such as targeted drug delivery, sensing and imaging. We draw inspiration from biological cell membranes and focus on their principal constituents. In this paper, we present our results on the role of molecular architecture on the interfacial, structural and dynamical properties of bio-inspired membranes. We focus on four lipid architectures with variations in the head group shape and the hydrocarbon tail length. Each lipid species is composed of a hydrophilic head group and two hydrophobic tails. In addition, we study a model of the Cholesterol molecule to understand the interfacial properties of a bilayer membrane composed of rigid, single-tail molecular species. We demonstrate the properties of the bilayer membranes to be determined by the molecular architecture and rigidity of the constituent species. Finally, we demonstrate the formation of a stable mixed bilayer membrane composed of Cholesterol and one of the phospholipid species. Our approach can be adopted to design multi-component bilayer membranes with tunable interfacial and mechanical properties. We use a Molecular Dynamics-based mesoscopic simulation technique called Dissipative Particle Dynamics that resolves the molecular details of the components through soft-sphere coarse-grained models and reproduces the hydrodynamic behavior of the system over extended time scales.

The energy transfer based fluorescent approach to detect the formation of silica supported phosphatidylcholine and phosphatidylserine containing bilayers

The work introduces the quenching of silica coated Tb(III) complexes by merocyanine 540 (MC540) and copper ions as a tool to reveal the adsorption of phosphatidylcholine (PC) and phosphatidylserine (PS) at various PS:PC ratio onto silica nanoparticles doped with Tb(III) complex. The binding of MC540 with PC-based bilayers and copper ions with PS-based ones are the basis of their use as organic and inorganic probes to sense PS:PC ratio in silica supported mixed bilayers. The enrichment of mixed bilayers by PS results in the displacement of MC540 anions, while it enhances the binding with copper ions. The displacement or binding of probe ions results in the diverse luminescence response of Tb(III)-centred luminescence. The reestablishment of the steady state and time resolved luminescence is observed, when MC540 anions are applied as probes. The use of copper ions as probes results in the opposite quenching effect. The developed route enables to discriminate the formation of the phospholipids bilayers onto silica surface from those in the bulk of solution under various concentration conditions.

Lipid Flip-Flop in Binary Membranes Composed of Phosphatidylserine and Phosphatidylcholine

The kinetics and thermodynamics of lipid flip-flop in bilayers composed of 1,2-dipalmitoyl-sn-glycero-3-phospho-L-serine (DPPS) and 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) were studied using sum-frequency vibrational spectroscopy. The kinetics of DSPC and DPPS flip-flop were examined as a function of temperature and bilayer composition. The rate of DSPC flip-flop did not exhibit any significant dependence on bilayer composition while the rate of DPPS flip-flop was inversely dependent on the mole fraction of DPPS. The transition-state thermodynamics for DSPC and DPPS lipids in these mixed bilayers were determined in order to identify the energetic impact of the phosphatidylserine headgroup on lipid flip-flop. The thermodynamics for the DSPC component remained statistically identical to bilayers composed entirely of DSPC. The activation energy for the DPPS component showed a linear correlation with the mole fraction of DPPS for all bilayer compositions. The enthalpy and entropy for DPPS flip-flop did not increase linearly with the fraction of DPPS but did directly correlate with the molecular area. The DPPS component also exhibited enthalpy-entropy compensation which suggests that lipid hydration may play a significant role in membrane dynamics.

MemBuilder: a web-based graphical interface to build heterogeneously mixed membrane bilayers for the GROMACS biomolecular simulation program

Molecular dynamics (MD) simulations have had a profound impact on studies of membrane proteins during past two decades, but the accuracy of MD simulations of membranes is limited by the quality of membrane models and the applied force fields. Membrane models used in MD simulations mostly contain one kind of lipid molecule. This is far from reality, for biological membranes always contain more than one kind of lipid molecule. Moreover, the lipid composition and their distribution are functionally important. As a result, there is a necessity to prepare more realistic lipid membranes containing different types of lipids at physiological concentrations. To automate and simplify the building process of heterogeneous lipid bilayers as well as providing molecular topologies for included lipids based on both united and all-atom force fields, we provided MemBuilder as a web-based graphical user interface. MemBuilder is a free web server available from www.membuilder.org.

Molecular Dynamics and NMR Spectroscopy Studies of E. coli Lipopolysaccharide Structure and Dynamics

Lipopolysaccharide (LPS), a component of Gram-negative bacterial outer membranes, comprises three regions: lipid A, core oligosaccharide, and O-antigen polysaccharide. Using the CHARMM36 lipid and carbohydrate force fields, we have constructed a model of an Escherichia coli R1 (core) O6 (antigen) LPS molecule. Several all-atom bilayers are built and simulated with lipid A only (LIPA) and varying lengths of 0 (LPS0), 5 (LPS5), and 10 (LPS10) O6 antigen repeating units; a single unit of O6 antigen contains five sugar residues. From 1H,1H-NOESY experiments, cross-relaxation rates are obtained from an O-antigen polysaccharide sample. Although some experimental deviations are due to spin-diffusion, the remaining effective proton-proton distances show generally very good agreement between NMR experiments and molecular dynamics simulations. The simulation results show that increasing the LPS molecular length has an impact on LPS structure and dynamics and also on LPS bilayer properties. Terminal residues in a LPS bilayer are more flexible and extended along the membrane normal. As the core and O-antigen are added, per-lipid area increases and lipid bilayer order decreases. In addition, results from mixed LPS0/5 and LPS0/10 bilayer simulations show that the LPS O-antigen conformations at a higher concentration of LPS5 and LPS10 are more orthogonal to the membrane and less flexible. The O-antigen concentration of mixed LPS bilayers does not have a significant effect on per-lipid area and hydrophobic thickness. Analysis of ion and water penetration shows that water molecules can penetrate inside the inner core region, and hydration is critical to maintain the integrity of the bilayer structure.

Lipid-II forms potential “landing terrain” for lantibiotics in simulated bacterial membrane

Bacterial cell wall is targeted by many antibiotics. Among them are lantibiotics, which realize their function via interaction with plasma membrane lipid-II molecule — a chemically conserved part of the cell wall synthesis pathway. To investigate structural and dynamic properties of this molecule, we have performed a series of nearly microsecond-long molecular dynamics simulations of lipid-II and some of its analogs in zwitterionic single component and charged mixed simulated phospholipid bilayers (the reference and the mimic of the bacterial plasma membrane, respectively). Extensive analysis revealed that lipid-II forms a unique “amphiphilic pattern” exclusively on the surface of the simulated bacterial membrane (and not in the reference one). We hypothesize that many lantibiotics exploit the conserved features of lipid-II along with characteristic modulation of the bacterial membrane as the “landing site”. This putative recognition mechanism opens new opportunities for studies on lantibiotics action and design of novel armament against resistant bacterial strains.

Prediction of Micelle/Water and Liposome/Water Partition Coefficients Based on Molecular Dynamics Simulations, COSMO-RS, and COSMOmic

Liposomes and micelles find various applications as potential solubilizers in extraction processes or in drug delivery systems. Thermodynamic and transport processes governing the interactions of different kinds of solutes in liposomes or micelles can be analyzed regarding the free energy profiles of the solutes in the system. However, free energy profiles in heterogeneous systems such as micelles are experimentally almost not accessible. Therefore, the development of predictive methods is desirable. Molecular dynamics (MD) simulations reliably simulate the structure and dynamics of lipid membranes and micelles, whereas COSMO-RS accurately reproduces solvation free energies in different solvents. For the first time, free energy profiles in micellar systems, as well as mixed lipid bilayers, are investigated, taking advantage of both methods: MD simulations and COSMO-RS, referred to as COSMOmic (Klamt, A.; Huniar, U.; Spycher, S.; Keldenich, J. COSMOmic: A Mechanistic Approach to the Calculation of Membrane-Water Partition Coefficients and Internal Distributions within Membranes and Micelles. J. Phys. Chem. B 2008, 112, 12148-12157). All-atom molecular dynamics simulations of the system SDS/water and CTAB/water have been applied in order to retrieve representative micelle structures for further analysis with COSMOmic. For the system CTAB/water, different surfactant concentrations were considered, which results in different micelle sizes. Free energy profiles of more than 200 solutes were predicted and validated by means of experimental partition coefficients. To our knowledge, these are the first quantitative predictions of micelle/water partition coefficients, which are based on whole free energy profiles from molecular methods. Further, the partitioning in lipid bilayer systems containing different hydrophobic tail groups (DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine), SOPC (stearoyl-oleoylphosphatidylcholine), DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine), and POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine)) as well as mixed bilayers was calculated. Experimental partition coefficients (log P) were reproduced with a root-mean-square error (RMSE) of 0.62. To determine the influence of cholesterol as an important component of cellular membranes, free energy profiles in the presence of cholesterol were calculated and shown to be in good agreement with experimental data.

Interaction of Human Synovial Phospholipase A2 with Mixed Lipid Bilayers: A Coarse-Grain and All-Atom Molecular Dynamics Simulation Study

Human secreted phospholipase A2s have been shown to promote inflammation in mammals by catalyzing the first step of the arachidonic acid pathway by breaking down phospholipids, producing fatty acids, including arachidonic acid. They bind to the membrane water interface to access their phospholipid substrates from the membrane. Their binding modes on membrane surfaces are regulated by diverse factors, including membrane charge, fluidity, and heterogeneity. The influence of these factors on the binding modes of the enzymes is not well understood. Here we have studied several human synovial phospholipase A2 (hs-PLA2)/mixed bilayer systems through a combined coarse-grain and all-atom molecular dynamics simulation. It was found that hydrophobic residues Leu2, Val3, Ala18, Leu19, Phe23, Gly30, and Phe63 that form the edge of the entrance of the hydrophobic binding pocket in hs-PLA2 tend to penetrate into the hydrophobic area of lipid bilayers, and more than half of the total amino acid residues make contact with the lipid headgroups. Each enzyme molecule forms 19-38 hydrogen bonds with the bilayer to which it binds, most of which are with the phosphate groups. Analysis of the root-mean-square deviation (rmsd) shows that residues Val30-Thr40, Tyr66-Gln80, and Lys107-Arg118 have relatively large rmsds during all-atom molecular dynamics simulations, in accordance with the observation of an enlarged entrance region of the hydrophobic binding pocket. The amino acid sequences forming the entrance of the binding pocket prefer to interact with lipid molecules that are more fluid or negatively charged, and the opening of the binding pocket would be larger when the lipid components are more fluid.

Dendritic Amphiphiles Strongly Affect the Biophysical Properties of DPPC Bilayer Membranes

Molecular dynamics (MD) simulations were used to gain insight on the molecular interactions in a model biological membrane comprised of a bilayer with DPPC (dipalmitoylphosphotidylcholine) and antimicrobial dendritic amphiphile molecules [RCONHC(CH(2)CH(2)COOH)(3), where R is the saturated hydrocarbon tail (R = n-C(n)H(2n+1)), to be abbreviated as 3CAmn]. This study analyzes different biophysical properties of the equilibrated mixed bilayers, at 300 and 325 K, to determine how the presence of the 3CAmn, in varying concentrations and tail lengths, affects the lipid bilayer. Lipid tail order parameter data, bilayer thickness trends, and qualitative lipid tail tilt observations suggest that a molar ratio of 0.2 3CAm19/DPPC is sufficient to induce a phase transition in the bilayer from gel to liquid crystalline at 300 K. These results also imply that the phase transition temperature of the mixed bilayer decreases upon incorporation of higher concentrations of 3CAm19. Hydrogen bonding takes place between the 3CAmn and DPPC at specific sites, as evidenced by the radial distribution function. Increased hydrogen bonding and the smaller headgroup size of the 3CAmn molecule result in a decrease in the total lateral area with higher concentrations of 3CAm19. Diffusion constants of 3CAmn varied with concentration and tail length; diffusion constants of DPPC and 3CAm19 increased with increasing 3CAm19 concentration at 300 K and shorter 3CAmn tails had higher diffusion constants at both temperatures. These computational studies provide a comprehensive understanding of the biophysical changes to model biological membranes by the association of 3CAmn.

How Cholesterol Tilt Modulates the Mechanical Properties of Saturated and Unsaturated Lipid Membranes

Although there have been great advances in understanding the effect of cholesterol on various properties of lipid membranes, its mechanistic role in determining the elasticity of bilayers at the molecular level is not fully resolved. Indeed, to date the molecular mechanisms that drive the experimentally detected differences in properties of saturated and unsaturated lipid bilayers that contain cholesterol remain unclear. By quantifying the cholesterol orientational degrees of freedom from atomistic molecular dynamics simulations of mixed lipid-cholesterol membranes, we address this question from the perspective of cholesterol tilt and splay. Following the fluctuations in orientations of cholesterol and of lipid molecules in simulations, we have extracted tilt and splay moduli both for cholesterol molecules and hydrocarbon lipid tails. This has further allowed us to estimate the contributions of these modes to the response of membranes to elastic deformations. We find that tilt and splay deformations importantly contribute to the overall elasticity of the mixed lipid membranes, and that they can account for the experimentally established differences between the liquid ordered sphingomyelin (SM)/cholesterol bilayers and fluid dioleoylphosphocholine (DOPC)/cholesterol bilayers. These findings underscore the importance of tilt and splay moduli, derived from simulations with relative computational ease, as useful metrics for quantitatively characterizing the mechanical properties of mixed lipid-cholesterol bilayers.