Intermonolayer friction and surface shear viscosity of lipid bilayer membranes

The flow behavior of lipid bilayer membranes is characterized by a surface viscosity for in-plane shear deformations, and an intermonolayer friction coefficient for slip between the two leaflets of the bilayer. Both properties have been studied for a variety of coarse-grained double-tailed model lipids, using equilibrium and nonequilibrium molecular dynamics simulations. For lipids with two identical tails, the surface shear viscosity rises rapidly with tail length, while the intermonolayer friction coefficient is less sensitive to the tail length. Interdigitation of lipid tails across the bilayer midsurface, as observed for lipids with two distinct tails, strongly enhances the intermonolayer friction coefficient, but hardly affects the surface shear viscosity. The simulation results are compared against the available experimental data.     

Coarse-grained molecular dynamics simulations of phase transitions in mixed lipid systems containing LPA, DOPA, and DOPE lipids

The mechanisms that mediate biomembrane shape transformations are of considerable interest in cell biology. Recent in vitro experiments show that the chemical transformation of minor membrane lipids can induce dramatic shape changes in biomembranes. Specifically, it was observed that the addition of DOPA to DOPE has no effect on the stability of the bilayer structure of the membrane. In contrast, the addition of LPA to DOPE stabilizes the bilayer phase of DOPE, increasing the temperature of a phase transition from the bilayer to the inverted hexagonal phase. This result suggests that the chemical conversion of DOPA to LPA is sufficient for triggering a dramatic change in the shape of biomembranes. The LPA/DOPA/DOPE mixture of lipids provides a simple model system for understanding the molecular events driving the shape change. In this work, we used coarse-grained molecular dynamics simulations to study the phase transitions of this lipid mixture. We show that despite the simplicity of the coarse-grained model, it reproduces the experimentally observed phase changes of: 1), pure LPA and DOPA with respect to changes in the concentration of cations; and 2), LPA/DOPE and DOPA/DOPE mixtures with respect to temperature. The good agreement between the model and experiments suggests that the computationally inexpensive coarse-grained approach can be used to infer macroscopic membrane properties. Furthermore, analysis of the shape of the lipid molecules demonstrates that the phase behavior of single-lipid systems is consistent with molecular packing theory. However, the phase stability of mixed lipid systems exhibits significant deviations from this theory, which suggests that the elastic energy of the lipids, neglected in the packing theory, plays an important role.     

Lipid Clustering Correlates with Membrane Curvature as Revealed by Molecular Simulations of Complex Lipid Bilayers

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP₂), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP₂, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.     

Peptide Aggregation and Pore Formation in a Lipid Bilayer: A Combined Coarse-Grained and All Atom Molecular Dynamics Study

We present a simulation study where different resolutions, namely coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, are used sequentially to combine the long timescale reachable by CG simulations with the high resolution of AA simulations, to describe the complete processes of peptide aggregation and pore formation by alamethicin peptides in a hydrated lipid bilayer. In the 1-μs CG simulations the peptides spontaneously aggregate in the lipid bilayer and exhibit occasional transitions between the membrane-spanning and the surface-bound configurations. One of the CG systems at t = 1 μs is reverted to an AA representation and subjected to AA simulation for 50 ns, during which water molecules penetrate the lipid bilayer through interactions with the peptide aggregates, and the membrane starts leaking water. During the AA simulation significant deviations from the α-helical structure of the peptides are observed, however, the size and arrangement of the clusters are not affected within the studied time frame. Solid-state NMR experiments designed to match closely the setup used in the molecular dynamics simulations provide strong support for our finding that alamethicin peptides adopt a diverse set of configurations in a lipid bilayer, which is in sharp contrast to the prevailing view of alamethicin oligomers formed by perfectly aligned helical alamethicin peptides in a lipid bilayer.     

Combined monte carlo and molecular dynamics simulation of fully hydrated dioleyl and palmitoyl-oleyl phosphatidylcholine lipid bilayers

We have applied a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer to hydrated bilayers of dioleyl-phosphatidylcholine (DOPC) and palmitoyl-oleyl phosphatidylcholine (POPC). The procedure consists of alternating molecular dynamics trajectory calculations in a constant surface tension and temperature ensemble with configurational bias Monte Carlo moves to different regions of the configuration space of the bilayer in a constant volume and temperature ensemble. The procedure is applied to bilayers of 128 molecules of POPC with 4628 water molecules, and 128 molecules of DOPC with 4825 water molecules. Progress toward equilibration is almost three times as fast in central processing unit (CPU) time compared with a purely molecular dynamics (MD) simulation. Equilibration is complete, as judged by the lack of energy drift in 200-ps runs of continuous MD. After the equilibrium state was reached, as determined by agreement between the simulation volume per lipid molecule with experiment, continuous MD was run in an ensemble in which the lateral area was restrained to fluctuate about a mean value and a pressure of 1 atm applied normal to the bilayer surface. Three separate continuous MD runs, 200 ps in duration each, separated by 10,000 CBMC steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of POPC and DOPC. Analysis of the hydration environment in the headgroups supports a mechanism by which unsaturation contributes to reduced transition temperatures. In this view, the relatively horizontal orientation of the unsaturated bond increases the area per lipid, resulting in increased water penetration between the headgroups. As a result the headgroup-headgroup interactions are attenuated and shielded, and this contributes to the lowered transition temperature.     

Transmembrane helices can induce domain formation in crowded model membranes

We studied compositionally heterogeneous multi-component model membranes comprised of saturated lipids, unsaturated lipids, cholesterol, and α-helical TM protein models using coarse-grained molecular dynamics simulations. Reducing the mismatch between the length of the saturated and unsaturated lipid tails reduced the driving force for segregation into liquid-ordered (l(o)) and liquid-disordered (l(d)) lipid domains. Cholesterol depletion had a similar effect, and binary lipid mixtures without cholesterol did not undergo large-scale phase separation under the simulation conditions. The phase-separating ternary dipalmitoyl-phosphatidylcholine (DPPC)/dilinoleoyl-PC (DLiPC)/cholesterol bilayer was found to segregate into l(o) and l(d) domains also in the presence of a high concentration of ΤΜ helices. The l(d) domain was highly crowded with TM helices (protein-to-lipid ratio ~1:5), slowing down lateral diffusion by a factor of 5-10 as compared to the dilute case, with anomalous (sub)-diffusion on the μs time scale. The membrane with the less strongly unsaturated palmitoyl-linoleoyl-PC instead of DLiPC, which in the absence of TM α-helices less strongly deviated from ideal mixing, could be brought closer to a miscibility critical point by introducing a high concentration of TM helices. Finally, the 7-TM protein bacteriorhodopsin was found to partition into the l(d) domains irrespective of hydrophobic matching. These results show that it is possible to directly study the lateral reorganization of lipids and proteins in compositionally heterogeneous and crowded model biomembranes with coarse-grained molecular dynamics simulations, a step toward simulations of realistic, compositionally complex cellular membranes. This article is part of a Special Issue entitled: Protein Folding in Membranes.     

1-Alkanols and membranes: A story of attraction

Although 1-alkanols have long been known to act as penetration enhancers and anesthetics, the mode of operation is not yet understood. In this study, long-time molecular dynamics simulations have been performed to investigate the effect of 1-alkanols of various carbon chain lengths onto the structure and dynamics of dimyristoylphosphatidylcholine bilayers. The simulations were complemented by microcalorimetry, continuous bleaching and film balance experiments. In the simulations, all investigated 1-alkanols assembled inside the lipid bilayer within tens of nanoseconds. Their hydroxyl groups bound preferentially to the lipid carbonyl group and the hydrocarbon chains stretched into the hydrophobic core of the bilayer. Both molecular dynamics simulations and experiments showed that all 1-alkanols drastically affected the bilayer properties. Insertion of long-chain 1-alkanols decreased the area per lipid while increasing the thickness of the bilayer and the order of the lipids. The bilayer elasticity was reduced and the diffusive motion of the lipids within the bilayer plane was suppressed. On the other hand, integration of ethanol into the bilayer enlarged the area per lipid. The bilayer became softer and lipid diffusion was enhanced.     

Exploring peptide-membrane interactions with coarse-grained MD simulations

The interaction of α-helical peptides with lipid bilayers is central to our understanding of the physicochemical principles of biological membrane organization and stability. Mutations that alter the position or orientation of an α-helix within a membrane, or that change the probability that the α-helix will insert into the membrane, can alter a range of membrane protein functions. We describe a comparative coarse-grained molecular dynamics simulation methodology, based on self-assembly of a lipid bilayer in the presence of an α-helical peptide, which allows us to model membrane transmembrane helix insertion. We validate this methodology against available experimental data for synthetic model peptides (WALP23 and LS3). Simulation-based estimates of apparent free energies of insertion into a bilayer of cystic fibrosis transmembrane regulator-derived helices correlate well with published data for translocon-mediated insertion. Comparison of values of the apparent free energy of insertion from self-assembly simulations with those from coarse-grained molecular dynamics potentials of mean force for model peptides, and with translocon-mediated insertion of cystic fibrosis transmembrane regulator-derived peptides suggests a nonequilibrium model of helix insertion into bilayers.     

1-Alkanols and membranes: a story of attraction

Although 1-alkanols have long been known to act as penetration enhancers and anesthetics, the mode of operation is not yet understood. In this study, long-time molecular dynamics simulations have been performed to investigate the effect of 1-alkanols of various carbon chain lengths onto the structure and dynamics of dimyristoylphosphatidylcholine bilayers. The simulations were complemented by microcalorimetry, continuous bleaching and film balance experiments. In the simulations, all investigated 1-alkanols assembled inside the lipid bilayer within tens of nanoseconds. Their hydroxyl groups bound preferentially to the lipid carbonyl group and the hydrocarbon chains stretched into the hydrophobic core of the bilayer. Both molecular dynamics simulations and experiments showed that all 1-alkanols drastically affected the bilayer properties. Insertion of long-chain 1-alkanols decreased the area per lipid while increasing the thickness of the bilayer and the order of the lipids. The bilayer elasticity was reduced and the diffusive motion of the lipids within the bilayer plane was suppressed. On the other hand, integration of ethanol into the bilayer enlarged the area per lipid. The bilayer became softer and lipid diffusion was enhanced.     

Effect of sodium chloride on a lipid bilayer

Electrostatic interactions govern structural and dynamical properties of membranes and can vary considerably with the composition of the aqueous buffer. We studied the influence of sodium chloride on a pure POPC lipid bilayer by fluorescence correlation spectroscopy experiments and molecular dynamics simulations. Increasing sodium chloride concentration was found to decrease the self-diffusion of POPC lipids within the bilayer. Self-diffusion coefficients calculated from the 100 ns simulations agree with those measured on a millisecond timescale, suggesting that most of the relaxation processes relevant for lipid diffusion are faster than the simulation timescale. As the dominant effect, the molecular dynamics simulations revealed a tight binding of sodium ions to the carbonyl oxygens of on average three lipids leading to larger complexes with reduced mobility. Additionally, the bilayer thickens by approximately 2 A, which increases the order parameter of the fatty acyl chains. Sodium binding alters the electrostatic potential, which is largely compensated by a changed polarization of the aqueous medium and a lipid dipole reorientation.     

Molecular dynamics simulation of lipid packing and mobility in bilayer membranes

A model of lipid bilayer membrane in water has been developed. Parameters have been selected that allow molecular dynamics simulation of lipid bilayers in the all-atom approximation. The calculated indices of packing and mobility of lipid molecules for the liquid crystalline state of the bilayer agree well with the experimental data. Based on the model of the liquid crystalline state of the membrane, a system in the gel-like state has been constructed. The gel-state model reproduces well the packing of lipids in real bilayers, whereas the mobility of molecules proves to be overestimated.     

Predicting the structure of apolipoprotein A-I in reconstituted high-density lipoprotein disks.

In reconstituted high-density lipoproteins, apolipoprotein A-I and phosphatidylcholines combine to form disks in which the amphipathic alpha-helices of apolipoprotein A-1 bind to the edge of a lipid bilayer core, shielding the hydrophic lipid tails from the aqueous environment. We have employed experimental data, sequence analysis, and molecular modeling to construct an atomic model of such a reconstituted high-density lipoprotein disk consisting of two apolipoprotein A-I proteins and 160 palmitoyloleoylphosphatidylcholine lipids. The initial globular domain (1-47) of apolipoprotein A-I was excluded from the model, which was hydrated with an 8-A shell of water molecules. Molecular dynamics and simulated annealing were used to test the stability of the model. Both head-to-tail and head-to-head forms of a reconstituted high-density lipoprotein were simulated. In our simulations the protein contained and adhered to the lipid bilayer while providing good coverage of the lipid tails.     

Combined Monte Carlo and molecular dynamics simulation of hydrated lipid-cholesterol lipid bilayers at low cholesterol concentration

We have applied a hybrid equilibration and sampling procedure for the atomic level simulation of a hydrated lipid bilayer to systems consisting of dipalmitoyl phosphatidylcholine (DPPC) and cholesterol, and palmitoyl-oleyl phosphatidylcholine (POPC) at low (approximately 6%) cholesterol concentration. The procedure is applied to bilayers of 94 molecules of DPPC, 6 molecules of cholesterol, and 3205 water molecules, and to bilayers of 120 molecules of POPC, 8 molecules of cholesterol, and 4268 water molecules, at a temperature of 325 K. After equilibration, three separate 400-ps continuous molecular dynamics runs, separated by 10,000 configurational bias Monte Carlo steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of DPPC and cholesterol, and of pure DPPC, and pure POPC. Certain properties of the bilayers are indistinguishable from cholesterol-free bilayers, including lateral diffusion and electron density. Other properties, most notably the order parameter profile, show the effect of cholesterol even at low concentrations.     

A Steroid in a Lipid Bilayer: Localization, Orientation, and Energetics

Steroid hormones are known to freely partition into lipid bilayers. As a case study, we investigated the behavior of the steroid hormone cortisone in a model lipid bilayer. First, we looked at energy barriers involved in the partitioning of a single molecule into a bilayer using umbrella sampling molecular dynamics simulations. A rather wide well of −4.5 kcal/mol was observed in the interfacial region between the lipid headgroup and tailgroup. Next, using two unconstrained molecular dynamics simulations with cortisone initially positioned at distinct locations within a bilayer, we studied the preferred location and orientation of the molecule. Finally, we observed how cortisone molecules could spontaneously insert and localize in a bilayer from bulk solution. The three independent approaches produced a converged picture of how cortisone behaves in a model lipid bilayer.     

Hydrogen-bonding propensities of sphingomyelin in solution and in a bilayer assembly: a molecular dynamics study

Sphingomyelin is enriched within lipid microdomains of the cell membrane termed lipid rafts. These microdomains play a part in regulating a variety of cellular events. Computer simulations of the hydrogen-bonding properties of sphingolipids, believed to be central to the organization of these domains, can delineate the possible molecular interactions that underlie this lipid structure. We have therefore used molecular dynamics simulations to unravel the hydrogen-bonding behavior of palmitoylsphingomyelin (PSM). A series of eight simulations of 3 ns each of a single PSM molecule in water showed that the sphingosine OH and NH groups can form hydrogen bonds with the phosphate oxygens of their own polar head, in agreement with NMR data. Simulations of PSM in a bilayer assembly were carried out for 8 ns with three different force field parameterizations. The major physico-chemical parameters of the simulated bilayer agree with those established experimentally. The sphingosine OH group was mainly involved in intramolecular hydrogen bonds, in contrast to the almost exclusive intermolecular hydrogen bonds formed by the amide NH moiety. During the bilayer simulations the intermolecular hydrogen bonds among lipids formed a dynamic network characterized by the presence of hydrogen-bonded lipid clusters of up to nine PSM molecules.     

The interaction of phospholipase A2 with a phospholipid bilayer: coarse-grained molecular dynamics simulations

A number of membrane-active enzymes act in a complex environment formed by the interface between a lipid bilayer and bulk water. Although x-ray diffraction studies yield structures of isolated enzyme molecules, a detailed characterization of their interactions with the interface requires a measure of how deeply such a membrane-associated protein penetrates into a lipid bilayer. Here, we apply coarse-grained (CG) molecular dynamics (MD) simulations to probe the interaction of porcine pancreatic phospholipase A2 (PLA2) with a lipid bilayer containing palmitoyl-oleoyl-phosphatidyl choline and palmitoyl-oleoyl-phosphatidyl glycerol molecules. We also used a configuration from a CG-MD trajectory to initiate two atomistic (AT) MD simulations. The results of the CG and AT simulations are evaluated by comparison with available experimental data. The membrane-binding surface of PLA2 consists of a patch of hydrophobic residues surrounded by polar and basic residues. We show this proposed footprint interacts preferentially with the anionic headgroups of the palmitoyl-oleoyl-phosphatidyl glycerol molecules. Thus, both electrostatic and hydrophobic interactions determine the location of PLA2 relative to the bilayer. From a general perspective, this study demonstrates that CG-MD simulations may be used to reveal the orientation and location of a membrane-surface-bound protein relative to a lipid bilayer, which may subsequently be refined by AT-MD simulations to probe more detailed interactions.     

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.     

Experimental validation of molecular dynamics simulations of lipid bilayers: a new approach

A novel protocol has been developed for comparing the structural properties of lipid bilayers determined by simulation with those determined by diffraction experiments, which makes it possible to test critically the ability of molecular dynamics simulations to reproduce experimental data. This model-independent method consists of analyzing data from molecular dynamics bilayer simulations in the same way as experimental data by determining the structure factors of the system and, via Fourier reconstruction, the overall transbilayer scattering-density profiles. Multi-nanosecond molecular dynamics simulations of a dioleoylphosphatidylcholine bilayer at 66% RH (5.4 waters/lipid) were performed in the constant pressure and temperature ensemble using the united-atom GROMACS and the all-atom CHARMM22/27 force fields with the GROMACS and NAMD software packages, respectively. The quality of the simulated bilayer structures was evaluated by comparing simulation with experimental results for bilayer thickness, area/lipid, individual molecular-component distributions, continuous and discrete structure factors, and overall scattering-density profiles. Neither the GROMACS nor the CHARMM22/27 simulations reproduced experimental data within experimental error. The widths of the simulated terminal methyl distributions showed a particularly strong disagreement with the experimentally observed distributions. A comparison of the older CHARMM22 with the newer CHARMM27 force fields shows that significant progress is being made in the development of atomic force fields for describing lipid bilayer systems empirically.     

Bilayer edge and curvature effects on partitioning of lipids by tail length: atomistic simulations

The partitioning of lipids among different microenvironments in a bilayer is of considerable relevance to characterization of composition variations in biomembranes. Atomistic simulation has been ill-suited to model equilibrated lipid mixtures because the time required for diffusive exchange of lipids among microenvironments exceeds typical submicrosecond molecular dynamics trajectories. A method to facilitate local composition fluctuations, using Monte Carlo mutations to change lipid structures within the semigrand-canonical ensemble (at a fixed difference in component chemical potentials, Δμ), was recently implemented to address this challenge. This technique was applied here to mixtures of dimyristoylphosphatidylcholine and a shorter-tail lipid, either symmetric (didecanoylphosphatidylcholine (DDPC)) or asymmetric (hexanoyl-myristoylphosphatidylcholine), arranged in two types of structure: bilayer ribbons and buckled bilayers. In ribbons, the shorter-tail component showed a clear enrichment at the highly curved rim, more so for hexanoyl-myristoylphosphatidylcholine than for DDPC. Results on buckled bilayers were variable. Overall, the DDPC content of buckled bilayers tended to exceed by several percent the DDPC content of flat ones simulated at the same Δμ, but only for mixtures with low overall DDPC content. Within the buckled bilayer structure, no correlation could be resolved between the sign or magnitude of the local curvature of a leaflet and the mean local lipid composition. Results are discussed in terms of packing constraints, surface area/volume ratios, and curvature elasticity.