Lipid tail chain asymmetry and the strength of membrane-induced interactions between membrane proteins

Many lipids are composed of asymmetric tail chains that differ by their molecular weight (MW) and/or degree of saturation. Previous studies found that membrane moduli vary with the degree of lipid tail asymmetry. However, to date little is known regarding the effect (if any) of tail asymmetry on the membrane-induced interactions between embedded proteins. In this paper we use a self-consistent field model to examine the effect of lipid tail asymmetry on membrane proteins. We first examine the case where the overall tail length (sum of both chains) is held constant, which implies that the membrane thickness remains constant as well, independent of tail asymmetry. We find that, in these systems, the membrane area stretch and bending moduli decrease with increasing chain asymmetry, thereby reducing the magnitude of the membrane-induced barrier to protein aggregation. Since in symmetric lipid bilayers the energy barrier is typically of order ∼ 1-2 times the thermal energy kT, the asymmetry-induced reduction in barrier height may increase the probability of protein aggregation significantly. In systems where one tail chain is held constant, increasing asymmetry involves changes in the bilayer thickness which are found to dominate any effect arising from the asymmetry.     

Molecular dynamics simulation of the evolution of hydrophobic defects in one monolayer of a phosphatidylcholine bilayer: relevance for membrane fusion mechanisms

The spontaneous formation of the phospholipid bilayer underlies the permeability barrier function of the biological membrane. Tears or defects that expose water to the acyl chains are spontaneously healed by lipid lateral diffusion. However, mechanical barriers, e.g., protein aggregates held in place, could sustain hydrophobic defects. Such defects have been postulated to occur in processes such as membrane fusion. This gives rise to a new question in bilayer structure: What do the lipids do in the absence of lipid lateral diffusion to minimize the free energy of a hydrophobic defect? As a first step to understand this rather fundamental question about bilayer structure, we performed molecular dynamic simulations of up to 10 ns of a planar bilayer from which lipids have been deleted randomly from one monolayer. In one set of simulations, approximately one-half of the lipids in the defect monolayer were restrained to form a mechanical barrier. In the second set, lipids were free to diffuse around. The question was simply whether the defects caused by removing a lipid would aggregate together, forming a large hydrophobic cavity, or whether the membrane would adjust in another way. When there are no mechanical barriers, the lipids in the defect monolayer simply spread out and thin with little effect on the other intact monolayer. In the presence of a mechanical barrier, the behavior of the lipids depends on the size of the defect. When 3 of 64 lipids are removed, the remaining lipids adjust the lower one-half of their chains, but the headgroup structure changes little and the intact monolayer is unaffected. When 6 to 12 lipids are removed, the defect monolayer thins, lipid disorder increases, and lipids from the intact monolayer move toward the defect monolayer. Whereas this is a highly simplified model of a fusion site, this engagement of the intact monolayer into the fusion defect is strikingly consistent with recent results for influenza hemagglutinin mediated fusion.     

Lipid tail chain asymmetry and the strength of membrane-induced interactions between membrane proteins

Many lipids are composed of asymmetric tail chains that differ by their molecular weight (MW) and/or degree of saturation. Previous studies found that membrane moduli vary with the degree of lipid tail asymmetry. However, to date little is known regarding the effect (if any) of tail asymmetry on the membrane-induced interactions between embedded proteins. In this paper we use a self-consistent field model to examine the effect of lipid tail asymmetry on membrane proteins. We first examine the case where the overall tail length (sum of both chains) is held constant, which implies that the membrane thickness remains constant as well, independent of tail asymmetry. We find that, in these systems, the membrane area stretch and bending moduli decrease with increasing chain asymmetry, thereby reducing the magnitude of the membrane-induced barrier to protein aggregation. Since in symmetric lipid bilayers the energy barrier is typically of order approximately 1-2 times the thermal energy kT, the asymmetry-induced reduction in barrier height may increase the probability of protein aggregation significantly. In systems where one tail chain is held constant, increasing asymmetry involves changes in the bilayer thickness which are found to dominate any effect arising from the asymmetry.     

Regulation of channel function due to physical energetic coupling with a lipid bilayer

Regulation of membrane protein functions due to hydrophobic coupling with a lipid bilayer has been investigated. An energy formula describing interactions between lipid bilayer and integral ion channels with different structures, which is based on the screened Coulomb interaction approximation, has been developed. Here the interaction energy is represented as being due to charge-based interactions between channel and lipid bilayer. The hydrophobic bilayer thickness channel length mismatch is found to induce channel destabilization exponentially while negative lipid curvature linearly. Experimental parameters related to channel dynamics are consistent with theoretical predictions. To measure comparable energy parameters directly in the system and to elucidate the mechanism at an atomistic level we performed molecular dynamics (MD) simulations of the ion channel forming peptide–lipid complexes. MD simulations indicate that peptides and lipids experience electrostatic and van der Waals interactions for short period of time when found within each other’s proximity. The energies from these two interactions are found to be similar to the energies derived theoretically using the screened Coulomb and the van der Waals interactions between peptides (in ion channel) and lipids (in lipid bilayer) due to mainly their charge properties. The results of in silico MD studies taken together with experimental observable parameters and theoretical energetic predictions suggest that the peptides induce ion channels inside lipid membranes due to peptide–lipid physical interactions. This study provides a new insight helping better understand of the underlying mechanisms of membrane protein functions in cell membrane leading to important biological implications.     

Artificial biomembrane morphology: a dissipative particle dynamics study

Artificial membranes mimicking biological structures are rapidly breaking new ground in the areas of medicine and soft-matter physics. In this endeavor, we use dissipative particle dynamics simulation to investigate the morphology and behavior of lipid-based biomembranes under conditions of varied lipid density and self-interaction. Our results show that a less-than-normal initial lipid density does not create the traditional membrane; but instead results in the formation of a ‘net’, or at very low densities, a series of disparate ‘clumps’ similar to the micelles formed by lipids in nature. When the initial lipid density is high, a membrane forms, but due to the large number of lipids, the naturally formed membrane would be larger than the simulation box, leading to ‘rippling’ behavior as the excess repulsive force of the membrane interior overcomes the bending energy of the membrane. Once the density reaches a certain point however, ‘bubbles’ appear inside the membrane, reducing the rippling behavior and eventually generating a relatively flat, but thick, structure with micelles of water inside the membrane itself. Our simulations also demonstrate that the interaction parameter between individual lipids plays a significant role in the formation and behavior of lipid membrane assemblies, creating similar structures as the initial lipid density distribution. This work provides a comprehensive approach to the intricacies of lipid membranes, and offers a guideline to design biological or polymeric membranes through self-assembly processes as well as develop novel cellular manipulation and destruction techniques.     

Molecular dynamics study of bipolar tetraether lipid membranes

Membranes composed of bipolar tetraether lipids have been studied by a series of 25-ns molecular dynamics simulations to understand the microscopic structure and dynamics as well as membrane area elasticity. By comparing macrocyclic and acyclic tetraether and diether archaeal lipids, the effect of tail linkage of the two phytanyl-chained lipids on the membrane properties is elucidated. Tetraether lipids show smaller molecular area and lateral mobility. For the latter, calculated diffusion coefficients are indeed one order-of-magnitude smaller than that of the diether lipid. These two tetraether membranes are alike in many physical properties except for membrane area elasticity. The macrocyclic tetraether membrane shows a higher elastic area expansion modulus than its acyclic counterpart by a factor of three. Free energy profiles of a water molecule crossing the membranes show no major difference in barrier height; however, a significant difference is observed near the membrane center due to the lack of the slip-plane in tetraether membranes.     

Dynamics of fusion pores connecting membranes of different tensions

The energetics underlying the expansion of fusion pores connecting biological or lipid bilayer membranes is elucidated. The energetics necessary to deform membranes as the pore enlarges, in some combination with the action of the fusion proteins, must determine pore growth. The dynamics of pore growth is considered for the case of two homogeneous fusing membranes under different tensions. It is rigorously shown that pore growth can be quantitatively described by treating the pore as a quasiparticle that moves in a medium with a viscosity determined by that of the membranes. Motion is subject to tension, bending, and viscous forces. Pore dynamics and lipid flow through the pore were calculated using Lagrange's equations, with dissipation caused by intra- and intermonolayer friction. These calculations show that the energy barrier that restrains pore enlargement depends only on the sum of the tensions; a difference in tension between the fusing membranes is irrelevant. In contrast, lipid flux through the fusion pore depends on the tension difference but is independent of the sum. Thus pore growth is not affected by tension-driven lipid flux from one membrane to the other. The calculations of the present study explain how increases in tension through osmotic swelling of vesicles cause enlargement of pores between the vesicles and planar bilayer membranes. In a similar fashion, swelling of secretory granules after fusion in biological systems could promote pore enlargement during exocytosis. The calculations also show that pore expansion can be caused by pore lengthening; lengthening may be facilitated by fusion proteins.     

Effect of lipid peroxidation on the properties of lipid bilayers: a molecular dynamics study

Lipid peroxidation plays an important role in cell membrane damage. We investigated the effect of lipid peroxidation on the properties of 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) lipid bilayers using molecular dynamics simulations. We focused on four main oxidation products of linoleic acid with either a hydroperoxide or an aldehyde group: 9-trans, cis-hydroperoxide linoleic acid, 13-trans, cis-hydroperoxide linoleic acid, 9-oxo-nonanoic acid, and 12-oxo-9-dodecenoic acid. These oxidized chains replaced the sn-2 linoleate chain. The properties of PLPC lipid bilayers were characterized as a function of the concentration of oxidized lipids, with concentrations from 2.8% to 50% for each oxidation product. The introduction of oxidized functional groups in the lipid tail leads to an important conformational change in the lipids: the oxidized tails bend toward the water phase and the oxygen atoms form hydrogen bonds with water and the polar lipid headgroup. This conformational change leads to an increase in the average area per lipid and, correspondingly, to a decrease of the bilayer thickness and the deuterium order parameters for the lipid tails, especially evident at high concentrations of oxidized lipid. Water defects are observed in the bilayers more frequently as the concentration of the oxidized lipids is increased. The changes in the structural properties of the bilayer and the water permeability are associated with the tendency of the oxidized lipid tails to bend toward the water interface. Our results suggest that one mechanism of cell membrane damage is the increase in membrane permeability due to the presence of oxidized lipids.     

Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.     

Lipid Tail Protrusion in Simulations Predicts Fusogenic Activity of Influenza Fusion Peptide Mutants and Conformational Models

Fusion peptides from influenza hemagglutinin act on membranes to promote membrane fusion, but the mechanism by which they do so remains unknown. Recent theoretical work has suggested that contact of protruding lipid tails may be an important feature of the transition state for membrane fusion. If this is so, then influenza fusion peptides would be expected to promote tail protrusion in proportion to the ability of the corresponding full-length hemagglutinin to drive lipid mixing in fusion assays. We have performed molecular dynamics simulations of influenza fusion peptides in lipid bilayers, comparing the X-31 influenza strain against a series of N-terminal mutants. As hypothesized, the probability of lipid tail protrusion correlates well with the lipid mixing rate induced by each mutant. This supports the conclusion that tail protrusion is important to the transition state for fusion. Furthermore, it suggests that tail protrusion can be used to examine how fusion peptides might interact with membranes to promote fusion. Previous models for native influenza fusion peptide structure in membranes include a kinked helix, a straight helix, and a helical hairpin. Our simulations visit each of these conformations. Thus, the free energy differences between each are likely low enough that specifics of the membrane environment and peptide construct may be sufficient to modulate the equilibrium between them. However, the kinked helix promotes lipid tail protrusion in our simulations much more strongly than the other two structures. We therefore predict that the kinked helix is the most fusogenic of these three conformations.     

How cholesterol is distributed between monolayers in asymmetric lipid membranes

The distribution of cholesterol in asymmetric lipid bilayers was studied by extensive coarse-grained molecular dynamics simulations. The effects of the lipid head group charge, acyl chain saturation, spontaneous membrane curvature and surface tension of the membrane were investigated. Four asymmetric bilayers containing DOPC, DOPS, DSPC or DSPS lipids were simulated on a time scale extended to tens of microseconds. We show that cholesterol strongly prefers anionic lipids to neutral and saturated lipid tails to unsaturated with a distribution ratio of ~0.7 in neutral/anionic bilayers and of ~0.4 in unsaturated/saturated bilayers. Multiple flip-flop transitions of cholesterol were observed directly, and their mean times ranged from 80 to 250?ns. It was shown that the distribution of cholesterol in the asymmetric membrane depends not only on the type of lipid, but also on the local membrane curvature and the surface tension. The membrane curvature enhances the influence of the lipid head groups on cholesterol distribution, while non-optimal surface tension caused by different areas per lipid in different monolayers increases the effect of the lipid tail saturation. It was clearly seen that the monolayers of asymmetric bilayers are interdependent. Mean distances from the bilayer center to cholesterol molecules depend not only on the type of the lipid in the considered monolayer but also on the composition of the opposite monolayer.     

Effect of Ceramide Tail Length on the Structure of Model Stratum Corneum Lipid Bilayers

Lipid bilayers composed of non-hydroxy sphingosine ceramide (CER NS), cholesterol (CHOL), and free fatty acids (FFAs), which are components of the human skin barrier, are studied via molecular dynamics simulations. Since mixtures of these lipids exist in dense gel phases with little molecular mobility at physiological conditions, care must be taken to ensure that the simulations become decorrelated from the initial conditions. Thus, we propose and validate an equilibration protocol based on simulated tempering, in which the simulation takes a random walk through temperature space, allowing the system to break out of metastable configurations and hence become decorrelated from its initial configuration. After validating the equilibration protocol, which we refer to as random-walk molecular dynamics, the effects of the lipid composition and ceramide tail length on bilayer properties are studied. Systems containing pure CER NS, CER NS + CHOL, and CER NS + CHOL + FFA, with the CER NS fatty acid tail length varied within each CER NS-CHOL-FFA composition, are simulated. The bilayer thickness is found to depend on the structure of the center of the bilayer, which arises as a result of the tail-length asymmetry between the lipids studied. The hydrogen bonding between the lipid headgroups and with water is found to change with the overall lipid composition, but is mostly independent of the CER fatty acid tail length. Subtle differences in the lateral packing of the lipid tails are also found as a function of CER tail length. Overall, these results provide insight into the experimentally observed trend of altered barrier properties in skin systems where there are more CERs with shorter tails present.     

Stability and membrane orientation of the fukutin transmembrane domain: a combined multiscale molecular dynamics and circular dichroism study

The N-terminal domain of fukutin-I has been implicated in the localization of the protein in the endoplasmic reticulum and Golgi Apparatus. It has been proposed to mediate this through its interaction with the thinner lipid bilayers found in these compartments. Here we have employed multiscale molecular dynamics simulations and circular dichroism spectroscopy to explore the structure, stability, and orientation of the short 36-residue N-terminus of fukutin-I (FK1TMD) in lipids with differing tail lengths. Our results show that FK1TMD adopts a stable helical conformation in phosphatidylcholine lipids when oriented with its principal axis perpendicular to the bilayer plane. The stability of the helix is largely insensitive to the lipid tail length, preventing hydrophobic mismatch by virtue of its mobility and ability to tilt within the lipid bilayers. This suggests that changes in FK1TMD tilt in response to bilayer properties may be implicated in the regulation of its trafficking. Coarse-grained simulations of the complex Golgi membrane suggest the N-terminal domain may induce the formation of microdomains in the surrounding membrane through its preferential interaction with 1,2-dipalmitoyl-sn-glycero-3-phosphatidylinositol 4,5-bisphosphate lipids.     

Atomic-scale structure and electrostatics of anionic palmitoyloleoylphosphatidylglycerol lipid bilayers with Na+ counterions

Anionic palmitoyloleoylphosphatidylglycerol (POPG) is one of the most abundant lipids in nature, yet its atomic-scale properties have not received significant attention. Here we report extensive 150-ns molecular dynamics simulations of a pure POPG lipid membrane with sodium counterions. It turns out that the average area per lipid of the POPG bilayer under physiological conditions is approximately 19% smaller than that of a bilayer built from its zwitterionic phosphatidylcholine analog, palmitoyloleoylphosphatidylcholine. This suggests that there are strong attractive interactions between anionic POPG lipids, which overcome the electrostatic repulsion between negative charges of PG headgroups. We demonstrate that interlipid counterion bridges and strong intra- and intermolecular hydrogen bonding play a key role in this seemingly counterintuitive behavior. In particular, the substantial strength and stability of ion-mediated binding between anionic lipid headgroups leads to complexation of PG molecules and ions and formation of large PG-ion clusters that act in a concerted manner. The ion-mediated binding seems to provide a possible molecular-level explanation for the low permeability of PG-containing bacterial membranes to organic solvents: highly polar interactions at the water/membrane interface are able to create a high free energy barrier for hydrophobic molecules such as benzene.     

Small scale membrane mechanics

Large scale changes to lipid bilayer shapes are well represented by the Helfrich model. However, there are membrane processes that take place at smaller length scales that this model cannot address. In this work, we present a one-dimensional continuum model that captures the mechanics of the lipid bilayer membrane at the length scale of the lipids themselves. The model is developed using the Cosserat theory of surfaces with lipid orientation, or ‘tilt’, as the fundamental degree of freedom. The Helfrich model can be recovered as a special case when the curvatures are small and the lipid tilt is everywhere zero. We use the tilt model to study local membrane deformations in response to a protein inclusion. Parameter estimates and boundary conditions are obtained from a coarse-grained molecular model using dissipative particle dynamics (DPD) to capture the same phenomenon. The continuum model is able to reproduce the membrane bending, stretch and lipid tilt as seen in the DPD model. The lipid tilt angle relaxes to the bulk tilt angle within 5–6 nm from the protein inclusion. Importantly, for large tilt gradients induced by the proteins, the tilt energy contribution is larger than the bending energy contribution. Thus, the continuum model of tilt accurately captures behaviors at length scales shorter than the membrane thickness.     

Molecular dynamics simulation of lipid reorientation at bilayer edges

Understanding cellular membrane processes is critical for the study of events such as viral entry, neurotransmitter exocytosis, and immune activation. Supported lipid bilayers are commonly used to model these membrane processes experimentally. Despite the relative simplicity of such a system, many important structural and dynamic parameters are not experimentally observable with current techniques. Computational approaches allow the development of a high-resolution model of bilayer processes. We have performed molecular dynamics simulations of dimyristoylphosphatidylcholine (DMPC) bilayers to model the creation of bilayer gaps—a common process in bilayer patterning—and to analyze their structure and dynamics. We propose a model for gap formation in which the bilayer edges form metastable micelle-like structures on a nanosecond timescale. Molecules near edges structurally resemble lipids in ungapped bilayers but undergo small-scale motions more rapidly. These data suggest that lipids may undergo rapid local rearrangements during membrane fusion, facilitating the formation of fusion intermediates thought key to the infection cycle of viruses such as influenza, Ebola, and HIV.     

Lycopene,but not zeaxanthin,serves as a skeleton for the formation of an orthorhombic organization of intercellular lipids within the lamellae in the stratum corneum: Molecular dynamics simulations of the hydrated ceramide NS bilayer model

Carotenoids play an important role in the protection of biomembranes against oxidative damage. Their function depends on the surroundings and the organization of the lipid membrane they are embedded in. Carotenoids are located parallel or perpendicular to the surface of the lipid bilayer. The influence of carotenoids on the organization of the lipid bilayer in the stratum corneum has not been thoroughly considered. Here, the orientation of the exemplary cutaneous carotenoids lycopene and zeaxanthin in a hydrated ceramide NS24 bilayer model and the influence of carotenoids on the lateral organization of the lipid bilayer model were studied by means of molecular dynamics simulations for 32 °C and 37 °C. The results confirm that lycopene is located parallel and zeaxanthin perpendicular to the surface of the lipid bilayer. The lycopene-loaded lipid bilayer appeared to have a strong orthorhombic organization, while zeaxanthin-loaded and pure lipid bilayers were organized in a disordered hexagonal-like and liquid-like state, respectively. The effect is stronger at 32 °C compared to 37 °C based on p-values. Therefore, it was assumed that carotenoids without hydroxyl polar groups in their structure facilitate the formation of the orthorhombic organization of lipids, which provides the skin barrier function. It was shown that the distance between carotenoid atoms matched the distance between atoms in the lipids, indicating that parallel located carotenoids without hydroxyl groups serve as a skeleton for lipid membranes inside the lamellae. The obtained results provide reasonable prediction of the overall qualitative properties of lipid model systems and show the importance of parallel-oriented carotenoids in the development and maintenance of the skin barrier function.     

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 (T_m). 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 T_m. 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.     

Bilayer conformation of fusion peptide of influenza virus hemagglutinin: a molecular dynamics simulation study

Unraveling the conformation of membrane-bound viral fusion peptides is essential for understanding how those peptides destabilize the bilayer topology of lipids that is important for virus-cell membrane fusion. Here, molecular dynamics (MD) simulations were performed to investigate the conformation of the 20 amino acids long fusion peptide of influenza hemagglutinin of strain X31 bound to a dimyristoyl phosphatidylcholine (DMPC) bilayer. The simulations revealed that the peptide adopts a kinked conformation, in agreement with the NMR structures of a related peptide in detergent micelles. The peptide is located at the amphipathic interface between the headgroups and hydrocarbon chains of the lipid by an energetically favorable arrangement: The hydrophobic side chains of the peptides are embedded into the hydrophobic region and the hydrophilic side chains are in the headgroup region. The N-terminus of the peptide is localized close to the amphipathic interface. The molecular dynamics simulations also revealed that the peptide affects the surrounding bilayer structure. The average hydrophobic thickness of the lipid phase close to the N-terminus is reduced in comparison with the average hydrophobic thickness of a pure dimyristoyl phosphatidylcholine bilayer.